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(390,389,388,386,406,405,390,380,381,381,382,384,394,393,388,381,399,396,378,368,369,373,374,379,385,385,395,387,400,390,365,350,365,374,367,375,382,380,378,363,375,366,341,326,338,345,336,342,347,360,360,334,347,336,305,289,303,308,294,299,306,313,321,287,296,283,248,235,241,244,237,241,251,259,264,229,237,228,197,182,182,176,172,175,185,195,206,175,185,174,140,130,133,130,131,136,149,155,161,131,145,134,93,87,86,80,79,91,108,105,112,78,87,74,32,25,26,25,22,36) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 390 394 385 382 347 306 251 185 149 108 NA [2,] 389 393 385 380 360 313 259 195 155 105 NA [3,] 388 388 395 378 360 321 264 206 161 112 NA [4,] 386 381 387 363 334 287 229 175 131 78 NA [5,] 406 399 400 375 347 296 237 185 145 87 NA [6,] 405 396 390 366 336 283 228 174 134 74 NA [7,] 390 378 365 341 305 248 197 140 93 32 NA [8,] 380 368 350 326 289 235 182 130 87 25 NA [9,] 381 369 365 338 303 241 182 133 86 26 NA [10,] 381 373 374 345 308 244 176 130 80 25 NA [11,] 382 374 367 336 294 237 172 131 79 22 NA [12,] 384 379 375 342 299 241 175 136 91 36 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -1 -1 0 -2 13 7 8 10 6 -3 NA [2,] -1 -5 10 -2 0 8 5 11 6 7 NA [3,] -2 -7 -8 -15 -26 -34 -35 -31 -30 -34 NA [4,] 20 18 13 12 13 9 8 10 14 9 NA [5,] -1 -3 -10 -9 -11 -13 -9 -11 -11 -13 NA [6,] -15 -18 -25 -25 -31 -35 -31 -34 -41 -42 NA [7,] -10 -10 -15 -15 -16 -13 -15 -10 -6 -7 NA [8,] 1 1 15 12 14 6 0 3 -1 1 NA [9,] 0 4 9 7 5 3 -6 -3 -6 -1 NA [10,] 1 1 -7 -9 -14 -7 -4 1 -1 -3 NA [11,] 2 5 8 6 5 4 3 5 12 14 NA [12,] 10 6 7 5 7 10 10 13 17 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/1feml1281011453.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/2feml1281011453.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/3feml1281011453.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/485361281011453.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,] 108.0 105.0 112.0 78.0 87.0 74.0 32.0 25 26 25 22.0 36 [2,] 185.0 195.0 206.0 175.0 185.0 174.0 140.0 130 133 130 131.0 136 [3,] 326.5 336.5 340.5 310.5 321.5 309.5 276.5 262 272 276 265.5 270 [4,] 385.0 385.0 388.0 381.0 399.0 390.0 365.0 350 365 373 367.0 375 [5,] 394.0 393.0 395.0 387.0 406.0 405.0 390.0 380 381 381 382.0 384 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 226.5720 241.5684 249.5655 207.5742 214.5771 201.5778 164.0810 152.0792 [2,] 426.4280 431.4316 431.4345 413.4258 428.4229 417.4222 388.9190 371.9208 [,9] [,10] [,11] [,12] [1,] 156.0836 154.5875 147.585 150.5861 [2,] 387.9164 397.4125 383.415 389.4139 $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(108, 185, 326.5, 385, 394, 105, 195, 336.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/585361281011453.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,] -3 -5.0 -35 8.0 -13.0 -42 -16.0 -1 -6.0 -14.0 2 5 [2,] -1 -1.0 -34 9.0 -11.0 -35 -15.0 1 -3.0 -7.0 4 7 [3,] 3 5.5 -28 12.5 -10.5 -31 -11.5 2 1.5 -3.5 5 10 [4,] 8 8.0 -8 14.0 -9.0 -25 -10.0 12 5.0 1.0 8 10 [5,] 13 11.0 -2 20.0 -9.0 -15 -6.0 15 9.0 1.0 14 13 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1.496759 1.003241 -40.99064 10.0018 -11.49928 -35.9964 -13.9982 -3.496039 [2,] 7.496759 9.996759 -15.00936 14.9982 -9.50072 -26.0036 -9.0018 7.496039 [,9] [,10] [,11] [,12] [1,] -2.497119 -7.497119 3.001441 8.42 [2,] 5.497119 0.497119 6.998559 11.58 $out [1] -1 -3 17 $group [1] 5 5 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-3, -1, 3, 8, 13, -5, -1, 5.5, 8, 11, -35, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6jx3r1281011453.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] [,9] [,10] [,11] [1,] 380.0 368.0 350.0 326.0 289 235.0 172.0 130 79.0 22.0 NA [2,] 381.5 373.5 366.0 339.5 301 241.0 179.0 132 86.5 25.5 NA [3,] 387.0 380.0 380.0 354.0 321 265.5 212.5 157 112.0 55.0 NA [4,] 390.0 393.5 388.5 376.5 347 301.0 244.0 185 147.0 96.0 NA [5,] 390.0 399.0 400.0 382.0 360 321.0 264.0 206 161.0 112.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 383.1231 370.8779 369.7376 337.1241 300.0191 238.1336 182.8531 132.8263 [2,] 390.8769 389.1221 390.2624 370.8759 341.9809 292.8664 242.1469 181.1737 [,9] [,10] [,11] [1,] 84.40554 22.84448 NA [2,] 139.59446 87.15552 NA $out [1] 406 405 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(380, 381.5, 387, 390, 390, 368, 373.5, 380, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7jx3r1281011453.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,] 237.2 262.0 243.5000 [2,] 243.0 271.0 252.7500 [3,] 262.0 293.0 269.5625 [4,] 288.7 324.0 294.1875 [5,] 297.3 340.5 303.0000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 241.1559 268.8263 250.6626 [2,] 282.8441 317.1737 288.4624 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(237.2, 243, 262, 288.7, 297.3, 262, 271, 293, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1feml1281011453.ps tmp/1feml1281011453.png",intern=TRUE)) character(0) > try(system("convert tmp/2feml1281011453.ps tmp/2feml1281011453.png",intern=TRUE)) character(0) > try(system("convert tmp/3feml1281011453.ps tmp/3feml1281011453.png",intern=TRUE)) character(0) > try(system("convert tmp/485361281011453.ps tmp/485361281011453.png",intern=TRUE)) character(0) > try(system("convert tmp/585361281011453.ps tmp/585361281011453.png",intern=TRUE)) character(0) > try(system("convert tmp/6jx3r1281011453.ps tmp/6jx3r1281011453.png",intern=TRUE)) character(0) > try(system("convert tmp/7jx3r1281011453.ps tmp/7jx3r1281011453.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.180 0.933 1.602