Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
105.2 105.2 105.6 105.6 106.2 106.3 106.4 106.9 107.2 107.3 107.3 107.4 107.55 107.87 108.37 108.38 107.92 108.03 108.14 108.3 108.64 108.66 109.04 109.03 109.03 109.54 109.75 109.83 109.65 109.82 109.95 110.12 110.15 110.2 109.99 110.14 110.14 110.81 110.97 110.99 109.73 109.81 110.02 110.18 110.21 110.25 110.36 110.51 110.64 110.95 111.18 111.19 111.69 111.7 111.83 111.77 111.73 112.01 111.86 112.04
Data Y:
151.7 121.3 133.0 119.6 122.2 117.4 106.7 87.5 81.0 110.3 87.0 55.7 146.0 137.5 138.5 135.6 107.3 99.0 91.4 68.4 82.6 98.4 71.3 47.6 130.8 113.6 125.7 113.6 97.1 104.4 91.8 75.1 89.2 110.2 78.4 68.4 122.8 129.7 159.1 139.0 102.2 113.6 81.5 77.4 87.6 101.2 87.2 64.9 133.1 118.0 135.9 125.7 108.0 128.3 84.7 86.4 92.2 95.8 92.3 54.3
Sample Range:
(leave blank to include all observations)
From:
To:
Chart options
Label y-axis:
Label x-axis:
R Code
n <- length(x) c <- array(NA,dim=c(401)) l <- array(NA,dim=c(401)) mx <- 0 mxli <- -999 for (i in 1:401) { l[i] <- (i-201)/100 if (l[i] != 0) { x1 <- (x^l[i] - 1) / l[i] } else { x1 <- log(x) } c[i] <- cor(x1,y) if (mx < abs(c[i])) { mx <- abs(c[i]) mxli <- l[i] } } c mx mxli if (mxli != 0) { x1 <- (x^mxli - 1) / mxli } else { x1 <- log(x) } r<-lm(y~x) se <- sqrt(var(r$residuals)) r1 <- lm(y~x1) se1 <- sqrt(var(r1$residuals)) bitmap(file='test1.png') plot(l,c,main='Box-Cox Linearity Plot',xlab='Lambda',ylab='correlation') grid() dev.off() bitmap(file='test2.png') plot(x,y,main='Linear Fit of Original Data',xlab='x',ylab='y') abline(r) grid() mtext(paste('Residual Standard Deviation = ',se)) dev.off() bitmap(file='test3.png') plot(x1,y,main='Linear Fit of Transformed Data',xlab='x',ylab='y') abline(r1) grid() mtext(paste('Residual Standard Deviation = ',se1)) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Box-Cox Linearity Plot',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'# observations x',header=TRUE) a<-table.element(a,n) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'maximum correlation',header=TRUE) a<-table.element(a,mx) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'optimal lambda(x)',header=TRUE) a<-table.element(a,mxli) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Residual SD (orginial)',header=TRUE) a<-table.element(a,se) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Residual SD (transformed)',header=TRUE) a<-table.element(a,se1) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation