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(100,96.21064363,96.31280765,107.1793443,114.9066592,92.56060184,114.9995356,107.1236185,117.7765394,107.3650971,106.2970187,114.5072908,98.0031578,103.0649206,100.2879168,104.6066685,111.1544534,104.9874617,109.9284852,111.5352466,132.4974459,100.3436426,123.0983561,114.2379493,104.569518,109.0833101,106.9843039,133.6769759,124.8537197,122.5132349,116.8013374,116.0118882,129.7575926,125.1973623,143.7912139,127.9465032,130.2962757,108.4424631,129.3675118,143.6797622,131.8844618,117.6186496,118.9560695,104.8202842,134.624315,140.401226,143.8005015,153.4317823,153.2924677,127.3149438,153.5525216,136.9276493,131.7730101,144.3391845,107.4208229,113.6249652,124.2221603,102.0618557,96.36853348,111.6838488) > 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] 60 > (np <- floor(n / par1)) [1] 5 > 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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 100.00000 98.00316 104.5695 130.2963 153.29247 NA [2,] 96.21064 103.06492 109.0833 108.4425 127.31494 NA [3,] 96.31281 100.28792 106.9843 129.3675 153.55252 NA [4,] 107.17934 104.60667 133.6770 143.6798 136.92765 NA [5,] 114.90666 111.15445 124.8537 131.8845 131.77301 NA [6,] 92.56060 104.98746 122.5132 117.6186 144.33918 NA [7,] 114.99954 109.92849 116.8013 118.9561 107.42082 NA [8,] 107.12362 111.53525 116.0119 104.8203 113.62497 NA [9,] 117.77654 132.49745 129.7576 134.6243 124.22216 NA [10,] 107.36510 100.34364 125.1974 140.4012 102.06186 NA [11,] 106.29702 123.09836 143.7912 143.8005 96.36853 NA [12,] 114.50729 114.23795 127.9465 153.4318 111.68385 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/1fcmb1257794714.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/2tf9k1257794714.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/3eopk1257794714.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/4g2ea1257794714.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] [1,] 98.00316 96.21064 96.3128 104.6067 111.1545 92.5606 107.4208 104.8203 [2,] 100.00000 103.06492 100.2879 107.1793 114.9067 104.9875 109.9285 107.1236 [3,] 104.56952 108.44246 106.9843 133.6770 124.8537 117.6186 114.9995 111.5352 [4,] 130.29628 109.08331 129.3675 136.9276 131.7730 122.5132 116.8013 113.6250 [5,] 153.29247 109.08331 153.5525 143.6798 131.8845 144.3392 118.9561 116.0119 [,9] [,10] [,11] [,12] [1,] 117.7765 100.3436 96.36853 111.6838 [2,] 124.2222 102.0619 106.29702 114.2379 [3,] 129.7576 107.3651 123.09836 114.5073 [4,] 132.4974 125.1974 143.79121 127.9465 [5,] 134.6243 140.4012 143.80050 127.9465 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 83.16225 104.1899 86.43674 112.6569 112.9360 105.2350 110.1432 106.9414 [2,] 125.97679 112.6950 127.53187 154.6971 136.7714 130.0023 119.8559 116.1291 [,9] [,10] [,11] [,12] [1,] 123.9103 91.0176 96.60505 104.8209 [2,] 135.6049 123.7126 149.59166 124.1937 $out [1] 127.3149 153.4318 $group [1] 2 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(98.0031578, 100, 104.569518, 130.2962757, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/59lgc1257794714.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] [1,] 92.5606 98.00316 104.5695 104.8203 96.36853 NA [2,] 98.1564 101.70428 112.5476 118.2874 109.55234 NA [3,] 107.1515 107.45797 123.6835 131.0904 125.76855 NA [4,] 114.7070 112.88660 128.8520 142.0405 140.63342 NA [5,] 117.7765 123.09836 143.7912 153.4318 153.55252 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 99.60265 102.3576 116.2469 120.2564 111.5923 NA [2,] 114.70031 112.5583 131.1200 141.9243 139.9448 NA $out [1] 132.4974 $group [1] 2 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(92.56060184, 98.156403825, 107.1514814, 114.706975, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/65hqv1257794714.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,] 4.604381 11.19160 6.018390 [2,] 8.184152 18.78889 7.574069 [3,] 17.348778 39.56534 17.196062 [4,] 20.500980 49.60528 29.413950 [5,] 23.969049 57.23971 37.494195 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 11.73099 25.50978 7.234746 [2,] 22.96657 53.62090 27.157378 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(4.60438069370768, 8.184152177133, 17.3487777622828, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1fcmb1257794714.ps tmp/1fcmb1257794714.png") > system("convert tmp/2tf9k1257794714.ps tmp/2tf9k1257794714.png") > system("convert tmp/3eopk1257794714.ps tmp/3eopk1257794714.png") > system("convert tmp/4g2ea1257794714.ps tmp/4g2ea1257794714.png") > system("convert tmp/59lgc1257794714.ps tmp/59lgc1257794714.png") > system("convert tmp/65hqv1257794714.ps tmp/65hqv1257794714.png") > > > proc.time() user system elapsed 1.071 0.821 1.492