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(100,87.14054095,112.0054296,112.312101,109.474134,104.9746116,100.4926851,104.2154743,120.1768388,112.1028355,108.1481575,116.802197,102.1699512,95.15358705,120.6707808,111.5234277,119.9669448,113.3697401,110.0717661,111.5567342,132.424212,107.900558,122.1626615,124.3992258,110.4450505,101.5874013,122.3203962,125.2582826,125.4411543,108.9902468,118.9243879,116.7242723,134.1724901,116.8530994,124.5732995,130.9914031,123.4239103,111.4536725,124.5135991,139.2589613,129.8596099,112.3460359,131.381655,133.0004776,134.3220552,144.2379719,134.1278719,150.1891559,140.722563,114.8389975,143.1973003,140.2738676,112.1248303,102.8951536,100.5090242,103.3513901,111.4134533,104.5887587,101.7840983,114.7007441,108.7426474) > 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] 61 > (np <- floor(n / par1)) [1] 5 > 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 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 100.00000 102.16995 110.4451 123.4239 140.7226 108.7426 [2,] 87.14054 95.15359 101.5874 111.4537 114.8390 NA [3,] 112.00543 120.67078 122.3204 124.5136 143.1973 NA [4,] 112.31210 111.52343 125.2583 139.2590 140.2739 NA [5,] 109.47413 119.96694 125.4412 129.8596 112.1248 NA [6,] 104.97461 113.36974 108.9902 112.3460 102.8952 NA [7,] 100.49269 110.07177 118.9244 131.3817 100.5090 NA [8,] 104.21547 111.55673 116.7243 133.0005 103.3514 NA [9,] 120.17684 132.42421 134.1725 134.3221 111.4135 NA [10,] 112.10284 107.90056 116.8531 144.2380 104.5888 NA [11,] 108.14816 122.16266 124.5733 134.1279 101.7841 NA [12,] 116.80220 124.39923 130.9914 150.1892 114.7007 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -12.8594590 -7.016364 -8.8576492 -11.970238 -25.883566 NA [2,] 24.8648886 25.517194 20.7329949 13.059927 28.358303 NA [3,] 0.3066714 -9.147353 2.9378864 14.745362 -2.923433 NA [4,] -2.8379670 8.443517 0.1828717 -9.399351 -28.149037 NA [5,] -4.4995224 -6.597205 -16.4509075 -17.513574 -9.229677 NA [6,] -4.4819265 -3.297974 9.9341411 19.035619 -2.386129 NA [7,] 3.7227892 1.484968 -2.2001156 1.618823 2.842366 NA [8,] 15.9613645 20.867478 17.4482178 1.321578 8.062063 NA [9,] -8.0740033 -24.523654 -17.3193907 9.915917 -6.824695 NA [10,] -3.9546780 14.262103 7.7202001 -10.110100 -2.804660 NA [11,] 8.6540395 2.236564 6.4181036 16.061284 12.916646 NA [12,] -14.6322458 -13.954175 -7.5674928 -9.466593 -5.958097 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/1f3tm1257954796.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/2sv2b1257954796.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/3pqf11257954796.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/4k3xw1257954796.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] [1,] 100.0000 87.14054 120.6708 111.5234 109.4741 102.8952 100.4927 103.3514 [2,] 102.1700 95.15359 120.6708 112.3121 112.1248 104.9746 100.5090 104.2155 [3,] 109.5938 101.58740 122.3204 125.2583 119.9669 108.9902 110.0718 111.5567 [4,] 123.4239 111.45367 124.5136 139.2590 125.4412 112.3460 118.9244 116.7243 [5,] 140.7226 114.83900 124.5136 140.2739 129.8596 113.3697 131.3817 133.0005 [,9] [,10] [,11] [,12] [1,] 111.4135 104.5888 101.7841 114.7007 [2,] 120.1768 107.9006 108.1482 116.8022 [3,] 132.4242 112.1028 122.1627 124.3992 [4,] 134.1725 116.8531 124.5733 130.9914 [5,] 134.3221 116.8531 134.1279 150.1892 $n [1] 6 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 95.88436 90.0698 119.6051 106.2177 110.5577 103.7816 97.05952 102.7180 [2,] 123.30334 113.1050 125.0357 144.2989 129.3762 114.1989 123.08402 120.3954 [,9] [,10] [,11] [,12] [1,] 122.5349 105.7770 110.5567 114.3732 [2,] 142.3135 118.4287 133.7686 134.4253 $out [1] 112.0054 143.1973 144.2380 $group [1] 3 3 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(100, 102.1699512, 109.59384895, 123.4239103, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5c36a1257954796.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] [1,] -12.859459 20.73299 -9.1473531 -9.3993514 -17.513574 -4.481927 1.484968 [2,] -12.859459 20.73299 -2.9234327 -9.3993514 -16.450907 -3.297974 1.484968 [3,] -11.970238 24.86489 0.3066714 -2.8379670 -9.229677 -2.386129 1.618823 [4,] -8.857649 25.51719 2.9378864 0.1828717 -6.597205 9.934141 2.842366 [5,] -7.016364 28.35830 2.9378864 8.4435171 -4.499522 19.035619 3.722789 [,8] [,9] [,10] [,11] [,12] [1,] 1.321578 -24.523654 -10.110100 2.236564 -14.632246 [2,] 8.062063 -17.319391 -3.954678 6.418104 -13.954175 [3,] 15.961365 -8.074003 -2.804660 8.654040 -9.466593 [4,] 17.448218 -6.824695 7.720200 12.916646 -7.567493 [5,] 20.867478 -6.824695 14.262103 16.061284 -5.958097 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -14.797907 21.48439 -3.834922 -9.608742 -16.192278 -11.73591 0.6596887 [2,] -9.142569 28.24539 4.448265 3.932808 -2.267075 6.96365 2.5779565 [,8] [,9] [,10] [,11] [,12] [1,] 9.329131 -15.4895291 -11.054100 4.062186 -13.979407 [2,] 22.593598 -0.6584775 5.444779 13.245893 -4.953779 $out [1] -25.883566 13.059927 14.745362 -28.149037 -2.200116 9.915917 $group [1] 1 2 3 4 7 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-12.85945905, -12.85945905, -11.9702378, -8.8576492, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6541h1257954796.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] [1,] 100.0000 95.15359 101.5874 111.4537 100.5090 108.7426 [2,] 102.3541 108.98616 113.5847 123.9688 103.1233 108.7426 [3,] 108.8111 112.46324 120.6224 132.1911 111.7691 108.7426 [4,] 112.2075 121.41672 125.3497 136.7905 127.5564 108.7426 [5,] 120.1768 132.42421 134.1725 150.1892 143.1973 108.7426 $n [1] 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 104.3169 106.7936 115.2563 126.3430 100.6250 108.7426 [2,] 113.3053 118.1329 125.9885 138.0392 122.9133 108.7426 $out [1] 87.14054 $group [1] 1 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(100, 102.3540797, 108.81114575, 112.20746825, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7z4hf1257954796.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,] 102.0348 101.5874 103.3036 [2,] 113.0228 109.8328 110.0933 [3,] 117.6479 116.0349 114.3688 [4,] 125.1334 123.3598 123.2445 [5,] 127.4165 132.4242 127.1747 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 112.1242 109.8651 108.3704 [2,] 123.1717 122.2046 120.3671 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(102.03483986, 113.02278668, 117.64793122, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > system("convert tmp/1f3tm1257954796.ps tmp/1f3tm1257954796.png") > system("convert tmp/2sv2b1257954796.ps tmp/2sv2b1257954796.png") > system("convert tmp/3pqf11257954796.ps tmp/3pqf11257954796.png") > system("convert tmp/4k3xw1257954796.ps tmp/4k3xw1257954796.png") > system("convert tmp/5c36a1257954796.ps tmp/5c36a1257954796.png") > system("convert tmp/6541h1257954796.ps tmp/6541h1257954796.png") > system("convert tmp/7z4hf1257954796.ps tmp/7z4hf1257954796.png") > > > proc.time() user system elapsed 1.272 0.998 19.769