R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(255 + ,280.2 + ,299.9 + ,339.2 + ,374.2 + ,393.5 + ,389.2 + ,381.7 + ,375.2 + ,369 + ,357.4 + ,352.1 + ,346.5 + ,342.9 + ,340.3 + ,328.3 + ,322.9 + ,314.3 + ,308.9 + ,294 + ,285.6 + ,281.2 + ,280.3 + ,278.8 + ,274.5 + ,270.4 + ,263.4 + ,259.9 + ,258 + ,262.7 + ,284.7 + ,311.3 + ,322.1 + ,327 + ,331.3 + ,333.3 + ,321.4 + ,327 + ,320 + ,314.7 + ,316.7 + ,314.4 + ,321.3 + ,318.2 + ,307.2 + ,301.3 + ,287.5 + ,277.7 + ,274.4 + ,258.8 + ,253.3 + ,251 + ,248.4 + ,249.5 + ,246.1 + ,244.5 + ,243.6 + ,244 + ,240.8 + ,249.8 + ,248 + ,259.4 + ,260.5 + ,260.8 + ,261.3 + ,259.5 + ,256.6 + ,257.9 + ,256.5 + ,254.2 + ,253.3 + ,253.8 + ,255.5 + ,257.1 + ,257.3 + ,253.2 + ,252.8 + ,252 + ,250.7 + ,252.2 + ,250 + ,251 + ,253.4 + ,251.2 + ,255.6 + ,261.1 + ,258.9 + ,259.9 + ,261.2 + ,264.7 + ,267.1 + ,266.4 + ,267.7 + ,268.6 + ,267.5 + ,268.5 + ,268.5 + ,270.5 + ,270.9 + ,270.1 + ,269.3 + ,269.8 + ,270.1 + ,264.9 + ,263.7 + ,264.8 + ,263.7 + ,255.9 + ,276.2 + ,360.1 + ,380.5 + ,373.7 + ,369.8 + ,366.6 + ,359.3 + ,345.8 + ,326.2 + ,324.5 + ,328.1 + ,327.5 + ,324.4 + ,316.5 + ,310.9 + ,301.5 + ,291.7 + ,290.4 + ,287.4 + ,277.7 + ,281.6 + ,288 + ,276 + ,272.9 + ,283 + ,283.3 + ,276.8 + ,284.5 + ,282.7 + ,281.2 + ,287.4 + ,283.1 + ,284 + ,285.5 + ,289.2 + ,292.5 + ,296.4 + ,305.2 + ,303.9 + ,311.5 + ,316.3 + ,316.7 + ,322.5 + ,317.1 + ,309.8 + ,303.8 + ,290.3 + ,293.7 + ,291.7 + ,296.5 + ,289.1 + ,288.5 + ,293.8 + ,297.7 + ,305.4 + ,302.7 + ,302.5 + ,303 + ,294.5 + ,294.1 + ,294.5 + ,297.1 + ,289.4 + ,292.4 + ,287.9 + ,286.6 + ,280.5 + ,272.4 + ,269.2 + ,270.6 + ,267.3 + ,262.5 + ,266.8 + ,268.8 + ,263.1 + ,261.2 + ,266 + ,262.5 + ,265.2 + ,261.3 + ,253.7 + ,249.2 + ,239.1 + ,236.4 + ,235.2 + ,245.2 + ,246.2 + ,247.7 + ,251.4 + ,253.3 + ,254.8 + ,250 + ,249.3 + ,241.5 + ,243.3 + ,248 + ,253 + ,252.9 + ,251.5 + ,251.6 + ,253.5 + ,259.8 + ,334.1 + ,448 + ,445.8 + ,445 + ,448.2 + ,438.2 + ,439.8 + ,423.4 + ,410.8 + ,408.4 + ,406.7 + ,405.9 + ,402.7 + ,405.1 + ,399.6 + ,386.5 + ,381.4 + ,375.2 + ,357.7 + ,359 + ,355 + ,352.7 + ,344.4 + ,343.8 + ,338 + ,339 + ,333.3 + ,334.4 + ,328.3 + ,330.7 + ,330 + ,331.6 + ,351.2 + ,389.4 + ,410.9 + ,442.8 + ,462.8 + ,466.9 + ,461.7 + ,439.2 + ,430.3 + ,416.1 + ,402.5 + ,397.3 + ,403.3 + ,395.9 + ,387.8 + ,378.6 + ,377.1 + ,370.4 + ,362 + ,350.3 + ,348.2 + ,344.6 + ,343.5 + ,342.8 + ,347.6 + ,346.6 + ,349.5 + ,342.1 + ,342 + ,342.8 + ,339.3 + ,348.2 + ,333.7 + ,334.7 + ,354 + ,367.7 + ,363.3 + ,358.4 + ,353.1 + ,343.1 + ,344.6 + ,344.4 + ,333.9 + ,331.7 + ,324.3 + ,321.2 + ,322.4 + ,321.7 + ,320.5 + ,312.8 + ,309.7 + ,315.6 + ,309.7 + ,304.6 + ,302.5 + ,301.5 + ,298.8 + ,291.3 + ,293.6 + ,294.6 + ,285.9 + ,297.6 + ,301.1 + ,293.8 + ,297.7 + ,292.9 + ,292.1 + ,287.2 + ,288.2 + ,283.8 + ,299.9 + ,292.4 + ,293.3 + ,300.8 + ,293.7 + ,293.1 + ,294.4 + ,292.1 + ,291.9 + ,282.5 + ,277.9 + ,287.5 + ,289.2 + ,285.6 + ,293.2 + ,290.8 + ,283.1 + ,275 + ,287.8 + ,287.8 + ,287.4 + ,284 + ,277.8 + ,277.6 + ,304.9 + ,294 + ,300.9 + ,324 + ,332.9 + ,341.6 + ,333.4 + ,348.2 + ,344.7 + ,344.7 + ,329.3 + ,323.5 + ,323.2 + ,317.4 + ,330.1 + ,329.2 + ,334.9 + ,315.8 + ,315.4 + ,319.6 + ,317.3 + ,313.8 + ,315.8 + ,311.3) > x <- c(87.28 + ,87.28 + ,87.09 + ,86.92 + ,87.59 + ,90.72 + ,90.69 + ,90.3 + ,89.55 + ,88.94 + ,88.41 + ,87.82 + ,87.07 + ,86.82 + ,86.4 + ,86.02 + ,85.66 + ,85.32 + ,85 + ,84.67 + ,83.94 + ,82.83 + ,81.95 + ,81.19 + ,80.48 + ,78.86 + ,69.47 + ,68.77 + ,70.06 + ,73.95 + ,75.8 + ,77.79 + ,81.57 + ,83.07 + ,84.34 + ,85.1 + ,85.25 + ,84.26 + ,83.63 + ,86.44 + ,85.3 + ,84.1 + ,83.36 + ,82.48 + ,81.58 + ,80.47 + ,79.34 + ,82.13 + ,81.69 + ,80.7 + ,79.88 + ,79.16 + ,78.38 + ,77.42 + ,76.47 + ,75.46 + ,74.48 + ,78.27 + ,80.7 + ,79.91 + ,78.75 + ,77.78 + ,81.14 + ,81.08 + ,80.03 + ,78.91 + ,78.01 + ,76.9 + ,75.97 + ,81.93 + ,80.27 + ,78.67 + ,77.42 + ,76.16 + ,74.7 + ,76.39 + ,76.04 + ,74.65 + ,73.29 + ,71.79 + ,74.39 + ,74.91 + ,74.54 + ,73.08 + ,72.75 + ,71.32 + ,70.38 + ,70.35 + ,70.01 + ,69.36 + ,67.77 + ,69.26 + ,69.8 + ,68.38 + ,67.62 + ,68.39 + ,66.95 + ,65.21 + ,66.64 + ,63.45 + ,60.66 + ,62.34 + ,60.32 + ,58.64 + ,60.46 + ,58.59 + ,61.87 + ,61.85 + ,67.44 + ,77.06 + ,91.74 + ,93.15 + ,94.15 + ,93.11 + ,91.51 + ,89.96 + ,88.16 + ,86.98 + ,88.03 + ,86.24 + ,84.65 + ,83.23 + ,81.7 + ,80.25 + ,78.8 + ,77.51 + ,76.2 + ,75.04 + ,74 + ,75.49 + ,77.14 + ,76.15 + ,76.27 + ,78.19 + ,76.49 + ,77.31 + ,76.65 + ,74.99 + ,73.51 + ,72.07 + ,70.59 + ,71.96 + ,76.29 + ,74.86 + ,74.93 + ,71.9 + ,71.01 + ,77.47 + ,75.78 + ,76.6 + ,76.07 + ,74.57 + ,73.02 + ,72.65 + ,73.16 + ,71.53 + ,69.78 + ,67.98 + ,69.96 + ,72.16 + ,70.47 + ,68.86 + ,67.37 + ,65.87 + ,72.16 + ,71.34 + ,69.93 + ,68.44 + ,67.16 + ,66.01 + ,67.25 + ,70.91 + ,69.75 + ,68.59 + ,67.48 + ,66.31 + ,64.81 + ,66.58 + ,65.97 + ,64.7 + ,64.7 + ,60.94 + ,59.08 + ,58.42 + ,57.77 + ,57.11 + ,53.31 + ,49.96 + ,49.4 + ,48.84 + ,48.3 + ,47.74 + ,47.24 + ,46.76 + ,46.29 + ,48.9 + ,49.23 + ,48.53 + ,48.03 + ,54.34 + ,53.79 + ,53.24 + ,52.96 + ,52.17 + ,51.7 + ,58.55 + ,78.2 + ,77.03 + ,76.19 + ,77.15 + ,75.87 + ,95.47 + ,109.67 + ,112.28 + ,112.01 + ,107.93 + ,105.96 + ,105.06 + ,102.98 + ,102.2 + ,105.23 + ,101.85 + ,99.89 + ,96.23 + ,94.76 + ,91.51 + ,91.63 + ,91.54 + ,85.23 + ,87.83 + ,87.38 + ,84.44 + ,85.19 + ,84.03 + ,86.73 + ,102.52 + ,104.45 + ,106.98 + ,107.02 + ,99.26 + ,94.45 + ,113.44 + ,157.33 + ,147.38 + ,171.89 + ,171.95 + ,132.71 + ,126.02 + ,121.18 + ,115.45 + ,110.48 + ,117.85 + ,117.63 + ,124.65 + ,109.59 + ,111.27 + ,99.78 + ,98.21 + ,99.2 + ,97.97 + ,89.55 + ,87.91 + ,93.34 + ,94.42 + ,93.2 + ,90.29 + ,91.46 + ,89.98 + ,88.35 + ,88.41 + ,82.44 + ,79.89 + ,75.69 + ,75.66 + ,84.5 + ,96.73 + ,87.48 + ,82.39 + ,83.48 + ,79.31 + ,78.16 + ,72.77 + ,72.45 + ,68.46 + ,67.62 + ,68.76 + ,70.07 + ,68.55 + ,65.3 + ,58.96 + ,59.17 + ,62.37 + ,66.28 + ,55.62 + ,55.23 + ,55.85 + ,56.75 + ,50.89 + ,53.88 + ,52.95 + ,55.08 + ,53.61 + ,58.78 + ,61.85 + ,55.91 + ,53.32 + ,46.41 + ,44.57 + ,50 + ,50 + ,53.36 + ,46.23 + ,50.45 + ,49.07 + ,45.85 + ,48.45 + ,49.96 + ,46.53 + ,50.51 + ,47.58 + ,48.05 + ,46.84 + ,47.67 + ,49.16 + ,55.54 + ,55.82 + ,58.22 + ,56.19 + ,57.77 + ,63.19 + ,54.76 + ,55.74 + ,62.54 + ,61.39 + ,69.6 + ,79.23 + ,80 + ,93.68 + ,107.63 + ,100.18 + ,97.3 + ,90.45 + ,80.64 + ,80.58 + ,75.82 + ,85.59 + ,89.35 + ,89.42 + ,104.73 + ,95.32 + ,89.27 + ,90.44 + ,86.97 + ,79.98 + ,81.22 + ,87.35 + ,83.64 + ,82.22 + ,94.4 + ,102.18) > #'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(x1,y) + if (mx < abs(c[i])) + { + mx <- abs(c[i]) + mxli <- l[i] + } + } > c [1] 0.5538905 0.5545455 0.5552009 0.5558566 0.5565126 0.5571690 0.5578258 [8] 0.5584828 0.5591402 0.5597978 0.5604558 0.5611140 0.5617725 0.5624313 [15] 0.5630904 0.5637497 0.5644092 0.5650690 0.5657290 0.5663893 0.5670497 [22] 0.5677104 0.5683713 0.5690323 0.5696935 0.5703549 0.5710164 0.5716781 [29] 0.5723400 0.5730019 0.5736640 0.5743262 0.5749885 0.5756509 0.5763134 [36] 0.5769759 0.5776385 0.5783012 0.5789639 0.5796267 0.5802895 0.5809522 [43] 0.5816150 0.5822778 0.5829406 0.5836033 0.5842660 0.5849287 0.5855913 [50] 0.5862538 0.5869163 0.5875786 0.5882409 0.5889030 0.5895650 0.5902269 [57] 0.5908886 0.5915502 0.5922116 0.5928728 0.5935339 0.5941947 0.5948553 [64] 0.5955157 0.5961758 0.5968357 0.5974953 0.5981546 0.5988137 0.5994724 [71] 0.6001309 0.6007890 0.6014467 0.6021042 0.6027612 0.6034179 0.6040742 [78] 0.6047300 0.6053855 0.6060405 0.6066951 0.6073492 0.6080029 0.6086561 [85] 0.6093088 0.6099609 0.6106126 0.6112637 0.6119142 0.6125642 0.6132136 [92] 0.6138624 0.6145106 0.6151582 0.6158052 0.6164514 0.6170971 0.6177420 [99] 0.6183862 0.6190298 0.6196726 0.6203146 0.6209559 0.6215965 0.6222362 [106] 0.6228751 0.6235133 0.6241506 0.6247870 0.6254226 0.6260573 0.6266910 [113] 0.6273239 0.6279559 0.6285869 0.6292169 0.6298460 0.6304741 0.6311011 [120] 0.6317272 0.6323521 0.6329761 0.6335989 0.6342207 0.6348413 0.6354608 [127] 0.6360792 0.6366964 0.6373124 0.6379272 0.6385408 0.6391532 0.6397643 [134] 0.6403742 0.6409827 0.6415900 0.6421959 0.6428005 0.6434038 0.6440056 [141] 0.6446061 0.6452051 0.6458028 0.6463989 0.6469936 0.6475868 0.6481785 [148] 0.6487687 0.6493573 0.6499444 0.6505298 0.6511137 0.6516960 0.6522766 [155] 0.6528555 0.6534328 0.6540083 0.6545822 0.6551543 0.6557246 0.6562932 [162] 0.6568600 0.6574249 0.6579880 0.6585493 0.6591086 0.6596661 0.6602216 [169] 0.6607752 0.6613269 0.6618765 0.6624242 0.6629698 0.6635134 0.6640549 [176] 0.6645944 0.6651317 0.6656669 0.6661999 0.6667308 0.6672595 0.6677860 [183] 0.6683102 0.6688322 0.6693519 0.6698693 0.6703844 0.6708971 0.6714075 [190] 0.6719155 0.6724211 0.6729243 0.6734250 0.6739232 0.6744190 0.6749122 [197] 0.6754029 0.6758910 0.6763765 0.6768595 0.6773398 0.6778175 0.6782925 [204] 0.6787648 0.6792344 0.6797012 0.6801653 0.6806266 0.6810851 0.6815408 [211] 0.6819936 0.6824436 0.6828906 0.6833348 0.6837760 0.6842142 0.6846495 [218] 0.6850817 0.6855110 0.6859371 0.6863602 0.6867802 0.6871971 0.6876108 [225] 0.6880214 0.6884288 0.6888329 0.6892339 0.6896316 0.6900260 0.6904171 [232] 0.6908048 0.6911893 0.6915703 0.6919480 0.6923223 0.6926931 0.6930605 [239] 0.6934244 0.6937847 0.6941416 0.6944949 0.6948447 0.6951908 0.6955333 [246] 0.6958722 0.6962075 0.6965390 0.6968669 0.6971910 0.6975114 0.6978280 [253] 0.6981408 0.6984498 0.6987550 0.6990563 0.6993538 0.6996473 0.6999369 [260] 0.7002226 0.7005043 0.7007821 0.7010558 0.7013255 0.7015912 0.7018528 [267] 0.7021102 0.7023636 0.7026129 0.7028580 0.7030989 0.7033357 0.7035682 [274] 0.7037965 0.7040205 0.7042403 0.7044558 0.7046669 0.7048738 0.7050763 [281] 0.7052744 0.7054681 0.7056574 0.7058423 0.7060227 0.7061987 0.7063702 [288] 0.7065372 0.7066996 0.7068575 0.7070109 0.7071597 0.7073039 0.7074435 [295] 0.7075785 0.7077088 0.7078344 0.7079554 0.7080717 0.7081833 0.7082901 [302] 0.7083922 0.7084896 0.7085821 0.7086699 0.7087529 0.7088310 0.7089043 [309] 0.7089727 0.7090363 0.7090950 0.7091488 0.7091977 0.7092417 0.7092807 [316] 0.7093148 0.7093439 0.7093681 0.7093872 0.7094014 0.7094105 0.7094146 [323] 0.7094137 0.7094077 0.7093967 0.7093806 0.7093594 0.7093331 0.7093017 [330] 0.7092651 0.7092235 0.7091767 0.7091248 0.7090677 0.7090055 0.7089381 [337] 0.7088655 0.7087877 0.7087047 0.7086165 0.7085231 0.7084245 0.7083207 [344] 0.7082116 0.7080973 0.7079778 0.7078530 0.7077229 0.7075876 0.7074470 [351] 0.7073012 0.7071500 0.7069936 0.7068320 0.7066650 0.7064927 0.7063152 [358] 0.7061323 0.7059442 0.7057508 0.7055520 0.7053480 0.7051387 0.7049240 [365] 0.7047041 0.7044788 0.7042483 0.7040124 0.7037713 0.7035248 0.7032731 [372] 0.7030160 0.7027537 0.7024860 0.7022131 0.7019349 0.7016514 0.7013626 [379] 0.7010685 0.7007692 0.7004646 0.7001547 0.6998396 0.6995192 0.6991936 [386] 0.6988628 0.6985267 0.6981854 0.6978388 0.6974871 0.6971301 0.6967680 [393] 0.6964007 0.6960281 0.6956505 0.6952676 0.6948796 0.6944865 0.6940882 [400] 0.6936848 0.6932763 > mx [1] 0.7094146 > mxli [1] 1.21 > 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/wessaorg/rcomp/tmp/1adgt1321542896.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2gdgi1321542896.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3f1c81321542896.ps",horizontal=F,onefile=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 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/4e9221321542896.tab") > > try(system("convert tmp/1adgt1321542896.ps tmp/1adgt1321542896.png",intern=TRUE)) character(0) > try(system("convert tmp/2gdgi1321542896.ps tmp/2gdgi1321542896.png",intern=TRUE)) character(0) > try(system("convert tmp/3f1c81321542896.ps tmp/3f1c81321542896.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.255 0.179 1.460