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(0.103629933239820 + ,0.0087703074390444 + ,0.0155637426660695 + ,0.0966084256018938 + ,0.0683025606568869 + ,-0.00653314961359896 + ,-0.0094367257568807 + ,0.00935003832683242 + ,0.118449494805129 + ,0.127456615819888 + ,0.00876382231245558 + ,-0.0129824429665246 + ,-0.00415418317902038 + ,0.0266421877779411 + ,0.194032448777481 + ,-0.00711070355878007 + ,0.0205688036309368 + ,-0.00442856389169322 + ,-0.00315893309607418 + ,0.0064652048591256 + ,0.131873973578195 + ,0.241305044412940 + ,-0.0310262850916416 + ,0.0245251063758616 + ,-0.0156222322911219 + ,-0.0206043854226294 + ,0.0326967759710612 + ,0.125988059619388 + ,-0.0435502067574021 + ,0.0123515713707718 + ,-0.00434737918082817 + ,-0.00353479514087951 + ,0.147562133549002 + ,0.0470693504162938 + ,0.105815471844778 + ,-0.0341342674195459 + ,-0.00414417660174138 + ,-0.0064652048591256 + ,0.0468257381816528 + ,-0.225104141948705 + ,3.12331458758308e-05 + ,0.0166685620235825 + ,0.00237689216950798 + ,0 + ,0.0357693020416008 + ,0.0102083822537082 + ,-0.0675670602819736 + ,0.0495991980182566 + ,0.0101167955535715 + ,0.0668105237905223 + ,-0.223001756937975 + ,0.0896858022691731 + ,-0.106531513400583 + ,-0.00678504840645644 + ,0.0134485427007007 + ,0.0300000000000011 + ,-0.0989691535751547 + ,-0.115464366100554 + ,0.0988679532745067 + ,-0.0709498589393718) > par1 = '12' > #'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!) > par1 <- as.numeric(par1) > (n <- length(x)) [1] 60 > (np <- floor(n / par1)) [1] 5 > arr <- array(NA,dim=c(par1,np+1)) > darr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > dx <- diff(x) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + darr[j,ari[j]] <- dx[i] + if (j == par1) j = 0 + } > ari [1] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.103629933 -0.004154183 -0.015622232 -4.144177e-03 0.010116796 NA [2,] 0.008770307 0.026642188 -0.020604385 -6.465205e-03 0.066810524 NA [3,] 0.015563743 0.194032449 0.032696776 4.682574e-02 -0.223001757 NA [4,] 0.096608426 -0.007110704 0.125988060 -2.251041e-01 0.089685802 NA [5,] 0.068302561 0.020568804 -0.043550207 3.123315e-05 -0.106531513 NA [6,] -0.006533150 -0.004428564 0.012351571 1.666856e-02 -0.006785048 NA [7,] -0.009436726 -0.003158933 -0.004347379 2.376892e-03 0.013448543 NA [8,] 0.009350038 0.006465205 -0.003534795 0.000000e+00 0.030000000 NA [9,] 0.118449495 0.131873974 0.147562134 3.576930e-02 -0.098969154 NA [10,] 0.127456616 0.241305044 0.047069350 1.020838e-02 -0.115464366 NA [11,] 0.008763822 -0.031026285 0.105815472 -6.756706e-02 0.098867953 NA [12,] -0.012982443 0.024525106 -0.034134267 4.959920e-02 -0.070949859 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.094859626 0.030796371 -0.004982153 -0.002321028 0.05669373 NA [2,] 0.006793435 0.167390261 0.053301161 0.053290943 -0.28981228 NA [3,] 0.081044683 -0.201143152 0.093291284 -0.271929880 0.31268756 NA [4,] -0.028305865 0.027679507 -0.169538266 0.225135375 -0.19621732 NA [5,] -0.074835710 -0.024997368 0.055901778 0.016637329 0.09974646 NA [6,] -0.002903576 0.001269631 -0.016698951 -0.014291670 0.02023359 NA [7,] 0.018786764 0.009624138 0.000812584 -0.002376892 0.01655146 NA [8,] 0.109099456 0.125408769 0.151096929 0.035769302 -0.12896915 NA [9,] 0.009007121 0.109431071 -0.100492783 -0.025560920 -0.01649521 NA [10,] -0.118692794 -0.272331330 0.058746121 -0.077775443 0.21433232 NA [11,] -0.021746265 0.055551391 -0.139949739 0.117166258 -0.16981781 NA [12,] 0.008828260 -0.040147339 0.029990091 -0.039482402 NA NA > arr.mean <- array(NA,dim=par1) > arr.median <- array(NA,dim=par1) > arr.midrange <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.mean[j] <- mean(arr[j,],na.rm=TRUE) + arr.median[j] <- median(arr[j,],na.rm=TRUE) + arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2 + } > overall.mean <- mean(x) > overall.median <- median(x) > overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2 > postscript(file="/var/www/html/rcomp/tmp/1977l1259702878.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.mean,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2p5ow1259702878.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.median,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3d94c1259702878.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.midrange,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4qjdj1259702878.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(arr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [1,] -0.015622232 -0.020604385 0.01556374 -0.007110704 -1.065315e-01 [2,] -0.004154183 -0.006465205 0.01556374 -0.007110704 -4.355021e-02 [3,] -0.004144177 0.008770307 0.03269678 0.089685802 3.123315e-05 [4,] 0.010116796 0.026642188 0.04682574 0.096608426 2.056880e-02 [5,] 0.010116796 0.066810524 0.04682574 0.125988060 6.830256e-02 [,6] [,7] [,8] [,9] [,10] [1,] -0.006785048 -0.009436726 -0.003534795 -0.09896915 -0.11546437 [2,] -0.006533150 -0.004347379 0.000000000 0.03576930 0.01020838 [3,] -0.004428564 -0.003158933 0.006465205 0.11844949 0.04706935 [4,] 0.012351571 0.002376892 0.009350038 0.13187397 0.12745662 [5,] 0.016668562 0.002376892 0.009350038 0.14756213 0.24130504 [,11] [,12] [1,] -0.067567060 -0.07094986 [2,] -0.031026285 -0.03413427 [3,] 0.008763822 -0.01298244 [4,] 0.098867953 0.02452511 [5,] 0.105815472 0.04959920 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.014228014 -0.01462329 0.01060713 0.01639813 -0.04527510 -0.017772460 [2,] 0.005939661 0.03216391 0.05478642 0.16297348 0.04533756 0.008915332 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.007910286 -0.0001415087 0.05054218 -0.03577796 -0.08301912 -0.05443101 [2,] 0.001592420 0.0130719184 0.18635681 0.12991666 0.10054676 0.02846612 $out [1] 0.10362993 0.19403245 -0.22300176 -0.22510414 0.01344854 0.03000000 $group [1] 1 3 3 4 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0156222322911219, -0.00415418317902038, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5504j1259702878.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(darr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.004982153 0.006793435 -0.27192988 -0.19621732 -0.07483571 -0.016698951 [2,] -0.004982153 0.006793435 -0.20114315 -0.16953827 -0.02499737 -0.014291670 [3,] -0.002321028 0.053290943 0.08104468 -0.02830586 0.01663733 -0.002903576 [4,] 0.030796371 0.053301161 0.09329128 0.02767951 0.05590178 0.001269631 [5,] 0.056693728 0.053301161 0.31268756 0.22513538 0.09974646 0.020233591 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.002376892 0.0357693 -0.025560920 -0.27233133 -0.16981781 -0.04014734 [2,] 0.000812584 0.0357693 -0.025560920 -0.11869279 -0.13994974 -0.03981487 [3,] 0.009624138 0.1090995 -0.016495213 -0.07777544 -0.02174627 -0.01532707 [4,] 0.016551457 0.1254088 0.009007121 0.05874612 0.05555139 0.01940918 [5,] 0.018786764 0.1510969 0.009007121 0.21433232 0.11716626 0.02999009 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.02760204 0.02042870 -0.1270019 -0.1676594 -0.04052580 -0.01389915 [2,] 0.02295999 0.08615319 0.2890913 0.1110477 0.07380046 0.00809200 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.001496910 0.04576044 -0.040920903 -0.20315333 -0.1598869 -0.06211407 [2,] 0.020745186 0.17243848 0.007930478 0.04760245 0.1163943 0.03145992 $out [1] -0.09485963 0.16739026 -0.28981228 -0.12896915 0.10943107 -0.10049278 $group [1] 1 2 2 8 9 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0049821531315075, -0.0049821531315075, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6qyr71259702878.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(arr) > names(z) <- c(1:np) > (boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.012982443 -0.031026285 -0.043550207 -0.006465205 -0.223001757 NA [2,] 0.001115336 -0.004291374 -0.018113309 -0.005304691 -0.102750333 NA [3,] 0.012456890 0.013517004 0.004408388 0.001204063 0.001665874 NA [4,] 0.100119179 0.079258081 0.076442411 0.026218932 0.048405262 NA [5,] 0.127456616 0.194032449 0.147562134 0.049599198 0.098867953 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.03269943 -0.02459046 -0.03871911 -0.01317407 -0.06727721 NA [2,] 0.05761321 0.05162447 0.04753589 0.01558220 0.07060896 NA $out [1] 0.24130504 -0.22510414 -0.06756706 $group [1] 2 4 4 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(-0.0129824429665246, 0.00111533634942831, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7xacm1259702878.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.mean,arr.median,arr.midrange)) > names(z) <- list('mean','median','midrange') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] -0.012235825 -0.012982443 -0.0114907016 [2,] 0.001015577 -0.003651555 0.0009619837 [3,] 0.014127038 0.007614514 0.0073817553 [4,] 0.020468004 0.039883063 0.0393348476 [5,] 0.022970780 0.089685802 0.0838216378 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 0.005254655 -0.01224192 -0.01012037 [2,] 0.022999420 0.02747095 0.02488388 $out [1] 0.06693715 0.06211501 0.11844949 $group [1] 1 1 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.0122358245448571, 0.00101557683190313, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1977l1259702878.ps tmp/1977l1259702878.png") > system("convert tmp/2p5ow1259702878.ps tmp/2p5ow1259702878.png") > system("convert tmp/3d94c1259702878.ps tmp/3d94c1259702878.png") > system("convert tmp/4qjdj1259702878.ps tmp/4qjdj1259702878.png") > system("convert tmp/5504j1259702878.ps tmp/5504j1259702878.png") > system("convert tmp/6qyr71259702878.ps tmp/6qyr71259702878.png") > system("convert tmp/7xacm1259702878.ps tmp/7xacm1259702878.png") > > > proc.time() user system elapsed 1.221 0.946 1.547