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(5714837,6015267,3415340,115612,615159,9714993,2414534,6614721,8414518,9513816,7713310,9212714,411919,3611400,710617,2410172,529569,269509,508895,308560,198397,977631,717793,57435,667628,57696,577284,677273,447200,327096,207113,616966,276765,866845,6725,576898,187085,966784,586883,476728,406853,217097,827154,526998,456977,427122,867075,157185,87059,657279,87566,67453,297553,227345,97406,187524,337651,397672,797737,717980,717784,897934,297932,108061,187927,887990,998136,288286,768212,988170,508491,438167) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 5714837 411919 667628 187085 87059 717784 NA [2,] 6015267 3611400 57696 966784 657279 897934 NA [3,] 3415340 710617 577284 586883 87566 297932 NA [4,] 115612 2410172 677273 476728 67453 108061 NA [5,] 615159 529569 447200 406853 297553 187927 NA [6,] 9714993 269509 327096 217097 227345 887990 NA [7,] 2414534 508895 207113 827154 97406 998136 NA [8,] 6614721 308560 616966 526998 187524 288286 NA [9,] 8414518 198397 276765 456977 337651 768212 NA [10,] 9513816 977631 866845 427122 397672 988170 NA [11,] 7713310 717793 6725 867075 797737 508491 NA [12,] 9212714 57435 576898 157185 717980 438167 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 300430 3199481 -609932 779699 570220 180150 NA [2,] -2599927 -2900783 519588 -379901 -569713 -600002 NA [3,] -3299728 1699555 99989 -110155 -20113 -189871 NA [4,] 499547 -1880603 -230073 -69875 230100 79866 NA [5,] 9099834 -260060 -120104 -189756 -70208 700063 NA [6,] -7300459 239386 -119983 610057 -129939 110146 NA [7,] 4200187 -200335 409853 -300156 90118 -709850 NA [8,] 1799797 -110163 -340201 -70021 150127 479926 NA [9,] 1099298 779234 590080 -29855 60021 219958 NA [10,] -1800506 -259838 -860120 439953 400065 -479679 NA [11,] 1499404 -660358 570173 -709890 -79757 -70324 NA [12,] -8800795 610193 -389813 -70126 -196 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/15s2y1272292561.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/2g2211272292561.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/3g2211272292561.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/4qt1m1272292561.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] [,6] [,7] [,8] [,9] [1,] 87059.0 57696 87566.0 67453 187927.0 217097.0 97406.0 187524 198397 [2,] 187085.0 657279 297932.0 108061 297553.0 227345.0 207113.0 288286 276765 [3,] 539773.5 932359 582083.5 296170 427026.5 298302.5 668024.5 417779 397314 [4,] 717784.0 3611400 710617.0 677273 529569.0 887990.0 998136.0 616966 768212 [5,] 717784.0 6015267 710617.0 677273 615159.0 887990.0 998136.0 616966 768212 [,10] [,11] [,12] [1,] 397672 6725 57435.0 [2,] 427122 508491 157185.0 [3,] 922238 757765 507532.5 [4,] 988170 867075 717980.0 [5,] 988170 867075 717980.0 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 197455.5 -973144.5 315888.3 -70990.12 277368.7 -127834.9 157789.1 [2,] 882091.5 2837862.5 848278.7 663330.12 576684.3 724439.9 1178259.9 [,8] [,9] [,10] [,11] [,12] [1,] 205769.8 80314.81 560343.9 526466.7 145801.6 [2,] 629788.2 714313.19 1284132.1 989063.3 869263.4 $out [1] 5714837 3415340 2410172 9714993 2414534 6614721 8414518 9513816 7713310 [10] 9212714 $group [1] 1 3 4 6 7 8 9 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(87059, 187085, 539773.5, 717784, 717784, 57696, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5qt1m1272292561.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] [,7] [,8] [1,] -609932 -2900783.0 -189871 -230073.0 -260060 -129939.0 -709850.0 -340201 [2,] 180150 -2599927.0 -189871 -230073.0 -189756 -129939.0 -300156.0 -110163 [3,] 435325 -584857.5 -65134 4995.5 -95156 -4918.5 -55108.5 40053 [4,] 779699 -379901.0 99989 230100.0 700063 239386.0 409853.0 479926 [5,] 779699 519588.0 99989 499547.0 700063 610057.0 409853.0 479926 [,9] [,10] [,11] [,12] [1,] -29855 -1800506.0 -709890.0 -389813 [2,] 60021 -860120.0 -660358.0 -389813 [3,] 405019 -369758.5 -75040.5 -70126 [4,] 779234 400065.0 570173.0 -196 [5,] 1099298 439953.0 1499404.0 -196 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 48596.53 -2016846 -252103.1 -291830.9 -669118 -243145.1 -513087.2 [2,] 822053.47 847131 121835.1 301821.9 478806 233308.1 402870.2 [,8] [,9] [,10] [,11] [,12] [1,] -340573.5 -58896.61 -1182618.5 -868772.7 -345428.4 [2,] 420679.5 868934.61 443101.5 718691.7 205176.4 $out [1] 3199481 -3299728 1699555 -1880603 9099834 -7300459 4200187 1799797 [9] -8800795 610193 $group [1] 1 3 3 4 5 6 7 8 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-609932, 180150, 435325, 779699, 779699, -2900783, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6qt1m1272292561.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] [,7] [1,] 115612 57435.0 6725 157185.0 67453.0 108061.0 NA [2,] 2914937 289034.5 241939 311975.0 92486.0 293109.0 NA [3,] 6314994 519232.0 512049 466852.5 262449.0 613137.5 NA [4,] 8813616 847712.0 642297 707018.5 527475.5 892962.0 NA [5,] 9714993 977631.0 866845 966784.0 797737.0 998136.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 3624567 264415.4 329443 286670.5 64047.37 339540.5 NA [2,] 9005421 774048.6 694655 647034.5 460850.63 886734.5 NA $out [1] 3611400 2410172 $group [1] 2 2 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(115612, 2914937, 6314994, 8813616, 9714993, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7cu3k1272292562.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,] 414043.5 296170.0 368542.8 [2,] 894071.7 407546.5 449492.1 [3,] 1582964.6 523653.0 492010.6 [4,] 1900367.4 712894.8 662126.5 [5,] 2195209.3 932359.0 761294.0 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1123986 384381.6 395026.6 [2,] 2041943 662924.4 588994.5 $out [1] 1833844 $group [1] 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(414043.5, 894071.666666667, 1582964.58333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/15s2y1272292561.ps tmp/15s2y1272292561.png",intern=TRUE)) character(0) > try(system("convert tmp/2g2211272292561.ps tmp/2g2211272292561.png",intern=TRUE)) character(0) > try(system("convert tmp/3g2211272292561.ps tmp/3g2211272292561.png",intern=TRUE)) character(0) > try(system("convert tmp/4qt1m1272292561.ps tmp/4qt1m1272292561.png",intern=TRUE)) character(0) > try(system("convert tmp/5qt1m1272292561.ps tmp/5qt1m1272292561.png",intern=TRUE)) character(0) > try(system("convert tmp/6qt1m1272292561.ps tmp/6qt1m1272292561.png",intern=TRUE)) character(0) > try(system("convert tmp/7cu3k1272292562.ps tmp/7cu3k1272292562.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.164 0.911 1.537