R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-redhat-linux-gnu (64-bit) 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(476,475,470,461,455,456,517,525,523,519,509,512,519,517,510,509,501,507,569,580,578,565,547,555,562,561,555,544,537,543,594,611,613,611,594,595,591,589,584,573,567,569,621,629,628,612,595,597,593,590,580,574,573,573,620,626,620,588,566,577,561,549,532,526,511,499,555,565,542,527,510,514,517,508,493,490,469,478,528,534,518,506,502,516,528,533,536,537,524,536,587,597,581,564,558,575,580,575,563,552,537,545,601,604,586,564,549,551,556) > 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] 109 > (np <- floor(n / par1)) [1] 9 > 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 9 9 9 9 9 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 476 519 562 591 593 561 517 528 580 556 [2,] 475 517 561 589 590 549 508 533 575 NA [3,] 470 510 555 584 580 532 493 536 563 NA [4,] 461 509 544 573 574 526 490 537 552 NA [5,] 455 501 537 567 573 511 469 524 537 NA [6,] 456 507 543 569 573 499 478 536 545 NA [7,] 517 569 594 621 620 555 528 587 601 NA [8,] 525 580 611 629 626 565 534 597 604 NA [9,] 523 578 613 628 620 542 518 581 586 NA [10,] 519 565 611 612 588 527 506 564 564 NA [11,] 509 547 594 595 566 510 502 558 549 NA [12,] 512 555 595 597 577 514 516 575 551 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/wessaorg/rcomp/tmp/1grdt1305035791.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2t5g81305035791.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3k8jp1305035791.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/436bq1305035791.ps",horizontal=F,onefile=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,] 476.0 475 470 461 455 456 517 525 518 506 502 512 [2,] 519.0 517 510 509 501 499 555 565 542 527 510 516 [3,] 558.5 549 536 537 524 536 587 597 581 564 549 555 [4,] 580.0 575 563 552 537 545 601 611 613 588 566 577 [5,] 593.0 590 584 574 573 573 621 629 628 612 595 597 $n [1] 10 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 528.022 518.4533 508.0867 514.3533 505.04 511.7733 562.7733 572.7733 [2,] 588.978 579.5467 563.9133 559.6467 542.96 560.2267 611.2267 621.2267 [,9] [,10] [,11] [,12] [1,] 543.6067 531.8733 519.5067 522.8733 [2,] 618.3933 596.1267 578.4933 587.1267 $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(476, 519, 558.5, 580, 593, 475, 517, 549, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/5pizc1305035791.ps",horizontal=F,onefile=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] [1,] 455.0 501.0 537.0 567.0 566.0 499.0 469.0 524.0 537.0 556 [2,] 465.5 509.5 549.5 578.5 573.5 512.5 491.5 534.5 550.0 556 [3,] 492.5 533.0 578.0 593.0 584.0 529.5 507.0 547.5 563.5 556 [4,] 518.0 567.0 603.0 616.5 606.5 552.0 517.5 578.0 583.0 556 [5,] 525.0 580.0 613.0 629.0 626.0 565.0 534.0 597.0 604.0 556 $n [1] 12 12 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 468.5544 506.7739 553.5983 575.668 568.9485 511.4838 495.1412 527.6594 [2,] 516.4456 559.2261 602.4017 610.332 599.0515 547.5162 518.8588 567.3406 [,9] [,10] [1,] 548.4485 556 [2,] 578.5515 556 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(455, 465.5, 492.5, 518, 525, 501, 509.5, 533, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/64i7c1305035791.ps",horizontal=F,onefile=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,] 35.05630 93.0 36.000 [2,] 37.58202 104.0 46.000 [3,] 38.09689 111.5 53.625 [4,] 39.68300 116.0 59.500 [5,] 41.02472 118.0 71.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 37.13862 106.0267 47.46756 [2,] 39.05517 116.9733 59.78244 $out [1] 34.12111 85.00000 $group [1] 1 2 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(35.0563039187597, 37.5820214057018, 38.0968927259926, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1grdt1305035791.ps tmp/1grdt1305035791.png",intern=TRUE)) character(0) > try(system("convert tmp/2t5g81305035791.ps tmp/2t5g81305035791.png",intern=TRUE)) character(0) > try(system("convert tmp/3k8jp1305035791.ps tmp/3k8jp1305035791.png",intern=TRUE)) character(0) > try(system("convert tmp/436bq1305035791.ps tmp/436bq1305035791.png",intern=TRUE)) character(0) > try(system("convert tmp/5pizc1305035791.ps tmp/5pizc1305035791.png",intern=TRUE)) character(0) > try(system("convert tmp/64i7c1305035791.ps tmp/64i7c1305035791.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 3.860 0.550 4.544