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(15136
+ ,16733
+ ,20016
+ ,17708
+ ,18019
+ ,19227
+ ,22893
+ ,23739
+ ,21133
+ ,22591
+ ,26786
+ ,29740
+ ,15028
+ ,17977
+ ,20008
+ ,21354
+ ,19498
+ ,22125
+ ,25817
+ ,28779
+ ,20960
+ ,22254
+ ,27392
+ ,29945
+ ,16933
+ ,17892
+ ,20533
+ ,23569
+ ,22417
+ ,22084
+ ,26580
+ ,27454
+ ,24081
+ ,23451
+ ,28991
+ ,31386
+ ,16896
+ ,20045
+ ,23471
+ ,21747
+ ,25621
+ ,23859
+ ,25500
+ ,30998
+ ,24475
+ ,23145
+ ,29701
+ ,34365
+ ,17556
+ ,22077
+ ,25702
+ ,22214
+ ,26886
+ ,23191
+ ,27831
+ ,35406
+ ,23195
+ ,25110
+ ,30009
+ ,36242
+ ,18450
+ ,21845
+ ,26488
+ ,22394
+ ,28057
+ ,25451
+ ,24872
+ ,33424
+ ,24052
+ ,28449
+ ,33533
+ ,37351
+ ,19969
+ ,21701
+ ,26249
+ ,24493
+ ,24603
+ ,26485
+ ,30723
+ ,34569
+ ,26689
+ ,26157
+ ,32064
+ ,38870
+ ,21337
+ ,19419
+ ,23166
+ ,28286
+ ,24570
+ ,24001
+ ,33151
+ ,24878
+ ,26804
+ ,28967
+ ,33311
+ ,40226
+ ,20504
+ ,23060
+ ,23562
+ ,27562
+ ,23940
+ ,24584
+ ,34303
+ ,25517
+ ,23494
+ ,29095
+ ,32903
+ ,34379
+ ,16991
+ ,21109
+ ,23740
+ ,25552
+ ,21752
+ ,20294
+ ,29009
+ ,25500
+ ,24166
+ ,26960
+ ,31222
+ ,38641
+ ,14672
+ ,17543
+ ,25453
+ ,32683
+ ,22449
+ ,22316
+ ,27595
+ ,25451
+ ,25421
+ ,25288
+ ,32568
+ ,35110
+ ,16052
+ ,22146
+ ,21198
+ ,19543
+ ,22084
+ ,23816
+ ,29961
+ ,26773
+ ,26635
+ ,26972
+ ,30207
+ ,38687
+ ,16974
+ ,21697
+ ,24179
+ ,23757
+ ,25013
+ ,24019
+ ,30345
+ ,24488
+ ,25156
+ ,25650
+ ,30923
+ ,37240
+ ,17466
+ ,19463
+ ,24352
+ ,26805
+ ,25236
+ ,24735
+ ,29356
+ ,31234
+ ,22724
+ ,28496
+ ,32857
+ ,37198
+ ,13652
+ ,22784
+ ,23565
+ ,26323
+ ,23779
+ ,27549
+ ,29660
+ ,23356)
> 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] 176
> (np <- floor(n / par1))
[1] 14
> 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,] 15136 15028 16933 16896 17556 18450 19969 21337 20504 16991 14672 16052
 [2,] 16733 17977 17892 20045 22077 21845 21701 19419 23060 21109 17543 22146
 [3,] 20016 20008 20533 23471 25702 26488 26249 23166 23562 23740 25453 21198
 [4,] 17708 21354 23569 21747 22214 22394 24493 28286 27562 25552 32683 19543
 [5,] 18019 19498 22417 25621 26886 28057 24603 24570 23940 21752 22449 22084
 [6,] 19227 22125 22084 23859 23191 25451 26485 24001 24584 20294 22316 23816
 [7,] 22893 25817 26580 25500 27831 24872 30723 33151 34303 29009 27595 29961
 [8,] 23739 28779 27454 30998 35406 33424 34569 24878 25517 25500 25451 26773
 [9,] 21133 20960 24081 24475 23195 24052 26689 26804 23494 24166 25421 26635
[10,] 22591 22254 23451 23145 25110 28449 26157 28967 29095 26960 25288 26972
[11,] 26786 27392 28991 29701 30009 33533 32064 33311 32903 31222 32568 30207
[12,] 29740 29945 31386 34365 36242 37351 38870 40226 34379 38641 35110 38687
      [,13] [,14]
 [1,] 16974 17466
 [2,] 21697 19463
 [3,] 24179 24352
 [4,] 23757 26805
 [5,] 25013 25236
 [6,] 24019 24735
 [7,] 30345 29356
 [8,] 24488 31234
 [9,] 25156 22724
[10,] 25650 28496
[11,] 30923 32857
[12,] 37240 37198
> 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] 21143.42 22594.75 23780.92 24985.25 26284.92 27030.50 27714.33 27343.00
 [9] 26908.58 25411.33 25545.75 25339.50 25786.75 26660.17
> arr.sd
 [1] 4260.626 4523.441 4300.259 4810.754 5479.476 5486.091 5416.724 5880.496
 [9] 4733.506 5700.807 6017.624 5952.324 5092.400 5591.494
> arr.range
 [1] 14604 14917 14453 17469 18686 18901 18901 20807 13875 21650 20438 22635
[13] 20266 19732
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
  -301.6728       0.2173  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
       -2.489          1.089  

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

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

Coefficients:
(Intercept)     arr.mean  
 -1439.2090       0.7783  

> postscript(file="/var/www/html/rcomp/tmp/1c4a31243855909.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/2a3ou1243855909.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 
> 
> #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/3tf6k1243855909.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/4mqu71243855909.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/5fx9i1243855909.tab") 
> 
> system("convert tmp/1c4a31243855909.ps tmp/1c4a31243855909.png")
> system("convert tmp/2a3ou1243855909.ps tmp/2a3ou1243855909.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.556   0.310   1.001