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(244.576,241.572,240.541,236.089,236.997,264.579,270.349,269.645,267.037,258.113,262.813,267.413,267.366,264.777,258.863,254.844,254.868,277.267,285.351,286.602,283.042,276.687,277.915,277.128,277.103,275.037,270.150,267.140,264.993,287.259,291.186,292.300,288.186,281.477,282.656,280.190,280.408,276.836,275.216,274.352,271.311,289.802,290.726,292.300,278.506,269.826,265.861,269.034,264.176,255.198,253.353,246.057,235.372,258.556,260.993,254.663,250.643,243.422,247.105,248.541,245.039,237.080,237.085,225.554,226.839,247.934,248.333,246.969,245.098,246.263,255.765,264.319,268.347,273.046,273.963,267.430,271.993,292.710,295.881) > #'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] -9.912032e-03 -9.866493e-03 -9.820938e-03 -9.775368e-03 -9.729782e-03 [6] -9.684180e-03 -9.638563e-03 -9.592930e-03 -9.547282e-03 -9.501618e-03 [11] -9.455938e-03 -9.410243e-03 -9.364533e-03 -9.318806e-03 -9.273064e-03 [16] -9.227307e-03 -9.181534e-03 -9.135745e-03 -9.089941e-03 -9.044121e-03 [21] -8.998286e-03 -8.952435e-03 -8.906568e-03 -8.860686e-03 -8.814789e-03 [26] -8.768875e-03 -8.722947e-03 -8.677002e-03 -8.631043e-03 -8.585067e-03 [31] -8.539076e-03 -8.493070e-03 -8.447048e-03 -8.401010e-03 -8.354957e-03 [36] -8.308888e-03 -8.262804e-03 -8.216704e-03 -8.170589e-03 -8.124458e-03 [41] -8.078312e-03 -8.032150e-03 -7.985972e-03 -7.939779e-03 -7.893571e-03 [46] -7.847347e-03 -7.801107e-03 -7.754852e-03 -7.708582e-03 -7.662296e-03 [51] -7.615994e-03 -7.569677e-03 -7.523345e-03 -7.476997e-03 -7.430633e-03 [56] -7.384254e-03 -7.337860e-03 -7.291450e-03 -7.245024e-03 -7.198583e-03 [61] -7.152127e-03 -7.105655e-03 -7.059167e-03 -7.012665e-03 -6.966146e-03 [66] -6.919612e-03 -6.873063e-03 -6.826499e-03 -6.779918e-03 -6.733323e-03 [71] -6.686712e-03 -6.640085e-03 -6.593443e-03 -6.546786e-03 -6.500113e-03 [76] -6.453425e-03 -6.406721e-03 -6.360002e-03 -6.313268e-03 -6.266518e-03 [81] -6.219752e-03 -6.172971e-03 -6.126175e-03 -6.079364e-03 -6.032537e-03 [86] -5.985694e-03 -5.938836e-03 -5.891963e-03 -5.845075e-03 -5.798171e-03 [91] -5.751251e-03 -5.704317e-03 -5.657367e-03 -5.610401e-03 -5.563420e-03 [96] -5.516424e-03 -5.469413e-03 -5.422386e-03 -5.375343e-03 -5.328286e-03 [101] -5.281213e-03 -5.234124e-03 -5.187021e-03 -5.139902e-03 -5.092767e-03 [106] -5.045618e-03 -4.998453e-03 -4.951273e-03 -4.904077e-03 -4.856866e-03 [111] -4.809640e-03 -4.762398e-03 -4.715142e-03 -4.667869e-03 -4.620582e-03 [116] -4.573279e-03 -4.525961e-03 -4.478628e-03 -4.431280e-03 -4.383916e-03 [121] -4.336537e-03 -4.289142e-03 -4.241733e-03 -4.194308e-03 -4.146868e-03 [126] -4.099413e-03 -4.051942e-03 -4.004456e-03 -3.956955e-03 -3.909439e-03 [131] -3.861907e-03 -3.814361e-03 -3.766799e-03 -3.719222e-03 -3.671629e-03 [136] -3.624022e-03 -3.576399e-03 -3.528761e-03 -3.481108e-03 -3.433440e-03 [141] -3.385756e-03 -3.338058e-03 -3.290344e-03 -3.242615e-03 -3.194871e-03 [146] -3.147112e-03 -3.099337e-03 -3.051548e-03 -3.003743e-03 -2.955923e-03 [151] -2.908088e-03 -2.860238e-03 -2.812373e-03 -2.764493e-03 -2.716597e-03 [156] -2.668687e-03 -2.620761e-03 -2.572821e-03 -2.524865e-03 -2.476894e-03 [161] -2.428908e-03 -2.380907e-03 -2.332891e-03 -2.284860e-03 -2.236814e-03 [166] -2.188752e-03 -2.140676e-03 -2.092585e-03 -2.044478e-03 -1.996357e-03 [171] -1.948221e-03 -1.900069e-03 -1.851903e-03 -1.803721e-03 -1.755525e-03 [176] -1.707313e-03 -1.659087e-03 -1.610845e-03 -1.562589e-03 -1.514318e-03 [181] -1.466031e-03 -1.417730e-03 -1.369414e-03 -1.321083e-03 -1.272736e-03 [186] -1.224375e-03 -1.175999e-03 -1.127608e-03 -1.079203e-03 -1.030782e-03 [191] -9.823461e-04 -9.338955e-04 -8.854300e-04 -8.369496e-04 -7.884543e-04 [196] -7.399442e-04 -6.914191e-04 -6.428793e-04 -5.943245e-04 -5.457549e-04 [201] -4.971704e-04 -4.485711e-04 -3.999570e-04 -3.513280e-04 -3.026842e-04 [206] -2.540256e-04 -2.053522e-04 -1.566640e-04 -1.079610e-04 -5.924315e-05 [211] -1.051054e-05 3.823684e-05 8.699899e-05 1.357759e-04 1.845676e-04 [216] 2.333740e-04 2.821952e-04 3.310311e-04 3.798818e-04 4.287472e-04 [221] 4.776273e-04 5.265221e-04 5.754316e-04 6.243559e-04 6.732948e-04 [226] 7.222484e-04 7.712167e-04 8.201997e-04 8.691973e-04 9.182096e-04 [231] 9.672365e-04 1.016278e-03 1.065334e-03 1.114405e-03 1.163491e-03 [236] 1.212591e-03 1.261706e-03 1.310835e-03 1.359979e-03 1.409137e-03 [241] 1.458310e-03 1.507498e-03 1.556700e-03 1.605917e-03 1.655148e-03 [246] 1.704394e-03 1.753655e-03 1.802930e-03 1.852219e-03 1.901523e-03 [251] 1.950842e-03 2.000175e-03 2.049522e-03 2.098884e-03 2.148261e-03 [256] 2.197652e-03 2.247057e-03 2.296477e-03 2.345912e-03 2.395360e-03 [261] 2.444824e-03 2.494301e-03 2.543793e-03 2.593300e-03 2.642821e-03 [266] 2.692356e-03 2.741906e-03 2.791470e-03 2.841048e-03 2.890641e-03 [271] 2.940248e-03 2.989870e-03 3.039506e-03 3.089156e-03 3.138820e-03 [276] 3.188499e-03 3.238192e-03 3.287900e-03 3.337622e-03 3.387358e-03 [281] 3.437108e-03 3.486873e-03 3.536652e-03 3.586445e-03 3.636252e-03 [286] 3.686074e-03 3.735910e-03 3.785760e-03 3.835624e-03 3.885502e-03 [291] 3.935395e-03 3.985302e-03 4.035223e-03 4.085158e-03 4.135108e-03 [296] 4.185071e-03 4.235049e-03 4.285041e-03 4.335047e-03 4.385067e-03 [301] 4.435101e-03 4.485150e-03 4.535212e-03 4.585288e-03 4.635379e-03 [306] 4.685484e-03 4.735602e-03 4.785735e-03 4.835882e-03 4.886043e-03 [311] 4.936218e-03 4.986407e-03 5.036610e-03 5.086827e-03 5.137058e-03 [316] 5.187302e-03 5.237561e-03 5.287834e-03 5.338121e-03 5.388422e-03 [321] 5.438737e-03 5.489065e-03 5.539408e-03 5.589764e-03 5.640135e-03 [326] 5.690519e-03 5.740917e-03 5.791329e-03 5.841755e-03 5.892195e-03 [331] 5.942648e-03 5.993116e-03 6.043597e-03 6.094092e-03 6.144601e-03 [336] 6.195124e-03 6.245660e-03 6.296211e-03 6.346775e-03 6.397353e-03 [341] 6.447944e-03 6.498550e-03 6.549169e-03 6.599801e-03 6.650448e-03 [346] 6.701108e-03 6.751782e-03 6.802470e-03 6.853171e-03 6.903886e-03 [351] 6.954615e-03 7.005357e-03 7.056113e-03 7.106883e-03 7.157666e-03 [356] 7.208463e-03 7.259273e-03 7.310097e-03 7.360935e-03 7.411786e-03 [361] 7.462651e-03 7.513529e-03 7.564421e-03 7.615326e-03 7.666245e-03 [366] 7.717177e-03 7.768123e-03 7.819083e-03 7.870056e-03 7.921042e-03 [371] 7.972042e-03 8.023055e-03 8.074082e-03 8.125122e-03 8.176176e-03 [376] 8.227243e-03 8.278323e-03 8.329417e-03 8.380524e-03 8.431645e-03 [381] 8.482779e-03 8.533926e-03 8.585087e-03 8.636261e-03 8.687448e-03 [386] 8.738649e-03 8.789863e-03 8.841090e-03 8.892330e-03 8.943584e-03 [391] 8.994851e-03 9.046131e-03 9.097425e-03 9.148732e-03 9.200052e-03 [396] 9.251385e-03 9.302731e-03 9.354091e-03 9.405463e-03 9.456849e-03 [401] 9.508248e-03 > mx [1] 0.009508248 > mxli [1] 2 > if (mxli != 0) + { + x1 <- (x^mxli - 1) / mxli + } else { + x1 <- log(x) + } > postscript(file="/var/www/html/rcomp/tmp/1f51n1259611015.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/2oeev1259611015.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/3ua1j1259611015.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/4j4vt1259611015.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/5nzgn1259611015.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/6pb1l1259611015.tab") > > system("convert tmp/1f51n1259611015.ps tmp/1f51n1259611015.png") > system("convert tmp/2oeev1259611015.ps tmp/2oeev1259611015.png") > system("convert tmp/3ua1j1259611015.ps tmp/3ua1j1259611015.png") > system("convert tmp/4j4vt1259611015.ps tmp/4j4vt1259611015.png") > system("convert tmp/5nzgn1259611015.ps tmp/5nzgn1259611015.png") > > > proc.time() user system elapsed 1.068 0.779 1.313