R version 2.6.0 (2007-10-03)
Copyright (C) 2007 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(0.00141557
+ ,0.001430964
+ ,0.00145722
+ ,0.001528139
+ ,0.00152196
+ ,0.001504125
+ ,0.001494319
+ ,0.001434645
+ ,0.00141161
+ ,0.001429132
+ ,0.001457646
+ ,0.001429132
+ ,0.001422965
+ ,0.001427508
+ ,0.001433416
+ ,0.001405819
+ ,0.001344975
+ ,0.001357137
+ ,0.001340863
+ ,0.001335916
+ ,0.001350801
+ ,0.001412697
+ ,0.001410723
+ ,0.001384623
+ ,0.001393066
+ ,0.001405038
+ ,0.00137956
+ ,0.001352134
+ ,0.001358754
+ ,0.001412203
+ ,0.001458392
+ ,0.001439899
+ ,0.001447082
+ ,0.001437112
+ ,0.001392779
+ ,0.001349118
+ ,0.001388781
+ ,0.00136137
+ ,0.001365636
+ ,0.001374079
+ ,0.001399891
+ ,0.001368101
+ ,0.001346822
+ ,0.001399987
+ ,0.001391634
+ ,0.001387077
+ ,0.001373709
+ ,0.001366456
+ ,0.001352757
+ ,0.001346822
+ ,0.001369016
+ ,0.001389351
+ ,0.001384058
+ ,0.001375375
+ ,0.001367095
+ ,0.001362003
+ ,0.001305986
+ ,0.001306559
+ ,0.001298027
+ ,0.001258193
+ ,0.001214471
+ ,0.001184023
+ ,0.001136008
+ ,0.001095042
+ ,0.001074779
+ ,0.001077089
+ ,0.001022973
+ ,0.000977405
+ ,0.001062808
+ ,0.001020932
+ ,0.001030587
+ ,0.001064356
+ ,0.001073603
+ ,0.001039518
+ ,0.001056033
+ ,0.00104279
+ ,0.001007661
+ ,0.001026488
+ ,0.001015493
+ ,0.001026075
+ ,0.001036457
+ ,0.001042076
+ ,0.001037867
+ ,0.001003734
+ ,0.00097871)
> 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] 85
> (np <- floor(n / par1))
[1] 7
> 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]
 [1,] 0.001415570 0.001422965 0.001393066 0.001388781 0.001352757 0.001214471
 [2,] 0.001430964 0.001427508 0.001405038 0.001361370 0.001346822 0.001184023
 [3,] 0.001457220 0.001433416 0.001379560 0.001365636 0.001369016 0.001136008
 [4,] 0.001528139 0.001405819 0.001352134 0.001374079 0.001389351 0.001095042
 [5,] 0.001521960 0.001344975 0.001358754 0.001399891 0.001384058 0.001074779
 [6,] 0.001504125 0.001357137 0.001412203 0.001368101 0.001375375 0.001077089
 [7,] 0.001494319 0.001340863 0.001458392 0.001346822 0.001367095 0.001022973
 [8,] 0.001434645 0.001335916 0.001439899 0.001399987 0.001362003 0.000977405
 [9,] 0.001411610 0.001350801 0.001447082 0.001391634 0.001305986 0.001062808
[10,] 0.001429132 0.001412697 0.001437112 0.001387077 0.001306559 0.001020932
[11,] 0.001457646 0.001410723 0.001392779 0.001373709 0.001298027 0.001030587
[12,] 0.001429132 0.001384623 0.001349118 0.001366456 0.001258193 0.001064356
             [,7]
 [1,] 0.001073603
 [2,] 0.001039518
 [3,] 0.001056033
 [4,] 0.001042790
 [5,] 0.001007661
 [6,] 0.001026488
 [7,] 0.001015493
 [8,] 0.001026075
 [9,] 0.001036457
[10,] 0.001042076
[11,] 0.001037867
[12,] 0.001003734
> 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] 0.001459538 0.001385620 0.001402095 0.001376962 0.001342937 0.001080039
[7] 0.001033983
> arr.sd
[1] 4.191413e-05 3.736231e-05 3.791928e-05 1.651653e-05 4.103587e-05
[6] 6.911840e-05 1.981572e-05
> arr.range
[1] 0.000454536 0.000393898 0.000402359 0.000357197 0.000381690 0.000187983
[7] 0.000058110
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
  6.036e-05   -1.749e-02  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
    -10.53735       -0.03833  

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

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

Coefficients:
(Intercept)     arr.mean  
 -0.0007516    0.8255703  

> postscript(file="/var/www/html/rcomp/tmp/1s3tc1196351589.ps",horizontal=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/2loxr1196351589.ps",horizontal=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 
> 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/37th51196351589.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/4g74r1196351589.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/5ugb51196351589.tab") 
> 
> system("convert tmp/1s3tc1196351589.ps tmp/1s3tc1196351589.png")
> system("convert tmp/2loxr1196351589.ps tmp/2loxr1196351589.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.748   0.292   0.997