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(131.6,132.05,132.4,132.57,133.02,133.47,133.66,133.96,134.19,134.93,134.9,135.05,135.16,135.23,135.15,135.12,137.29,137.41,137.44,137.62,137.78,137.98,138.06,138.16,138.28,138.33,138.43,138.44,138.41,138.55,138.64,138.72,138.9,139.02,139.04,139.15,139.3,140.73,141.84,141.95,142.1,142.36,142.58,142.75,142.85,143.03,143.19,143.62,143.89,144.69,147.51,147.78,148.04,148.21,148.29,148.34,148.33,148.38,148.37,148.37) > 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,] 131.60 135.16 138.28 139.30 143.89 NA [2,] 132.05 135.23 138.33 140.73 144.69 NA [3,] 132.40 135.15 138.43 141.84 147.51 NA [4,] 132.57 135.12 138.44 141.95 147.78 NA [5,] 133.02 137.29 138.41 142.10 148.04 NA [6,] 133.47 137.41 138.55 142.36 148.21 NA [7,] 133.66 137.44 138.64 142.58 148.29 NA [8,] 133.96 137.62 138.72 142.75 148.34 NA [9,] 134.19 137.78 138.90 142.85 148.33 NA [10,] 134.93 137.98 139.02 143.03 148.38 NA [11,] 134.90 138.06 139.04 143.19 148.37 NA [12,] 135.05 138.16 139.15 143.62 148.37 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.45 0.07 0.05 1.43 0.80 NA [2,] 0.35 -0.08 0.10 1.11 2.82 NA [3,] 0.17 -0.03 0.01 0.11 0.27 NA [4,] 0.45 2.17 -0.03 0.15 0.26 NA [5,] 0.45 0.12 0.14 0.26 0.17 NA [6,] 0.19 0.03 0.09 0.22 0.08 NA [7,] 0.30 0.18 0.08 0.17 0.05 NA [8,] 0.23 0.16 0.18 0.10 -0.01 NA [9,] 0.74 0.20 0.12 0.18 0.05 NA [10,] -0.03 0.08 0.02 0.16 -0.01 NA [11,] 0.15 0.10 0.11 0.43 0.00 NA [12,] 0.11 0.12 0.15 0.27 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/17tko1258140703.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/2qt861258140703.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/34lpa1258140703.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/493fd1258140703.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,] 131.60 132.05 132.40 132.57 133.02 133.47 133.66 133.96 134.19 134.93 [2,] 135.16 135.23 135.15 135.12 137.29 137.41 137.44 137.62 137.78 137.98 [3,] 138.28 138.33 138.43 138.44 138.41 138.55 138.64 138.72 138.90 139.02 [4,] 139.30 140.73 141.84 141.95 142.10 142.36 142.58 142.75 142.85 143.03 [5,] 143.89 144.69 147.51 147.78 148.04 148.21 148.29 148.34 148.33 148.38 [,11] [,12] [1,] 134.90 135.05 [2,] 138.06 138.16 [3,] 139.04 139.15 [4,] 143.19 143.62 [5,] 148.37 148.37 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 135.3547 134.4437 133.7029 133.6139 135.0113 135.0523 135.0081 135.0952 [2,] 141.2053 142.2163 143.1571 143.2661 141.8087 142.0477 142.2719 142.3448 [,9] [,10] [,11] [,12] [1,] 135.3176 135.4517 135.4152 135.2920 [2,] 142.4824 142.5883 142.6648 143.0080 $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(131.6, 135.16, 138.28, 139.3, 143.89, 132.05, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5a75u1258140703.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,] 0.05 -0.08 -0.03 -0.03 0.12 0.03 0.05 -0.01 0.05 -0.03 0.10 0.110 [2,] 0.07 0.10 0.01 0.15 0.14 0.08 0.08 0.10 0.12 -0.01 0.10 0.115 [3,] 0.45 0.35 0.11 0.26 0.17 0.09 0.17 0.16 0.18 0.02 0.11 0.135 [4,] 0.80 1.11 0.17 0.45 0.26 0.19 0.18 0.18 0.20 0.08 0.15 0.210 [5,] 1.43 1.11 0.27 0.45 0.26 0.22 0.30 0.23 0.20 0.16 0.15 0.270 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.06581616 -0.3636635 -0.003055597 0.04802076 0.0852083 0.01227428 [2,] 0.96581616 1.0636635 0.223055597 0.47197924 0.2547917 0.16772572 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.09934025 0.1034722 0.1234722 -0.04359377 0.07467013 0.05995 [2,] 0.24065975 0.2165278 0.2365278 0.08359377 0.14532987 0.21005 $out [1] 2.82 2.17 0.45 0.74 0.43 0.00 $group [1] 2 4 5 9 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0.0500000000000114, 0.0699999999999932, 0.450000000000017, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6vp251258140703.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,] 131.600 135.120 138.280 140.730 147.510 NA [2,] 132.485 135.195 138.420 141.895 147.645 NA [3,] 133.565 137.425 138.595 142.470 148.250 NA [4,] 134.545 137.880 138.960 142.940 148.355 NA [5,] 135.050 138.160 139.150 143.620 148.380 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 132.6254 136.2004 138.3487 141.9934 147.9262 NA [2,] 134.5046 138.6496 138.8413 142.9466 148.5738 NA $out [1] 139.30 143.89 144.69 $group [1] 4 5 5 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(131.6, 132.485, 133.565, 134.545, 135.05, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7ojsx1258140703.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,] 137.646 138.280 137.2300 [2,] 139.119 138.420 138.5150 [3,] 140.061 138.595 139.9475 [4,] 140.539 138.960 140.4100 [5,] 140.870 139.150 140.8900 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 139.4133 138.3487 139.0832 [2,] 140.7087 138.8413 140.8118 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(137.646, 139.119, 140.061, 140.539, 140.87, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/17tko1258140703.ps tmp/17tko1258140703.png") > system("convert tmp/2qt861258140703.ps tmp/2qt861258140703.png") > system("convert tmp/34lpa1258140703.ps tmp/34lpa1258140703.png") > system("convert tmp/493fd1258140703.ps tmp/493fd1258140703.png") > system("convert tmp/5a75u1258140703.ps tmp/5a75u1258140703.png") > system("convert tmp/6vp251258140703.ps tmp/6vp251258140703.png") > system("convert tmp/7ojsx1258140703.ps tmp/7ojsx1258140703.png") > > > proc.time() user system elapsed 1.241 0.972 1.491