R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. 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(333,332,331,329,349,348,333,323,324,324,325,327,329,333,329,333,355,358,338,326,320,322,322,324,326,330,331,333,364,363,341,327,313,321,312,312,312,314,312,319,356,351,329,313,298,303,278,275,276,276,273,287,320,313,281,266,258,259,237,231,237,236,229,243,271,262,227,208,212,222,200,193,204,203,190,209,240,234,210,195,202,204,180,169,178,181,163,174,194,187,160,143,151,154,141,127,134,138,120,129,151,152,124,99,104,109,96,87,94,89,63,76,100,104,80,55,60,71,62,61) > 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,] 333 329 326 312 276 237 204 178 134 94 NA [2,] 332 333 330 314 276 236 203 181 138 89 NA [3,] 331 329 331 312 273 229 190 163 120 63 NA [4,] 329 333 333 319 287 243 209 174 129 76 NA [5,] 349 355 364 356 320 271 240 194 151 100 NA [6,] 348 358 363 351 313 262 234 187 152 104 NA [7,] 333 338 341 329 281 227 210 160 124 80 NA [8,] 323 326 327 313 266 208 195 143 99 55 NA [9,] 324 320 313 298 258 212 202 151 104 60 NA [10,] 324 322 321 303 259 222 204 154 109 71 NA [11,] 325 322 312 278 237 200 180 141 96 62 NA [12,] 327 324 312 275 231 193 169 127 87 61 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/15om71281956430.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/freestat/rcomp/tmp/25om71281956430.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/freestat/rcomp/tmp/35om71281956430.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/freestat/rcomp/tmp/4gyma1281956430.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,] 94.0 89 63 76 100.0 104.0 80 55 60 71.0 62.0 61 [2,] 178.0 181 163 174 194.0 187.0 160 143 151 154.0 141.0 127 [3,] 256.5 256 251 265 295.5 287.5 254 237 235 240.5 218.5 212 [4,] 326.0 330 329 329 355.0 351.0 333 323 313 321.0 312.0 312 [5,] 333.0 333 331 333 364.0 363.0 341 327 324 324.0 325.0 327 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 182.5533 181.5537 168.0598 187.5558 215.058 205.5591 167.5623 147.0648 [2,] 330.4467 330.4463 333.9402 342.4442 375.942 369.4409 340.4377 326.9352 [,9] [,10] [,11] [,12] [1,] 154.0583 157.0601 133.0616 119.5666 [2,] 315.9417 323.9399 303.9384 304.4334 $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(94, 178, 256.5, 326, 333, 89, 181, 256, 330, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5gyma1281956430.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,] 323.0 320.0 312.0 275.0 231.0 193 169.0 127.0 87.0 55.0 NA [2,] 324.5 323.0 317.0 300.5 258.5 210 192.5 147.0 101.5 61.5 NA [3,] 330.0 329.0 328.5 312.5 274.5 228 203.5 161.5 122.0 73.5 NA [4,] 333.0 335.5 337.0 324.0 284.0 240 209.5 179.5 136.0 91.5 NA [5,] 333.0 338.0 364.0 356.0 320.0 271 234.0 194.0 152.0 104.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,] 326.1231 323.2987 319.3779 301.7815 262.8693 214.3168 195.7462 146.6765 [2,] 333.8769 334.7013 337.6221 323.2185 286.1307 241.6832 211.2538 176.3235 [,9] [,10] [,11] [1,] 106.2643 59.8168 NA [2,] 137.7357 87.1832 NA $out [1] 349 348 355 358 240 $group [1] 1 1 2 2 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(323, 324.5, 330, 333, 333, 320, 323, 329, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6gyma1281956430.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,] 88.16248 244 138.000 [2,] 93.18038 255 144.625 [3,] 95.13923 262 150.750 [4,] 96.95788 265 157.250 [5,] 101.28096 272 165.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 93.41628 257.4389 144.9917 [2,] 96.86217 266.5611 156.5083 $out [1] 86.98282 239.00000 $group [1] 1 2 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(88.1624762709296, 93.180378464695, 95.139225334457, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/15om71281956430.ps tmp/15om71281956430.png",intern=TRUE)) character(0) > try(system("convert tmp/25om71281956430.ps tmp/25om71281956430.png",intern=TRUE)) character(0) > try(system("convert tmp/35om71281956430.ps tmp/35om71281956430.png",intern=TRUE)) character(0) > try(system("convert tmp/4gyma1281956430.ps tmp/4gyma1281956430.png",intern=TRUE)) character(0) > try(system("convert tmp/5gyma1281956430.ps tmp/5gyma1281956430.png",intern=TRUE)) character(0) > try(system("convert tmp/6gyma1281956430.ps tmp/6gyma1281956430.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.740 1.280 2.044