R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(18004,17537,20366,22782,19169,13807,29743,25591,29096,26482,22405,27044,17970,18730,19684,19785,18479,10698,31956,29506,34506,27165,26736,23691,18157,17328,18205,20995,17382,9367,31124,26551,30651,25859,25100,25778,20418,18688,20424,24776,19814,12738,31566,30111,30019,31934,25826,26835,20205,17789,20520,22518,15572,11509,25447,24090,27786,26195,20516,22759,19028,16971,20036,22485,18730) > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > 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 [1] 2.040749e-02 2.037303e-02 2.033821e-02 2.030304e-02 2.026750e-02 [6] 2.023160e-02 2.019535e-02 2.015875e-02 2.012177e-02 2.008444e-02 [11] 2.004674e-02 2.000867e-02 1.997024e-02 1.993145e-02 1.989228e-02 [16] 1.985275e-02 1.981284e-02 1.977256e-02 1.973191e-02 1.969088e-02 [21] 1.964948e-02 1.960770e-02 1.956555e-02 1.952302e-02 1.948011e-02 [26] 1.943682e-02 1.939315e-02 1.934910e-02 1.930466e-02 1.925984e-02 [31] 1.921464e-02 1.916905e-02 1.912307e-02 1.907671e-02 1.902996e-02 [36] 1.898282e-02 1.893530e-02 1.888738e-02 1.883907e-02 1.879037e-02 [41] 1.874128e-02 1.869180e-02 1.864192e-02 1.859165e-02 1.854098e-02 [46] 1.848992e-02 1.843847e-02 1.838662e-02 1.833437e-02 1.828172e-02 [51] 1.822868e-02 1.817524e-02 1.812140e-02 1.806716e-02 1.801252e-02 [56] 1.795748e-02 1.790205e-02 1.784621e-02 1.778997e-02 1.773333e-02 [61] 1.767629e-02 1.761885e-02 1.756101e-02 1.750277e-02 1.744412e-02 [66] 1.738507e-02 1.732562e-02 1.726577e-02 1.720552e-02 1.714486e-02 [71] 1.708380e-02 1.702234e-02 1.696048e-02 1.689821e-02 1.683554e-02 [76] 1.677247e-02 1.670900e-02 1.664513e-02 1.658085e-02 1.651617e-02 [81] 1.645109e-02 1.638561e-02 1.631973e-02 1.625345e-02 1.618677e-02 [86] 1.611969e-02 1.605221e-02 1.598432e-02 1.591604e-02 1.584737e-02 [91] 1.577829e-02 1.570881e-02 1.563894e-02 1.556867e-02 1.549801e-02 [96] 1.542695e-02 1.535549e-02 1.528364e-02 1.521140e-02 1.513876e-02 [101] 1.506573e-02 1.499231e-02 1.491849e-02 1.484429e-02 1.476969e-02 [106] 1.469471e-02 1.461934e-02 1.454358e-02 1.446743e-02 1.439090e-02 [111] 1.431399e-02 1.423668e-02 1.415900e-02 1.408093e-02 1.400249e-02 [116] 1.392366e-02 1.384445e-02 1.376486e-02 1.368490e-02 1.360456e-02 [121] 1.352384e-02 1.344275e-02 1.336129e-02 1.327946e-02 1.319725e-02 [126] 1.311468e-02 1.303173e-02 1.294842e-02 1.286475e-02 1.278071e-02 [131] 1.269630e-02 1.261154e-02 1.252641e-02 1.244092e-02 1.235507e-02 [136] 1.226887e-02 1.218232e-02 1.209540e-02 1.200814e-02 1.192052e-02 [141] 1.183256e-02 1.174424e-02 1.165558e-02 1.156658e-02 1.147723e-02 [146] 1.138753e-02 1.129750e-02 1.120713e-02 1.111642e-02 1.102537e-02 [151] 1.093399e-02 1.084228e-02 1.075023e-02 1.065786e-02 1.056516e-02 [156] 1.047213e-02 1.037877e-02 1.028510e-02 1.019110e-02 1.009678e-02 [161] 1.000215e-02 9.907201e-03 9.811937e-03 9.716361e-03 9.620474e-03 [166] 9.524278e-03 9.427776e-03 9.330968e-03 9.233858e-03 9.136446e-03 [171] 9.038734e-03 8.940725e-03 8.842420e-03 8.743822e-03 8.644931e-03 [176] 8.545751e-03 8.446282e-03 8.346528e-03 8.246489e-03 8.146168e-03 [181] 8.045567e-03 7.944687e-03 7.843531e-03 7.742101e-03 7.640399e-03 [186] 7.538427e-03 7.436186e-03 7.333679e-03 7.230908e-03 7.127874e-03 [191] 7.024581e-03 6.921029e-03 6.817222e-03 6.713161e-03 6.608847e-03 [196] 6.504284e-03 6.399474e-03 6.294417e-03 6.189118e-03 6.083576e-03 [201] 5.977796e-03 5.871778e-03 5.765526e-03 5.659040e-03 5.552323e-03 [206] 5.445378e-03 5.338206e-03 5.230810e-03 5.123191e-03 5.015352e-03 [211] 4.907296e-03 4.799023e-03 4.690536e-03 4.581838e-03 4.472931e-03 [216] 4.363816e-03 4.254496e-03 4.144973e-03 4.035249e-03 3.925327e-03 [221] 3.815208e-03 3.704894e-03 3.594389e-03 3.483693e-03 3.372810e-03 [226] 3.261741e-03 3.150488e-03 3.039054e-03 2.927441e-03 2.815650e-03 [231] 2.703685e-03 2.591547e-03 2.479238e-03 2.366761e-03 2.254118e-03 [236] 2.141310e-03 2.028340e-03 1.915211e-03 1.801924e-03 1.688481e-03 [241] 1.574885e-03 1.461138e-03 1.347241e-03 1.233198e-03 1.119009e-03 [246] 1.004678e-03 8.902066e-04 7.755964e-04 6.608497e-04 5.459689e-04 [251] 4.309558e-04 3.158127e-04 2.005416e-04 8.514474e-05 -3.037593e-05 [256] -1.460183e-04 -2.617802e-04 -3.776596e-04 -4.936545e-04 -6.097627e-04 [261] -7.259822e-04 -8.423110e-04 -9.587469e-04 -1.075288e-03 -1.191932e-03 [266] -1.308677e-03 -1.425521e-03 -1.542463e-03 -1.659499e-03 -1.776628e-03 [271] -1.893848e-03 -2.011157e-03 -2.128553e-03 -2.246033e-03 -2.363597e-03 [276] -2.481242e-03 -2.598966e-03 -2.716766e-03 -2.834642e-03 -2.952591e-03 [281] -3.070611e-03 -3.188701e-03 -3.306858e-03 -3.425080e-03 -3.543365e-03 [286] -3.661712e-03 -3.780119e-03 -3.898584e-03 -4.017104e-03 -4.135678e-03 [291] -4.254305e-03 -4.372981e-03 -4.491707e-03 -4.610478e-03 -4.729294e-03 [296] -4.848154e-03 -4.967054e-03 -5.085993e-03 -5.204970e-03 -5.323982e-03 [301] -5.443029e-03 -5.562107e-03 -5.681215e-03 -5.800352e-03 -5.919516e-03 [306] -6.038704e-03 -6.157916e-03 -6.277149e-03 -6.396402e-03 -6.515673e-03 [311] -6.634960e-03 -6.754262e-03 -6.873577e-03 -6.992902e-03 -7.112238e-03 [316] -7.231581e-03 -7.350930e-03 -7.470284e-03 -7.589641e-03 -7.708999e-03 [321] -7.828356e-03 -7.947712e-03 -8.067064e-03 -8.186411e-03 -8.305751e-03 [326] -8.425082e-03 -8.544404e-03 -8.663714e-03 -8.783010e-03 -8.902293e-03 [331] -9.021559e-03 -9.140807e-03 -9.260036e-03 -9.379244e-03 -9.498430e-03 [336] -9.617592e-03 -9.736729e-03 -9.855839e-03 -9.974921e-03 -1.009397e-02 [341] -1.021299e-02 -1.033198e-02 -1.045094e-02 -1.056986e-02 -1.068874e-02 [346] -1.080758e-02 -1.092639e-02 -1.104515e-02 -1.116387e-02 -1.128255e-02 [351] -1.140118e-02 -1.151976e-02 -1.163830e-02 -1.175679e-02 -1.187522e-02 [356] -1.199361e-02 -1.211194e-02 -1.223022e-02 -1.234844e-02 -1.246660e-02 [361] -1.258471e-02 -1.270276e-02 -1.282074e-02 -1.293866e-02 -1.305652e-02 [366] -1.317432e-02 -1.329205e-02 -1.340971e-02 -1.352730e-02 -1.364482e-02 [371] -1.376228e-02 -1.387966e-02 -1.399697e-02 -1.411420e-02 -1.423136e-02 [376] -1.434844e-02 -1.446544e-02 -1.458237e-02 -1.469921e-02 -1.481597e-02 [381] -1.493265e-02 -1.504925e-02 -1.516576e-02 -1.528219e-02 -1.539853e-02 [386] -1.551478e-02 -1.563095e-02 -1.574702e-02 -1.586300e-02 -1.597889e-02 [391] -1.609469e-02 -1.621040e-02 -1.632600e-02 -1.644152e-02 -1.655693e-02 [396] -1.667225e-02 -1.678747e-02 -1.690259e-02 -1.701760e-02 -1.713252e-02 [401] -1.724733e-02 > mx [1] 0.02040749 > mxli [1] -2 > if (mxli != 0) + { + x1 <- (x^mxli - 1) / mxli + } else { + x1 <- log(x) + } > postscript(file="/var/www/html/rcomp/tmp/1oft51261329554.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(l,c,main='Box-Cox Normality Plot',xlab='Lambda',ylab='correlation') > mtext(paste('Optimal Lambda =',mxli)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/28yfc1261329554.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(x,main='Histogram of Original Data',xlab='X',ylab='frequency') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3qnxf1261329554.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(x1,main='Histogram of Transformed Data',xlab='X',ylab='frequency') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/40t471261329554.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(x) > qqline(x) > grid() > mtext('Original Data') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5pof91261329554.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(x1) > qqline(x1) > grid() > mtext('Transformed Data') > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/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="/var/www/html/rcomp/tmp/6532q1261329554.tab") > > try(system("convert tmp/1oft51261329554.ps tmp/1oft51261329554.png",intern=TRUE)) character(0) > try(system("convert tmp/28yfc1261329554.ps tmp/28yfc1261329554.png",intern=TRUE)) character(0) > try(system("convert tmp/3qnxf1261329554.ps tmp/3qnxf1261329554.png",intern=TRUE)) character(0) > try(system("convert tmp/40t471261329554.ps tmp/40t471261329554.png",intern=TRUE)) character(0) > try(system("convert tmp/5pof91261329554.ps tmp/5pof91261329554.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.053 0.759 1.260