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(115.65
+ ,116.00
+ ,115.92
+ ,116.10
+ ,116.44
+ ,116.65
+ ,117.45
+ ,117.58
+ ,117.43
+ ,117.24
+ ,117.25
+ ,117.29
+ ,117.83
+ ,118.22
+ ,118.11
+ ,118.23
+ ,118.15
+ ,118.23
+ ,119.03
+ ,119.38
+ ,118.97
+ ,118.78
+ ,118.97
+ ,118.94
+ ,119.86
+ ,120.09
+ ,120.13
+ ,120.15
+ ,119.90
+ ,120.00
+ ,120.84
+ ,121.17
+ ,120.81
+ ,121.00
+ ,121.12
+ ,121.29
+ ,122.09
+ ,121.88
+ ,121.31
+ ,121.33
+ ,121.45
+ ,121.67
+ ,122.78
+ ,122.84
+ ,122.34
+ ,122.37
+ ,122.72
+ ,122.68
+ ,122.78
+ ,123.08
+ ,122.92
+ ,123.51
+ ,124.18
+ ,124.05
+ ,124.36
+ ,123.87
+ ,123.84
+ ,123.85
+ ,123.83
+ ,123.84
+ ,124.27
+ ,124.56
+ ,124.57
+ ,124.87
+ ,125.08
+ ,124.86
+ ,124.89
+ ,124.58
+ ,124.83
+ ,124.97
+ ,125.19
+ ,125.42
+ ,125.74
+ ,126.07
+ ,126.35
+ ,126.69
+ ,126.85
+ ,127.12
+ ,127.43
+ ,127.49
+ ,128.05
+ ,127.85
+ ,128.35
+ ,128.29
+ ,128.38
+ ,128.80
+ ,129.18
+ ,130.14
+ ,130.77
+ ,131.19
+ ,131.32
+ ,131.41
+ ,131.61
+ ,131.69
+ ,131.94
+ ,131.70
+ ,132.54
+ ,132.74
+ ,133.02
+ ,132.76
+ ,133.05
+ ,132.74
+ ,133.16
+ ,133.10
+ ,133.37
+ ,133.15
+ ,133.18
+ ,133.29
+ ,133.76
+ ,134.51
+ ,134.82
+ ,134.71
+ ,134.52
+ ,134.86
+ ,135.11
+ ,135.28
+ ,135.61
+ ,135.22
+ ,135.47
+ ,135.42
+ ,135.85
+ ,136.27
+ ,136.30
+ ,136.85
+ ,137.05
+ ,137.03
+ ,137.45
+ ,137.49
+ ,137.55
+ ,138.04
+ ,138.03
+ ,137.75
+ ,138.27
+ ,138.99
+ ,139.74
+ ,139.70
+ ,139.97
+ ,140.21
+ ,140.78
+ ,140.80
+ ,140.64
+ ,140.42
+ ,140.85
+ ,140.96
+ ,141.04
+ ,141.71
+ ,141.60
+ ,142.11
+ ,142.59
+ ,142.56
+ ,143.00
+ ,143.18
+ ,143.15
+ ,143.10
+ ,143.45
+ ,143.59
+ ,143.92
+ ,144.66
+ ,144.34
+ ,144.82
+ ,144.49
+ ,144.41
+ ,144.99
+ ,144.95
+ ,145.00
+ ,145.66
+ ,146.68
+ ,147.38
+ ,147.94
+ ,149.12
+ ,149.95
+ ,150.19
+ ,151.16
+ ,151.74
+ ,152.56
+ ,152.09
+ ,152.46
+ ,152.66
+ ,152.38
+ ,152.59
+ ,152.88
+ ,153.29
+ ,152.35
+ ,152.49
+ ,152.20
+ ,151.57
+ ,151.55
+ ,151.79
+ ,151.52
+ ,151.76
+ ,151.92
+ ,152.20
+ ,152.75
+ ,153.49
+ ,153.78
+ ,154.10
+ ,154.62
+ ,154.65
+ ,154.81
+ ,154.92
+ ,155.40
+ ,155.63
+ ,155.76)
> 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] 203
> (np <- floor(n / par1))
[1] 16
> 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]
 [1,] 115.65 117.83 119.86 122.09 122.78 124.27 125.74 128.38 132.54 133.76
 [2,] 116.00 118.22 120.09 121.88 123.08 124.56 126.07 128.80 132.74 134.51
 [3,] 115.92 118.11 120.13 121.31 122.92 124.57 126.35 129.18 133.02 134.82
 [4,] 116.10 118.23 120.15 121.33 123.51 124.87 126.69 130.14 132.76 134.71
 [5,] 116.44 118.15 119.90 121.45 124.18 125.08 126.85 130.77 133.05 134.52
 [6,] 116.65 118.23 120.00 121.67 124.05 124.86 127.12 131.19 132.74 134.86
 [7,] 117.45 119.03 120.84 122.78 124.36 124.89 127.43 131.32 133.16 135.11
 [8,] 117.58 119.38 121.17 122.84 123.87 124.58 127.49 131.41 133.10 135.28
 [9,] 117.43 118.97 120.81 122.34 123.84 124.83 128.05 131.61 133.37 135.61
[10,] 117.24 118.78 121.00 122.37 123.85 124.97 127.85 131.69 133.15 135.22
[11,] 117.25 118.97 121.12 122.72 123.83 125.19 128.35 131.94 133.18 135.47
[12,] 117.29 118.94 121.29 122.68 123.84 125.42 128.29 131.70 133.29 135.42
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]
 [1,] 135.85 138.27 141.04 143.92 147.94 152.88
 [2,] 136.27 138.99 141.71 144.66 149.12 153.29
 [3,] 136.30 139.74 141.60 144.34 149.95 152.35
 [4,] 136.85 139.70 142.11 144.82 150.19 152.49
 [5,] 137.05 139.97 142.59 144.49 151.16 152.20
 [6,] 137.03 140.21 142.56 144.41 151.74 151.57
 [7,] 137.45 140.78 143.00 144.99 152.56 151.55
 [8,] 137.49 140.80 143.18 144.95 152.09 151.79
 [9,] 137.55 140.64 143.15 145.00 152.46 151.52
[10,] 138.04 140.42 143.10 145.66 152.66 151.76
[11,] 138.03 140.85 143.45 146.68 152.38 151.92
[12,] 137.75 140.96 143.59 147.38 152.59 152.20
> 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] 116.7500 118.5700 120.5300 122.1217 123.6758 124.8408 127.1900 130.6775
 [9] 133.0083 134.9408 137.1383 140.1108 142.5900 145.1083 151.2367 152.1267
> arr.sd
 [1] 0.7014919 0.4915652 0.5523010 0.5822032 0.5012976 0.3134981 0.8692421
 [8] 1.2464358 0.2554793 0.5231454 0.7129814 0.8282891 0.8100393 1.0064504
[15] 1.5811120 0.5571410
> arr.range
 [1] 1.93 1.55 1.43 1.53 1.58 1.15 2.61 3.56 0.83 1.85 2.19 2.69 2.55 3.46 4.72
[16] 1.77
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
   -1.17031      0.01427  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
      -11.542          2.276  

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

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

Coefficients:
(Intercept)     arr.mean  
   -4.51419      0.05075  

> postscript(file="/var/www/html/rcomp/tmp/16oqz1291463379.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/26oqz1291463379.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/3sp651291463379.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/4dp5s1291463379.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/5yq3y1291463379.tab") 
> 
> try(system("convert tmp/16oqz1291463379.ps tmp/16oqz1291463379.png",intern=TRUE))
character(0)
> try(system("convert tmp/26oqz1291463379.ps tmp/26oqz1291463379.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.524   0.307   1.118