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(126.51,131.02,136.51,138.04,132.92,129.61,122.96,124.04,121.29,124.56,118.53,113.14,114.15,122.17,129.23,131.19,129.12,128.28,126.83,138.13,140.52,146.83,135.14,131.84,125.7,128.98,133.25,136.76,133.24,128.54,121.08,120.23,119.08,125.75,126.89,126.6,121.89,123.44,126.46,129.49,127.78,125.29,119.02,119.96,122.86,131.89,132.73,135.01,136.71,142.73,144.43,144.93,138.75,130.22,122.19,128.4,140.43,153.5,149.33,142.97) > 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] 60 > (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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 126.51 114.15 125.70 121.89 136.71 NA [2,] 131.02 122.17 128.98 123.44 142.73 NA [3,] 136.51 129.23 133.25 126.46 144.43 NA [4,] 138.04 131.19 136.76 129.49 144.93 NA [5,] 132.92 129.12 133.24 127.78 138.75 NA [6,] 129.61 128.28 128.54 125.29 130.22 NA [7,] 122.96 126.83 121.08 119.02 122.19 NA [8,] 124.04 138.13 120.23 119.96 128.40 NA [9,] 121.29 140.52 119.08 122.86 140.43 NA [10,] 124.56 146.83 125.75 131.89 153.50 NA [11,] 118.53 135.14 126.89 132.73 149.33 NA [12,] 113.14 131.84 126.60 135.01 142.97 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 4.51 8.02 3.28 1.55 6.02 NA [2,] 5.49 7.06 4.27 3.02 1.70 NA [3,] 1.53 1.96 3.51 3.03 0.50 NA [4,] -5.12 -2.07 -3.52 -1.71 -6.18 NA [5,] -3.31 -0.84 -4.70 -2.49 -8.53 NA [6,] -6.65 -1.45 -7.46 -6.27 -8.03 NA [7,] 1.08 11.30 -0.85 0.94 6.21 NA [8,] -2.75 2.39 -1.15 2.90 12.03 NA [9,] 3.27 6.31 6.67 9.03 13.07 NA [10,] -6.03 -11.69 1.14 0.84 -4.17 NA [11,] -5.39 -3.30 -0.29 2.28 -6.36 NA [12,] 1.01 -6.14 -4.71 1.70 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/1ph9z1258141217.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/2jz7y1258141217.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/354yw1258141217.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/4gp4h1258141217.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] [1,] 121.89 122.17 126.46 129.49 127.78 128.28 119.02 119.96 119.08 124.56 [2,] 121.89 123.44 129.23 131.19 129.12 128.28 121.08 120.23 121.29 125.75 [3,] 125.70 128.98 133.25 136.76 132.92 128.54 122.19 124.04 122.86 131.89 [4,] 126.51 131.02 136.51 138.04 133.24 129.61 122.96 128.40 140.43 146.83 [5,] 126.51 131.02 144.43 144.93 138.75 130.22 122.96 138.13 140.52 153.50 [,11] [,12] [1,] 118.53 126.60 [2,] 126.89 126.60 [3,] 132.73 131.84 [4,] 135.14 135.01 [5,] 135.14 142.97 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 122.4355 123.624 128.1060 131.9198 130.0088 127.6002 120.8616 118.2671 [2,] 128.9645 134.336 138.3940 141.6002 135.8312 129.4798 123.5184 129.8129 [,9] [,10] [,11] [,12] [1,] 109.3357 116.9949 126.9006 125.8975 [2,] 136.3843 146.7851 138.5594 137.7825 $out [1] 114.15 136.71 142.73 125.29 126.83 149.33 113.14 $group [1] 1 1 2 6 7 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(121.89, 121.89, 125.7, 126.51, 126.51, 122.17, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/51dqh1258141217.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,] 1.55 1.70 0.50 -6.18 -4.70 -8.03 -0.85 -2.75 3.27 -11.69 -6.36 -6.140 [2,] 3.28 3.02 1.53 -5.12 -4.70 -7.46 0.94 -1.15 6.31 -6.03 -5.39 -5.425 [3,] 4.51 4.27 1.96 -3.52 -3.31 -6.65 1.08 2.39 6.67 -4.17 -3.30 -1.850 [4,] 6.02 5.49 3.03 -2.07 -2.49 -6.27 6.21 2.90 9.03 0.84 -0.29 1.355 [5,] 8.02 7.06 3.51 -1.71 -0.84 -6.27 11.30 2.90 13.07 1.14 2.28 1.700 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 2.573923 2.524704 0.9001038 -5.675122 -4.871580 -7.490851 -2.643769 [2,] 6.446077 6.015296 3.0198962 -1.364878 -1.748420 -5.809149 4.803769 [,8] [,9] [,10] [,11] [,12] [1,] -0.4717198 4.748055 -9.0243247 -6.9036472 -7.2062 [2,] 5.2517198 8.591945 0.6843247 0.3036472 3.5062 $out [1] -8.53 -1.45 12.03 $group [1] 5 6 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1.55000000000000, 3.27999999999999, 4.51000000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6deth1258141217.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,] 113.140 114.150 119.080 119.020 122.190 NA [2,] 122.125 127.555 123.390 122.375 133.465 NA [3,] 125.535 130.210 126.745 125.875 141.580 NA [4,] 131.970 136.635 131.110 130.690 144.680 NA [5,] 138.040 146.830 136.760 135.010 153.500 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 121.0446 126.0686 123.2239 122.0825 136.4648 NA [2,] 130.0254 134.3514 130.2661 129.6675 146.6952 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(113.14, 122.125, 125.535, 131.97, 138.04, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7u0501258141217.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,] 122.416 122.190 122.0200 [2,] 127.270 124.870 125.7725 [3,] 129.790 130.410 130.8325 [4,] 133.250 132.825 132.0250 [5,] 136.506 136.760 136.2900 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 127.0625 126.7817 127.9807 [2,] 132.5175 134.0383 133.6843 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(122.416, 127.27, 129.79, 133.25, 136.506, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1ph9z1258141217.ps tmp/1ph9z1258141217.png") > system("convert tmp/2jz7y1258141217.ps tmp/2jz7y1258141217.png") > system("convert tmp/354yw1258141217.ps tmp/354yw1258141217.png") > system("convert tmp/4gp4h1258141217.ps tmp/4gp4h1258141217.png") > system("convert tmp/51dqh1258141217.ps tmp/51dqh1258141217.png") > system("convert tmp/6deth1258141217.ps tmp/6deth1258141217.png") > system("convert tmp/7u0501258141217.ps tmp/7u0501258141217.png") > > > proc.time() user system elapsed 1.232 0.975 1.522