R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(276986
+ ,260633
+ ,291551
+ ,275383
+ ,275302
+ ,231693
+ ,238829
+ ,274215
+ ,277808
+ ,299060
+ ,286629
+ ,232313
+ ,294053
+ ,267510
+ ,309739
+ ,280733
+ ,287298
+ ,235672
+ ,256449
+ ,288997
+ ,290789
+ ,321898
+ ,291834
+ ,241380
+ ,295469
+ ,258200
+ ,306102
+ ,281480
+ ,283101
+ ,237414
+ ,274834
+ ,299340
+ ,300383
+ ,340862
+ ,318794
+ ,265740
+ ,322656
+ ,281563
+ ,323461
+ ,312579
+ ,310784
+ ,262785
+ ,273754
+ ,320036
+ ,310336
+ ,342206
+ ,320052
+ ,265582
+ ,326988
+ ,300713
+ ,346414
+ ,317325
+ ,326208
+ ,270657
+ ,278158
+ ,324584
+ ,321801
+ ,343542
+ ,354040
+ ,278179
+ ,330246
+ ,307344
+ ,375874
+ ,335309
+ ,339271
+ ,280264
+ ,293689
+ ,341161
+ ,345097
+ ,368712
+ ,369403
+ ,288384
+ ,340981
+ ,319072
+ ,374214
+ ,344529
+ ,337271
+ ,281016
+ ,282224
+ ,320984
+ ,325426
+ ,366276
+ ,380296
+ ,300727
+ ,359326
+ ,327610
+ ,383563
+ ,352405
+ ,329351
+ ,294486
+ ,333454
+ ,334339
+ ,358000
+ ,396057
+ ,386976
+ ,307155
+ ,363909
+ ,344700
+ ,397561
+ ,376791
+ ,337085
+ ,299252
+ ,323136
+ ,329091
+ ,346991
+ ,461999
+ ,436533
+ ,360372
+ ,415467
+ ,382110
+ ,432197
+ ,424254
+ ,386728
+ ,354508
+ ,375765
+ ,367986
+ ,402378
+ ,426516
+ ,433313
+ ,338461
+ ,416834
+ ,381099
+ ,445673
+ ,412408
+ ,393997
+ ,348241
+ ,380134
+ ,373688
+ ,393588
+ ,434192
+ ,430731
+ ,344468
+ ,411891
+ ,370497
+ ,437305
+ ,411270
+ ,385495
+ ,341273
+ ,384217
+ ,373223
+ ,415771
+ ,448634
+ ,454341
+ ,350297
+ ,419104
+ ,398027
+ ,456059
+ ,430052
+ ,399757
+ ,362731
+ ,384896
+ ,385349
+ ,432289
+ ,468891
+ ,442702
+ ,370178
+ ,439400
+ ,393900
+ ,468700
+ ,438800
+ ,430100
+ ,366300
+ ,391000
+ ,380900
+ ,431400
+ ,465400
+ ,471500
+ ,387500
+ ,446400
+ ,421500
+ ,504800
+ ,492071
+ ,421253
+ ,396682
+ ,428000
+ ,421900
+ ,465600
+ ,525793
+ ,499855
+ ,435287
+ ,479499
+ ,473027
+ ,554410
+ ,489574
+ ,462157
+ ,420331)
> 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] 186
> (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]
 [1,] 276986 294053 295469 322656 326988 330246 340981 359326 363909 415467
 [2,] 260633 267510 258200 281563 300713 307344 319072 327610 344700 382110
 [3,] 291551 309739 306102 323461 346414 375874 374214 383563 397561 432197
 [4,] 275383 280733 281480 312579 317325 335309 344529 352405 376791 424254
 [5,] 275302 287298 283101 310784 326208 339271 337271 329351 337085 386728
 [6,] 231693 235672 237414 262785 270657 280264 281016 294486 299252 354508
 [7,] 238829 256449 274834 273754 278158 293689 282224 333454 323136 375765
 [8,] 274215 288997 299340 320036 324584 341161 320984 334339 329091 367986
 [9,] 277808 290789 300383 310336 321801 345097 325426 358000 346991 402378
[10,] 299060 321898 340862 342206 343542 368712 366276 396057 461999 426516
[11,] 286629 291834 318794 320052 354040 369403 380296 386976 436533 433313
[12,] 232313 241380 265740 265582 278179 288384 300727 307155 360372 338461
       [,11]  [,12]  [,13]  [,14]  [,15]
 [1,] 416834 411891 419104 439400 446400
 [2,] 381099 370497 398027 393900 421500
 [3,] 445673 437305 456059 468700 504800
 [4,] 412408 411270 430052 438800 492071
 [5,] 393997 385495 399757 430100 421253
 [6,] 348241 341273 362731 366300 396682
 [7,] 380134 384217 384896 391000 428000
 [8,] 373688 373223 385349 380900 421900
 [9,] 393588 415771 432289 431400 465600
[10,] 434192 448634 468891 465400 525793
[11,] 430731 454341 442702 471500 499855
[12,] 344468 350297 370178 387500 435287
> 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] 268366.8 280529.3 288476.6 303816.2 315717.4 331229.5 331084.7 346893.5
 [9] 364785.0 394973.6 396254.4 398684.5 412502.9 422075.0 454928.4
> arr.sd
 [1] 22733.99 25923.47 27940.98 26030.60 27986.71 32531.84 32934.45 31651.70
 [9] 47316.72 32114.09 32686.77 37139.08 34207.78 36910.41 41494.24
> arr.range
 [1]  67367  86226 103448  79421  83383  95610  99280 101571 162747  94852
[11] 101205 113068 106160 105200 129111
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
  3.515e+03    8.227e-02  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
      -1.5515         0.9344  

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

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

Coefficients:
(Intercept)     arr.mean  
  2.399e+04    2.201e-01  

> postscript(file="/var/wessaorg/rcomp/tmp/1qzpl1322940356.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/wessaorg/rcomp/tmp/2x38l1322940356.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/3j92d1322940356.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/wessaorg/rcomp/tmp/44d931322940356.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/wessaorg/rcomp/tmp/5phb01322940356.tab") 
> 
> try(system("convert tmp/1qzpl1322940356.ps tmp/1qzpl1322940356.png",intern=TRUE))
character(0)
> try(system("convert tmp/2x38l1322940356.ps tmp/2x38l1322940356.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.668   0.124   0.826