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(167.16,179.84,174.44,180.35,193.17,195.16,202.43,189.91,195.98,212.09,205.81,204.31,196.07,199.98,199.1,198.31,195.72,223.04,238.41,259.73,326.54,335.15,321.81,368.62,369.59,425,439.72,362.23,328.76,348.55,328.18,329.34,295.55,237.38,226.85,220.14,239.36,224.69,230.98,233.47,256.7,253.41,224.95,210.37,191.09,198.85,211.04,206.25,201.19,194.37,191.08,192.87,181.61,157.67,196.14,246.35,271.9) > 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] 57 > (np <- floor(n / par1)) [1] 4 > 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 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] 167.16 196.07 369.59 239.36 201.19 [2,] 179.84 199.98 425.00 224.69 194.37 [3,] 174.44 199.10 439.72 230.98 191.08 [4,] 180.35 198.31 362.23 233.47 192.87 [5,] 193.17 195.72 328.76 256.70 181.61 [6,] 195.16 223.04 348.55 253.41 157.67 [7,] 202.43 238.41 328.18 224.95 196.14 [8,] 189.91 259.73 329.34 210.37 246.35 [9,] 195.98 326.54 295.55 191.09 271.90 [10,] 212.09 335.15 237.38 198.85 NA [11,] 205.81 321.81 226.85 211.04 NA [12,] 204.31 368.62 220.14 206.25 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 12.68 3.91 55.41 -14.67 -6.82 [2,] -5.40 -0.88 14.72 6.29 -3.29 [3,] 5.91 -0.79 -77.49 2.49 1.79 [4,] 12.82 -2.59 -33.47 23.23 -11.26 [5,] 1.99 27.32 19.79 -3.29 -23.94 [6,] 7.27 15.37 -20.37 -28.46 38.47 [7,] -12.52 21.32 1.16 -14.58 50.21 [8,] 6.07 66.81 -33.79 -19.28 25.55 [9,] 16.11 8.61 -58.17 7.76 NA [10,] -6.28 -13.34 -10.53 12.19 NA [11,] -1.50 46.81 -6.71 -4.79 NA [12,] -8.24 0.97 19.22 -5.06 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/1yyy61289932370.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/2yyy61289932370.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/39pyr1289932370.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/49pyr1289932370.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,] 167.16 179.84 174.44 180.35 181.61 157.67 196.14 189.91 191.09 198.850 [2,] 196.07 194.37 191.08 192.87 193.17 195.16 202.43 210.37 195.98 205.470 [3,] 201.19 199.98 199.10 198.31 195.72 223.04 224.95 246.35 271.90 224.735 [4,] 239.36 224.69 230.98 233.47 256.70 253.41 238.41 259.73 295.55 286.265 [5,] 239.36 224.69 230.98 233.47 328.76 253.41 238.41 329.34 326.54 335.150 [,11] [,12] [1,] 205.810 204.310 [2,] 208.425 205.280 [3,] 218.945 213.195 [4,] 274.330 294.380 [5,] 321.810 368.620 $n [1] 5 5 5 5 5 5 5 5 5 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 170.6014 178.5560 170.9068 169.6221 150.8299 181.8807 199.5266 211.4723 [2,] 231.7786 221.4040 227.2932 226.9979 240.6101 264.1993 250.3734 281.2277 [,9] [,10] [,11] [,12] [1,] 201.5441 160.9070 166.8801 142.806 [2,] 342.2559 288.5630 271.0100 283.584 $out [1] 369.59 425.00 439.72 362.23 348.55 328.18 $group [1] 1 2 3 4 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(167.16, 196.07, 201.19, 239.36, 239.36, 179.84, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5jyfu1289932370.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,] -14.67 -5.40 -0.79 -33.47 -23.94 -28.46 -14.58 -33.79 -58.170 -13.340 [2,] -6.82 -3.29 -0.79 -11.26 -3.29 -20.37 -12.52 -19.28 -25.205 -11.935 [3,] 3.91 -0.88 1.79 -2.59 1.99 7.27 1.16 6.07 8.185 -8.405 [4,] 12.68 6.29 2.49 12.82 19.79 15.37 21.32 25.55 12.360 2.955 [5,] 12.68 14.72 5.91 23.23 27.32 38.47 50.21 66.81 16.110 12.190 [,11] [,12] [1,] -6.710 -8.240 [2,] -5.750 -6.650 [3,] -3.145 -2.045 [4,] 22.655 10.095 [5,] 46.810 19.220 $n [1] 5 5 5 5 5 5 5 5 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -9.86865 -7.649204 -0.5276397 -19.60487 -14.31827 -17.98379 -22.75126 [2,] 17.68865 5.889204 4.1076397 14.42487 18.29827 32.52379 25.07126 [,8] [,9] [,10] [,11] [,12] [1,] -25.60677 -21.49135 -20.1681 -25.58495 -15.27355 [2,] 37.74677 37.86135 3.3581 19.29495 11.18355 $out [1] 55.41 -77.49 $group [1] 1 3 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-14.6700000000000, -6.81999999999999, 3.91, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6jyfu1289932370.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] [1,] 167.160 195.720 220.140 191.090 181.61 [2,] 180.095 198.705 266.465 208.310 191.08 [3,] 194.165 230.725 329.050 224.820 194.37 [4,] 203.370 324.175 365.910 236.415 201.19 [5,] 212.090 368.620 439.720 256.700 201.19 $n [1] 12 12 12 12 9 $conf [,1] [,2] [,3] [,4] [,5] [1,] 183.5491 173.4973 283.6925 212.0011 189.0454 [2,] 204.7809 287.9527 374.4075 237.6389 199.6946 $out [1] 157.67 246.35 271.90 $group [1] 5 5 5 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(167.16, 180.095, 194.165, 203.37, 212.09, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7jyfu1289932370.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,] 231.1920 195.7200 209.5300 [2,] 235.1200 199.5400 215.4425 [3,] 243.0768 216.0700 224.6100 [4,] 247.1020 224.8425 233.2812 [5,] 256.2120 246.3500 245.7650 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 237.6117 204.5294 216.4736 [2,] 248.5418 227.6106 232.7464 $out [1] 271.9 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(231.192, 235.12, 243.07675, 247.102, 256.212, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1yyy61289932370.ps tmp/1yyy61289932370.png",intern=TRUE)) character(0) > try(system("convert tmp/2yyy61289932370.ps tmp/2yyy61289932370.png",intern=TRUE)) character(0) > try(system("convert tmp/39pyr1289932370.ps tmp/39pyr1289932370.png",intern=TRUE)) character(0) > try(system("convert tmp/49pyr1289932370.ps tmp/49pyr1289932370.png",intern=TRUE)) character(0) > try(system("convert tmp/5jyfu1289932370.ps tmp/5jyfu1289932370.png",intern=TRUE)) character(0) > try(system("convert tmp/6jyfu1289932370.ps tmp/6jyfu1289932370.png",intern=TRUE)) character(0) > try(system("convert tmp/7jyfu1289932370.ps tmp/7jyfu1289932370.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.165 0.928 1.849