R version 2.8.0 (2008-10-20) Copyright (C) 2008 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(97.7,101.5,119.6,108.1,117.8,125.5,89.2,92.3,104.6,122.8,96.0,94.6,93.3,101.1,114.2,104.7,113.3,118.2,83.6,73.9,99.5,97.7,103.0,106.3,92.2,101.8,122.8,111.8,106.3,121.5,81.9,85.4,110.9,117.3,106.3,105.5,101.3,105.9,126.3,111.9,108.9,127.2,94.2,85.7,116.2,107.2,110.6,112.0,104.5,112.0,132.8,110.8,128.7,136.8,94.9,88.8,123.2,125.3,122.7,125.7,116.3,118.7,142.0,127.9,131.9,152.3,110.8,99.1,135.0,133.2,131.0,133.9,119.9,136.9,148.9,145.1,142.4,159.6,120.7,109.0,142.0) > 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] 81 > (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] 7 7 7 7 7 7 7 7 7 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 97.7 93.3 92.2 101.3 104.5 116.3 119.9 [2,] 101.5 101.1 101.8 105.9 112.0 118.7 136.9 [3,] 119.6 114.2 122.8 126.3 132.8 142.0 148.9 [4,] 108.1 104.7 111.8 111.9 110.8 127.9 145.1 [5,] 117.8 113.3 106.3 108.9 128.7 131.9 142.4 [6,] 125.5 118.2 121.5 127.2 136.8 152.3 159.6 [7,] 89.2 83.6 81.9 94.2 94.9 110.8 120.7 [8,] 92.3 73.9 85.4 85.7 88.8 99.1 109.0 [9,] 104.6 99.5 110.9 116.2 123.2 135.0 142.0 [10,] 122.8 97.7 117.3 107.2 125.3 133.2 NA [11,] 96.0 103.0 106.3 110.6 122.7 131.0 NA [12,] 94.6 106.3 105.5 112.0 125.7 133.9 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 3.8 7.8 9.6 4.6 7.5 2.4 17.0 [2,] 18.1 13.1 21.0 20.4 20.8 23.3 12.0 [3,] -11.5 -9.5 -11.0 -14.4 -22.0 -14.1 -3.8 [4,] 9.7 8.6 -5.5 -3.0 17.9 4.0 -2.7 [5,] 7.7 4.9 15.2 18.3 8.1 20.4 17.2 [6,] -36.3 -34.6 -39.6 -33.0 -41.9 -41.5 -38.9 [7,] 3.1 -9.7 3.5 -8.5 -6.1 -11.7 -11.7 [8,] 12.3 25.6 25.5 30.5 34.4 35.9 33.0 [9,] 18.2 -1.8 6.4 -9.0 2.1 -1.8 NA [10,] -26.8 5.3 -11.0 3.4 -2.6 -2.2 NA [11,] -1.4 3.3 -0.8 1.4 3.0 2.9 NA [12,] -1.3 -14.1 -4.2 -7.5 -9.4 -14.0 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/18f561228911650.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/2o4xy1228911650.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/32ymt1228911650.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/4pt3r1228911650.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] [,10] [,11] [1,] 92.2 101.10 114.2 104.70 106.3 118.20 81.90 73.90 99.50 97.70 96.00 [2,] 95.5 101.65 121.2 109.45 111.1 123.50 86.40 85.55 107.75 107.20 103.00 [3,] 101.3 105.90 126.3 111.80 117.8 127.20 94.20 88.80 116.20 120.05 108.45 [4,] 110.4 115.35 137.4 119.90 130.3 144.55 102.85 95.70 129.10 125.30 122.70 [5,] 119.9 118.70 148.9 127.90 142.4 159.60 120.70 109.00 142.00 133.20 131.00 [,12] [1,] 94.60 [2,] 105.50 [3,] 109.15 [4,] 125.70 [5,] 133.90 $n [1] 7 7 7 7 7 7 7 7 7 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 92.40196 97.71858 116.6256 105.5594 106.3341 114.6293 84.37633 82.73858 [2,] 110.19804 114.08142 135.9744 118.0406 129.2659 139.7707 104.02367 94.86142 [,9] [,10] [,11] [,12] [1,] 103.4501 108.3749 95.74286 96.12035 [2,] 128.9499 131.7251 121.15714 122.17965 $out [1] 136.9 145.1 $group [1] 2 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(92.2, 95.5, 101.3, 110.4, 119.9, 101.1, 101.65, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5dmc71228911650.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] [,11] [,12] [1,] 2.4 12.0 -14.40 -5.50 4.90 -41.90 -11.7 25.50 -9.00 -26.8 -1.40 -14.10 [2,] 4.2 15.6 -14.25 -2.85 7.90 -40.55 -10.7 25.55 -1.80 -11.0 -0.80 -14.00 [3,] 7.5 20.4 -11.50 4.00 15.20 -38.90 -8.5 30.50 0.15 -2.4 2.15 -8.45 [4,] 8.7 20.9 -10.25 9.15 17.75 -35.45 -1.5 33.70 6.40 3.4 3.00 -4.20 [5,] 9.6 23.3 -9.50 17.90 20.40 -33.00 3.5 35.90 18.20 5.3 3.30 -1.30 $n [1] 7 7 7 7 7 7 7 7 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 4.812673 17.23493 -13.888735 -3.166206 9.317739 -41.94564 -13.994092 [2,] 10.187327 23.56507 -9.111265 11.166206 21.082261 -35.85436 -3.005908 [,8] [,9] [,10] [,11] [,12] [1,] 25.63295 -5.139265 -11.688465 -0.3011227 -14.771317 [2,] 35.36705 5.439265 6.888465 4.6011227 -2.128683 $out [1] 17.0 -22.0 -3.8 12.3 $group [1] 1 3 3 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(2.40000000000001, 4.2, 7.5, 8.7, 9.6, 12, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6d8l51228911650.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,] 89.20 83.60 81.90 94.20 88.80 99.10 109.0 [2,] 95.30 95.50 97.00 103.60 107.65 117.50 120.7 [3,] 103.05 102.05 106.30 109.75 122.95 131.45 142.0 [4,] 118.70 109.80 114.55 114.10 127.20 134.45 145.1 [5,] 125.50 118.20 122.80 127.20 136.80 152.30 159.6 $n [1] 12 12 12 12 12 12 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 92.3771 95.52767 98.29533 104.9609 114.0331 123.719 129.1493 [2,] 113.7229 108.57233 114.30467 114.5391 131.8669 139.181 154.8507 $out [1] 73.9 85.7 $group [1] 2 4 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(89.2, 95.3, 103.05, 118.7, 125.5, 83.6, 95.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7j9jr1228911650.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,] 90.6000 88.800 90.6250 [2,] 107.3643 103.600 105.7250 [3,] 115.0929 110.475 114.3313 [4,] 120.0500 118.925 119.5625 [5,] 134.4429 127.200 134.0250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 109.3068 103.4852 108.0199 [2,] 120.8789 117.4648 120.6426 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(90.6, 107.364285714286, 115.092857142857, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/18f561228911650.ps tmp/18f561228911650.png") > system("convert tmp/2o4xy1228911650.ps tmp/2o4xy1228911650.png") > system("convert tmp/32ymt1228911650.ps tmp/32ymt1228911650.png") > system("convert tmp/4pt3r1228911650.ps tmp/4pt3r1228911650.png") > system("convert tmp/5dmc71228911650.ps tmp/5dmc71228911650.png") > system("convert tmp/6d8l51228911650.ps tmp/6d8l51228911650.png") > system("convert tmp/7j9jr1228911650.ps tmp/7j9jr1228911650.png") > > > proc.time() user system elapsed 1.201 0.923 1.522