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(-5.57188623643403e-05 + ,0.000673127001655781 + ,0.00253297316798237 + ,-0.00312808800265659 + ,6.06736685254191e-05 + ,0.00111459533092098 + ,-0.000554109520648384 + ,-0.00105302564390300 + ,0.00104708036627266 + ,-0.000154081652574897 + ,0.00384608131631652 + ,-0.00119915735301702 + ,-0.000654000489824518 + ,0.00147251443605623 + ,-0.000700006059835063 + ,-0.00061908177055739 + ,0.000837030387252926 + ,0.00112904739358168 + ,-0.000635971834993874 + ,0.000438787250423991 + ,0.00336388391216243 + ,0.000921403408594523 + ,-0.000666757721707773 + ,-0.00102031957840900 + ,-0.000796646941257639 + ,-0.00332057539851693 + ,-0.000182532249361493 + ,0.000146444976543607 + ,0.000563712163999504 + ,-0.00101555676532265 + ,-0.000213138748455123 + ,-0.00116972320289465 + ,0.000718163786614911 + ,0.00147451936379633 + ,-0.00186871515930147 + ,-0.00152438264067532 + ,-0.00079919301702642 + ,0.00214361864863578 + ,0.00246081566765619 + ,-0.00105955030987362 + ,-0.000841465941827341 + ,-0.000224431304780554 + ,-0.000524348696489676 + ,0.00168130096323857 + ,-0.00179203973258881 + ,8.57375031507516e-05 + ,-0.000771092454595768 + ,-0.000903583373710958 + ,0.000148807028873512 + ,-0.001557405462672 + ,0.000344209067865560 + ,-0.00101577409715098 + ,0.000968280070164763 + ,-0.000691907060623672 + ,-0.000647161091857317 + ,0.000306078767940221 + ,-0.00143204637452955 + ,0.00171652552326860 + ,0.00100568755956631 + ,-0.00162283353151299 + ,-5.03611946386278e-05 + ,-0.00139199227766638 + ,0.0008205486056255 + ,0.00106165242946242 + ,-0.00149168951466355 + ,0.000113318551053342 + ,-0.000798705188071938 + ,0.000725657525939686 + ,0.000996463288562474 + ,0.000551870084508547 + ,-0.000750413138897793 + ,-0.00107155858159102 + ,0.000778446826870656 + ,-0.000420986168065270 + ,-0.00123476222127950 + ,5.58572848751406e-05 + ,0.000674648736605934 + ,-0.00070571283485043 + ,0.000562433086011177 + ,0.00155795278362076 + ,0.000246475697694884 + ,-0.000138379316298359 + ,-0.000372926968554567 + ,0.00161826258007914 + ,-0.00118416871764147 + ,-1.79064625009521e-05 + ,-0.000530448371326978 + ,0.000513060612174883 + ,-0.000488611530751738 + ,-0.000344435665802584 + ,0.000842398928486991 + ,-0.00078884664263527 + ,-0.00030481479842588 + ,-0.00072415539117241 + ,0.000502998571283987 + ,-0.000437715789622652 + ,-0.000201547985612967 + ,-0.000674925273024028 + ,-0.000139783804754848 + ,-0.000849121145158096 + ,6.60224785191873e-05 + ,0.000113764899993581 + ,0.00159428581141906 + ,0.00113767761833825 + ,0.00355745661277584 + ,0.00266114068025193 + ,0.00242685751739785 + ,-0.00123747069021010 + ,0.000630850445637065 + ,-0.000959554038128756 + ,-0.000666732374209195 + ,-0.00142693285098644 + ,-0.00158022058363810 + ,0.00143352595627639 + ,-0.00138297408790195 + ,0.00105401647418059) > 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] 116 > (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] 10 10 10 10 10 10 10 10 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [1,] -5.571886e-05 -0.0006540005 -0.0007966469 -0.0007991930 0.0001488070 [2,] 6.731270e-04 0.0014725144 -0.0033205754 0.0021436186 -0.0015574055 [3,] 2.532973e-03 -0.0007000061 -0.0001825322 0.0024608157 0.0003442091 [4,] -3.128088e-03 -0.0006190818 0.0001464450 -0.0010595503 -0.0010157741 [5,] 6.067367e-05 0.0008370304 0.0005637122 -0.0008414659 0.0009682801 [6,] 1.114595e-03 0.0011290474 -0.0010155568 -0.0002244313 -0.0006919071 [7,] -5.541095e-04 -0.0006359718 -0.0002131387 -0.0005243487 -0.0006471611 [8,] -1.053026e-03 0.0004387873 -0.0011697232 0.0016813010 0.0003060788 [9,] 1.047080e-03 0.0033638839 0.0007181638 -0.0017920397 -0.0014320464 [10,] -1.540817e-04 0.0009214034 0.0014745194 0.0000857375 0.0017165255 [11,] 3.846081e-03 -0.0006667577 -0.0018687152 -0.0007710925 0.0010056876 [12,] -1.199157e-03 -0.0010203196 -0.0015243826 -0.0009035834 -0.0016228335 [,6] [,7] [,8] [,9] [,10] [1,] -5.036119e-05 7.784468e-04 -1.184169e-03 -2.015480e-04 0.0006308504 [2,] -1.391992e-03 -4.209862e-04 -1.790646e-05 -6.749253e-04 -0.0009595540 [3,] 8.205486e-04 -1.234762e-03 -5.304484e-04 -1.397838e-04 -0.0006667324 [4,] 1.061652e-03 5.585728e-05 5.130606e-04 -8.491211e-04 -0.0014269329 [5,] -1.491690e-03 6.746487e-04 -4.886115e-04 6.602248e-05 -0.0015802206 [6,] 1.133186e-04 -7.057128e-04 -3.444357e-04 1.137649e-04 0.0014335260 [7,] -7.987052e-04 5.624331e-04 8.423989e-04 1.594286e-03 -0.0013829741 [8,] 7.256575e-04 1.557953e-03 -7.888466e-04 1.137678e-03 0.0010540165 [9,] 9.964633e-04 2.464757e-04 -3.048148e-04 3.557457e-03 NA [10,] 5.518701e-04 -1.383793e-04 -7.241554e-04 2.661141e-03 NA [11,] -7.504131e-04 -3.729270e-04 5.029986e-04 2.426858e-03 NA [12,] -1.071559e-03 1.618263e-03 -4.377158e-04 -1.237471e-03 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.0007288459 2.126515e-03 -0.0025239285 0.0029428117 -1.706212e-03 [2,] 0.0018598462 -2.172520e-03 0.0031380431 0.0003171970 1.901615e-03 [3,] -0.0056610612 8.092429e-05 0.0003289772 -0.0035203660 -1.359983e-03 [4,] 0.0031887617 1.456112e-03 0.0004172672 0.0002180844 1.984054e-03 [5,] 0.0010539217 2.920170e-04 -0.0015792689 0.0006170346 -1.660187e-03 [6,] -0.0016687049 -1.765019e-03 0.0008024180 -0.0002999174 4.474597e-05 [7,] -0.0004989161 1.074759e-03 -0.0009565845 0.0022056497 9.532399e-04 [8,] 0.0021001060 2.925097e-03 0.0018878870 -0.0034733407 -1.738125e-03 [9,] -0.0012011620 -2.442481e-03 0.0007563556 0.0018777772 3.148572e-03 [10,] 0.0040001630 -1.588161e-03 -0.0033432345 -0.0008568300 -7.108380e-04 [11,] -0.0050452387 -3.535619e-04 0.0003443325 -0.0001324909 -2.628521e-03 [12,] 0.0005451569 2.236726e-04 0.0007251896 0.0010523904 1.572472e-03 [,6] [,7] [,8] [,9] [,10] [1,] -0.0013416311 -0.0011994330 0.0011662623 -4.733773e-04 -0.0015904045 [2,] 0.0022125409 -0.0008137761 -0.0005125419 5.351415e-04 0.0002928217 [3,] 0.0002411038 0.0012906195 0.0010435090 -7.093373e-04 -0.0007602005 [4,] -0.0025533419 0.0006187915 -0.0010016721 9.151436e-04 -0.0001532877 [5,] 0.0016050081 -0.0013803616 0.0001441759 4.774242e-05 0.0030137465 [6,] -0.0009120237 0.0012681459 0.0011868346 1.480521e-03 -0.0028165000 [7,] 0.0015243627 0.0009955197 -0.0016312456 -4.566082e-04 0.0024369906 [8,] 0.0002708058 -0.0013114771 0.0004840318 2.419779e-03 NA [9,] -0.0004445932 -0.0003848550 -0.0004193406 -8.963159e-04 NA [10,] -0.0013022832 -0.0002345477 0.0012271540 -2.342832e-04 NA [11,] -0.0003211454 0.0019911895 -0.0009407144 -3.664328e-03 NA [12,] 0.0018500054 -0.0028024313 0.0002361678 1.868321e-03 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/1jh8j1259954830.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/2ibfx1259954830.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/3o2dr1259954830.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/4v2u61259954830.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.0011841687 -0.0033205754 -0.0012347622 -0.0014269329 -1.580221e-03 [2,] -0.0007966469 -0.0013919923 -0.0006667324 -0.0010595503 -8.414659e-04 [3,] -0.0001286334 -0.0005479557 -0.0001611580 -0.0007341015 6.334807e-05 [4,] 0.0001488070 0.0006731270 0.0008205486 0.0001464450 6.746487e-04 [5,] 0.0007784468 0.0021436186 0.0025329732 0.0010616524 9.682801e-04 [,6] [,7] [,8] [,9] [,10] [1,] -1.015557e-03 -0.0013829741 -0.0011697232 -0.0017920397 -0.0007241554 [2,] -6.919071e-04 -0.0006471611 -0.0007888466 -0.0003048148 -0.0001383793 [3,] -5.555638e-05 -0.0005392291 0.0005822224 0.0007181638 0.0005518701 [4,] 1.114595e-03 0.0005624331 0.0011376776 0.0010470804 0.0014745194 [5,] 1.433526e-03 0.0015942858 0.0016813010 0.0010470804 0.0026611407 [,11] [,12] [1,] -0.0018687152 -0.0016228335 [2,] -0.0007504131 -0.0012374707 [3,] -0.0003729270 -0.0010715586 [4,] 0.0010056876 -0.0009035834 [5,] 0.0024268575 -0.0004377158 $n [1] 10 10 10 10 10 10 10 10 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.0006010199 -0.0015797716 -0.0009042629 -0.0013366648 -0.0006941633 [2,] 0.0003437531 0.0004838602 0.0005819468 -0.0001315381 0.0008208594 [,6] [,7] [,8] [,9] [,10] [1,] -0.0009581570 -1.143591e-03 -0.0003803459 6.165667e-06 -0.0002975899 [2,] 0.0008470442 6.513237e-05 0.0015447907 1.430162e-03 0.0014013301 [,11] [,12] [1,] -0.0012978067 -0.0012474059 [2,] 0.0005519527 -0.0008957113 $out [1] -0.003128088 0.003363884 0.003557457 0.003846081 0.001618263 $group [1] 4 9 9 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.00118416871764147, -0.000796646941257639, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5e3bc1259954830.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.0025239285 -0.0021725205 -0.0035203660 -0.0025533419 -0.0016601871 [2,] -0.0015904045 -0.0005125419 -0.0013599832 -0.0001532877 -0.0013803616 [3,] -0.0008364051 0.0004261692 -0.0003142065 0.0005180293 0.0002180964 [4,] 0.0011662623 0.0019016145 0.0003289772 0.0014561122 0.0010539217 [5,] 0.0029428117 0.0031380431 0.0012906195 0.0031887617 0.0030137465 [,6] [,7] [,8] [,9] [,10] [1,] -0.0028165000 -0.0016312456 -0.0034733407 -0.0024424805 -0.0015881611 [2,] -0.0016687049 -0.0004989161 -0.0013114771 -0.0008963159 -0.0013022832 [3,] -0.0001275857 0.0009743798 0.0004840318 -0.0004193406 -0.0007108380 [4,] 0.0011868346 0.0015243627 0.0021001060 0.0007563556 -0.0002342832 [5,] 0.0014805209 0.0024369906 0.0029250967 0.0031485719 0.0012271540 [,11] [,12] [1,] -0.0050452387 0.0002236726 [2,] -0.0026285211 0.0002361678 [3,] -0.0003535619 0.0007251896 [4,] -0.0001324909 0.0015724723 [5,] 0.0019911895 0.0018683211 $n [1] 10 10 10 10 10 10 10 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.0022137458 -0.0007800396 -0.0011580785 -0.0002860910 -0.0009981685 [2,] 0.0005409355 0.0016323781 0.0005296654 0.0013221497 0.0014343614 [,6] [,7] [,8] [,9] [,10] [1,] -0.001554327 -0.000036531 -0.001312735 -0.0012897476 -0.0012733180 [2,] 0.001299156 0.001985291 0.002280799 0.0004510664 -0.0001483579 [,11] [,12] [1,] -0.001668138 2.140257e-05 [2,] 0.000961014 1.428977e-03 $out [1] -0.005661061 0.004000163 -0.003343235 -0.002802431 $group [1] 3 10 10 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.00252392845725929, -0.00159040448376582, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/63i2s1259954830.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.128088e-03 -1.020320e-03 -0.0033205754 -0.0017920397 -0.0016228335 [2,] -8.035676e-04 -6.603791e-04 -0.0013470529 -0.0008725247 -0.0012239102 [3,] 2.477403e-06 -9.014726e-05 -0.0005048928 -0.0006477206 -0.0002491770 [4,] 1.080838e-03 1.025225e-03 0.0003550786 0.0008835192 0.0006562446 [5,] 3.846081e-03 3.363884e-03 0.0014745194 0.0024608157 0.0017165255 [,6] [,7] [,8] [,9] [,10] [1,] -1.491690e-03 -0.0012347622 -0.0011841687 -1.237471e-03 -0.0015802206 [2,] -9.351319e-04 -0.0003969566 -0.0006273019 -4.382366e-04 -0.0014049535 [3,] 3.147868e-05 0.0001511665 -0.0003910757 8.989369e-05 -0.0008131432 [4,] 7.731031e-04 0.0007265478 0.0002425461 2.010572e-03 0.0008424335 [5,] 1.061652e-03 0.0016182626 0.0008423989 3.557457e-03 0.0014335260 $n [1] 12 12 12 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.0008570126 -0.0008589628 -0.0012812464 -0.0014486640 -0.0011067283 [2,] 0.0008619674 0.0006786683 0.0002714608 0.0001532228 0.0006083742 [,6] [,7] [,8] [,9] [,10] [1,] -0.0007476587 -0.0003612714 -7.878192e-04 -0.001027024 -0.0020685658 [2,] 0.0008106161 0.0006636044 5.667755e-06 0.001206812 0.0004422794 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(-0.00312808800265659, -0.000803567582275692, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/77lx11259954830.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,] -8.220843e-04 -0.0010715586 -1.070527e-03 [2,] -3.118809e-04 -0.0005435924 -3.613711e-04 [3,] -1.547058e-05 -0.0001448957 -9.383849e-06 [4,] 3.807004e-04 0.0003076091 2.152489e-04 [5,] 7.111803e-04 0.0007181638 6.680700e-04 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.0003313616 -0.0005331344 -0.0002723841 [2,] 0.0003004204 0.0002433430 0.0002536164 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.000822084328741102, -0.000311880895060564, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > system("convert tmp/1jh8j1259954830.ps tmp/1jh8j1259954830.png") > system("convert tmp/2ibfx1259954830.ps tmp/2ibfx1259954830.png") > system("convert tmp/3o2dr1259954830.ps tmp/3o2dr1259954830.png") > system("convert tmp/4v2u61259954830.ps tmp/4v2u61259954830.png") > system("convert tmp/5e3bc1259954830.ps tmp/5e3bc1259954830.png") > system("convert tmp/63i2s1259954830.ps tmp/63i2s1259954830.png") > system("convert tmp/77lx11259954830.ps tmp/77lx11259954830.png") > > > proc.time() user system elapsed 1.202 0.928 1.460