R version 2.7.2 (2008-08-25) 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(299.63,305.945,382.252,348.846,335.367,373.617,312.612,312.232,337.161,331.476,350.103,345.127,297.256,295.979,361.007,321.803,354.937,349.432,290.979,349.576,327.625,349.377,336.777,339.134,323.321,318.86,373.583,333.03,408.556,414.646,291.514,348.857,349.368,375.765,364.136,349.53,348.167,332.856,360.551,346.969,392.815,372.02,371.027,342.672,367.343,390.786,343.785,362.6,349.468,340.624,369.536,407.782,392.239) > 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] 53 > (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 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] 299.630 297.256 323.321 348.167 349.468 [2,] 305.945 295.979 318.860 332.856 340.624 [3,] 382.252 361.007 373.583 360.551 369.536 [4,] 348.846 321.803 333.030 346.969 407.782 [5,] 335.367 354.937 408.556 392.815 392.239 [6,] 373.617 349.432 414.646 372.020 NA [7,] 312.612 290.979 291.514 371.027 NA [8,] 312.232 349.576 348.857 342.672 NA [9,] 337.161 327.625 349.368 367.343 NA [10,] 331.476 349.377 375.765 390.786 NA [11,] 350.103 336.777 364.136 343.785 NA [12,] 345.127 339.134 349.530 362.600 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 6.315 -1.277 -4.461 -15.311 -8.844 [2,] 76.307 65.028 54.723 27.695 28.912 [3,] -33.406 -39.204 -40.553 -13.582 38.246 [4,] -13.479 33.134 75.526 45.846 -15.543 [5,] 38.250 -5.505 6.090 -20.795 NA [6,] -61.005 -58.453 -123.132 -0.993 NA [7,] -0.380 58.597 57.343 -28.355 NA [8,] 24.929 -21.951 0.511 24.671 NA [9,] -5.685 21.752 26.397 23.443 NA [10,] 18.627 -12.600 -11.629 -47.001 NA [11,] -4.976 2.357 -14.606 18.815 NA [12,] -47.871 -15.813 -1.363 -13.132 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/1632r1225549441.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/2vrca1225549441.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/3c3ka1225549441.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/4gof21225549441.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] [1,] 297.256 295.979 360.551 321.803 335.367 349.4320 290.9790 312.2320 [2,] 299.630 305.945 361.007 333.030 354.937 360.7260 291.2465 327.4520 [3,] 323.321 318.860 369.536 346.969 392.239 372.8185 302.0630 345.7645 [4,] 348.167 332.856 373.583 348.846 392.815 394.1315 341.8195 349.2165 [5,] 349.468 340.624 382.252 348.846 408.556 414.6460 371.0270 349.5760 [,9] [,10] [,11] [,12] [1,] 327.6250 331.4760 336.7770 339.1340 [2,] 332.3930 340.4265 340.2810 342.1305 [3,] 343.2645 362.5710 346.9440 347.3285 [4,] 358.3555 383.2755 357.1195 356.0650 [5,] 367.3430 390.7860 364.1360 362.6000 $n [1] 5 5 5 5 5 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 289.0249 299.8448 360.6498 335.7935 365.4745 346.4282 262.1103 328.5705 [2,] 357.6171 337.8752 378.4222 358.1445 419.0035 399.2088 342.0157 362.9585 [,9] [,10] [,11] [,12] [1,] 322.7541 328.7203 333.6416 336.3202 [2,] 363.7749 396.4217 360.2464 358.3368 $out [1] 407.782 $group [1] 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(297.256, 299.63, 323.321, 348.167, 349.468, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5pjxl1225549441.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] [1,] -15.311 27.695 -40.553 -15.543 -20.7950 -123.1320 -28.3550 -21.951 -5.6850 [2,] -8.844 28.912 -39.204 -13.479 -13.1500 -92.0685 -14.3675 -10.720 8.0335 [3,] -4.461 54.723 -33.406 33.134 0.2925 -59.7290 28.4815 12.591 22.5975 [4,] -1.277 65.028 -13.582 45.846 22.1700 -29.7230 57.9700 24.800 24.9200 [5,] 6.315 76.307 -13.582 75.526 38.2500 -0.9930 58.5970 24.929 26.3970 [,10] [,11] [,12] [1,] -47.0010 -14.6060 -47.8710 [2,] -29.8005 -9.7910 -31.8420 [3,] -12.1145 -1.3095 -14.4725 [4,] 3.4990 10.5860 -7.2475 [5,] 18.6270 18.8150 -1.3630 $n [1] 5 5 5 5 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -9.8078231 29.20353 -51.51044 -8.784896 -27.6103 -108.98194 -28.66512 [2,] 0.8858231 80.24247 -15.30156 75.052896 28.1953 -10.47606 85.62813 [,8] [,9] [,10] [,11] [,12] [1,] -15.4698 9.257165 -38.42111 -17.40733 -33.902155 [2,] 40.6518 35.937835 14.19211 14.78833 4.957155 $out [1] 38.246 $group [1] 3 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-15.3110000000000, -8.844, -4.46100000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6i6vz1225549441.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,] 299.6300 290.9790 291.5140 332.8560 340.624 [2,] 312.4220 309.5295 328.1755 345.3770 349.468 [3,] 336.2640 337.9555 349.4490 361.5755 369.536 [4,] 349.4745 349.5040 374.6740 371.5235 392.239 [5,] 382.2520 361.0070 414.6460 392.8150 407.782 $n [1] 12 12 12 12 5 $conf [,1] [,2] [,3] [,4] [,5] [1,] 319.3641 319.7229 328.2407 349.6499 339.3141 [2,] 353.1639 356.1881 370.6573 373.5011 399.7579 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(299.63, 312.422, 336.264, 349.4745, 382.252, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/712il1225549441.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,] 316.5330 302.0630 309.2980 [2,] 330.9513 333.2928 332.4183 [3,] 348.8990 346.9565 346.0708 [4,] 365.6184 366.0535 364.7530 [5,] 377.4288 392.2390 375.1236 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 333.0871 332.0141 331.3227 [2,] 364.7109 361.8989 360.8188 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(316.533, 330.951325, 348.899, 365.6184, 377.42875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1632r1225549441.ps tmp/1632r1225549441.png") > system("convert tmp/2vrca1225549441.ps tmp/2vrca1225549441.png") > system("convert tmp/3c3ka1225549441.ps tmp/3c3ka1225549441.png") > system("convert tmp/4gof21225549441.ps tmp/4gof21225549441.png") > system("convert tmp/5pjxl1225549441.ps tmp/5pjxl1225549441.png") > system("convert tmp/6i6vz1225549441.ps tmp/6i6vz1225549441.png") > system("convert tmp/712il1225549441.ps tmp/712il1225549441.png") > > > proc.time() user system elapsed 1.416 0.949 1.845