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(78.33,78.21,78.94,77.94,77.31,75.75,77.73,77.90,77.45,77.46,77.97,77.23,76.56,76.70,76.51,76.03,76.69,76.38,76.80,76.63,77.17,78.63,78.89,76.94,77.50,79.27,79.77,78.62,78.60,77.88,78.71,79.27,80.12,81.12,81.48,82.81,82.39,82.41,82.20,81.99,81.61,83.51,84.05,82.99,83.54,84.44,84.24,83.88,84.17,84.59,84.76,85.14,85.22,84.77,84.50,84.56,83.79,83.96,84.80,84.89,84.78,84.80,84.44,84.65,84.22,84.08,85.29,85.00,84.63,84.92,84.61,84.50,84.29,84.50,84.41,84.71,84.21,83.86,84.40,83.71,84.42,85.26,85.08,85.65,85.74,85.89,86.08,85.49,85.97,85.84,86.72,85.42,83.87,85.45,85.35,84.27,83.13,83.79,83.70,83.76,83.47,83.78,84.83,84.43,84.90,85.36,85.49,85.29)
> par1 = '12'
> par1 <- 7
> (n <- length(x))
[1] 108
> (np <- floor(n / par1))
[1] 15
> arr <- array(NA,dim=c(par1,np))
> j <- 0
> k <- 1
> for (i in 1:(np*par1))
+ {
+ j = j + 1
+ arr[j,k] <- x[i]
+ if (j == par1) {
+ j = 0
+ k=k+1
+ }
+ }
> arr
      [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
[1,] 78.33 77.90 76.51 78.63 78.60 82.81 84.05 84.59 83.79 84.65 84.61 83.86
[2,] 78.21 77.45 76.03 78.89 77.88 82.39 82.99 84.76 83.96 84.22 84.50 84.40
[3,] 78.94 77.46 76.69 76.94 78.71 82.41 83.54 85.14 84.80 84.08 84.29 83.71
[4,] 77.94 77.97 76.38 77.50 79.27 82.20 84.44 85.22 84.89 85.29 84.50 84.42
[5,] 77.31 77.23 76.80 79.27 80.12 81.99 84.24 84.77 84.78 85.00 84.41 85.26
[6,] 75.75 76.56 76.63 79.77 81.12 81.61 83.88 84.50 84.80 84.63 84.71 85.08
[7,] 77.73 76.70 77.17 78.62 81.48 83.51 84.17 84.56 84.44 84.92 84.21 85.65
     [,13] [,14] [,15]
[1,] 85.74 85.42 83.70
[2,] 85.89 83.87 83.76
[3,] 86.08 85.45 83.47
[4,] 85.49 85.35 83.78
[5,] 85.97 84.27 84.83
[6,] 85.84 83.13 84.43
[7,] 86.72 83.79 84.90
> arr.mean <- array(NA,dim=np)
> arr.sd <- array(NA,dim=np)
> arr.range <- array(NA,dim=np)
> for (j in 1:np)
+ {
+ arr.mean[j] <- mean(arr[,j],na.rm=TRUE)
+ arr.sd[j] <- sd(arr[,j],na.rm=TRUE)
+ arr.range[j] <- max(arr[,j],na.rm=TRUE) - min(arr[,j],na.rm=TRUE)
+ }
> arr.mean
 [1] 77.74429 77.32429 76.60143 78.51714 79.59714 82.41714 83.90143 84.79143
 [9] 84.49429 84.68429 84.46143 84.62571 85.96143 84.46857 84.12429
> arr.sd
 [1] 1.0160357 0.5424284 0.3547971 0.9849486 1.3524387 0.6097189 0.4934041
 [8] 0.2845130 0.4489936 0.4295291 0.1742057 0.7275497 0.3831200 0.9394401
[15] 0.5848891
> arr.range
 [1] 3.19 1.41 1.14 2.83 3.60 1.90 1.45 0.72 1.10 1.21 0.50 1.94 1.23 2.32 1.43
> (lm1 <- lm(arr.sd~arr.mean))

Call:
lm(formula = arr.sd ~ arr.mean)

Coefficients:
(Intercept)     arr.mean  
    3.95132     -0.04048  

> (lnlm1 <- lm(log(arr.sd)~log(arr.mean)))

Call:
lm(formula = log(arr.sd) ~ log(arr.mean))

Coefficients:
  (Intercept)  log(arr.mean)  
       21.914         -5.107  

> (lm2 <- lm(arr.range~arr.mean))

Call:
lm(formula = arr.range ~ arr.mean)

Coefficients:
(Intercept)     arr.mean  
    12.5329      -0.1313  

> postscript(file="/var/www/html/rcomp/tmp/1nix81292945185.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.mean,arr.sd,main='Standard Deviation-Mean Plot',xlab='mean',ylab='standard deviation')
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/2g9wb1292945185.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.mean,arr.range,main='Range-Mean Plot',xlab='mean',ylab='range')
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Standard Deviation-Mean Plot',4,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Section',header=TRUE)
> a<-table.element(a,'Mean',header=TRUE)
> a<-table.element(a,'Standard Deviation',header=TRUE)
> a<-table.element(a,'Range',header=TRUE)
> a<-table.row.end(a)
> for (j in 1:np) {
+ a<-table.row.start(a)
+ a<-table.element(a,j,header=TRUE)
+ a<-table.element(a,arr.mean[j])
+ a<-table.element(a,arr.sd[j] )
+ a<-table.element(a,arr.range[j] )
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/31acz1292945185.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Regression: S.E.(k) = alpha + beta * Mean(k)',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'alpha',header=TRUE)
> a<-table.element(a,lm1$coefficients[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'beta',header=TRUE)
> a<-table.element(a,lm1$coefficients[[2]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,summary(lm1)$coefficients[2,2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'T-STAT',header=TRUE)
> a<-table.element(a,summary(lm1)$coefficients[2,3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value',header=TRUE)
> a<-table.element(a,summary(lm1)$coefficients[2,4])
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/4nab51292945185.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Regression: ln S.E.(k) = alpha + beta * ln Mean(k)',2,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'alpha',header=TRUE)
> a<-table.element(a,lnlm1$coefficients[[1]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'beta',header=TRUE)
> a<-table.element(a,lnlm1$coefficients[[2]])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,summary(lnlm1)$coefficients[2,2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'T-STAT',header=TRUE)
> a<-table.element(a,summary(lnlm1)$coefficients[2,3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value',header=TRUE)
> a<-table.element(a,summary(lnlm1)$coefficients[2,4])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Lambda',header=TRUE)
> a<-table.element(a,1-lnlm1$coefficients[[2]])
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/5qbra1292945185.tab") 
> 
> try(system("convert tmp/1nix81292945185.ps tmp/1nix81292945185.png",intern=TRUE))
character(0)
> try(system("convert tmp/2g9wb1292945185.ps tmp/2g9wb1292945185.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.507   0.279   1.246