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(132.1,125,127.1,101.5,85.7,79.3,70.9,77.1,83.9,96.2,111.7,127.2,143.6,134.9,135.6,105.3,86.4,74.6,67.6,73.4,78.5,98.2,118.6,136.9,137.9,115.6,119.3,98.5,84.3,73.5,60.7,69.5,77.9,113.9,126.3,135.1,130.5,113.1,110,90.8,85.4,72.5,64.7,67.2,77.9,105.2,107.2,120.7,121.3,107.9,105.6,81.3,71.7,64.8,57.3,61.9,70.1,88.8,106.8,110.7,114.1,108,111.5,86.8,78.4,68,57.3,65.3,73.3,88.6,101.3,122.9,126.6,114.1,124.7) > 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] 75 > (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] 7 7 7 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 132.1 143.6 137.9 130.5 121.3 114.1 126.6 [2,] 125.0 134.9 115.6 113.1 107.9 108.0 114.1 [3,] 127.1 135.6 119.3 110.0 105.6 111.5 124.7 [4,] 101.5 105.3 98.5 90.8 81.3 86.8 NA [5,] 85.7 86.4 84.3 85.4 71.7 78.4 NA [6,] 79.3 74.6 73.5 72.5 64.8 68.0 NA [7,] 70.9 67.6 60.7 64.7 57.3 57.3 NA [8,] 77.1 73.4 69.5 67.2 61.9 65.3 NA [9,] 83.9 78.5 77.9 77.9 70.1 73.3 NA [10,] 96.2 98.2 113.9 105.2 88.8 88.6 NA [11,] 111.7 118.6 126.3 107.2 106.8 101.3 NA [12,] 127.2 136.9 135.1 120.7 110.7 122.9 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/1c59y1291149020.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/2nerj1291149020.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/3nerj1291149020.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/4y5841291149020.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,] 114.10 107.90 105.60 81.30 71.70 64.8 57.3 61.90 70.1 88.6 101.30 110.70 [2,] 123.95 110.55 110.75 86.80 78.40 68.0 57.3 65.30 73.3 88.8 106.80 120.70 [3,] 130.50 114.10 119.30 94.65 84.85 73.0 62.7 68.35 77.9 97.2 109.45 125.05 [4,] 135.00 120.30 125.90 101.50 85.70 74.6 67.6 73.40 78.5 105.2 118.60 135.10 [5,] 143.60 134.90 135.60 105.30 86.40 79.3 70.9 77.10 83.9 113.9 126.30 136.90 $n [1] 7 7 7 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 123.9011 108.2775 110.2527 85.16803 80.14126 68.74279 56.05617 63.12524 [2,] 137.0989 119.9225 128.3473 104.13197 89.55874 77.25721 69.34383 73.57476 [,9] [,10] [,11] [,12] [1,] 74.54583 86.62147 101.8386 115.7615 [2,] 81.25417 107.77853 117.0614 134.3385 $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(114.1, 123.95, 130.5, 135, 143.6, 107.9, 110.55, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5y5841291149020.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,] 70.90 67.60 60.7 64.70 57.30 57.30 114.10 [2,] 81.60 76.55 75.7 75.20 67.45 70.65 119.40 [3,] 98.85 101.75 106.2 98.00 85.05 87.70 124.70 [4,] 126.05 135.25 122.8 111.55 107.35 109.75 125.65 [5,] 132.10 143.60 137.9 130.50 121.30 122.90 126.60 $n [1] 12 12 12 12 12 12 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 78.57606 74.97654 84.71737 81.42052 66.85134 69.86623 118.9987 [2,] 119.12394 128.52346 127.68263 114.57948 103.24866 105.53377 130.4013 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(70.9, 81.6, 98.85, 126.05, 132.1, 67.6, 76.55, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6y5841291149020.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.748754 13.6 3.9000 [2,] 5.560645 14.6 6.2000 [3,] 9.166221 24.5 9.8625 [4,] 9.786089 26.6 12.3375 [5,] 10.727446 30.0 15.1500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 7.238967 19.02672 7.063145 [2,] 11.093475 29.97328 12.661855 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(4.74875422260057, 5.56064482224489, 9.16622111050557, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1c59y1291149020.ps tmp/1c59y1291149020.png",intern=TRUE)) character(0) > try(system("convert tmp/2nerj1291149020.ps tmp/2nerj1291149020.png",intern=TRUE)) character(0) > try(system("convert tmp/3nerj1291149020.ps tmp/3nerj1291149020.png",intern=TRUE)) character(0) > try(system("convert tmp/4y5841291149020.ps tmp/4y5841291149020.png",intern=TRUE)) character(0) > try(system("convert tmp/5y5841291149020.ps tmp/5y5841291149020.png",intern=TRUE)) character(0) > try(system("convert tmp/6y5841291149020.ps tmp/6y5841291149020.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.690 1.239 1.911