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.000212923517308801 + ,-0.000366405397205687 + ,-0.00557947052088672 + ,5.39349385913969e-05 + ,0.00152371966762789 + ,0.00175850683193595 + ,3.77672323594971e-05 + ,0.00278145948179106 + ,-0.00301556173149557 + ,0.00109580938363319 + ,-0.000724132225897542 + ,-0.000973482661448984 + ,-0.00171178063930695 + ,0.000125644167827909 + ,-0.00610402063012358 + ,-0.000405313640416258 + ,7.26579064577097e-05 + ,-0.00169069400551828 + ,0.00113463103249704 + ,-0.00502021989243185 + ,0.00290671911265335 + ,0.0026126296789493 + ,-0.00247365301824070 + ,-0.0041985012668281 + ,0.00220019470292948 + ,0.00403711355294332 + ,0.00334851349840243 + ,-0.00063160538606634 + ,-0.000495111342862893 + ,0.0019798998189597 + ,0.00156938199908718 + ,-0.00196786552076488 + ,0.00453550691249892 + ,0.000238202506318942 + ,-0.00340957657161807 + ,-0.00110204572510288 + ,-0.00571315621593679 + ,-0.000483401501108857 + ,-0.00345934859091221 + ,0.00157507010324105 + ,-0.00404358222050192 + ,-0.00189791760312953 + ,0.00308855952516814 + ,-0.000418917431675471 + ,0.00131829063658993 + ,-0.000343162709266577 + ,0.00313721513408243 + ,0.00281275625682837) > 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] 48 > (np <- floor(n / par1)) [1] 4 > 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] 4 4 4 4 4 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -2.129235e-04 -0.0017117806 0.0022001947 -0.0057131562 NA [2,] -3.664054e-04 0.0001256442 0.0040371136 -0.0004834015 NA [3,] -5.579471e-03 -0.0061040206 0.0033485135 -0.0034593486 NA [4,] 5.393494e-05 -0.0004053136 -0.0006316054 0.0015750701 NA [5,] 1.523720e-03 0.0000726579 -0.0004951113 -0.0040435822 NA [6,] 1.758507e-03 -0.0016906940 0.0019798998 -0.0018979176 NA [7,] 3.776723e-05 0.0011346310 0.0015693820 0.0030885595 NA [8,] 2.781459e-03 -0.0050202199 -0.0019678655 -0.0004189174 NA [9,] -3.015562e-03 0.0029067191 0.0045355069 0.0013182906 NA [10,] 1.095809e-03 0.0026126297 0.0002382025 -0.0003431627 NA [11,] -7.241322e-04 -0.0024736530 -0.0034095766 0.0031372151 NA [12,] -9.734827e-04 -0.0041985013 -0.0011020457 0.0028127563 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] -0.0001534819 0.0018374248 0.0018369189 0.0052297547 NA [2,] -0.0052130651 -0.0062296648 -0.0006886001 -0.0029759471 NA [3,] 0.0056334055 0.0056987070 -0.0039801189 0.0050344187 NA [4,] 0.0014697847 0.0004779715 0.0001364940 -0.0056186523 NA [5,] 0.0002347872 -0.0017633519 0.0024750112 0.0021456646 NA [6,] -0.0017207396 0.0028253250 -0.0004105178 0.0049864771 NA [7,] 0.0027436922 -0.0061548509 -0.0035372475 -0.0035074770 NA [8,] -0.0057970212 0.0079269390 0.0065033724 0.0017372081 NA [9,] 0.0041113711 -0.0002940894 -0.0042973044 -0.0016614533 NA [10,] -0.0018199416 -0.0050862827 -0.0036477791 0.0034803778 NA [11,] -0.0002493504 -0.0017248482 0.0023075308 -0.0003244589 NA [12,] -0.0007382980 0.0063986960 -0.0046111105 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/16ijk1260289704.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/2v7u11260289704.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/3uf531260289704.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/4jlvs1260289704.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.0057131562 -0.0004834015 -6.104021e-03 -0.0006316054 -0.0040435822 [2,] -0.0037124684 -0.0004249034 -5.841746e-03 -0.0005184595 -0.0022693468 [3,] -0.0009623521 -0.0001203806 -4.519410e-03 -0.0001756894 -0.0002112267 [4,] 0.0009936356 0.0020813789 -5.541755e-05 0.0008145025 0.0007981888 [5,] 0.0022001947 0.0040371136 3.348513e-03 0.0015750701 0.0015237197 [,6] [,7] [,8] [,9] [,10] [1,] -1.897918e-03 3.776723e-05 -0.005020220 -0.0030155617 -0.0003431627 [2,] -1.794306e-03 5.861991e-04 -0.003494043 -0.0008486355 -0.0000524801 [3,] 3.390641e-05 1.352007e-03 -0.001193391 0.0021125049 0.0006670059 [4,] 1.869203e-03 2.328971e-03 0.001181271 0.0037211130 0.0018542195 [5,] 1.979900e-03 3.088560e-03 0.002781459 0.0045355069 0.0026126297 [,11] [,12] [1,] -0.003409577 -0.0041985013 [2,] -0.002941615 -0.0026502735 [3,] -0.001598893 -0.0010377642 [4,] 0.001206541 0.0009196368 [5,] 0.003137215 0.0028127563 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.004680174 -0.002100344 -9.090609e-03 -0.0012287294 -0.002634580 [2,] 0.002755470 0.001859582 5.178959e-05 0.0008773507 0.002212126 [,6] [,7] [,8] [,9] [,10] [1,] -0.002860266 -2.478307e-05 -0.004886889 -0.001497596 -0.0008392868 [2,] 0.002928079 2.728796e-03 0.002500106 0.005722606 0.0021732987 [,11] [,12] [1,] -0.004875936 -0.003857993 [2,] 0.001678151 0.001782465 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.00571315621593679, -0.00371246842762187, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5xvfl1260289704.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.0001534819 -0.006229665 -0.0039801189 -0.0056186523 -0.0017633519 [2,] 0.0008417185 -0.005721365 0.0005271499 -0.0027410791 -0.0007642824 [3,] 0.0018371718 -0.004094506 0.0053339121 0.0003072328 0.0011902259 [4,] 0.0035335898 -0.001832274 0.0056660562 0.0009738781 0.0023103379 [5,] 0.0052297547 -0.000688600 0.0056987070 0.0014697847 0.0024750112 [,6] [,7] [,8] [,9] [,10] [1,] -0.001720740 -0.0061548509 -0.005797021 -0.0042973044 -0.0050862827 [2,] -0.001065629 -0.0048460492 -0.002029907 -0.0029793789 -0.0043670309 [3,] 0.001207404 -0.0035223622 0.004120290 -0.0009777714 -0.0027338603 [4,] 0.003905901 -0.0003818924 0.007215156 0.0019086408 0.0008302181 [5,] 0.004986477 0.0027436922 0.007926939 0.0041113711 0.0034803778 [,11] [,12] [1,] -0.0017248482 -0.004611110 [2,] -0.0010246536 -0.002674704 [3,] -0.0002869047 -0.000738298 [4,] 0.0010290902 0.002830199 [5,] 0.0023075308 0.006398696 $n [1] 4 4 4 4 4 4 4 4 4 4 4 3 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.0002894065 -0.007166888 0.001274176 -0.002627583 -0.001238724 [2,] 0.0039637501 -0.001022124 0.009393648 0.003242049 0.003619176 [,6] [,7] [,8] [,9] [,10] [1,] -0.002720105 -7.049046e-03 -0.003183309 -0.004839307 -0.006839687 [2,] 0.005134912 4.321688e-06 0.011423889 0.002883764 0.001371966 [,11] [,12] [1,] -0.001909362 -0.005759945 [2,] 0.001335553 0.004283349 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.000153481879896886, 0.000841718485058477, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6dok61260289704.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,] -3.015562e-03 -0.0061040206 -0.0034095766 -0.0057131562 NA [2,] -8.488074e-04 -0.0033360771 -0.0008668256 -0.0026786331 NA [3,] -8.757814e-05 -0.0010480038 0.0009037923 -0.0003810401 NA [4,] 1.309765e-03 0.0006301376 0.0027743541 0.0021939132 NA [5,] 2.781459e-03 0.0029067191 0.0045355069 0.0031372151 NA $n [1] 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.001072117 -0.0028570210 -0.0007569742 -0.002603441 NA [2,] 0.000896961 0.0007610133 0.0025645587 0.001841361 NA $out [1] -0.00557947 $group [1] 1 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.00301556173149557, -0.000848807443673263, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7nrmb1260289704.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.0029485816 -0.0015988926 -0.0037339956 [2,] -0.0010119613 -0.0011155778 -0.0011678865 [3,] -0.0003490651 -0.0001934580 -0.0002436184 [4,] 0.0008645537 0.0003504562 0.0005689332 [5,] 0.0014575849 0.0021125049 0.0017743718 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.001204956 -0.0008621260 -0.0010357935 [2,] 0.000506826 0.0004752099 0.0005485567 $out [1] -0.00451941 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.00294858156088002, -0.00101196125559438, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/16ijk1260289704.ps tmp/16ijk1260289704.png") > system("convert tmp/2v7u11260289704.ps tmp/2v7u11260289704.png") > system("convert tmp/3uf531260289704.ps tmp/3uf531260289704.png") > system("convert tmp/4jlvs1260289704.ps tmp/4jlvs1260289704.png") > system("convert tmp/5xvfl1260289704.ps tmp/5xvfl1260289704.png") > system("convert tmp/6dok61260289704.ps tmp/6dok61260289704.png") > system("convert tmp/7nrmb1260289704.ps tmp/7nrmb1260289704.png") > > > proc.time() user system elapsed 1.220 0.991 4.763