R version 2.7.0 (2008-04-22) 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. 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(284.4,212.8,226.9,308.4,262,227.9,236.1,320.4,271.9,232.8,237,313.4,261.4,226.8,249.9,314.3,286.1,226.5,260.4,311.4,294.7,232.6,257.2,339.2,279.1,249.8,269.8,345.7,293.8,254.7,277.5,363.4,313.4,272.8,300.1,369.5,330.8,287.8,305.9,386.1,335.2,288,308.3,402.3,352.8,316.1,324.9,404.8,393,318.9,327,442.3,383.1,331.6,361.4,445.9,386.6,357.2,373.6,466.2,409.6,369.8,378.6,487,419.2,376.7,392.8,506.1,458.4,387.4,426.9,565,464.8,444.5,449.5,556.1,499.6,451.9,434.9,553.8,510,432.9,453.2,547.6,485.8,452.6,456.6,565.7,514.8,464.3,430.9,588.3,503.1,442.6,448,554.5,504.5,427.3,473.1,526.2,547.5,440.2,468.7,574.5,492.6,432.6,479.8,575.7,474.6,405.3,434.6,535.1,452.6,429.5,417.2,551.8,464,416.6,422.9,553.6,458.6,427.6,429.2,534.2,481.7,416,440.2,538.7,473.8,439.9,446.8,597.5,467.2,439.4,447.4,568.5,485.9,442.1,430.5,600,464.5,423.6,437,574,443,410,420,532,432,420,411,512,449,382) > 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] 154 > (np <- floor(n / par1)) [1] 12 > 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] 13 13 13 13 13 13 13 13 13 13 12 12 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 284.4 261.4 279.1 330.8 393.0 409.6 464.8 485.8 504.5 474.6 458.6 467.2 [2,] 212.8 226.8 249.8 287.8 318.9 369.8 444.5 452.6 427.3 405.3 427.6 439.4 [3,] 226.9 249.9 269.8 305.9 327.0 378.6 449.5 456.6 473.1 434.6 429.2 447.4 [4,] 308.4 314.3 345.7 386.1 442.3 487.0 556.1 565.7 526.2 535.1 534.2 568.5 [5,] 262.0 286.1 293.8 335.2 383.1 419.2 499.6 514.8 547.5 452.6 481.7 485.9 [6,] 227.9 226.5 254.7 288.0 331.6 376.7 451.9 464.3 440.2 429.5 416.0 442.1 [7,] 236.1 260.4 277.5 308.3 361.4 392.8 434.9 430.9 468.7 417.2 440.2 430.5 [8,] 320.4 311.4 363.4 402.3 445.9 506.1 553.8 588.3 574.5 551.8 538.7 600.0 [9,] 271.9 294.7 313.4 352.8 386.6 458.4 510.0 503.1 492.6 464.0 473.8 464.5 [10,] 232.8 232.6 272.8 316.1 357.2 387.4 432.9 442.6 432.6 416.6 439.9 423.6 [11,] 237.0 257.2 300.1 324.9 373.6 426.9 453.2 448.0 479.8 422.9 446.8 437.0 [12,] 313.4 339.2 369.5 404.8 466.2 565.0 547.6 554.5 575.7 553.6 597.5 574.0 [,13] [1,] 443 [2,] 410 [3,] 420 [4,] 532 [5,] 432 [6,] 420 [7,] 411 [8,] 512 [9,] 449 [10,] 382 [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -71.6 -34.6 -29.3 -43.0 -74.1 -39.8 -20.3 -33.2 -77.2 -69.3 -31.0 [2,] 14.1 23.1 20.0 18.1 8.1 8.8 5.0 4.0 45.8 29.3 1.6 [3,] 81.5 64.4 75.9 80.2 115.3 108.4 106.6 109.1 53.1 100.5 105.0 [4,] -46.4 -28.2 -51.9 -50.9 -59.2 -67.8 -56.5 -50.9 21.3 -82.5 -52.5 [5,] -34.1 -59.6 -39.1 -47.2 -51.5 -42.5 -47.7 -50.5 -107.3 -23.1 -65.7 [6,] 8.2 33.9 22.8 20.3 29.8 16.1 -17.0 -33.4 28.5 -12.3 24.2 [7,] 84.3 51.0 85.9 94.0 84.5 113.3 118.9 157.4 105.8 134.6 98.5 [8,] -48.5 -16.7 -50.0 -49.5 -59.3 -47.7 -43.8 -85.2 -81.9 -87.8 -64.9 [9,] -39.1 -62.1 -40.6 -36.7 -29.4 -71.0 -77.1 -60.5 -60.0 -47.4 -33.9 [10,] 4.2 24.6 27.3 8.8 16.4 39.5 20.3 5.4 47.2 6.3 6.9 [11,] 76.4 82.0 69.4 79.9 92.6 138.1 94.4 106.5 95.9 130.7 150.7 [12,] -52.0 -60.1 -38.7 -11.8 -56.6 -100.2 -61.8 -50.0 -101.1 -95.0 -130.3 [,12] [,13] [1,] -27.8 -33 [2,] 8.0 10 [3,] 121.1 112 [4,] -82.6 -100 [5,] -43.8 -12 [6,] -11.6 -9 [7,] 169.5 101 [8,] -135.5 -63 [9,] -40.9 -67 [10,] 13.4 NA [11,] 137.0 NA [12,] -131.0 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/1vgvi1228415753.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/2x9c71228415753.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/3t1nb1228415753.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/4b6f71228415753.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,] 261.4 212.8 226.9 308.4 262.0 226.5 236.1 311.4 271.9 232.6 237.0 313.40 [2,] 330.8 287.8 305.9 386.1 335.2 288.0 308.3 402.3 352.8 316.1 312.5 387.15 [3,] 443.0 405.3 420.0 526.2 432.0 416.0 411.0 512.0 458.4 387.4 424.9 550.60 [4,] 467.2 427.6 447.4 535.1 485.9 440.2 430.9 553.8 473.8 432.6 447.4 569.50 [5,] 504.5 452.6 473.1 568.5 547.5 464.3 468.7 600.0 510.0 442.6 479.8 597.50 $n [1] 13 13 13 13 13 13 13 13 13 13 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 383.2277 344.0378 357.9928 460.9062 365.9613 349.3040 357.2751 445.6107 [2,] 502.7723 466.5622 482.0072 591.4938 498.0387 482.6960 464.7249 578.3893 [,9] [,10] [,11] [,12] [1,] 405.3762 336.3482 363.3712 467.4289 [2,] 511.4238 438.4518 486.4288 633.7711 $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(261.4, 330.8, 443, 467.2, 504.5, 212.8, 287.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5zi8f1228415753.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,] -77.2 1.6 53.1 -82.6 -65.7 -33.4 51.0 -87.8 -77.1 4.20 69.40 -131.00 [2,] -69.3 8.0 80.2 -67.8 -51.5 -11.6 85.9 -81.9 -62.1 6.60 80.95 -100.65 [3,] -34.6 10.0 105.0 -52.5 -47.2 16.1 101.0 -59.3 -47.4 14.90 95.15 -60.95 [4,] -31.0 20.0 109.1 -50.9 -39.1 24.2 118.9 -48.5 -39.1 25.95 133.85 -51.00 [5,] -20.3 29.3 121.1 -28.2 -23.1 33.9 157.4 -16.7 -29.4 47.20 150.70 -11.80 $n [1] 13 13 13 13 13 13 13 13 13 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -51.38356 4.741442 92.33564 -59.9058 -52.63384 0.4119691 86.53897 [2,] -17.81644 15.258558 117.66436 -45.0942 -41.76616 31.7880309 115.46103 [,8] [,9] [,10] [,11] [,12] [1,] -73.93632 -57.4789 6.074335 71.02195 -83.5957 [2,] -44.66368 -37.3211 23.725665 119.27805 -38.3043 $out [1] 45.8 21.3 -100.0 -107.3 -12.0 169.5 -135.5 $group [1] 2 4 4 5 5 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-77.2, -69.3, -34.6, -31, -20.3, 1.59999999999997, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6may41228415753.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,] 212.80 226.50 249.80 287.80 318.90 369.8 432.9 430.90 427.30 405.30 416.00 [2,] 230.35 241.25 271.30 307.10 344.40 383.0 447.0 450.30 454.45 420.05 434.55 [3,] 249.50 260.90 286.45 327.85 378.35 414.4 459.0 475.05 486.20 443.60 452.70 [4,] 296.40 303.05 329.55 369.45 417.65 472.7 528.8 534.65 536.85 504.85 507.95 [5,] 320.40 339.20 369.50 404.80 466.20 565.0 556.1 588.30 575.70 553.60 597.50 [,12] [,13] [1,] 423.60 382 [2,] 438.20 411 [3,] 455.95 426 [4,] 527.20 449 [5,] 600.00 449 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 219.3742 232.7126 259.8818 299.4117 344.9402 373.4872 421.6905 436.5774 [2,] 279.6258 289.0874 313.0182 356.2883 411.7598 455.3128 496.3095 513.5226 [,9] [,10] [,11] [,12] [,13] [1,] 448.6168 404.9222 419.2218 415.3565 407.0137 [2,] 523.7832 482.2778 486.1782 496.5435 444.9863 $out [1] 532 512 $group [1] 13 13 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" NA Warning message: In bxp(list(stats = c(212.8, 230.35, 249.5, 296.4, 320.4, 226.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7rs5o1228415753.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,] 359.4308 387.40 357.7000 [2,] 370.6885 413.50 371.9750 [3,] 394.1596 428.45 390.9500 [4,] 443.7077 485.20 436.9500 [5,] 488.4167 550.60 481.6125 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 360.8551 395.7471 361.3145 [2,] 427.4642 461.1529 420.5855 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(359.430769230769, 370.688461538462, 394.159615384615, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1vgvi1228415753.ps tmp/1vgvi1228415753.png") > system("convert tmp/2x9c71228415753.ps tmp/2x9c71228415753.png") > system("convert tmp/3t1nb1228415753.ps tmp/3t1nb1228415753.png") > system("convert tmp/4b6f71228415753.ps tmp/4b6f71228415753.png") > system("convert tmp/5zi8f1228415753.ps tmp/5zi8f1228415753.png") > system("convert tmp/6may41228415753.ps tmp/6may41228415753.png") > system("convert tmp/7rs5o1228415753.ps tmp/7rs5o1228415753.png") > > > proc.time() user system elapsed 3.522 1.717 3.769