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(376,375,374,372,370,369,370,372,373,373,374,376,371,374,369,363,357,366,362,366,361,362,358,363,360,360,348,345,332,333,323,327,332,337,336,337,343,337,326,321,309,302,293,287,292,292,289,302,310,295,276,264,257,243,227,226,226,229,224,240,244,226,208,199,193,180,167,164,166,173,169,191,193,166,143,147,139,129,115,108,106,116,108,135) > 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,] 376 371 360 343 310 244 193 NA [2,] 375 374 360 337 295 226 166 NA [3,] 374 369 348 326 276 208 143 NA [4,] 372 363 345 321 264 199 147 NA [5,] 370 357 332 309 257 193 139 NA [6,] 369 366 333 302 243 180 129 NA [7,] 370 362 323 293 227 167 115 NA [8,] 372 366 327 287 226 164 108 NA [9,] 373 361 332 292 226 166 106 NA [10,] 373 362 337 292 229 173 116 NA [11,] 374 358 336 289 224 169 108 NA [12,] 376 363 337 302 240 191 135 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 3 0 -6 -15 -18 -27 NA [2,] -1 -5 -12 -11 -19 -18 -23 NA [3,] -2 -6 -3 -5 -12 -9 4 NA [4,] -2 -6 -13 -12 -7 -6 -8 NA [5,] -1 9 1 -7 -14 -13 -10 NA [6,] 1 -4 -10 -9 -16 -13 -14 NA [7,] 2 4 4 -6 -1 -3 -7 NA [8,] 1 -5 5 5 0 2 -2 NA [9,] 0 1 5 0 3 7 10 NA [10,] 1 -4 -1 -3 -5 -4 -8 NA [11,] 2 5 1 13 16 22 27 NA [12,] -5 -3 6 8 4 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/1e4sy1281785195.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/2e4sy1281785195.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/3pvr11281785195.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/4pvr11281785195.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,] 193.0 166.0 143.0 147.0 139.0 129.0 115.0 108.0 106.0 116.0 108.0 135.0 [2,] 277.0 260.5 242.0 231.5 225.0 211.5 197.0 195.0 196.0 201.0 196.5 215.5 [3,] 343.0 337.0 326.0 321.0 309.0 302.0 293.0 287.0 292.0 292.0 289.0 302.0 [4,] 365.5 367.0 358.5 354.0 344.5 349.5 342.5 346.5 346.5 349.5 347.0 350.0 [5,] 376.0 375.0 374.0 372.0 370.0 369.0 370.0 372.0 373.0 373.0 374.0 376.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,] 290.1492 273.3999 256.4281 247.845 237.6365 219.5886 206.1097 196.5266 [2,] 395.8508 400.6001 395.5719 394.155 380.3635 384.4114 379.8903 377.4734 [,9] [,10] [,11] [,12] [1,] 202.1238 203.3182 199.1238 221.6788 [2,] 381.8762 380.6818 378.8762 382.3212 $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(193, 277, 343, 365.5, 376, 166, 260.5, 337, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5pvr11281785195.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,] -27.0 -23.0 -12.0 -13 -14.0 -16.0 -7.0 -5.0 0.0 -8.0 1.0 -5 [2,] -16.5 -18.5 -7.5 -10 -11.5 -13.5 -4.5 -1.0 0.5 -4.5 3.5 -3 [3,] -6.0 -12.0 -5.0 -7 -7.0 -10.0 -1.0 1.0 3.0 -4.0 13.0 3 [4,] -0.5 -8.0 -2.5 -6 0.0 -6.5 3.0 3.5 6.0 -2.0 19.0 6 [5,] 3.0 -1.0 4.0 -2 9.0 1.0 4.0 5.0 10.0 1.0 27.0 8 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -15.554942 -18.270431 -7.985919 -9.388735 -13.8676145 -14.180287 -5.478879 [2,] 3.554942 -5.729569 -2.014081 -4.611265 -0.1323855 -5.819713 3.478879 [,8] [,9] [,10] [,11] [,12] [1,] -1.687327 -0.2845113 -5.492960 3.74365 -2.805291 [2,] 3.687327 6.2845113 -2.507040 22.25635 8.805291 $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(-27, -16.5, -6, -0.5, 3, -23, -18.5, -12, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6h48m1281785195.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,] 369.0 357.0 323.0 287.0 224.0 164.0 106.0 NA [2,] 371.0 361.5 332.0 292.0 226.5 168.0 111.5 NA [3,] 373.0 363.0 336.5 302.0 241.5 185.5 132.0 NA [4,] 374.5 367.5 346.5 323.5 270.0 203.5 145.0 NA [5,] 376.0 374.0 360.0 343.0 310.0 244.0 193.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 371.4036 360.2634 329.8865 287.6326 221.6594 169.3082 116.7204 NA [2,] 374.5964 365.7366 343.1135 316.3674 261.3406 201.6918 147.2796 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(369, 371, 373, 374.5, 376, 357, 361.5, 363, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7h48m1281785195.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,] 264.2857 287.0 269.750 [2,] 265.3571 292.0 271.500 [3,] 276.1429 302.0 281.625 [4,] 289.6429 323.5 296.500 [5,] 313.8571 343.0 321.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 265.0660 287.6326 270.2223 [2,] 287.2197 316.3674 293.0277 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(264.285714285714, 265.357142857143, 276.142857142857, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1e4sy1281785195.ps tmp/1e4sy1281785195.png",intern=TRUE)) character(0) > try(system("convert tmp/2e4sy1281785195.ps tmp/2e4sy1281785195.png",intern=TRUE)) character(0) > try(system("convert tmp/3pvr11281785195.ps tmp/3pvr11281785195.png",intern=TRUE)) character(0) > try(system("convert tmp/4pvr11281785195.ps tmp/4pvr11281785195.png",intern=TRUE)) character(0) > try(system("convert tmp/5pvr11281785195.ps tmp/5pvr11281785195.png",intern=TRUE)) character(0) > try(system("convert tmp/6h48m1281785195.ps tmp/6h48m1281785195.png",intern=TRUE)) character(0) > try(system("convert tmp/7h48m1281785195.ps tmp/7h48m1281785195.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.984 1.514 2.218