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.

  Natural language support but running in an English locale

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(11974
+ ,10106
+ ,12069
+ ,11412
+ ,11180
+ ,10508
+ ,11288
+ ,10928
+ ,10199
+ ,11030
+ ,11234
+ ,13747
+ ,13912
+ ,12376
+ ,12264
+ ,11675
+ ,11271
+ ,10672
+ ,10933
+ ,10379
+ ,10187
+ ,10747
+ ,10970
+ ,12175
+ ,14200
+ ,11676
+ ,11258
+ ,10872
+ ,11148
+ ,10690
+ ,10684
+ ,11658
+ ,10178
+ ,10981
+ ,10773
+ ,11665
+ ,11359
+ ,10716
+ ,12928
+ ,12317
+ ,11641
+ ,10459
+ ,10953
+ ,10703
+ ,10703
+ ,11101
+ ,11334
+ ,13268
+ ,13145
+ ,12334
+ ,13153
+ ,11289
+ ,11374
+ ,10914
+ ,11299
+ ,11284
+ ,10694
+ ,11077
+ ,11104
+ ,12820
+ ,14915
+ ,11773
+ ,11608
+ ,11468
+ ,11511
+ ,11200
+ ,11164
+ ,10960
+ ,10667
+ ,11556
+ ,11372
+ ,12333
+ ,13102
+ ,11115
+ ,12572
+ ,11557
+ ,12059
+ ,11420
+ ,11185
+ ,11113
+ ,10706
+ ,11523
+ ,11391
+ ,12634
+ ,13469
+ ,11735
+ ,13281
+ ,11968
+ ,11623
+ ,11084
+ ,11509
+ ,11134
+ ,10438
+ ,11530
+ ,11491
+ ,13093
+ ,13106
+ ,11305
+ ,13113
+ ,12203
+ ,11309
+ ,11088
+ ,11234
+ ,11619
+ ,10942
+ ,11445
+ ,11291
+ ,13281
+ ,13726
+ ,11300
+ ,11983
+ ,11092
+ ,11093
+ ,10692
+ ,10786
+ ,11166
+ ,10553
+ ,11103
+ ,10969
+ ,12090
+ ,12544
+ ,12264
+ ,13783
+ ,11214
+ ,11453
+ ,10883
+ ,10381
+ ,10348
+ ,10024
+ ,10805
+ ,10796
+ ,11907
+ ,12261
+ ,11377
+ ,12689
+ ,11474
+ ,10992
+ ,10764
+ ,12164
+ ,10409
+ ,10398
+ ,10349
+ ,10865
+ ,11630
+ ,12221
+ ,10884
+ ,12019
+ ,11021
+ ,10799
+ ,10423
+ ,10484
+ ,10450
+ ,9906
+ ,11049
+ ,11281
+ ,12485
+ ,12849
+ ,11380
+ ,12079
+ ,11366
+ ,11328
+ ,10444
+ ,10854
+ ,10434
+ ,10137
+ ,10992
+ ,10906
+ ,12367
+ ,14371
+ ,11695
+ ,11546
+ ,10922
+ ,10670
+ ,10254
+ ,10573
+ ,10239
+ ,10253
+ ,11176
+ ,10719
+ ,11817)
> 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] 180
> (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,] 11974 13912 14200 11359 13145 14915 13102 13469 13106 13726 12544 12261
 [2,] 10106 12376 11676 10716 12334 11773 11115 11735 11305 11300 12264 11377
 [3,] 12069 12264 11258 12928 13153 11608 12572 13281 13113 11983 13783 12689
 [4,] 11412 11675 10872 12317 11289 11468 11557 11968 12203 11092 11214 11474
 [5,] 11180 11271 11148 11641 11374 11511 12059 11623 11309 11093 11453 10992
 [6,] 10508 10672 10690 10459 10914 11200 11420 11084 11088 10692 10883 10764
 [7,] 11288 10933 10684 10953 11299 11164 11185 11509 11234 10786 10381 12164
 [8,] 10928 10379 11658 10703 11284 10960 11113 11134 11619 11166 10348 10409
 [9,] 10199 10187 10178 10703 10694 10667 10706 10438 10942 10553 10024 10398
[10,] 11030 10747 10981 11101 11077 11556 11523 11530 11445 11103 10805 10349
[11,] 11234 10970 10773 11334 11104 11372 11391 11491 11291 10969 10796 10865
[12,] 13747 12175 11665 13268 12820 12333 12634 13093 13281 12090 11907 11630
      [,13] [,14] [,15]
 [1,] 12221 12849 14371
 [2,] 10884 11380 11695
 [3,] 12019 12079 11546
 [4,] 11021 11366 10922
 [5,] 10799 11328 10670
 [6,] 10423 10444 10254
 [7,] 10484 10854 10573
 [8,] 10450 10434 10239
 [9,]  9906 10137 10253
[10,] 11049 10992 11176
[11,] 11281 10906 10719
[12,] 12485 12367 11817
> 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] 11306.25 11463.42 11315.25 11456.83 11707.25 11710.58 11698.08 11862.92
 [9] 11828.00 11379.42 11366.83 11281.00 11085.17 11261.33 11186.25
> arr.sd
 [1]  978.9051 1067.1559 1018.3747  919.7144  896.3152 1092.0878  731.9274
 [8]  940.5687  865.6791  870.9811 1088.7752  787.5960  790.9469  821.4514
[15] 1147.9982
> arr.range
 [1] 3641 3725 4022 2809 2459 4248 2396 3031 2339 3173 3759 2340 2579 2712 4132
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
 1277.19404     -0.02990  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
       9.3252        -0.2668  

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

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

Coefficients:
(Intercept)     arr.mean  
 10126.6490      -0.6081  

> postscript(file="/var/www/html/freestat/rcomp/tmp/1a4wp1292934948.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/freestat/rcomp/tmp/2a4wp1292934948.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/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/3v5ud1292934948.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/freestat/rcomp/tmp/4zotj1292934948.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/freestat/rcomp/tmp/5dg9s1292934948.tab") 
> 
> try(system("convert tmp/1a4wp1292934948.ps tmp/1a4wp1292934948.png",intern=TRUE))
character(0)
> try(system("convert tmp/2a4wp1292934948.ps tmp/2a4wp1292934948.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.862   0.424   1.091