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.00542916130943278 + ,0.0445893441667106 + ,-0.0137888879399847 + ,0.0453152005372048 + ,0.00473315051318341 + ,-0.00629302867546508 + ,-0.070117568265751 + ,-0.0217843333183126 + ,-0.0095666262457008 + ,0.0165637840451361 + ,0.00824993070407797 + ,-0.0831327333629231 + ,0.0144161451837631 + ,0.034621459194465 + ,0.0184381942383567 + ,-0.074932085622151 + ,-0.0264502766059977 + ,-0.0112128979582415 + ,-0.10436727046515 + ,0.0715680268198404 + ,0.121060489346963 + ,-0.0324214631834054 + ,-0.00529431524629609 + ,-0.00146846269678377 + ,-0.00570151995165197 + ,-0.025987421858386 + ,-0.0685516945712228 + ,-0.0733842752561877 + ,-0.0163445778492365 + ,-0.0110134858421602 + ,-0.087660645529663 + ,0.109058909700739 + ,0.0339807510912559 + ,0.00300957055239715 + ,-0.117182246964976 + ,0.0294137488267335 + ,0.0893484543689374 + ,0.0298847434899963 + ,0.0611729025233987 + ,0.00324519869485862 + ,0.0212178770497218 + ,0.00660658211635031 + ,0.00372487320385017 + ,-0.00838101340628686 + ,0.0170033508018199 + ,0.0621313219493351 + ,0.0197401442742924 + ,-0.0417164562197347 + ,0.0228259755720078 + ,-0.042255817474836 + ,0.0309137428469852 + ,-0.00918278678351927 + ,0.0146293075798733 + ,0.0297847213520409 + ,0.00986110949240185 + ,0.0614197587407885 + ,0.0126952510044700 + ,0.00800151002534761 + ,0.0341797308392325 + ,-0.0180464496214532 + ,-0.0153479524622885 + ,0.0463170445239729 + ,0.0195281991212937 + ,0.0371269245007414 + ,-0.00744692508954433 + ,0.0462493884790173 + ,-0.0395877264208058 + ,0.0382572916842234 + ,0.0374023801364467 + ,0.039717646637445 + ,-0.0103197596326754 + ,0.00344718875977555 + ,0.00885007204463633 + ,-0.0791417702355047 + ,0.0158234735164099 + ,0.016538541536172 + ,0.0144056817515094 + ,0.00282613662406027 + ,0.0534415346690993 + ,-0.0114429989198319 + ,0.0380793061367662 + ,0.0110221503256857 + ,-0.0190504147856999 + ,0.0414994449048327 + ,0.00816603652061595 + ,0.041019836581154 + ,-0.0171632971999998 + ,-0.0239258971120320 + ,-0.0219604522663097 + ,0.00948947211086537 + ,0.039949080719774 + ,-0.0982011343440696 + ,0.0127754891139968 + ,-0.134456574642575 + ,0.0790854375582093 + ,-0.0486237536443106 + ,0.0893991774345677 + ,0.0305540121228036 + ,-0.137690614650960 + ,-0.0322141353773482 + ,0.03442439486507 + ,-0.164911125871352 + ,-0.183934800698021 + ,0.0234325246699871 + ,-0.0233244578414444 + ,0.0283600031799320 + ,-0.0273434825220074 + ,0.0285743363772815) > 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] 108 > (np <- floor(n / par1)) [1] 9 > 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] 9 9 9 9 9 9 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [1,] 0.005429161 0.014416145 -0.005701520 0.089348454 0.022825976 [2,] 0.044589344 0.034621459 -0.025987422 0.029884743 -0.042255817 [3,] -0.013788888 0.018438194 -0.068551695 0.061172903 0.030913743 [4,] 0.045315201 -0.074932086 -0.073384275 0.003245199 -0.009182787 [5,] 0.004733151 -0.026450277 -0.016344578 0.021217877 0.014629308 [6,] -0.006293029 -0.011212898 -0.011013486 0.006606582 0.029784721 [7,] -0.070117568 -0.104367270 -0.087660646 0.003724873 0.009861109 [8,] -0.021784333 0.071568027 0.109058910 -0.008381013 0.061419759 [9,] -0.009566626 0.121060489 0.033980751 0.017003351 0.012695251 [10,] 0.016563784 -0.032421463 0.003009571 0.062131322 0.008001510 [11,] 0.008249931 -0.005294315 -0.117182247 0.019740144 0.034179731 [12,] -0.083132733 -0.001468463 0.029413749 -0.041716456 -0.018046450 [,6] [,7] [,8] [,9] [,10] [1,] -0.015347952 0.008850072 0.008166037 0.08939918 NA [2,] 0.046317045 -0.079141770 0.041019837 0.03055401 NA [3,] 0.019528199 0.015823474 -0.017163297 -0.13769061 NA [4,] 0.037126925 0.016538542 -0.023925897 -0.03221414 NA [5,] -0.007446925 0.014405682 -0.021960452 0.03442439 NA [6,] 0.046249388 0.002826137 0.009489472 -0.16491113 NA [7,] -0.039587726 0.053441535 0.039949081 -0.18393480 NA [8,] 0.038257292 -0.011442999 -0.098201134 0.02343252 NA [9,] 0.037402380 0.038079306 0.012775489 -0.02332446 NA [10,] 0.039717647 0.011022150 -0.134456575 0.02836000 NA [11,] -0.010319760 -0.019050415 0.079085438 -0.02734348 NA [12,] 0.003447189 0.041499445 -0.048623754 0.02857434 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.039160183 0.020205314 -0.020285902 -0.059463711 -0.065081793 [2,] -0.058378232 -0.016183265 -0.042564273 0.031288159 0.073169560 [3,] 0.059104088 -0.093370280 -0.004832581 -0.057927704 -0.040096530 [4,] -0.040582050 0.048481809 0.057039697 0.017972678 0.023812094 [5,] -0.011026179 0.015237379 0.005331092 -0.014611295 0.015155414 [6,] -0.063824540 -0.093154373 -0.076647160 -0.002881709 -0.019923612 [7,] 0.048333235 0.175935297 0.196719555 -0.012105887 0.051558649 [8,] 0.012217707 0.049492463 -0.075078159 0.025384364 -0.048724508 [9,] 0.026130410 -0.153481953 -0.030971181 0.045127971 -0.004693741 [10,] -0.008313853 0.027127148 -0.120191818 -0.042391178 0.026178221 [11,] -0.091382664 0.003825853 0.146595996 -0.061456600 -0.052226180 [12,] 0.097548879 -0.004233057 0.059934706 0.064542432 0.002698497 [,6] [,7] [,8] [,9] [,10] [1,] 0.0616649970 -0.087991842 0.032853800 -0.05884517 NA [2,] -0.0267888454 0.094965244 -0.058183134 -0.16824463 NA [3,] 0.0175987254 0.000715068 -0.006762600 0.10547648 NA [4,] -0.0445738496 -0.002132860 0.001965445 0.06663853 NA [5,] 0.0536963136 -0.011579545 0.031449924 -0.19933552 NA [6,] -0.0858371149 0.050615398 0.030459609 -0.01902367 NA [7,] 0.0778450181 -0.064884534 -0.138150215 0.20736733 NA [8,] -0.0008549115 0.049522305 0.110976623 -0.04675698 NA [9,] 0.0023152665 -0.027057156 -0.147232064 0.05168446 NA [10,] -0.0500374063 -0.030072565 0.213542012 -0.05570349 NA [11,] 0.0137669484 0.060549860 -0.127709191 0.05591782 NA [12,] 0.0054028833 -0.033333408 0.138022931 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/1mvkg1259862265.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/2ca1k1259862265.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/3sgaw1259862265.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/425501259862265.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.015347952 -0.07914177 -0.06855169 -0.074932086 -0.026450277 [2,] 0.005429161 -0.02598742 -0.01716330 -0.032214135 -0.016344578 [3,] 0.008850072 0.03055401 0.01582347 -0.009182787 0.004733151 [4,] 0.022825976 0.04101984 0.01952820 0.016538542 0.014629308 [5,] 0.022825976 0.04631704 0.06117290 0.045315201 0.034424395 [,6] [,7] [,8] [,9] [,10] [,11] [1,] -0.011212898 -0.18393480 -0.09820113 -0.02332446 -0.032421463 -0.027343483 [2,] -0.011013486 -0.08766065 -0.01144300 0.01269525 0.003009571 -0.019050415 [3,] 0.002826137 -0.03958773 0.02343252 0.01700335 0.011022150 -0.005294315 [4,] 0.009489472 0.00986111 0.06141976 0.03740238 0.028360003 0.019740144 [5,] 0.029784721 0.05344153 0.10905891 0.03807931 0.062131322 0.034179731 [,12] [1,] -0.083132733 [2,] -0.041716456 [3,] -0.001468463 [4,] 0.028574336 [5,] 0.041499445 $n [1] 9 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.0003122501 -0.004736477 -0.003500715 -0.03485920 -0.01157976 [2,] 0.0180123942 0.065844502 0.035147662 0.01649362 0.02104606 [,6] [,7] [,8] [,9] [,10] [,11] [1,] -0.007972088 -0.09094918 -0.01494186 0.003990929 -0.002329078 -0.02572401 [2,] 0.013624361 0.01177373 0.06180691 0.030015772 0.024373378 0.01513538 [,12] [1,] -0.03848828 [2,] 0.03555135 $out [1] 0.08934845 0.08939918 -0.13769061 0.04624939 -0.16491113 0.12106049 [7] -0.13445657 -0.11718225 0.07908544 $group [1] 1 1 3 6 6 9 10 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0153479524622885, 0.00542916130943278, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5mbnz1259862265.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] [1,] -0.08799184 -0.16824463 -0.093370280 -0.044573850 -0.014611295 [2,] -0.05946371 -0.05818313 -0.040096530 -0.002132860 -0.011579545 [3,] -0.02028590 -0.02678885 -0.004832581 0.017972678 0.005331092 [4,] 0.03285380 0.03128816 0.017598725 0.048481809 0.015237379 [5,] 0.06166500 0.09496524 0.059104088 0.066638530 0.053696314 [,6] [,7] [,8] [,9] [,10] [,11] [1,] -0.093154373 -0.13815022 -0.07507816 -0.030971181 -0.12019182 -0.127709191 [2,] -0.076647160 -0.01210589 -0.04675698 -0.030971181 -0.05003741 -0.061456600 [3,] -0.019923612 0.05155865 0.01221771 -0.004693741 -0.03007257 0.003825853 [4,] -0.002881709 0.17593530 0.04949246 0.026130410 0.02617822 0.055917819 [5,] 0.050615398 0.20736733 0.11097662 0.051684461 0.02712715 0.146595996 [,12] [1,] -0.03333341 [2,] -0.00076728 [3,] 0.03266879 [4,] 0.08104566 [5,] 0.13802293 $n [1] 9 9 9 9 9 9 9 9 9 9 9 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.06890646 -0.07391039 -0.03521875 -0.00868438 -0.008792488 -0.05877342 [2,] 0.02833465 0.02033270 0.02555359 0.04462974 0.019454672 0.01892619 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.04747637 -0.03847367 -0.03476725 -0.07021280 -0.05799134 -0.01303309 [2,] 0.15059367 0.06290908 0.02537976 0.01006767 0.06564305 0.07837067 $out [1] 0.1054765 -0.1993355 -0.1534820 -0.1472321 0.2135420 $group [1] 3 5 9 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.087991842280141, -0.0594637108789411, -0.0202859019067340, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6oss81259862265.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] [1,] -0.0217843333 -0.104367270 -0.11718225 -0.041716456 -0.0422558175 [2,] -0.0177866106 -0.029435870 -0.07096798 0.003485036 -0.0005906384 [3,] -0.0007799391 -0.003381389 -0.01367903 0.018371748 0.0136622793 [4,] 0.0124068574 0.026529827 0.01621166 0.045528823 0.0303492321 [5,] 0.0453152005 0.071568027 0.10905891 0.089348454 0.0614197587 [,6] [,7] [,8] [,9] [,10] [1,] -0.039587726 -0.019050415 -0.09820113 -1.839348e-01 NA [2,] -0.008883342 -0.004308431 -0.03627483 -8.495238e-02 NA [3,] 0.028327562 0.012713916 -0.00449863 5.403341e-05 NA [4,] 0.038987469 0.027308924 0.02636228 2.956417e-02 NA [5,] 0.046317045 0.053441535 0.07908544 8.939918e-02 NA $n [1] 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.01455138 -0.02890772 -0.05344225 -0.000804706 -0.0004496033 0.006493363 [2,] 0.01299150 0.02214494 0.02608419 0.037548201 0.0277741619 0.050161760 [,7] [,8] [,9] [,10] [1,] -0.001706972 -0.03306784 -0.05217773 NA [2,] 0.027134804 0.02407058 0.05228580 NA $out [1] -0.07011757 -0.08313273 0.12106049 -0.07914177 -0.13445657 $group [1] 1 1 2 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(-0.0217843333183126, -0.0177866106291487, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7iad81259862265.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.042076824 -0.009182787 -0.0078377969 [2,] -0.010544012 -0.003381389 -0.0037143475 [3,] -0.002000390 0.006791611 0.0007636579 [4,] 0.013529359 0.016413412 0.0149061777 [5,] 0.026678437 0.030554012 0.0250488156 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.012980417 -0.002236930 -0.007729289 [2,] 0.008979636 0.015820153 0.009256604 $out [1] -0.03958773 -0.03889977 $group [1] 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.0420768236993628, -0.0105440122100338, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1mvkg1259862265.ps tmp/1mvkg1259862265.png") > system("convert tmp/2ca1k1259862265.ps tmp/2ca1k1259862265.png") > system("convert tmp/3sgaw1259862265.ps tmp/3sgaw1259862265.png") > system("convert tmp/425501259862265.ps tmp/425501259862265.png") > system("convert tmp/5mbnz1259862265.ps tmp/5mbnz1259862265.png") > system("convert tmp/6oss81259862265.ps tmp/6oss81259862265.png") > system("convert tmp/7iad81259862265.ps tmp/7iad81259862265.png") > > > proc.time() user system elapsed 1.202 0.953 1.439