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(84.9,81.9,95.9,81,89.2,102.5,89.8,88.8,83.2,90.2,100.4,187.1,87.6,85.4,86.1,86.7,89.1,103.7,86.9,85.2,80.8,91.2,102.8,182.5,80.9,83.1,88.3,86.6,93,105.3,93.8,86.4,87,96.7,100.5,196.7,86.8,88.2,93.8,85,90.4,115.9,94.9,87.7,91.7,95.9,106.8,204.5,90.2,90.5,93.2,97.8,99.4,120,108.2,98.5,104.3,102.9,111.1,188.1,93.8,94.5,112.4,102.5,115.8,136.5,122.1,110.6,116.4,112.6,121.5,199.3) > 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] 72 > (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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 84.9 87.6 80.9 86.8 90.2 93.8 NA [2,] 81.9 85.4 83.1 88.2 90.5 94.5 NA [3,] 95.9 86.1 88.3 93.8 93.2 112.4 NA [4,] 81.0 86.7 86.6 85.0 97.8 102.5 NA [5,] 89.2 89.1 93.0 90.4 99.4 115.8 NA [6,] 102.5 103.7 105.3 115.9 120.0 136.5 NA [7,] 89.8 86.9 93.8 94.9 108.2 122.1 NA [8,] 88.8 85.2 86.4 87.7 98.5 110.6 NA [9,] 83.2 80.8 87.0 91.7 104.3 116.4 NA [10,] 90.2 91.2 96.7 95.9 102.9 112.6 NA [11,] 100.4 102.8 100.5 106.8 111.1 121.5 NA [12,] 187.1 182.5 196.7 204.5 188.1 199.3 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/182oc1291665697.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/282oc1291665697.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/382oc1291665697.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/4jbof1291665697.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,] 80.9 81.9 86.1 81.00 89.1 102.5 86.90 85.20 80.80 90.2 100.4 182.5 [2,] 84.9 83.1 88.3 85.00 89.2 103.7 89.80 86.40 83.20 91.2 100.5 187.1 [3,] 87.2 86.8 93.5 86.65 91.7 110.6 94.35 88.25 89.35 96.3 104.8 192.4 [4,] 90.2 90.5 95.9 97.80 99.4 120.0 108.20 98.50 104.30 102.9 111.1 199.3 [5,] 93.8 94.5 95.9 102.50 99.4 136.5 122.10 110.60 116.40 112.6 121.5 204.5 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 83.78133 82.02676 88.59775 78.39359 85.12067 100.0860 82.4814 80.44511 [2,] 90.61867 91.57324 98.40225 94.90641 98.27933 121.1140 106.2186 96.05489 [,9] [,10] [,11] [,12] [1,] 75.73982 88.75312 97.96266 184.5306 [2,] 102.96018 103.84688 111.63734 200.2694 $out [1] 112.4 115.8 $group [1] 3 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(80.9, 84.9, 87.2, 90.2, 93.8, 81.9, 83.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5jbof1291665697.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,] 81.00 80.80 80.90 85.00 90.20 93.80 NA [2,] 84.05 85.75 86.50 87.95 95.50 106.55 NA [3,] 89.50 87.25 90.65 92.75 101.15 114.20 NA [4,] 98.15 97.00 98.60 101.35 109.65 121.80 NA [5,] 102.50 103.70 105.30 115.90 120.00 136.50 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 83.0689 82.1188 85.13111 86.63817 94.69609 107.2444 NA [2,] 95.9311 92.3812 96.16889 98.86183 107.60391 121.1556 NA $out [1] 187.1 182.5 196.7 204.5 188.1 199.3 $group [1] 1 2 3 4 5 6 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(81, 84.05, 89.5, 98.15, 102.5, 80.8, 85.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6bkn01291665697.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,] 4.420256 12.60 4.1750 [2,] 8.220563 21.30 7.2750 [3,] 8.860636 23.90 9.1625 [4,] 11.741967 30.35 12.6875 [5,] 13.805506 35.60 17.0000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 7.25450 19.77223 6.693822 [2,] 10.46677 28.02777 11.631178 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(4.42025640281948, 8.22056316683924, 8.8606364408691, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/182oc1291665697.ps tmp/182oc1291665697.png",intern=TRUE)) character(0) > try(system("convert tmp/282oc1291665697.ps tmp/282oc1291665697.png",intern=TRUE)) character(0) > try(system("convert tmp/382oc1291665697.ps tmp/382oc1291665697.png",intern=TRUE)) character(0) > try(system("convert tmp/4jbof1291665697.ps tmp/4jbof1291665697.png",intern=TRUE)) character(0) > try(system("convert tmp/5jbof1291665697.ps tmp/5jbof1291665697.png",intern=TRUE)) character(0) > try(system("convert tmp/6bkn01291665697.ps tmp/6bkn01291665697.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.674 1.221 1.860