Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
100 87.14054095 112.0054296 112.312101 109.474134 104.9746116 100.4926851 104.2154743 120.1768388 112.1028355 108.1481575 116.802197 102.1699512 95.15358705 120.6707808 111.5234277 119.9669448 113.3697401 110.0717661 111.5567342 132.424212 107.900558 122.1626615 124.3992258 110.4450505 101.5874013 122.3203962 125.2582826 125.4411543 108.9902468 118.9243879 116.7242723 134.1724901 116.8530994 124.5732995 130.9914031 123.4239103 111.4536725 124.5135991 139.2589613 129.8596099 112.3460359 131.381655 133.0004776 134.3220552 144.2379719 134.1278719 150.1891559 140.722563 114.8389975 143.1973003 140.2738676 112.1248303 102.8951536 100.5090242 103.3513901 111.4134533 104.5887587 101.7840983 114.7007441 108.7426474
Type of transformation
Full Box-Cox transform
Simple Box-Cox transform
Minimum lambda
-2
-8
-7
-6
-5
-4
-3
-2
-1
Maximum lambda
2
1
2
3
4
5
6
7
8
Constant term to be added before analysis is performed
(?)
Display table with original and transformed data?
No
Yes
Chart options
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(qnorm(ppoints(x), mean=0, sd=1),x1) if (mx < c[i]) { mx <- c[i] mxli <- l[i] } } c mx mxli if (mxli != 0) { x1 <- (x^mxli - 1) / mxli } else { x1 <- log(x) } bitmap(file='test1.png') plot(l,c,main='Box-Cox Normality Plot',xlab='Lambda',ylab='correlation') mtext(paste('Optimal Lambda =',mxli)) grid() dev.off() bitmap(file='test2.png') hist(x,main='Histogram of Original Data',xlab='X',ylab='frequency') grid() dev.off() bitmap(file='test3.png') hist(x1,main='Histogram of Transformed Data',xlab='X',ylab='frequency') grid() dev.off() bitmap(file='test4.png') qqnorm(x) qqline(x) grid() mtext('Original Data') dev.off() bitmap(file='test5.png') qqnorm(x1) qqline(x1) grid() mtext('Transformed Data') dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Box-Cox Normality 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',header=TRUE) a<-table.element(a,mxli) 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