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. Natural language support but running in an English locale 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(225,224,223,221,219,218,219,221,222,222,223,225,226,229,235,229,231,229,226,232,234,230,232,231,237,241,256,255,264,259,253,258,265,258,257,255,255,255,275,276,278,268,253,257,255,253,245,248,246,243,260,262,262,251,236,238,239,243,233,238,232,224,238,236,231,209,179,179,165,174,163,166,164,152,163,167,157,138,111,110,91,100,88,89) > 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] 84 > (np <- floor(n / par1)) [1] 7 > 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 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 225 226 237 255 246 232 164 NA [2,] 224 229 241 255 243 224 152 NA [3,] 223 235 256 275 260 238 163 NA [4,] 221 229 255 276 262 236 167 NA [5,] 219 231 264 278 262 231 157 NA [6,] 218 229 259 268 251 209 138 NA [7,] 219 226 253 253 236 179 111 NA [8,] 221 232 258 257 238 179 110 NA [9,] 222 234 265 255 239 165 91 NA [10,] 222 230 258 253 243 174 100 NA [11,] 223 232 257 245 233 163 88 NA [12,] 225 231 255 248 238 166 89 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 3 4 0 -3 -8 -12 NA [2,] -1 6 15 20 17 14 11 NA [3,] -2 -6 -1 1 2 -2 4 NA [4,] -2 2 9 2 0 -5 -10 NA [5,] -1 -2 -5 -10 -11 -22 -19 NA [6,] 1 -3 -6 -15 -15 -30 -27 NA [7,] 2 6 5 4 2 0 -1 NA [8,] 1 2 7 -2 1 -14 -19 NA [9,] 0 -4 -7 -2 4 9 9 NA [10,] 1 2 -1 -8 -10 -11 -12 NA [11,] 2 -1 -2 3 5 3 1 NA [12,] 1 6 0 -2 -6 -2 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/freestat/rcomp/tmp/19mml1282249270.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/freestat/rcomp/tmp/2kel61282249270.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/freestat/rcomp/tmp/3kel61282249270.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/freestat/rcomp/tmp/4dnkr1282249270.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] [,12] [1,] 225.0 224 223 221.0 219 209.0 179.0 179.0 165.0 174 163 166.0 [2,] 225.5 224 229 225.0 225 213.5 199.0 200.0 193.5 198 193 195.5 [3,] 232.0 229 238 236.0 231 229.0 226.0 232.0 234.0 230 232 231.0 [4,] 241.5 242 258 258.5 263 255.0 244.5 247.5 247.0 248 239 243.0 [5,] 255.0 255 275 276.0 278 268.0 253.0 258.0 265.0 258 257 255.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 222.4451 218.2507 220.6817 215.9943 208.307 204.2169 198.8281 203.6338 [2,] 241.5549 239.7493 255.3183 256.0057 253.693 253.7831 253.1719 260.3662 [,9] [,10] [,11] [,12] [1,] 202.0507 200.1408 204.5295 202.6338 [2,] 265.9493 259.8592 259.4705 259.3662 $out [1] 164 152 163 167 157 138 111 110 91 100 88 89 $group [1] 1 2 3 4 5 6 7 8 9 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(225, 225.5, 232, 241.5, 255, 224, 224, 229, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5dnkr1282249270.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,] -12.0 -1.0 -6.0 -10.0 -22.0 -30.0 -1.0 -19.0 -7.0 -12.0 -2 -6 [2,] -5.5 8.5 -2.0 -3.5 -15.0 -21.0 1.0 -8.0 -3.0 -10.5 0 -2 [3,] -1.0 14.0 -1.0 0.0 -10.0 -15.0 2.0 1.0 0.0 -8.0 2 -1 [4,] 1.5 16.0 1.5 2.0 -3.5 -4.5 4.5 1.5 6.5 0.0 3 1 [5,] 4.0 20.0 4.0 9.0 -1.0 1.0 6.0 7.0 9.0 2.0 5 1 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -5.180287 9.521121 -3.090144 -3.284511 -16.867614 -24.853534 -0.09014354 [2,] 3.180287 18.478879 1.090144 3.284511 -3.132386 -5.146466 4.09014354 [,8] [,9] [,10] [,11] [,12] [1,] -4.673247 -5.673247 -14.270431 0.2084484 -2.9350969 [2,] 6.673247 5.673247 -1.729569 3.7915516 0.9350969 $out [1] 6 $group [1] 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-12, -5.5, -1, 1.5, 4, -1, 8.5, 14, 16, 20, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6nwkc1282249270.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] [,8] [1,] 218.0 226.0 253.0 245.0 233.0 163.0 88.0 NA [2,] 220.0 229.0 254.0 253.0 238.0 170.0 95.5 NA [3,] 222.0 230.5 256.5 255.0 243.0 194.0 124.5 NA [4,] 223.5 232.0 258.5 271.5 255.5 231.5 160.0 NA [5,] 225.0 235.0 265.0 278.0 262.0 238.0 167.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 220.4036 229.1317 254.4475 246.5620 235.0181 165.9494 95.08112 NA [2,] 223.5964 231.8683 258.5525 263.4380 250.9819 222.0506 153.91888 NA $out [1] 237 241 $group [1] 3 3 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(218, 220, 222, 223.5, 225, 226, 229, 230.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7nwkc1282249270.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,] 205.8571 226.0 216.000 [2,] 210.5714 229.5 221.000 [3,] 218.7857 231.5 228.375 [4,] 230.5000 233.0 238.000 [5,] 235.7143 238.0 244.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 209.6962 229.9036 220.6212 [2,] 227.8753 233.0964 236.1288 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(205.857142857143, 210.571428571429, 218.785714285714, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/19mml1282249270.ps tmp/19mml1282249270.png",intern=TRUE)) character(0) > try(system("convert tmp/2kel61282249270.ps tmp/2kel61282249270.png",intern=TRUE)) character(0) > try(system("convert tmp/3kel61282249270.ps tmp/3kel61282249270.png",intern=TRUE)) character(0) > try(system("convert tmp/4dnkr1282249270.ps tmp/4dnkr1282249270.png",intern=TRUE)) character(0) > try(system("convert tmp/5dnkr1282249270.ps tmp/5dnkr1282249270.png",intern=TRUE)) character(0) > try(system("convert tmp/6nwkc1282249270.ps tmp/6nwkc1282249270.png",intern=TRUE)) character(0) > try(system("convert tmp/7nwkc1282249270.ps tmp/7nwkc1282249270.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.003 1.514 2.307