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.00450750869087398 + ,-0.0149657506674553 + ,-0.0220859233394083 + ,-0.0148439243168503 + ,0.0676027781490227 + ,0.0264092299234277 + ,-0.0116301716341945 + ,0.0136476387387367 + ,-0.0197271433713688 + ,0.0303927147155863 + ,0.000995418935294206 + ,-0.00874063700973733 + ,-0.0231462237863890 + ,-0.00377238720197872 + ,0.00120873614259887 + ,-0.0176359614054821 + ,0.00250200103515695 + ,-0.0282901697726092 + ,-0.00423879675659204 + ,0.0230884796585779 + ,-0.00859471068972903 + ,0.0198693493392492 + ,-0.0093442264177522 + ,-0.0258281935447543 + ,0.0198614750918193 + ,-0.0122547342511849 + ,-0.00589049815189684 + ,0.0298004421068244 + ,0.0434323456332011 + ,-0.0223164871309812 + ,0.0154280097077713 + ,0.0156548149573578 + ,0.0365282768534479 + ,0.0241846907151588 + ,-0.0322215315404314 + ,-0.0141831668172313 + ,0.0254232968852305 + ,-0.0193736352949016 + ,-0.048536332661388 + ,-0.0766717425093701 + ,0.0281974074828454 + ,0.00917312366714304 + ,0.0616193814776786 + ,-0.0503369625485599 + ,-0.0070958483480044 + ,0.0252260185757131 + ,-0.0147649614724515 + ,0.0297338279951919 + ,0.0427693034085115 + ,0.0334986515495981 + ,-0.00431066374073932 + ,0.0483508046074947 + ,-0.0289633986573534 + ,0.0411760286974008 + ,0.00490436151035459 + ,0.00572516695084428 + ,0.0020827414405662 + ,0.0152702987490803 + ,0.0226514137977885 + ,0.0169257233733748 + ,0.00183213307281401 + ,0.0147224276771016 + ,-0.00123533731905361 + ,0.000655329694370428 + ,-0.00854895738497048 + ,-0.00289709050455367 + ,0.0236354130433993 + ,0.023075931765511 + ,0.0289109594323294 + ,-0.0376782888435601 + ,0.00381606701300569 + ,0.0202008973428546 + ,-0.00412021969633953 + ,0.0418280922509324 + ,0.067769119091173 + ,-0.0525211864056023 + ,0.0247371173448568 + ,0.00923061642650255 + ,0.0112743669564025 + ,-0.000125546542039254 + ,-0.0227570917899331 + ,0.00893231402208485 + ,-0.0388414437336418 + ,-0.0884921733798761 + ,-0.0686977274158266) > 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] 85 > (np <- floor(n / par1)) [1] 7 > 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] 8 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [1,] 0.004507509 -0.023146224 0.019861475 0.025423297 0.042769303 [2,] -0.014965751 -0.003772387 -0.012254734 -0.019373635 0.033498652 [3,] -0.022085923 0.001208736 -0.005890498 -0.048536333 -0.004310664 [4,] -0.014843924 -0.017635961 0.029800442 -0.076671743 0.048350805 [5,] 0.067602778 0.002502001 0.043432346 0.028197407 -0.028963399 [6,] 0.026409230 -0.028290170 -0.022316487 0.009173124 0.041176029 [7,] -0.011630172 -0.004238797 0.015428010 0.061619381 0.004904362 [8,] 0.013647639 0.023088480 0.015654815 -0.050336963 0.005725167 [9,] -0.019727143 -0.008594711 0.036528277 -0.007095848 0.002082741 [10,] 0.030392715 0.019869349 0.024184691 0.025226019 0.015270299 [11,] 0.000995419 -0.009344226 -0.032221532 -0.014764961 0.022651414 [12,] -0.008740637 -0.025828194 -0.014183167 0.029733828 0.016925723 [,6] [,7] [,8] [1,] 0.0018321331 -0.0041202197 -0.06869773 [2,] 0.0147224277 0.0418280923 NA [3,] -0.0012353373 0.0677691191 NA [4,] 0.0006553297 -0.0525211864 NA [5,] -0.0085489574 0.0247371173 NA [6,] -0.0028970905 0.0092306164 NA [7,] 0.0236354130 0.0112743670 NA [8,] 0.0230759318 -0.0001255465 NA [9,] 0.0289109594 -0.0227570918 NA [10,] -0.0376782888 0.0089323140 NA [11,] 0.0038160670 -0.0388414437 NA [12,] 0.0202008973 -0.0884921734 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] -0.019473259 0.019373837 -0.0321162093 -0.04479693 -0.0092706519 [2,] -0.007120173 0.004981123 0.0063642361 -0.02916270 -0.0378093153 [3,] 0.007241999 -0.018844698 0.0356909403 -0.02813541 0.0526614683 [4,] 0.082446702 0.020137962 0.0136319035 0.10486915 -0.0773142033 [5,] -0.041193548 -0.030792171 -0.0657488328 -0.01902428 0.0701394274 [6,] -0.038039402 0.024051373 0.0377444968 0.05244626 -0.0362716672 [7,] 0.025277810 0.027327276 0.0002268052 -0.11195634 0.0008208054 [8,] -0.033374782 -0.031683190 0.0208734619 0.04324111 -0.0036424255 [9,] 0.050119858 0.028464060 -0.0123435861 0.03232187 0.0131875573 [10,] -0.029397296 -0.029213576 -0.0564062223 -0.03999098 0.0073811150 [11,] -0.009736056 -0.016483967 0.0180383647 0.04449879 -0.0057256904 [12,] -0.014405587 0.045689669 0.0396064637 0.01303548 -0.0150935903 [,6] [,7] [,8] [1,] 0.0128902946 0.045948312 NA [2,] -0.0159577650 0.025941027 NA [3,] 0.0018906670 -0.120290305 NA [4,] -0.0092042871 0.077258304 NA [5,] 0.0056518669 -0.015506501 NA [6,] 0.0265325035 0.002043751 NA [7,] -0.0005594813 -0.011399913 NA [8,] 0.0058350277 -0.022631545 NA [9,] -0.0665892483 0.031689406 NA [10,] 0.0414943559 -0.047773758 NA [11,] 0.0163848303 -0.049650730 NA [12,] -0.0243211170 0.019794446 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/1olrz1260099989.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/2hgot1260099990.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/3dq6c1260099990.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/4fw091260099990.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.023146224 -0.019373635 -2.208592e-02 -0.07667174 -0.028963399 [2,] -0.013633222 -0.013610242 -1.398821e-02 -0.03507857 -0.003023478 [3,] 0.003169821 -0.003772387 -4.310664e-03 -0.01484392 0.024737117 [4,] 0.022642386 0.024110540 -1.330059e-05 0.01522789 0.035814877 [5,] 0.042769303 0.041828092 1.208736e-03 0.04835080 0.067602778 [,6] [,7] [,8] [,9] [,10] [1,] -0.028290170 -0.0116301716 -0.0001255465 -0.022757092 0.008932314 [2,] -0.012606789 0.0003327824 0.0027998102 -0.014160927 0.012101306 [3,] 0.009173124 0.0112743670 0.0136476387 -0.007095848 0.019869349 [4,] 0.017819923 0.0195317114 0.0193653734 0.015496850 0.024705355 [5,] 0.041176029 0.0236354130 0.0230884797 0.036528277 0.030392715 [,11] [,12] [1,] -0.038841444 -0.025828194 [2,] -0.023493247 -0.020005680 [3,] -0.009344226 -0.008740637 [4,] 0.002405743 0.018563310 [5,] 0.022651414 0.029733828 $n [1] 8 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.01709425 -0.02629863 -0.012656255 -0.04488613 0.001543478 -0.008997218 [2,] 0.02343390 0.01875386 0.004034927 0.01519828 0.047930756 0.027343465 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.0001909237 0.003754952 -0.02480699 0.01234242 -0.024810685 -0.03177342 [2,] 0.0227396576 0.023540326 0.01061530 0.02739628 0.006122232 0.01429214 $out [1] -0.06869773 -0.04853633 0.06776912 0.06161938 -0.05033696 -0.03767829 [7] -0.08849217 $group [1] 1 3 3 7 8 10 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.023146223786389, -0.0136332217413643, 0.00316982088184399, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5l27g1260099990.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.044796932 -0.037809315 -0.028135410 -0.077314203 -0.065748833 [2,] -0.025794734 -0.022560231 -0.023490054 0.002213808 -0.035992860 [3,] -0.009270652 -0.007120173 0.001890667 0.020137962 -0.019024284 [4,] 0.016132066 0.005672680 0.021466470 0.079852503 -0.004927317 [5,] 0.045948312 0.025941027 0.052661468 0.104869150 0.005651867 [,6] [,7] [,8] [,9] [,10] [1,] -0.03803940 -0.0113999135 -0.033374782 -0.0123435861 -0.056406222 [2,] -0.01711396 -0.0059796974 -0.027157368 0.0004219856 -0.043882369 [3,] 0.02405137 0.0002268052 -0.003642426 0.0284640600 -0.029397296 [4,] 0.03213850 0.0130493079 0.013354245 0.0320056364 -0.010916230 [5,] 0.05244626 0.0273272764 0.043241114 0.0501198581 0.007381115 [,11] [,12] [1,] -0.04965073 -0.02432112 [2,] -0.01311001 -0.01474959 [3,] -0.00572569 0.01303548 [4,] 0.01721160 0.02970045 [5,] 0.04449879 0.04568967 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.03430866 -0.023980412 -0.02495664 -0.02622661 -0.037576125 -0.005361401 [2,] 0.01576736 0.009740066 0.02873798 0.06650254 -0.000472443 0.053464147 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.01113701 -0.02783531 0.009602813 -0.04908414 -0.02383327 -0.01350937 [2,] 0.01159062 0.02055046 0.047325307 -0.00971045 0.01238189 0.03958032 $out [1] -0.12029031 0.07013943 -0.11195634 -0.06658925 0.04149436 $group [1] 3 5 7 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0447969321801321, -0.0257947343506667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/61ni81260099990.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.022085923 -0.028290170 -0.03222153 -0.076671743 -0.028963399 [2,] -0.014904837 -0.020391093 -0.01321895 -0.033954984 0.003493551 [3,] -0.003872609 -0.006416754 0.01554141 0.001038638 0.016098011 [4,] 0.020028434 0.001855369 0.02699257 0.026810352 0.037337340 [5,] 0.067602778 0.023088480 0.04343235 0.061619381 0.048350805 [,6] [,7] [,8] [1,] -0.008548957 -0.088492173 -0.06869773 [2,] -0.002066214 -0.030799268 -0.06869773 [3,] 0.002824100 0.004403384 -0.06869773 [4,] 0.021638415 0.018005742 -0.06869773 [5,] 0.028910959 0.067769119 -0.06869773 $n [1] 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.01980591 -0.016563514 -0.002799330 -0.02667684 0.0006616319 -0.00798774 [2,] 0.01206069 0.003730007 0.033882155 0.02875412 0.0315343902 0.01363594 [,7] [,8] [1,] -0.01785691 -0.06869773 [2,] 0.02666368 -0.06869773 $out [1] -0.03767829 $group [1] 6 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(-0.0220859233394083, -0.0149048374921528, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/71w2v1260099990.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.011838034 -0.0148439243 -0.010543752 [2,] -0.005770726 -0.0079182427 -0.003860970 [3,] 0.002862622 -0.0003012832 0.003928358 [4,] 0.008991412 0.0124610028 0.010507419 [5,] 0.018422756 0.0247371173 0.018403331 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.003870488 -0.009596394 -0.002625161 [2,] 0.009595732 0.008993828 0.010481877 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.0118380340326593, -0.00577072595692162, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1olrz1260099989.ps tmp/1olrz1260099989.png") > system("convert tmp/2hgot1260099990.ps tmp/2hgot1260099990.png") > system("convert tmp/3dq6c1260099990.ps tmp/3dq6c1260099990.png") > system("convert tmp/4fw091260099990.ps tmp/4fw091260099990.png") > system("convert tmp/5l27g1260099990.ps tmp/5l27g1260099990.png") > system("convert tmp/61ni81260099990.ps tmp/61ni81260099990.png") > system("convert tmp/71w2v1260099990.ps tmp/71w2v1260099990.png") > > > proc.time() user system elapsed 1.195 0.932 1.797