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(89507,87562,85209,82360,79054,79069,107551,115759,115585,110260,103444,102303,101397,97994,94044,91159,87239,89235,118647,125620,125154,117529,109459,108483,107137,104699,100804,96066,91971,93228,120144,127233,127166,118194,109940,106683,102834,99882,96666,92540,88744,89321,115870,122401,122030,113802,105791,103076,98658,96945,92497,90687,88796,90015,113228,118711,117460,106556,97347,92657,93118,89037,83570,81693,75956,73993,97088,102394,96549,89727,82336,82653,82303,79596,74472,73562,66618,69029,89899,93774,90305,83799,80320,82497) > 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] 84 > (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] 7 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 89507 101397 107137 102834 98658 93118 82303 NA [2,] 87562 97994 104699 99882 96945 89037 79596 NA [3,] 85209 94044 100804 96666 92497 83570 74472 NA [4,] 82360 91159 96066 92540 90687 81693 73562 NA [5,] 79054 87239 91971 88744 88796 75956 66618 NA [6,] 79069 89235 93228 89321 90015 73993 69029 NA [7,] 107551 118647 120144 115870 113228 97088 89899 NA [8,] 115759 125620 127233 122401 118711 102394 93774 NA [9,] 115585 125154 127166 122030 117460 96549 90305 NA [10,] 110260 117529 118194 113802 106556 89727 83799 NA [11,] 103444 109459 109940 105791 97347 82336 80320 NA [12,] 102303 108483 106683 103076 92657 82653 82497 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1945 -3403 -2438 -2952 -1713 -4081 -2707 NA [2,] -2353 -3950 -3895 -3216 -4448 -5467 -5124 NA [3,] -2849 -2885 -4738 -4126 -1810 -1877 -910 NA [4,] -3306 -3920 -4095 -3796 -1891 -5737 -6944 NA [5,] 15 1996 1257 577 1219 -1963 2411 NA [6,] 28482 29412 26916 26549 23213 23095 20870 NA [7,] 8208 6973 7089 6531 5483 5306 3875 NA [8,] -174 -466 -67 -371 -1251 -5845 -3469 NA [9,] -5325 -7625 -8972 -8228 -10904 -6822 -6506 NA [10,] -6816 -8070 -8254 -8011 -9209 -7391 -3479 NA [11,] -1141 -976 -3257 -2715 -4690 317 2177 NA [12,] -906 -1346 -3849 -4418 461 -350 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/1aubd1259783318.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/2jdqi1259783318.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/3ln6i1259783318.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/4t0nz1259783318.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,] 82303.0 79596.0 74472.0 73562.0 66618 69029 89899.0 93774.0 90305 [2,] 91312.5 88299.5 84389.5 82026.5 77505 76531 102319.5 109076.5 106067 [3,] 98658.0 96945.0 92497.0 90687.0 87239 89235 113228.0 118711.0 117460 [4,] 102115.5 98938.0 95355.0 91849.5 88770 89668 117258.5 124010.5 123592 [5,] 107137.0 104699.0 100804.0 96066.0 91971 93228 120144.0 127233.0 127166 [,10] [,11] [,12] [1,] 83799.0 80320.0 82497.0 [2,] 98141.5 89841.5 87655.0 [3,] 110260.0 103444.0 102303.0 [4,] 115665.5 107625.0 104879.5 [5,] 118194.0 109940.0 108483.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 92206.62 90591.86 85948.58 84820.86 80511.72 81389.8 104306.7 109792.7 [2,] 105109.38 103298.14 99045.42 96553.14 93966.28 97080.2 122149.3 127629.3 [,9] [,10] [,11] [,12] [1,] 106994.4 99794.95 92823.98 92016.8 [2,] 127925.6 120725.05 114064.02 112589.2 $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(82303, 91312.5, 98658, 102115.5, 107137, 79596, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5spbx1259783318.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] [,9] [,10] [1,] -4081.0 -5467.0 -4738.0 -6944 15.0 20870 3875.0 -3469.0 -10904 -9209.0 [2,] -3177.5 -4786.0 -3505.5 -4916 296.0 23154 5394.5 -2360.0 -8600 -8162.0 [3,] -2707.0 -3950.0 -2849.0 -3920 1219.0 26549 6531.0 -466.0 -7625 -8011.0 [4,] -2191.5 -3555.5 -1843.5 -3551 1626.5 27699 7031.0 -272.5 -6664 -7103.5 [5,] -1713.0 -2353.0 -910.0 -1891 2411.0 29412 8208.0 -67.0 -5325 -6816.0 [,11] [,12] [1,] -4690.0 -4418 [2,] -2986.0 -3849 [3,] -1141.0 -1126 [4,] -329.5 -350 [5,] 2177.0 461 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3295.823 -4684.835 -3841.520 -4735.156 424.4469 23834.8 5553.709 [2,] -2118.177 -3215.165 -1856.480 -3104.844 2013.5531 29263.2 7508.291 [,8] [,9] [,10] [,11] [,12] [1,] -1712.6213 -8781.148 -8643.119 -2727.4189 -3382.968 [2,] 780.6213 -6468.852 -7378.881 445.4189 1130.968 $out [1] -1963 -5845 -3479 $group [1] 5 8 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-4081, -3177.5, -2707, -2191.5, -1713, -5467, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6mjj71259783318.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] [,8] [1,] 79054.0 87239.0 91971 88744 88796 73993.0 66618.0 NA [2,] 83784.5 92601.5 98435 94603 91592 82014.5 74017.0 NA [3,] 95905.0 104940.0 106910 102955 97146 86303.5 81311.5 NA [4,] 108905.5 118088.0 119169 114836 109892 94833.5 86849.0 NA [5,] 115759.0 125620.0 127233 122401 118711 102394.0 93774.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 84447.14 93315.44 97453.08 93726.6 88799.25 80456.67 75458.74 NA [2,] 107362.86 116564.56 116366.92 112183.4 105492.75 92150.33 87164.26 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(79054, 83784.5, 95905, 108905.5, 115759, 87239, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7d1h91259783318.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,] 82625.43 87239.0 83099.50 [2,] 88237.79 91592.0 88405.12 [3,] 96664.71 100480.5 96490.62 [4,] 107306.71 111744.0 108346.25 [5,] 115127.43 118711.0 116543.50 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 87967.25 91289.04 87395.34 [2,] 105362.18 109671.96 105585.91 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(82625.4285714286, 88237.7857142857, 96664.7142857143, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1aubd1259783318.ps tmp/1aubd1259783318.png") > system("convert tmp/2jdqi1259783318.ps tmp/2jdqi1259783318.png") > system("convert tmp/3ln6i1259783318.ps tmp/3ln6i1259783318.png") > system("convert tmp/4t0nz1259783318.ps tmp/4t0nz1259783318.png") > system("convert tmp/5spbx1259783318.ps tmp/5spbx1259783318.png") > system("convert tmp/6mjj71259783318.ps tmp/6mjj71259783318.png") > system("convert tmp/7d1h91259783318.ps tmp/7d1h91259783318.png") > > > proc.time() user system elapsed 1.161 0.929 1.463