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(97.9,96.6,96.2,95.8,96,96.2,96.5,95.4,96.2,96.8,97,97.1,96.6,95.8,96.2,96,95.9,95.8,96,96.2,96.8,96.9,97.2,98.9,98.9,98.8,99,99.6,99.9,100.1,99.9,100.5,100.7,101.3,101.3,101.7,103,104.7,104.5,104.6,104.2,104.5,104.8,104.7,104.6,104.5,104.7,104.8,106.2,105.6,104.5,104.2,104.1,103.3,103.8,103.3,103.7,103.5,103.7,103.4,102.8) > 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] 61 > (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] 6 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 97.9 96.6 98.9 103.0 106.2 102.8 [2,] 96.6 95.8 98.8 104.7 105.6 NA [3,] 96.2 96.2 99.0 104.5 104.5 NA [4,] 95.8 96.0 99.6 104.6 104.2 NA [5,] 96.0 95.9 99.9 104.2 104.1 NA [6,] 96.2 95.8 100.1 104.5 103.3 NA [7,] 96.5 96.0 99.9 104.8 103.8 NA [8,] 95.4 96.2 100.5 104.7 103.3 NA [9,] 96.2 96.8 100.7 104.6 103.7 NA [10,] 96.8 96.9 101.3 104.5 103.5 NA [11,] 97.0 97.2 101.3 104.7 103.7 NA [12,] 97.1 98.9 101.7 104.8 103.4 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/1d98p1257943509.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/2xsjq1257943509.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/3p1pa1257943509.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/4jhtb1257943509.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,] 96.60 95.8 96.2 95.8 95.9 95.8 96.0 95.4 96.2 96.8 97.0 97.1 [2,] 97.90 96.6 96.2 96.0 96.0 96.2 96.5 96.2 96.8 96.9 97.2 98.9 [3,] 100.85 98.8 99.0 99.6 99.9 100.1 99.9 100.5 100.7 101.3 101.3 101.7 [4,] 103.00 104.7 104.5 104.2 104.1 103.3 103.8 103.3 103.7 103.5 103.7 103.4 [5,] 106.20 105.6 104.5 104.6 104.2 104.5 104.8 104.7 104.6 104.5 104.7 104.8 $n [1] 6 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 97.56034 93.07656 93.13524 93.8059 94.17656 95.08316 94.74184 [2,] 104.13966 104.52344 104.86476 105.3941 105.62344 105.11684 105.05816 [,8] [,9] [,10] [,11] [,12] [1,] 95.48316 95.82448 96.63646 96.70712 98.52031 [2,] 105.51684 105.57552 105.96354 105.89288 104.87969 $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(96.6, 97.9, 100.85, 103, 106.2, 95.8, 96.6, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5bfat1257943509.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,] 95.40 95.80 98.8 104.2 103.30 102.8 [2,] 96.10 95.95 99.3 104.5 103.45 102.8 [3,] 96.35 96.20 100.0 104.6 103.75 102.8 [4,] 96.90 96.85 101.0 104.7 104.35 102.8 [5,] 97.90 97.20 101.7 104.8 105.60 102.8 $n [1] 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 95.98511 95.7895 99.22462 104.5088 103.3395 102.8 [2,] 96.71489 96.6105 100.77538 104.6912 104.1605 102.8 $out [1] 98.9 103.0 106.2 $group [1] 2 4 5 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(95.4, 96.1, 96.35, 96.9, 97.9, 95.8, 95.95, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6t18c1257943509.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,] 3.168122 7.70 4.50 [2,] 3.645409 8.00 6.55 [3,] 4.011820 8.55 7.10 [4,] 4.170909 9.05 8.10 [5,] 4.572745 9.80 8.30 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 3.772136 8.071088 6.393035 [2,] 4.251505 9.028912 7.806965 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(3.16812247238013, 3.64540873104429, 4.01182036776363, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1d98p1257943509.ps tmp/1d98p1257943509.png") > system("convert tmp/2xsjq1257943509.ps tmp/2xsjq1257943509.png") > system("convert tmp/3p1pa1257943509.ps tmp/3p1pa1257943509.png") > system("convert tmp/4jhtb1257943509.ps tmp/4jhtb1257943509.png") > system("convert tmp/5bfat1257943509.ps tmp/5bfat1257943509.png") > system("convert tmp/6t18c1257943509.ps tmp/6t18c1257943509.png") > > > proc.time() user system elapsed 1.047 0.807 2.100