R version 2.6.0 (2007-10-03) Copyright (C) 2007 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. > y <- c(72.50,59.40,85.70,88.20,62.80,87.00,79.20,112.00,79.20,132.10,40.10,69.00,59.40,73.80,57.40,81.10,46.60,41.40,71.20,67.90,72.00,145.50,39.70,51.90,73.70,70.90,60.80,61.00,54.50,39.10,66.60,58.50,59.80,80.90,37.30,44.60,48.70,54.00,49.50,61.60,35.00,35.70,51.30,49.00,41.50,72.50,42.10,44.10,45.10,50.30,40.90,47.20,36.90,40.90,38.30,46.30,28.40,78.40,36.80,50.70,42.80) > x <- c(109.20,88.60,94.30,98.30,86.40,80.60,104.10,108.20,93.40,71.90,94.10,94.90,96.40,91.10,84.40,86.40,88.00,75.10,109.70,103.00,82.10,68.00,96.40,94.30,90.00,88.00,76.10,82.50,81.40,66.50,97.20,94.10,80.70,70.50,87.80,89.50,99.60,84.20,75.10,92.00,80.80,73.10,99.80,90.00,83.10,72.40,78.80,87.30,91.00,80.10,73.60,86.40,74.50,71.20,92.40,81.50,85.30,69.90,84.20,90.70,100.30) > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: Wessa P., (2007), Box-Cox Linearity Plot (v1.0.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_boxcoxlin.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description > 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 [1] -4.157716e-02 -4.130689e-02 -4.103649e-02 -4.076597e-02 -4.049532e-02 [6] -4.022455e-02 -3.995366e-02 -3.968265e-02 -3.941151e-02 -3.914025e-02 [11] -3.886887e-02 -3.859736e-02 -3.832574e-02 -3.805400e-02 -3.778214e-02 [16] -3.751016e-02 -3.723806e-02 -3.696584e-02 -3.669351e-02 -3.642106e-02 [21] -3.614850e-02 -3.587582e-02 -3.560302e-02 -3.533011e-02 -3.505709e-02 [26] -3.478395e-02 -3.451070e-02 -3.423734e-02 -3.396387e-02 -3.369029e-02 [31] -3.341659e-02 -3.314279e-02 -3.286887e-02 -3.259485e-02 -3.232072e-02 [36] -3.204648e-02 -3.177213e-02 -3.149768e-02 -3.122312e-02 -3.094846e-02 [41] -3.067369e-02 -3.039881e-02 -3.012384e-02 -2.984876e-02 -2.957357e-02 [46] -2.929829e-02 -2.902290e-02 -2.874741e-02 -2.847182e-02 -2.819613e-02 [51] -2.792034e-02 -2.764446e-02 -2.736847e-02 -2.709239e-02 -2.681621e-02 [56] -2.653993e-02 -2.626356e-02 -2.598709e-02 -2.571053e-02 -2.543387e-02 [61] -2.515712e-02 -2.488027e-02 -2.460334e-02 -2.432631e-02 -2.404919e-02 [66] -2.377198e-02 -2.349468e-02 -2.321729e-02 -2.293981e-02 -2.266224e-02 [71] -2.238459e-02 -2.210685e-02 -2.182902e-02 -2.155110e-02 -2.127310e-02 [76] -2.099502e-02 -2.071685e-02 -2.043859e-02 -2.016026e-02 -1.988184e-02 [81] -1.960334e-02 -1.932476e-02 -1.904609e-02 -1.876735e-02 -1.848853e-02 [86] -1.820963e-02 -1.793065e-02 -1.765159e-02 -1.737246e-02 -1.709324e-02 [91] -1.681396e-02 -1.653460e-02 -1.625516e-02 -1.597565e-02 -1.569606e-02 [96] -1.541641e-02 -1.513668e-02 -1.485687e-02 -1.457700e-02 -1.429706e-02 [101] -1.401705e-02 -1.373696e-02 -1.345681e-02 -1.317660e-02 -1.289631e-02 [106] -1.261596e-02 -1.233554e-02 -1.205505e-02 -1.177451e-02 -1.149389e-02 [111] -1.121321e-02 -1.093247e-02 -1.065167e-02 -1.037081e-02 -1.008988e-02 [116] -9.808896e-03 -9.527850e-03 -9.246745e-03 -8.965581e-03 -8.684358e-03 [121] -8.403079e-03 -8.121742e-03 -7.840350e-03 -7.558902e-03 -7.277399e-03 [126] -6.995841e-03 -6.714230e-03 -6.432566e-03 -6.150850e-03 -5.869081e-03 [131] -5.587262e-03 -5.305392e-03 -5.023472e-03 -4.741503e-03 -4.459486e-03 [136] -4.177420e-03 -3.895307e-03 -3.613147e-03 -3.330942e-03 -3.048691e-03 [141] -2.766395e-03 -2.484055e-03 -2.201671e-03 -1.919245e-03 -1.636776e-03 [146] -1.354266e-03 -1.071715e-03 -7.891238e-04 -5.064928e-04 -2.238228e-04 [151] 5.888553e-05 3.416316e-04 6.244147e-04 9.072343e-04 1.190090e-03 [156] 1.472980e-03 1.755905e-03 2.038864e-03 2.321856e-03 2.604880e-03 [161] 2.887936e-03 3.171024e-03 3.454142e-03 3.737290e-03 4.020467e-03 [166] 4.303673e-03 4.586907e-03 4.870168e-03 5.153456e-03 5.436770e-03 [171] 5.720109e-03 6.003473e-03 6.286861e-03 6.570273e-03 6.853707e-03 [176] 7.137163e-03 7.420640e-03 7.704139e-03 7.987657e-03 8.271195e-03 [181] 8.554752e-03 8.838327e-03 9.121919e-03 9.405528e-03 9.689153e-03 [186] 9.972794e-03 1.025645e-02 1.054012e-02 1.082380e-02 1.110750e-02 [191] 1.139121e-02 1.167493e-02 1.195866e-02 1.224240e-02 1.252615e-02 [196] 1.280991e-02 1.309367e-02 1.337745e-02 1.366123e-02 1.394502e-02 [201] 1.422881e-02 1.451261e-02 1.479641e-02 1.508022e-02 1.536403e-02 [206] 1.564784e-02 1.593165e-02 1.621547e-02 1.649928e-02 1.678309e-02 [211] 1.706690e-02 1.735072e-02 1.763452e-02 1.791833e-02 1.820213e-02 [216] 1.848593e-02 1.876972e-02 1.905350e-02 1.933728e-02 1.962105e-02 [221] 1.990482e-02 2.018857e-02 2.047232e-02 2.075606e-02 2.103978e-02 [226] 2.132350e-02 2.160720e-02 2.189089e-02 2.217457e-02 2.245823e-02 [231] 2.274188e-02 2.302552e-02 2.330913e-02 2.359274e-02 2.387632e-02 [236] 2.415989e-02 2.444344e-02 2.472696e-02 2.501047e-02 2.529396e-02 [241] 2.557743e-02 2.586087e-02 2.614430e-02 2.642770e-02 2.671107e-02 [246] 2.699442e-02 2.727775e-02 2.756105e-02 2.784432e-02 2.812757e-02 [251] 2.841079e-02 2.869398e-02 2.897714e-02 2.926027e-02 2.954337e-02 [256] 2.982643e-02 3.010947e-02 3.039247e-02 3.067545e-02 3.095838e-02 [261] 3.124128e-02 3.152415e-02 3.180698e-02 3.208977e-02 3.237253e-02 [266] 3.265525e-02 3.293793e-02 3.322057e-02 3.350317e-02 3.378573e-02 [271] 3.406825e-02 3.435073e-02 3.463316e-02 3.491555e-02 3.519790e-02 [276] 3.548020e-02 3.576246e-02 3.604467e-02 3.632683e-02 3.660895e-02 [281] 3.689102e-02 3.717304e-02 3.745501e-02 3.773693e-02 3.801880e-02 [286] 3.830062e-02 3.858238e-02 3.886410e-02 3.914576e-02 3.942736e-02 [291] 3.970891e-02 3.999041e-02 4.027185e-02 4.055324e-02 4.083456e-02 [296] 4.111583e-02 4.139704e-02 4.167819e-02 4.195928e-02 4.224031e-02 [301] 4.252128e-02 4.280219e-02 4.308303e-02 4.336381e-02 4.364453e-02 [306] 4.392518e-02 4.420577e-02 4.448629e-02 4.476674e-02 4.504713e-02 [311] 4.532745e-02 4.560770e-02 4.588788e-02 4.616799e-02 4.644803e-02 [316] 4.672800e-02 4.700790e-02 4.728773e-02 4.756748e-02 4.784716e-02 [321] 4.812676e-02 4.840629e-02 4.868575e-02 4.896512e-02 4.924442e-02 [326] 4.952364e-02 4.980279e-02 5.008185e-02 5.036084e-02 5.063974e-02 [331] 5.091857e-02 5.119731e-02 5.147597e-02 5.175455e-02 5.203304e-02 [336] 5.231145e-02 5.258977e-02 5.286801e-02 5.314617e-02 5.342423e-02 [341] 5.370221e-02 5.398010e-02 5.425790e-02 5.453562e-02 5.481324e-02 [346] 5.509077e-02 5.536821e-02 5.564556e-02 5.592282e-02 5.619998e-02 [351] 5.647705e-02 5.675403e-02 5.703091e-02 5.730769e-02 5.758438e-02 [356] 5.786097e-02 5.813746e-02 5.841386e-02 5.869016e-02 5.896635e-02 [361] 5.924245e-02 5.951844e-02 5.979434e-02 6.007013e-02 6.034582e-02 [366] 6.062141e-02 6.089689e-02 6.117227e-02 6.144754e-02 6.172271e-02 [371] 6.199777e-02 6.227272e-02 6.254757e-02 6.282231e-02 6.309694e-02 [376] 6.337146e-02 6.364587e-02 6.392017e-02 6.419435e-02 6.446843e-02 [381] 6.474239e-02 6.501624e-02 6.528998e-02 6.556360e-02 6.583711e-02 [386] 6.611050e-02 6.638377e-02 6.665693e-02 6.692997e-02 6.720289e-02 [391] 6.747570e-02 6.774838e-02 6.802095e-02 6.829339e-02 6.856571e-02 [396] 6.883791e-02 6.910999e-02 6.938195e-02 6.965378e-02 6.992549e-02 [401] 7.019708e-02 > mx [1] 0.07019708 > mxli [1] 2 > 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)) > postscript(file="/var/www/html/rcomp/tmp/1erbh1194302862.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(l,c,main='Box-Cox Linearity Plot',xlab='Lambda',ylab='correlation') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/20z1v1194302862.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x,y,main='Linear Fit of Original Data',xlab='x',ylab='y') > abline(r) > grid() > mtext(paste('Residual Standard Deviation = ',se)) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3pb1t1194302862.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x1,y,main='Linear Fit of Transformed Data',xlab='x',ylab='y') > abline(r1) > grid() > mtext(paste('Residual Standard Deviation = ',se1)) > dev.off() null device 1 > load(file='/var/www/html/rcomp/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="/var/www/html/rcomp/tmp/4yabe1194302862.tab") > > system("convert tmp/1erbh1194302862.ps tmp/1erbh1194302862.png") > system("convert tmp/20z1v1194302862.ps tmp/20z1v1194302862.png") > system("convert tmp/3pb1t1194302862.ps tmp/3pb1t1194302862.png") > > > proc.time() user system elapsed 1.059 0.519 1.256