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(81.28,69.39,67.63,51.25,103.97,133.83,162.37,172.91,163.01,151.50,111.73,88.58,74.29,63.98,61.18,76.48,107.98,124.97,145.57,140.20,143.84,138.80,104.06,74.70,60.18,55.16,35.62,56.18,85.44,114.08,133.64,67.14,95.58,89.37,75.24,69.18,54.49,57.50,62.16,76.67,110.04,127.38,156.47,167.56,153.54,124.08,100.97,79.17,68.13,61.77,54.31,60.30,84.18,104.05,114.66,105.55,96.61,70.94,63.91,58.61,44.53,49.58,57.39,76.76,104.57,125.41,143.11,136.35,135.15,131.70,96.87,70.63,66.29,63.49,62.97,66.43,101.49,127.69,133.21,158.72,148.61,134.31,100.99,75.16,59.74,52.87,52.07,57.38,79.43,101.40,120.19,134.38,135.97,113.83,84.38,70.28,65.96,56.36,49.57,68.33,90.32,117.06,134.69,131.67,129.25,118.77,88.44,76.79,75.28,73.89,76.24,88.58,105.83,115.84,127.76,131.75,119.63,93.38,75.55,51.79)
> 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]
 [1,]  81.28  74.29  60.18  54.49  68.13  44.53  66.29  59.74  65.96  75.28
 [2,]  69.39  63.98  55.16  57.50  61.77  49.58  63.49  52.87  56.36  73.89
 [3,]  67.63  61.18  35.62  62.16  54.31  57.39  62.97  52.07  49.57  76.24
 [4,]  51.25  76.48  56.18  76.67  60.30  76.76  66.43  57.38  68.33  88.58
 [5,] 103.97 107.98  85.44 110.04  84.18 104.57 101.49  79.43  90.32 105.83
 [6,] 133.83 124.97 114.08 127.38 104.05 125.41 127.69 101.40 117.06 115.84
 [7,] 162.37 145.57 133.64 156.47 114.66 143.11 133.21 120.19 134.69 127.76
 [8,] 172.91 140.20  67.14 167.56 105.55 136.35 158.72 134.38 131.67 131.75
 [9,] 163.01 143.84  95.58 153.54  96.61 135.15 148.61 135.97 129.25 119.63
[10,] 151.50 138.80  89.37 124.08  70.94 131.70 134.31 113.83 118.77  93.38
[11,] 111.73 104.06  75.24 100.97  63.91  96.87 100.99  84.38  88.44  75.55
[12,]  88.58  74.70  69.18  79.17  58.61  70.63  75.16  70.28  76.79  51.79
      [,11]
 [1,]    NA
 [2,]    NA
 [3,]    NA
 [4,]    NA
 [5,]    NA
 [6,]    NA
 [7,]    NA
 [8,]    NA
 [9,]    NA
[10,]    NA
[11,]    NA
[12,]    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/1tgsq1274027370.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/2tgsq1274027370.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/3tgsq1274027370.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/4zm2m1274027370.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,] 44.530 49.580 49.570 51.25  79.43 101.400 114.660 105.550  95.58  70.940
[2,] 59.740 55.160 52.070 57.38  85.44 114.080 127.760 131.670 119.63  93.380
[3,] 66.125 59.635 59.285 67.38 102.73 121.015 134.165 135.365 135.56 121.425
[4,] 74.290 63.980 62.970 76.67 105.83 127.380 145.570 158.720 148.61 134.310
[5,] 81.280 73.890 76.240 88.58 110.04 133.830 162.370 172.910 163.01 151.500
       [,11]  [,12]
[1,]  63.910 58.610
[2,]  75.550 69.180
[3,]  92.655 72.665
[4,] 100.990 76.790
[5,] 111.730 79.170

$n
 [1] 10 10 10 10 10 10 10 10 10 10 10 10

$conf
         [,1]     [,2]     [,3]     [,4]      [,5]     [,6]     [,7]     [,8]
[1,] 58.85524 55.22818 53.83893 57.74195  92.54234 114.3698 125.2664 121.8497
[2,] 73.39476 64.04182 64.73107 77.01805 112.91766 127.6602 143.0636 148.8803
         [,9]    [,10]     [,11]    [,12]
[1,] 121.0804 100.9747  79.94416 68.86274
[2,] 150.0396 141.8753 105.36584 76.46726

$out
[1] 35.62 67.14 88.58 51.79

$group
[1]  3  8 12 12

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

Warning message:
In bxp(list(stats = c(44.53, 59.74, 66.125, 74.29, 81.28, 49.58,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5zm2m1274027370.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]
[1,]  51.250  61.180  35.620  54.490  54.310  44.530  62.97  52.070  49.570
[2,]  75.335  74.495  58.180  69.415  61.035  64.010  66.36  58.560  67.145
[3,] 107.850 106.020  72.210 105.505  69.535 100.720 101.24  81.905  89.380
[4,] 156.935 139.500  92.475 140.460 100.330 133.425 133.76 117.010 124.010
[5,] 172.910 145.570 133.640 167.560 114.660 143.110 158.72 135.970 134.690
       [,10] [,11]
[1,]  51.790    NA
[2,]  75.415    NA
[3,]  90.980    NA
[4,] 117.735    NA
[5,] 131.750    NA

$n
 [1] 12 12 12 12 12 12 12 12 12 12  0

$conf
          [,1]      [,2]     [,3]     [,4]     [,5]      [,6]      [,7]
[1,]  70.63169  76.37078 56.56782  73.1009 51.61229  69.05935  70.49841
[2,] 145.06831 135.66922 87.85218 137.9091 87.45771 132.38065 131.98159
          [,8]      [,9]     [,10] [,11]
[1,]  55.24556  63.44349  71.67756    NA
[2,] 108.56444 115.31651 110.28244    NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(51.25, 75.335, 107.85, 156.935, 172.91, 61.18,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/69v171274027370.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,]  7.549865 24.310  6.92750
[2,] 10.705514 34.590 11.25250
[3,] 11.513018 38.975 17.17250
[4,] 18.948056 57.625 22.81375
[5,] 30.879190 80.560 35.16500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,]  7.75354 28.46858 11.89934
[2,] 15.27250 49.48142 22.44566

$out
[1] 105.77

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(7.54986453152932, 10.7055142217566, 11.5130184120251,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1tgsq1274027370.ps tmp/1tgsq1274027370.png",intern=TRUE))
character(0)
> try(system("convert tmp/2tgsq1274027370.ps tmp/2tgsq1274027370.png",intern=TRUE))
character(0)
> try(system("convert tmp/3tgsq1274027370.ps tmp/3tgsq1274027370.png",intern=TRUE))
character(0)
> try(system("convert tmp/4zm2m1274027370.ps tmp/4zm2m1274027370.png",intern=TRUE))
character(0)
> try(system("convert tmp/5zm2m1274027370.ps tmp/5zm2m1274027370.png",intern=TRUE))
character(0)
> try(system("convert tmp/69v171274027370.ps tmp/69v171274027370.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.031   0.790   1.566