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(130,127,122,117,112,113,149,157,157,147,137,132,125,123,117,114,111,112,144,150,149,134,123,116,117,111,105,102,95,93,124,130,124,115,106,105,105,101,95,93,84,87,116,120,117,109,105,107,109,109,108,107,99,103,131,137,135,124,118,121,121,118,113,107,100,102,130,136,133,120,112,109,110) > 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] 73 > (np <- floor(n / par1)) [1] 6 > 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 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 130 125 117 105 109 121 110 [2,] 127 123 111 101 109 118 NA [3,] 122 117 105 95 108 113 NA [4,] 117 114 102 93 107 107 NA [5,] 112 111 95 84 99 100 NA [6,] 113 112 93 87 103 102 NA [7,] 149 144 124 116 131 130 NA [8,] 157 150 130 120 137 136 NA [9,] 157 149 124 117 135 133 NA [10,] 147 134 115 109 124 120 NA [11,] 137 123 106 105 118 112 NA [12,] 132 116 105 107 121 109 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/1gx3r1305553135.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/226lz1305553135.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/3yg7n1305553135.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/4qk5j1305553135.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,] 105.0 101.0 95.0 93 84.0 87.0 116.0 120.0 117 109 105 105.0 [2,] 109.5 109.0 105.0 102 95.0 93.0 124.0 130.0 124 115 106 107.0 [3,] 117.0 114.5 110.5 107 99.5 102.5 130.5 136.5 134 122 115 112.5 [4,] 123.0 123.0 117.0 114 111.0 112.0 144.0 150.0 149 134 123 121.0 [5,] 130.0 127.0 122.0 117 112.0 113.0 149.0 157.0 157 147 137 132.0 $n [1] 7 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 108.9380 105.4695 102.7596 99.25961 89.17948 90.24439 117.5994 123.5994 [2,] 125.0620 123.5305 118.2404 114.74039 109.82052 114.75561 143.4006 149.4006 [,9] [,10] [,11] [,12] [1,] 117.8742 109.7444 104.0345 103.4695 [2,] 150.1258 134.2556 125.9655 121.5305 $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(105, 109.5, 117, 123, 130, 101, 109, 114.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5v7ri1305553135.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] [1,] 112.0 111 93.0 84.0 99.0 100.0 110 [2,] 119.5 115 103.5 94.0 107.5 108.0 110 [3,] 131.0 123 108.5 105.0 113.5 115.5 110 [4,] 148.0 139 120.5 112.5 127.5 125.5 110 [5,] 157.0 150 130.0 120.0 137.0 136.0 110 $n [1] 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 118.0010 112.0534 100.7462 96.56203 104.3779 107.5181 110 [2,] 143.9990 133.9466 116.2538 113.43797 122.6221 123.4819 110 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(112, 119.5, 131, 148, 157, 111, 115, 123, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/634jj1305553135.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,] 8.594572 24.0 9.000 [2,] 9.595896 26.0 12.250 [3,] 10.363774 27.5 13.875 [4,] 12.851848 35.0 15.250 [5,] 14.998889 40.0 19.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 8.878712 23.39504 12.50668 [2,] 11.848836 31.60496 15.24332 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(8.5945719303911, 9.59589573116497, 10.3637742185841, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1gx3r1305553135.ps tmp/1gx3r1305553135.png",intern=TRUE)) character(0) > try(system("convert tmp/226lz1305553135.ps tmp/226lz1305553135.png",intern=TRUE)) character(0) > try(system("convert tmp/3yg7n1305553135.ps tmp/3yg7n1305553135.png",intern=TRUE)) character(0) > try(system("convert tmp/4qk5j1305553135.ps tmp/4qk5j1305553135.png",intern=TRUE)) character(0) > try(system("convert tmp/5v7ri1305553135.ps tmp/5v7ri1305553135.png",intern=TRUE)) character(0) > try(system("convert tmp/634jj1305553135.ps tmp/634jj1305553135.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.627 0.198 1.945