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.

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(153.3
+ ,154.5
+ ,155.2
+ ,156.9
+ ,157
+ ,157.4
+ ,157.2
+ ,157.5
+ ,158
+ ,158.5
+ ,159
+ ,159.3
+ ,160
+ ,160.8
+ ,161.9
+ ,162.5
+ ,162.7
+ ,162.8
+ ,162.9
+ ,163
+ ,164
+ ,164.7
+ ,164.8
+ ,164.9
+ ,165
+ ,165.8
+ ,166.1
+ ,167.2
+ ,167.7
+ ,168.3
+ ,168.6
+ ,168.9
+ ,169.1
+ ,169.5
+ ,169.6
+ ,169.7
+ ,169.8
+ ,170.4
+ ,170.9
+ ,171.9
+ ,171.9
+ ,172
+ ,172
+ ,172.4
+ ,173
+ ,173.7
+ ,173.8
+ ,173.8
+ ,173.9
+ ,174.6
+ ,175
+ ,175.9
+ ,176
+ ,175.1
+ ,175.6
+ ,175.9
+ ,176.7
+ ,176.1
+ ,176.1
+ ,176.2
+ ,176.3
+ ,177.8
+ ,178.5
+ ,179.4
+ ,179.5
+ ,179.6
+ ,179.7
+ ,179.7
+ ,179.8
+ ,179.9
+ ,180.2
+ ,180.4
+ ,180.4
+ ,181.3
+ ,181.9
+ ,182.5
+ ,182.7
+ ,183.1
+ ,183.6
+ ,183.7
+ ,183.8
+ ,183.9
+ ,184.1
+ ,184.4
+ ,184.5
+ ,185.9
+ ,186.6
+ ,187.6
+ ,187.8
+ ,187.9
+ ,188
+ ,188.3
+ ,188.4
+ ,188.5
+ ,188.5
+ ,188.6
+ ,188.6
+ ,189.4
+ ,190
+ ,191.9
+ ,192.5
+ ,193
+ ,193.5
+ ,193.9
+ ,194.2
+ ,194.9
+ ,194.9
+ ,194.9
+ ,194.9
+ ,195.5
+ ,196
+ ,196.2
+ ,196.2
+ ,196.2
+ ,196.2
+ ,197
+ ,197.7
+ ,198
+ ,198.2
+ ,198.5
+ ,198.6
+ ,199.5
+ ,200
+ ,201.3
+ ,202.2
+ ,202.9
+ ,203.5
+ ,203.5
+ ,204
+ ,204.1
+ ,204.3
+ ,204.5
+ ,204.8
+ ,205.1
+ ,205.7
+ ,206.5
+ ,206.9
+ ,207.1
+ ,207.8
+ ,208
+ ,208.5
+ ,208.6
+ ,209
+ ,209.1
+ ,209.7
+ ,209.8
+ ,209.9
+ ,210
+ ,210.8
+ ,211.4
+ ,211.7
+ ,212
+ ,212.2
+ ,212.4
+ ,212.9
+ ,213.4
+ ,213.7
+ ,214
+ ,214.3
+ ,214.8
+ ,215
+ ,215.9
+ ,216.4
+ ,216.9
+ ,217.2
+ ,217.5
+ ,217.9
+ ,218.1
+ ,218.6
+ ,218.9
+ ,219.3
+ ,220.4
+ ,220.9
+ ,221
+ ,221.8
+ ,222
+ ,222.2
+ ,222.5
+ ,222.9
+ ,223.1
+ ,223.4
+ ,224
+ ,225.1
+ ,225.5
+ ,225.9
+ ,226.3
+ ,226.5
+ ,227
+ ,227.3
+ ,227.8
+ ,228.1
+ ,228.4
+ ,228.5
+ ,228.8
+ ,229
+ ,229.1
+ ,229.3
+ ,229.6
+ ,229.9
+ ,230
+ ,230.2
+ ,230.8
+ ,231
+ ,231.7
+ ,231.9
+ ,233
+ ,235.1
+ ,236
+ ,236.9
+ ,237.1
+ ,237.5
+ ,238.2
+ ,238.9
+ ,239.1
+ ,240
+ ,240.2
+ ,240.5
+ ,240.7
+ ,241.1
+ ,241.4
+ ,242.2
+ ,242.9
+ ,243.2
+ ,243.9)
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = '12'
> ylab = ''
> xlab = ''
> main = ''
> #'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] 224
> (np <- floor(n / par1))
[1] 18
> 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] 19 19 19 19 19 19 19 19 18 18 18 18
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 153.3 160.0 165.0 169.8 173.9 176.3 180.4 184.5 188.6 194.9 198.6 204.8
 [2,] 154.5 160.8 165.8 170.4 174.6 177.8 181.3 185.9 189.4 195.5 199.5 205.1
 [3,] 155.2 161.9 166.1 170.9 175.0 178.5 181.9 186.6 190.0 196.0 200.0 205.7
 [4,] 156.9 162.5 167.2 171.9 175.9 179.4 182.5 187.6 191.9 196.2 201.3 206.5
 [5,] 157.0 162.7 167.7 171.9 176.0 179.5 182.7 187.8 192.5 196.2 202.2 206.9
 [6,] 157.4 162.8 168.3 172.0 175.1 179.6 183.1 187.9 193.0 196.2 202.9 207.1
 [7,] 157.2 162.9 168.6 172.0 175.6 179.7 183.6 188.0 193.5 196.2 203.5 207.8
 [8,] 157.5 163.0 168.9 172.4 175.9 179.7 183.7 188.3 193.9 197.0 203.5 208.0
 [9,] 158.0 164.0 169.1 173.0 176.7 179.8 183.8 188.4 194.2 197.7 204.0 208.5
[10,] 158.5 164.7 169.5 173.7 176.1 179.9 183.9 188.5 194.9 198.0 204.1 208.6
[11,] 159.0 164.8 169.6 173.8 176.1 180.2 184.1 188.5 194.9 198.2 204.3 209.0
[12,] 159.3 164.9 169.7 173.8 176.2 180.4 184.4 188.6 194.9 198.5 204.5 209.1
      [,13] [,14] [,15] [,16] [,17] [,18] [,19]
 [1,] 209.7 213.7 218.6 223.4 228.5 231.9 240.5
 [2,] 209.8 214.0 218.9 224.0 228.8 233.0 240.7
 [3,] 209.9 214.3 219.3 225.1 229.0 235.1 241.1
 [4,] 210.0 214.8 220.4 225.5 229.1 236.0 241.4
 [5,] 210.8 215.0 220.9 225.9 229.3 236.9 242.2
 [6,] 211.4 215.9 221.0 226.3 229.6 237.1 242.9
 [7,] 211.7 216.4 221.8 226.5 229.9 237.5 243.2
 [8,] 212.0 216.9 222.0 227.0 230.0 238.2 243.9
 [9,] 212.2 217.2 222.2 227.3 230.2 238.9    NA
[10,] 212.4 217.5 222.5 227.8 230.8 239.1    NA
[11,] 212.9 217.9 222.9 228.1 231.0 240.0    NA
[12,] 213.4 218.1 223.1 228.4 231.7 240.2    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/12l1p1257935697.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/freestat/rcomp/tmp/2r0x01257935697.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/freestat/rcomp/tmp/3pvw01257935697.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/freestat/rcomp/tmp/4fz961257935697.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]
[1,] 153.30 154.50 155.20 156.90 157.00 157.40 157.20 157.50 158.00 158.50
[2,] 175.10 176.20 176.75 177.65 177.75 177.35 177.65 177.80 176.70 176.10
[3,] 194.90 195.50 196.00 196.20 196.20 196.20 196.20 197.00 195.95 196.45
[4,] 216.15 216.45 216.80 217.60 217.95 218.45 219.10 219.45 217.20 217.50
[5,] 240.50 240.70 241.10 241.40 242.20 242.90 243.20 243.90 238.90 239.10
      [,11] [,12]
[1,] 159.00 159.3
[2,] 176.10 176.2
[3,] 196.55 196.7
[4,] 217.90 218.1
[5,] 240.00 240.2

$n
 [1] 19 19 19 19 19 19 19 19 18 18 18 18

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 180.0203 180.9103 181.4828 181.7190 181.6284 181.3022 181.1753 181.9028
[2,] 209.7797 210.0897 210.5172 210.6810 210.7716 211.0978 211.2247 212.0972
         [,9]    [,10]    [,11]    [,12]
[1,] 180.8674 181.0322 180.9833 181.0960
[2,] 211.0326 211.8678 212.1167 212.3040

$out
numeric(0)

$group
numeric(0)

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12"

> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5b3zr1257935697.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]
[1,] 153.30 160.00 165.00 169.80 173.90 177.80 180.40 185.90 188.60 194.90
[2,] 156.05 162.20 166.65 171.40 175.05 178.95 182.20 187.10 190.95 196.10
[3,] 157.30 162.85 168.45 172.00 175.90 179.65 183.35 187.95 193.25 196.20
[4,] 158.25 164.35 169.30 173.35 176.10 179.85 183.85 188.45 194.55 197.85
[5,] 159.30 164.90 169.70 173.80 176.70 180.40 184.40 188.60 194.90 198.50
      [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]
[1,] 198.60 204.80 209.70 213.70 218.60 223.40 228.50 231.90 240.50
[2,] 200.65 206.10 209.95 214.55 219.85 225.30 229.05 235.55 240.90
[3,] 203.20 207.45 211.55 216.15 221.40 226.40 229.75 237.30 241.80
[4,] 204.05 208.55 212.30 217.35 222.35 227.55 230.50 239.00 243.05
[5,] 204.50 209.10 213.40 218.10 223.10 228.40 231.70 240.20 243.90

$n
 [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12  8

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 156.2966 161.8694 167.2413 171.1106 175.4211 179.2395 182.5974 187.3343
[2,] 158.3034 163.8306 169.6587 172.8894 176.3789 180.0605 184.1026 188.5657
        [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 191.608 195.4018 201.6492 206.3325 210.4781 214.8729 220.2597 225.3738
[2,] 194.892 196.9982 204.7508 208.5675 212.6219 217.4271 222.5403 227.4262
        [,17]    [,18]   [,19]
[1,] 229.0886 235.7264 240.599
[2,] 230.4114 238.8736 243.001

$out
[1] 176.3 184.5

$group
[1] 6 8

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" NA  

Warning message:
In bxp(list(stats = c(153.3, 156.05, 157.3, 158.25, 159.3, 160,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6h1rv1257935697.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,] 24.44461 80.60 38.475
[2,] 24.60424 80.95 39.600
[3,] 25.85319 85.35 40.125
[4,] 26.02424 86.10 41.075
[5,] 26.17652 87.20 41.650

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 25.20552 83.00105 39.45224
[2,] 26.50086 87.69895 40.79776

$out
numeric(0)

$group
numeric(0)

$names
[1] "S.D."  "Range" "IQR"  

Warning message:
In bxp(list(stats = c(24.44461326442, 24.6042397314727, 25.8531884867133,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/12l1p1257935697.ps tmp/12l1p1257935697.png")
> system("convert tmp/2r0x01257935697.ps tmp/2r0x01257935697.png")
> system("convert tmp/3pvw01257935697.ps tmp/3pvw01257935697.png")
> system("convert tmp/4fz961257935697.ps tmp/4fz961257935697.png")
> system("convert tmp/5b3zr1257935697.ps tmp/5b3zr1257935697.png")
> system("convert tmp/6h1rv1257935697.ps tmp/6h1rv1257935697.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.695   1.278   2.054