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(244,243,242,240,260,259,244,234,235,235,236,238,240,233,233,233,247,251,233,226,233,233,227,225,223,212,206,202,223,221,212,205,204,200,195,193,196,180,174,164,192,189,181,167,166,164,167,164,161,141,134,111,147,144,142,140,143,137,140,130,129,112,101,74,104,103,100,98,99,91,92,94,93,76,64,32,62,69,69,68,68,59,66,73,70,57,48,22,64,74,67,61,61,52,54,69,69,53,50,22,69,78,74,63,67,59,60,80,77,58,54,32,78,86,84,78,72,64,62,72) > 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,] 244 240 223 196 161 129 93 70 69 77 NA [2,] 243 233 212 180 141 112 76 57 53 58 NA [3,] 242 233 206 174 134 101 64 48 50 54 NA [4,] 240 233 202 164 111 74 32 22 22 32 NA [5,] 260 247 223 192 147 104 62 64 69 78 NA [6,] 259 251 221 189 144 103 69 74 78 86 NA [7,] 244 233 212 181 142 100 69 67 74 84 NA [8,] 234 226 205 167 140 98 68 61 63 78 NA [9,] 235 233 204 166 143 99 68 61 67 72 NA [10,] 235 233 200 164 137 91 59 52 59 64 NA [11,] 236 227 195 167 140 92 66 54 60 62 NA [12,] 238 225 193 164 130 94 73 69 80 72 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/19w5h1281904795.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/29w5h1281904795.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/3wkxv1281904795.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/4wkxv1281904795.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,] 69 53.0 48.0 22.0 62.0 69.0 67 61 61 52 54 69 [2,] 77 58.0 54.0 32.0 69.0 78.0 74 68 68 59 62 73 [3,] 145 126.5 117.5 92.5 125.5 123.5 121 119 121 114 116 112 [4,] 223 212.0 206.0 202.0 223.0 221.0 212 205 204 200 195 193 [5,] 244 243.0 242.0 240.0 260.0 259.0 244 234 235 235 236 238 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 72.05258 49.55546 41.55474 7.561222 48.55546 52.0515 52.0497 [2,] 217.94742 203.44454 193.44526 177.438778 202.44454 194.9485 189.9503 [,8] [,9] [,10] [,11] [,12] [1,] 50.54934 53.04898 43.55078 49.5479 52.04322 [2,] 187.45066 188.95102 184.44922 182.4521 171.95678 $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(69, 77, 145, 223, 244, 53, 58, 126.5, 212, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5wkxv1281904795.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,] 234.0 225.0 193.0 164.0 130.0 91.0 59 48 50.0 54 NA [2,] 235.5 230.0 201.0 165.0 135.5 93.0 63 53 56.0 60 NA [3,] 241.0 233.0 205.5 170.5 140.5 99.5 68 61 65.0 72 NA [4,] 244.0 236.5 216.5 185.0 143.5 103.5 71 68 71.5 78 NA [5,] 244.0 240.0 223.0 196.0 147.0 112.0 76 74 80.0 86 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 237.1231 230.0353 198.4303 161.3779 136.8511 94.71088 64.35115 54.1584 [2,] 244.8769 235.9647 212.5697 179.6221 144.1489 104.28912 71.64885 67.8416 [,9] [,10] [,11] [1,] 57.93035 63.79008 NA [2,] 72.06965 80.20992 NA $out [1] 260 259 247 251 161 111 129 74 93 32 22 22 32 $group [1] 1 1 2 2 5 5 6 6 7 7 8 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(234, 235.5, 241, 244, 244, 225, 230, 233, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6pbey1281904795.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,] 66.49946 169.0 111.000 [2,] 71.11617 174.5 125.250 [3,] 73.38630 182.5 131.875 [4,] 77.24342 192.0 141.500 [5,] 79.80560 218.0 160.500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 70.59162 174.5181 124.4633 [2,] 76.18098 190.4819 139.2867 $out [1] 89.60382 $group [1] 1 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(66.4994569735556, 71.1161697889638, 73.3862964871296, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/19w5h1281904795.ps tmp/19w5h1281904795.png",intern=TRUE)) character(0) > try(system("convert tmp/29w5h1281904795.ps tmp/29w5h1281904795.png",intern=TRUE)) character(0) > try(system("convert tmp/3wkxv1281904795.ps tmp/3wkxv1281904795.png",intern=TRUE)) character(0) > try(system("convert tmp/4wkxv1281904795.ps tmp/4wkxv1281904795.png",intern=TRUE)) character(0) > try(system("convert tmp/5wkxv1281904795.ps tmp/5wkxv1281904795.png",intern=TRUE)) character(0) > try(system("convert tmp/6pbey1281904795.ps tmp/6pbey1281904795.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.061 0.795 1.405