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(210,209,208,206,226,225,210,200,201,201,202,204,197,196,187,196,221,218,200,191,194,192,199,196,182,178,169,177,207,213,191,182,188,189,194,195,171,165,156,170,201,208,189,175,184,187,193,199,179,188,171,182,212,216,192,182,183,183,187,190,167,167,158,171,201,208,181,169,173,180,181,192,169,168,156,161,195,208,176,164,170,175,170,175,148,151,143,139,166,186,149,142,138,137,130,138,118,113,99,93,125,146,109,97,97,94,92,103,78,72,57,40,70,89,53,46,43,38,29,34) > 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,] 210 197 182 171 179 167 169 148 118 78 NA [2,] 209 196 178 165 188 167 168 151 113 72 NA [3,] 208 187 169 156 171 158 156 143 99 57 NA [4,] 206 196 177 170 182 171 161 139 93 40 NA [5,] 226 221 207 201 212 201 195 166 125 70 NA [6,] 225 218 213 208 216 208 208 186 146 89 NA [7,] 210 200 191 189 192 181 176 149 109 53 NA [8,] 200 191 182 175 182 169 164 142 97 46 NA [9,] 201 194 188 184 183 173 170 138 97 43 NA [10,] 201 192 189 187 183 180 175 137 94 38 NA [11,] 202 199 194 193 187 181 170 130 92 29 NA [12,] 204 196 195 199 190 192 175 138 103 34 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/1buyq1281709465.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/2buyq1281709465.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/343gt1281709465.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/443gt1281709465.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,] 118 113.0 143 93.0 125 146 109 97 97 94.0 92 103 [2,] 148 151.0 143 139.0 166 186 149 142 138 137.0 130 138 [3,] 170 167.5 157 170.5 201 208 185 172 178 181.5 184 191 [4,] 182 188.0 171 182.0 212 216 192 182 188 189.0 194 196 [5,] 210 209.0 208 206.0 226 225 210 200 201 201.0 202 204 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 153.0122 149.0133 143.0101 149.0155 178.0166 193.0108 163.5155 152.0144 [2,] 186.9878 185.9867 170.9899 191.9845 223.9834 222.9892 206.4845 191.9856 [,9] [,10] [,11] [,12] [1,] 153.018 155.5187 152.0230 162.0209 [2,] 202.982 207.4813 215.9770 219.9791 $out [1] 78 72 99 57 40 70 89 53 46 43 38 29 34 $group [1] 1 2 3 3 4 5 6 7 8 9 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(118, 148, 170, 182, 210, 113, 151, 167.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/543gt1281709465.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,] 200.0 187.0 169.0 156.0 171 158.0 156.0 130.0 92.0 29.0 NA [2,] 201.5 193.0 180.0 170.5 182 168.0 166.0 138.0 95.5 39.0 NA [3,] 207.0 196.0 188.5 185.5 185 176.5 170.0 142.5 101.0 49.5 NA [4,] 210.0 199.5 194.5 196.0 191 186.5 175.5 150.0 115.5 71.0 NA [5,] 210.0 200.0 213.0 208.0 192 208.0 176.0 166.0 125.0 89.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,] 203.1231 193.0353 181.8865 173.8693 180.8950 168.0620 165.667 137.0267 [2,] 210.8769 198.9647 195.1135 197.1307 189.1050 184.9380 174.333 147.9733 [,9] [,10] [,11] [1,] 91.87787 34.90459 NA [2,] 110.12213 64.09541 NA $out [1] 226 225 221 218 212 216 195 208 186 146 $group [1] 1 1 2 2 5 5 7 7 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(200, 201.5, 207, 210, 210, 187, 193, 196, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/67n0r1281709466.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,] 38.81709 132.0 23.750 [2,] 42.99564 144.0 29.750 [3,] 49.13902 156.5 36.125 [4,] 51.87468 164.5 41.500 [5,] 57.18595 173.0 53.750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 45.08923 147.1498 30.76575 [2,] 53.18882 165.8502 41.48425 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(38.8170924550856, 42.9956355242446, 49.1390233489582, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1buyq1281709465.ps tmp/1buyq1281709465.png",intern=TRUE)) character(0) > try(system("convert tmp/2buyq1281709465.ps tmp/2buyq1281709465.png",intern=TRUE)) character(0) > try(system("convert tmp/343gt1281709465.ps tmp/343gt1281709465.png",intern=TRUE)) character(0) > try(system("convert tmp/443gt1281709465.ps tmp/443gt1281709465.png",intern=TRUE)) character(0) > try(system("convert tmp/543gt1281709465.ps tmp/543gt1281709465.png",intern=TRUE)) character(0) > try(system("convert tmp/67n0r1281709466.ps tmp/67n0r1281709466.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.048 0.797 1.494