R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-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,557,561,549,532,526,511,499,555,565,542,527,510,514,517,508,493,490,469,478,528,534,518,506,502,516) > 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] 84 > (np <- floor(n / par1)) [1] 7 > 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] 7 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 476 519 562 591 593 561 517 NA [2,] 475 517 561 589 590 549 508 NA [3,] 470 510 555 584 580 532 493 NA [4,] 461 509 544 573 574 526 490 NA [5,] 455 501 537 567 573 511 469 NA [6,] 456 507 543 569 573 499 478 NA [7,] 517 569 594 621 620 555 528 NA [8,] 525 580 611 629 626 565 534 NA [9,] 523 578 613 628 620 542 518 NA [10,] 519 565 611 612 588 527 506 NA [11,] 509 547 594 595 566 510 502 NA [12,] 512 555 595 597 557 514 516 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/freestat/rcomp/tmp/1yncu1304594058.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/html/freestat/rcomp/tmp/2u9cu1304594058.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/html/freestat/rcomp/tmp/3qj3o1304594058.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/html/freestat/rcomp/tmp/4vhp61304594058.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.0 470.0 461.0 455 456.0 517.0 525.0 518.0 506.0 502.0 512 [2,] 518.0 512.5 501.5 499.5 485 488.5 541.5 549.5 532.5 523.0 509.5 515 [3,] 561.0 549.0 532.0 526.0 511 507.0 569.0 580.0 578.0 565.0 547.0 555 [4,] 576.5 575.0 567.5 558.5 552 556.0 607.0 618.5 616.5 599.5 580.0 576 [5,] 593.0 590.0 584.0 574.0 573 573.0 621.0 629.0 628.0 612.0 595.0 597 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 526.0647 511.676 492.5859 490.7662 470.9887 466.6901 529.8845 538.7943 [2,] 595.9353 586.324 571.4141 561.2338 551.0113 547.3099 608.1155 621.2057 [,9] [,10] [,11] [,12] [1,] 527.8366 519.3154 504.8985 518.5718 [2,] 628.1634 610.6846 589.1015 591.4282 $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, 518, 561, 576.5, 593, 475, 512.5, 549, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/52tuy1304594058.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] [1,] 455.0 501.0 537.0 567.0 557.0 499.0 469.0 NA [2,] 465.5 509.5 549.5 578.5 573.0 512.5 491.5 NA [3,] 492.5 533.0 578.0 593.0 584.0 529.5 507.0 NA [4,] 518.0 567.0 603.0 616.5 606.5 552.0 517.5 NA [5,] 525.0 580.0 613.0 629.0 626.0 565.0 534.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 468.5544 506.7739 553.5983 575.668 568.7204 511.4838 495.1412 NA [2,] 516.4456 559.2261 602.4017 610.332 599.2796 547.5162 518.8588 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" 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/html/freestat/rcomp/tmp/60yuu1304594058.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,] 40.22200 93.0 58.50 [2,] 41.92352 104.0 61.75 [3,] 43.17903 111.5 66.50 [4,] 44.68982 116.0 69.75 [5,] 47.13052 118.0 76.50 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 41.91730 106.0267 62.85115 [2,] 44.44076 116.9733 70.14885 $out [1] 36.97232 85.00000 84.00000 $group [1] 1 2 3 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(40.222002981074, 41.9235175270543, 43.1790289495123, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1yncu1304594058.ps tmp/1yncu1304594058.png",intern=TRUE)) character(0) > try(system("convert tmp/2u9cu1304594058.ps tmp/2u9cu1304594058.png",intern=TRUE)) character(0) > try(system("convert tmp/3qj3o1304594058.ps tmp/3qj3o1304594058.png",intern=TRUE)) character(0) > try(system("convert tmp/4vhp61304594058.ps tmp/4vhp61304594058.png",intern=TRUE)) character(0) > try(system("convert tmp/52tuy1304594058.ps tmp/52tuy1304594058.png",intern=TRUE)) character(0) > try(system("convert tmp/60yuu1304594058.ps tmp/60yuu1304594058.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.612 0.184 1.830