R version 2.12.1 (2010-12-16)
Copyright (C) 2010 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(1483509
+ ,8036554
+ ,4623093
+ ,5528662
+ ,4221032
+ ,8061847
+ ,7640066
+ ,2935533
+ ,8161548
+ ,2543967
+ ,13163450
+ ,3348436
+ ,3997440
+ ,2322911
+ ,2019457
+ ,3047748
+ ,5728767
+ ,2605173
+ ,5646743
+ ,13121544
+ ,3453409
+ ,1878333
+ ,4247362
+ ,23022552
+ ,7646203
+ ,9016602
+ ,3606568
+ ,3173510
+ ,17568772
+ ,10805045
+ ,31056269
+ ,15623385
+ ,6663443
+ ,35435745
+ ,2823250
+ ,5197089
+ ,4120632
+ ,8832767
+ ,3695374
+ ,8385805
+ ,3777904
+ ,5199532
+ ,5297275
+ ,14847382
+ ,5900158
+ ,4416718
+ ,3926429
+ ,4876884
+ ,2795297
+ ,3385527
+ ,3877941
+ ,3556729
+ ,4982836
+ ,2976325
+ ,2295026
+ ,2218752
+ ,4146062
+ ,3302091
+ ,3864505
+ ,5454794
+ ,1749836
+ ,6684048
+ ,2809918
+ ,4092664
+ ,5070470
+ ,9814477
+ ,6665318
+ ,3912554
+ ,6188129
+ ,3627991
+ ,3308767
+ ,3820332
+ ,4932979
+ ,5567917
+ ,5020814
+ ,3803273
+ ,3999984
+ ,4883104
+ ,13731747
+ ,47531824
+ ,8415570
+ ,22178158
+ ,61211654
+ ,18223748
+ ,17678085
+ ,49299580
+ ,25899948
+ ,34121754
+ ,9859231
+ ,29740892
+ ,21085212
+ ,43003866
+ ,59549247
+ ,18026465
+ ,4680597
+ ,5564728
+ ,11792347
+ ,10371624
+ ,3728446
+ ,5732978
+ ,4067638
+ ,2395508
+ ,5018801
+ ,22068888
+ ,7678580
+ ,15510095
+ ,6471239
+ ,14349204
+ ,35151574
+ ,8210488
+ ,5022664
+ ,13996871
+ ,12822431
+ ,14011552
+ ,20260980
+ ,23718976
+ ,45833049
+ ,30688420
+ ,16576062
+ ,14844405
+ ,16728286
+ ,43477680
+ ,57497427
+ ,24233726
+ ,24921208
+ ,9516725
+ ,27977239
+ ,21632046
+ ,22956809
+ ,9704324
+ ,19871149
+ ,5553842
+ ,5667858
+ ,4348188
+ ,10025042
+ ,10639796
+ ,8639184
+ ,10764378
+ ,12097733
+ ,3988414
+ ,4607102
+ ,7126895
+ ,6009625
+ ,21533237
+ ,5986771
+ ,5455310
+ ,1822874
+ ,3374062
+ ,2920748
+ ,2295942
+ ,6809829
+ ,3318281
+ ,13784645
+ ,7366577
+ ,1628637
+ ,4258976
+ ,7159779
+ ,8098401
+ ,6894240
+ ,3771246
+ ,3249726
+ ,3147380
+ ,4063037
+ ,9621916
+ ,5890158
+ ,2142901
+ ,3145007
+ ,1562168
+ ,3303103
+ ,5886910
+ ,3454270
+ ,6995348
+ ,6487869
+ ,12091976
+ ,3934625
+ ,3999749
+ ,3613526
+ ,4271706
+ ,4253390
+ ,5551591
+ ,4663041
+ ,2104104
+ ,5385399
+ ,6205877
+ ,7529500
+ ,17222705
+ ,6230913
+ ,6508275
+ ,4518884
+ ,4234991
+ ,5625388
+ ,5810139
+ ,6942187
+ ,3711188
+ ,4261281
+ ,1989945
+ ,5033342
+ ,7239565
+ ,11058795
+ ,7384772
+ ,3884771
+ ,3239201
+ ,2316403
+ ,4034947
+ ,3245271
+ ,2387251
+ ,2174886
+ ,3436080
+ ,3738956
+ ,1884730
+ ,1509144
+ ,42728366
+ ,3446317
+ ,4600683
+ ,2953615
+ ,3570060
+ ,2130208
+ ,2442943
+ ,4892020
+ ,3222192
+ ,3121617
+ ,3665542
+ ,5519432
+ ,4113468
+ ,1714614
+ ,3651985
+ ,2419548
+ ,2378854
+ ,2303949
+ ,2555534
+ ,1713005
+ ,1705960
+ ,6115046
+ ,3951044
+ ,3785568
+ ,4670530
+ ,2265100
+ ,1105643
+ ,2814152
+ ,3728673
+ ,2038949
+ ,2402919
+ ,2348814
+ ,2797822
+ ,902505
+ ,1331319
+ ,4204238
+ ,2212485
+ ,6797382
+ ,4532324
+ ,1778808
+ ,1890720
+ ,5463736
+ ,11368931
+ ,2040164
+ ,4276399
+ ,3714445
+ ,2068168
+ ,1003842
+ ,2858535
+ ,2355484
+ ,2719262
+ ,1897741
+ ,3945185
+ ,3799916
+ ,1017654
+ ,3052241
+ ,3932970
+ ,3598151
+ ,2296005
+ ,2202018
+ ,2461777
+ ,2452042
+ ,2185142
+ ,11968502
+ ,20395972
+ ,21756900
+ ,30024300
+ ,10811344
+ ,1819202
+ ,1276885
+ ,2946701
+ ,3587459
+ ,2832691
+ ,6674805
+ ,3868362
+ ,4302909
+ ,23265229
+ ,22348002
+ ,11883953
+ ,6634979
+ ,2935493
+ ,3425669
+ ,1171611
+ ,6875879
+ ,19451908
+ ,13885933
+ ,7643317
+ ,10797966
+ ,7297445
+ ,8739736
+ ,12455537
+ ,24291181
+ ,4215150
+ ,28652176
+ ,6851172
+ ,3746871
+ ,7327861
+ ,16829710
+ ,13778594
+ ,6463717
+ ,8956867
+ ,21204915
+ ,16115855
+ ,2536113
+ ,16645717
+ ,17003730
+ ,15969006
+ ,31020427
+ ,23798897
+ ,20770321
+ ,44410402
+ ,27037491
+ ,29627771
+ ,18189792
+ ,4654610
+ ,12307201
+ ,15300578
+ ,10623864
+ ,6880178
+ ,29947357
+ ,18611399
+ ,42432604
+ ,20208278
+ ,14004392
+ ,25737765
+ ,16735738
+ ,22450825
+ ,6880840
+ ,8510379
+ ,8182481
+ ,10948683
+ ,4805277
+ ,2589229
+ ,5658407
+ ,12862611
+ ,5666188
+ ,6875556
+ ,7098766
+ ,36083309
+ ,10200330
+ ,7784976)
> 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] 352
> (np <- floor(n / par1))
[1] 29
> 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]
 [1,]  1483509  3997440  7646203  4120632 2795297 1749836  4932979 17678085
 [2,]  8036554  2322911  9016602  8832767 3385527 6684048  5567917 49299580
 [3,]  4623093  2019457  3606568  3695374 3877941 2809918  5020814 25899948
 [4,]  5528662  3047748  3173510  8385805 3556729 4092664  3803273 34121754
 [5,]  4221032  5728767 17568772  3777904 4982836 5070470  3999984  9859231
 [6,]  8061847  2605173 10805045  5199532 2976325 9814477  4883104 29740892
 [7,]  7640066  5646743 31056269  5297275 2295026 6665318 13731747 21085212
 [8,]  2935533 13121544 15623385 14847382 2218752 3912554 47531824 43003866
 [9,]  8161548  3453409  6663443  5900158 4146062 6188129  8415570 59549247
[10,]  2543967  1878333 35435745  4416718 3302091 3627991 22178158 18026465
[11,] 13163450  4247362  2823250  3926429 3864505 3308767 61211654  4680597
[12,]  3348436 23022552  5197089  4876884 5454794 3820332 18223748  5564728
          [,9]    [,10]    [,11]    [,12]    [,13]   [,14]    [,15]    [,16]
 [1,] 11792347 35151574 16728286  5667858  5986771 7159779  3303103  4663041
 [2,] 10371624  8210488 43477680  4348188  5455310 8098401  5886910  2104104
 [3,]  3728446  5022664 57497427 10025042  1822874 6894240  3454270  5385399
 [4,]  5732978 13996871 24233726 10639796  3374062 3771246  6995348  6205877
 [5,]  4067638 12822431 24921208  8639184  2920748 3249726  6487869  7529500
 [6,]  2395508 14011552  9516725 10764378  2295942 3147380 12091976 17222705
 [7,]  5018801 20260980 27977239 12097733  6809829 4063037  3934625  6230913
 [8,] 22068888 23718976 21632046  3988414  3318281 9621916  3999749  6508275
 [9,]  7678580 45833049 22956809  4607102 13784645 5890158  3613526  4518884
[10,] 15510095 30688420  9704324  7126895  7366577 2142901  4271706  4234991
[11,]  6471239 16576062 19871149  6009625  1628637 3145007  4253390  5625388
[12,] 14349204 14844405  5553842 21533237  4258976 1562168  5551591  5810139
         [,17]    [,18]   [,19]   [,20]   [,21]    [,22]    [,23]    [,24]
 [1,]  6942187  3245271 2130208 2303949 2038949  5463736  3799916 21756900
 [2,]  3711188  2387251 2442943 2555534 2402919 11368931  1017654 30024300
 [3,]  4261281  2174886 4892020 1713005 2348814  2040164  3052241 10811344
 [4,]  1989945  3436080 3222192 1705960 2797822  4276399  3932970  1819202
 [5,]  5033342  3738956 3121617 6115046  902505  3714445  3598151  1276885
 [6,]  7239565  1884730 3665542 3951044 1331319  2068168  2296005  2946701
 [7,] 11058795  1509144 5519432 3785568 4204238  1003842  2202018  3587459
 [8,]  7384772 42728366 4113468 4670530 2212485  2858535  2461777  2832691
 [9,]  3884771  3446317 1714614 2265100 6797382  2355484  2452042  6674805
[10,]  3239201  4600683 3651985 1105643 4532324  2719262  2185142  3868362
[11,]  2316403  2953615 2419548 2814152 1778808  1897741 11968502  4302909
[12,]  4034947  3570060 2378854 3728673 1890720  3945185 20395972 23265229
         [,25]    [,26]    [,27]    [,28]    [,29]
 [1,] 22348002  8739736 21204915 18189792 16735738
 [2,] 11883953 12455537 16115855  4654610 22450825
 [3,]  6634979 24291181  2536113 12307201  6880840
 [4,]  2935493  4215150 16645717 15300578  8510379
 [5,]  3425669 28652176 17003730 10623864  8182481
 [6,]  1171611  6851172 15969006  6880178 10948683
 [7,]  6875879  3746871 31020427 29947357  4805277
 [8,] 19451908  7327861 23798897 18611399  2589229
 [9,] 13885933 16829710 20770321 42432604  5658407
[10,]  7643317 13778594 44410402 20208278 12862611
[11,] 10797966  6463717 27037491 14004392  5666188
[12,]  7297445  8956867 29627771 25737765  6875556
> 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]  5812308  5924287 12384657  6106405  3571324  4812042 16625064 26542467
 [9]  9098779 20094789 23672538  8787288  4918554  4895497  5320339  6336601
[17]  5091366  6306280  3272702  3059517  2769857  3642658  4946866  9430566
[25]  9529346 11859048 22178387 18241502  9347184
> arr.sd
 [1]  3312518  6184841 10834927  3231245   980351  2200853 18847812 17382379
 [9]  5891084 11870449 14609790  4887134  3388000  2561647  2463648  3694117
[17]  2608189 11503252  1153993  1435915  1647415  2725125  5606468  9900871
[25]  6490909  7876458 10338023 10513921  5632079
> arr.range
 [1] 11679941 21144219 32612495 11152008  3236042  8064641 57408381 54868650
 [9] 19673380 40810385 51943585 17544823 12156008  8059748  8788873 15118601
[17]  9068850 41219222  3804818  5009403  5894877 10365089 19378318 28747415
[25] 21176391 24905305 41874289 37777994 19861596
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
  578809.48         0.63  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
       -2.459          1.127  

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

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

Coefficients:
(Intercept)     arr.mean  
  2.115e+06    2.120e+00  

> postscript(file="/var/www/rcomp/tmp/1bijn1322934120.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/rcomp/tmp/2jnv11322934120.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/rcomp/tmp/37bz01322934120.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/rcomp/tmp/4iwmb1322934120.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/rcomp/tmp/5o2nc1322934120.tab") 
> 
> try(system("convert tmp/1bijn1322934120.ps tmp/1bijn1322934120.png",intern=TRUE))
character(0)
> try(system("convert tmp/2jnv11322934120.ps tmp/2jnv11322934120.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.900   0.152   1.077