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(556,555,554,552,572,571,556,546,547,547,548,550,555,549,555,550,566,573,543,535,542,541,535,536,548,546,548,548,561,563,527,527,541,534,522,527,539,533,532,519,538,542,503,502,522,511,492,500,509,511,505,493,518,518,474,471,483,461,439,446,461,449,441,424,447,448,404,403,411,386,359,370,385,369,368,352,378,383,334,323,330,303,275,284,301,281,284,272,297,300,240,236,247,218,192,201,223,197,195,175,197,204,142,142,151,127,100,114,139,112,123,108,132,140,76,71,81,57,38,46) > 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)) > ari <- array(0,dim=par1) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[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,] 556 555 548 539 509 461 385 301 223 139 NA [2,] 555 549 546 533 511 449 369 281 197 112 NA [3,] 554 555 548 532 505 441 368 284 195 123 NA [4,] 552 550 548 519 493 424 352 272 175 108 NA [5,] 572 566 561 538 518 447 378 297 197 132 NA [6,] 571 573 563 542 518 448 383 300 204 140 NA [7,] 556 543 527 503 474 404 334 240 142 76 NA [8,] 546 535 527 502 471 403 323 236 142 71 NA [9,] 547 542 541 522 483 411 330 247 151 81 NA [10,] 547 541 534 511 461 386 303 218 127 57 NA [11,] 548 535 522 492 439 359 275 192 100 38 NA [12,] 550 536 527 500 446 370 284 201 114 46 NA > arr.sd <- array(NA,dim=par1) > arr.range <- array(NA,dim=par1) > arr.iqr <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) + arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) + arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) + } > overall.sd <- sqrt(var(x)) > overall.range <- max(x) - min(x) > overall.iqr <- quantile(x,0.75) - quantile(x,0.25) > postscript(file="/var/www/html/rcomp/tmp/1cubi1281951901.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.sd,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2cubi1281951901.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.range,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3mms31281951901.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.iqr,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4mms31281951901.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,] 139 112 123 108.0 132.0 140 76 71 81 57.0 38 46 [2,] 301 281 284 272.0 297.0 300 240 236 247 218.0 192 201 [3,] 485 480 473 458.5 482.5 483 439 437 447 423.5 399 408 [4,] 548 546 548 548.0 561.0 563 527 527 541 534.0 522 527 [5,] 556 555 555 552.0 572.0 573 556 546 547 547.0 548 550 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 361.5890 347.5954 341.0951 320.5994 350.5951 351.5947 295.6034 291.6048 [2,] 608.4110 612.4046 604.9049 596.4006 614.4049 614.4053 582.3966 582.3952 [,9] [,10] [,11] [,12] [1,] 300.1059 265.6138 234.1188 245.1174 [2,] 593.8941 581.3862 563.8812 570.8826 $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(139, 301, 485, 548, 556, 112, 281, 480, 546, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5mms31281951901.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,] 546.0 535.0 522.0 492.0 439 359.0 275.0 192.0 100.0 38.0 NA [2,] 547.5 538.5 527.0 502.5 466 394.5 313.0 227.0 134.5 64.0 NA [3,] 553.0 546.0 543.5 520.5 488 417.5 343.0 259.5 163.0 94.5 NA [4,] 556.0 555.0 548.0 535.5 510 447.5 373.5 290.5 197.0 127.5 NA [5,] 556.0 573.0 563.0 542.0 518 461.0 385.0 301.0 223.0 140.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,] 549.1231 538.4742 533.9218 505.4485 467.9313 393.3263 315.4055 230.5372 [2,] 556.8769 553.5258 553.0782 535.5515 508.0687 441.6737 370.5945 288.4628 [,9] [,10] [,11] [1,] 134.4933 65.53722 NA [2,] 191.5067 123.46278 NA $out [1] 572 571 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(546, 547.5, 553, 556, 556, 535, 538.5, 546, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6fv9o1281951901.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) > names(z) <- list('S.D.','Range','IQR') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 152.7519 417.0 223.750 [2,] 161.7068 436.5 238.500 [3,] 169.6247 455.0 253.125 [4,] 178.6756 485.0 278.750 [5,] 189.0362 510.0 301.750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 161.8852 432.8788 234.7667 [2,] 177.3643 477.1212 271.4833 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(152.751941249712, 161.706815235661, 169.624722069753, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1cubi1281951901.ps tmp/1cubi1281951901.png",intern=TRUE)) character(0) > try(system("convert tmp/2cubi1281951901.ps tmp/2cubi1281951901.png",intern=TRUE)) character(0) > try(system("convert tmp/3mms31281951901.ps tmp/3mms31281951901.png",intern=TRUE)) character(0) > try(system("convert tmp/4mms31281951901.ps tmp/4mms31281951901.png",intern=TRUE)) character(0) > try(system("convert tmp/5mms31281951901.ps tmp/5mms31281951901.png",intern=TRUE)) character(0) > try(system("convert tmp/6fv9o1281951901.ps tmp/6fv9o1281951901.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.115 0.825 1.635