Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
107.11 107.57 107.81 108.75 109.43 109.62 109.54 109.53 109.84 109.67 109.79 109.56 110.22 110.40 110.69 110.72 110.89 110.58 110.94 110.91 111.22 111.09 111.00 111.06 111.55 112.32 112.64 112.36 112.04 112.37 112.59 112.89 113.22 112.85 113.06 112.99 113.32 113.74 113.91 114.52 114.96 114.91 115.30 115.44 115.52 116.08 115.94 115.56 115.88 116.66 117.41 117.68 117.85 118.21 118.92 119.03 119.17 118.95 118.92 118.90
Data Y:
107.56 107.70 107.67 107.67 107.72 108.35 108.25 108.26 108.31 108.33 108.36 108.36 108.97 109.62 109.60 109.64 109.65 109.64 109.93 109.81 109.77 110.10 110.40 110.50 111.89 112.10 111.92 112.15 112.16 112.17 112.32 112.38 112.34 113.14 113.18 113.21 113.76 113.99 113.95 113.93 114.01 114.10 114.11 114.10 114.12 114.68 114.71 114.73 115.81 116.01 116.12 116.49 116.51 116.60 117.01 117.01 117.12 117.22 118.38 118.80
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