R version 2.6.2 (2008-02-08)
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.

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(1.327
+ ,1.3231
+ ,1.3106
+ ,1.3084
+ ,1.3006
+ ,1.3018
+ ,1.2988
+ ,1.2984
+ ,1.2893
+ ,1.2941
+ ,1.2953
+ ,1.2909
+ ,1.2922
+ ,1.2958
+ ,1.2936
+ ,1.304
+ ,1.3005
+ ,1.2978
+ ,1.2901
+ ,1.2921
+ ,1.2972
+ ,1.2954
+ ,1.302
+ ,1.302
+ ,1.2925
+ ,1.2955
+ ,1.2987
+ ,1.2991
+ ,1.3007
+ ,1.2956
+ ,1.3022
+ ,1.3082
+ ,1.3137
+ ,1.3119
+ ,1.3132
+ ,1.3145
+ ,1.3145
+ ,1.3106
+ ,1.3134
+ ,1.316
+ ,1.323
+ ,1.3211
+ ,1.3225
+ ,1.3163
+ ,1.3083
+ ,1.31
+ ,1.3135
+ ,1.3152
+ ,1.3155
+ ,1.3156
+ ,1.3218
+ ,1.3183
+ ,1.3226
+ ,1.3325
+ ,1.3303
+ ,1.3296
+ ,1.3293
+ ,1.3351
+ ,1.3327
+ ,1.3265
+ ,1.3347
+ ,1.3348
+ ,1.3352
+ ,1.3318
+ ,1.3366
+ ,1.3358
+ ,1.3352
+ ,1.3373
+ ,1.3426
+ ,1.3418
+ ,1.3467
+ ,1.3532
+ ,1.355
+ ,1.3549
+ ,1.3577
+ ,1.3601
+ ,1.3606
+ ,1.3557
+ ,1.3582
+ ,1.3649
+ ,1.3596
+ ,1.3643
+ ,1.3605
+ ,1.3588
+ ,1.3613
+ ,1.3561
+ ,1.3615
+ ,1.3558
+ ,1.3535
+ ,1.3527
+ ,1.3486
+ ,1.3549
+ ,1.3538
+ ,1.3574
+ ,1.3516
+ ,1.3477
+ ,1.3444
+ ,1.3454
+ ,1.349
+ ,1.3448
+ ,1.3441
+ ,1.3453
+ ,1.3509
+ ,1.342
+ ,1.3453
+ ,1.3436
+ ,1.3482
+ ,1.3532
+ ,1.3513
+ ,1.347
+ ,1.3349
+ ,1.3355
+ ,1.3345
+ ,1.3287
+ ,1.3304
+ ,1.3314
+ ,1.3404
+ ,1.3403
+ ,1.3427
+ ,1.3397
+ ,1.3441
+ ,1.3461
+ ,1.346
+ ,1.3438
+ ,1.3467
+ ,1.3505
+ ,1.3588
+ ,1.3601
+ ,1.3618
+ ,1.364
+ ,1.3596
+ ,1.3621
+ ,1.3666
+ ,1.3753
+ ,1.3788
+ ,1.3782
+ ,1.3781
+ ,1.3771
+ ,1.3779
+ ,1.382
+ ,1.3803
+ ,1.3821
+ ,1.3833
+ ,1.3743
+ ,1.3722
+ ,1.3651
+ ,1.3659
+ ,1.3707
+ ,1.3663
+ ,1.3664
+ ,1.3694
+ ,1.3818
+ ,1.3794
+ ,1.3794
+ ,1.3729
+ ,1.365
+ ,1.3651
+ ,1.3591
+ ,1.3476
+ ,1.3405
+ ,1.3454
+ ,1.3476
+ ,1.3508
+ ,1.3493
+ ,1.3574
+ ,1.3615
+ ,1.3658
+ ,1.3664
+ ,1.3631
+ ,1.361
+ ,1.3705
+ ,1.3632
+ ,1.358
+ ,1.3588
+ ,1.3669
+ ,1.3696
+ ,1.3795
+ ,1.3824
+ ,1.3885
+ ,1.3897
+ ,1.386
+ ,1.3877
+ ,1.3867
+ ,1.3975
+ ,1.403
+ ,1.4049
+ ,1.4113
+ ,1.4106
+ ,1.4127
+ ,1.418
+ ,1.4179
+ ,1.4232
+ ,1.4165
+ ,1.4195
+ ,1.4109
+ ,1.4136
+ ,1.4089
+ ,1.4037
+ ,1.4146
+ ,1.4199
+ ,1.4173
+ ,1.4226
+ ,1.415
+ ,1.42
+ ,1.4299
+ ,1.4288
+ ,1.4166
+ ,1.4254
+ ,1.423
+ ,1.4309
+ ,1.4384
+ ,1.4391
+ ,1.4407
+ ,1.4447
+ ,1.4423
+ ,1.4479
+ ,1.4488
+ ,1.4547
+ ,1.4722
+ ,1.4666
+ ,1.4683
+ ,1.4579
+ ,1.4607
+ ,1.47
+ ,1.4639
+ ,1.4651
+ ,1.4654
+ ,1.4785
+ ,1.4814
+ ,1.4829
+ ,1.4809
+ ,1.4845
+ ,1.4874
+ ,1.4747
+ ,1.4738
+ ,1.4761
+ ,1.4666
+ ,1.4741
+ ,1.472
+ ,1.4554
+ ,1.4649
+ ,1.4718
+ ,1.4672
+ ,1.4675
+ ,1.4683
+ ,1.4509
+ ,1.4393
+ ,1.4416
+ ,1.4385
+ ,1.4349
+ ,1.438
+ ,1.4398
+ ,1.4516
+ ,1.4692
+ ,1.4721)
> par1 = '4'
> #'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] 255
> (np <- floor(n / par1))
[1] 63
> 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,] 1.3270 1.3006 1.2893 1.2922 1.3005 1.2972 1.2925 1.3007 1.3137 1.3145
[2,] 1.3231 1.3018 1.2941 1.2958 1.2978 1.2954 1.2955 1.2956 1.3119 1.3106
[3,] 1.3106 1.2988 1.2953 1.2936 1.2901 1.3020 1.2987 1.3022 1.3132 1.3134
[4,] 1.3084 1.2984 1.2909 1.3040 1.2921 1.3020 1.2991 1.3082 1.3145 1.3160
      [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
[1,] 1.3230 1.3083 1.3155 1.3226 1.3293 1.3347 1.3366 1.3426 1.3550 1.3606
[2,] 1.3211 1.3100 1.3156 1.3325 1.3351 1.3348 1.3358 1.3418 1.3549 1.3557
[3,] 1.3225 1.3135 1.3218 1.3303 1.3327 1.3352 1.3352 1.3467 1.3577 1.3582
[4,] 1.3163 1.3152 1.3183 1.3296 1.3265 1.3318 1.3373 1.3532 1.3601 1.3649
      [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27]  [,28]  [,29]  [,30]
[1,] 1.3596 1.3613 1.3535 1.3538 1.3444 1.3441 1.3453 1.3513 1.3345 1.3404
[2,] 1.3643 1.3561 1.3527 1.3574 1.3454 1.3453 1.3436 1.3470 1.3287 1.3403
[3,] 1.3605 1.3615 1.3486 1.3516 1.3490 1.3509 1.3482 1.3349 1.3304 1.3427
[4,] 1.3588 1.3558 1.3549 1.3477 1.3448 1.3420 1.3532 1.3355 1.3314 1.3397
      [,31]  [,32]  [,33]  [,34]  [,35]  [,36]  [,37]  [,38]  [,39]  [,40]
[1,] 1.3441 1.3467 1.3618 1.3666 1.3781 1.3803 1.3722 1.3663 1.3794 1.3651
[2,] 1.3461 1.3505 1.3640 1.3753 1.3771 1.3821 1.3651 1.3664 1.3794 1.3591
[3,] 1.3460 1.3588 1.3596 1.3788 1.3779 1.3833 1.3659 1.3694 1.3729 1.3476
[4,] 1.3438 1.3601 1.3621 1.3782 1.3820 1.3743 1.3707 1.3818 1.3650 1.3405
      [,41]  [,42]  [,43]  [,44]  [,45]  [,46]  [,47]  [,48]  [,49]  [,50]
[1,] 1.3454 1.3574 1.3631 1.3580 1.3795 1.3860 1.4030 1.4127 1.4165 1.4089
[2,] 1.3476 1.3615 1.3610 1.3588 1.3824 1.3877 1.4049 1.4180 1.4195 1.4037
[3,] 1.3508 1.3658 1.3705 1.3669 1.3885 1.3867 1.4113 1.4179 1.4109 1.4146
[4,] 1.3493 1.3664 1.3632 1.3696 1.3897 1.3975 1.4106 1.4232 1.4136 1.4199
      [,51]  [,52]  [,53]  [,54]  [,55]  [,56]  [,57]  [,58]  [,59]  [,60]
[1,] 1.4173 1.4299 1.4230 1.4407 1.4488 1.4683 1.4639 1.4814 1.4874 1.4666
[2,] 1.4226 1.4288 1.4309 1.4447 1.4547 1.4579 1.4651 1.4829 1.4747 1.4741
[3,] 1.4150 1.4166 1.4384 1.4423 1.4722 1.4607 1.4654 1.4809 1.4738 1.4720
[4,] 1.4200 1.4254 1.4391 1.4479 1.4666 1.4700 1.4785 1.4845 1.4761 1.4554
      [,61]  [,62]  [,63]
[1,] 1.4649 1.4683 1.4385
[2,] 1.4718 1.4509 1.4349
[3,] 1.4672 1.4393 1.4380
[4,] 1.4675 1.4416 1.4398
> 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] 1.317275 1.299900 1.292400 1.296400 1.295125 1.299150 1.296450 1.301675
 [9] 1.313325 1.313625 1.320725 1.311750 1.317800 1.328750 1.330900 1.334125
[17] 1.336225 1.346075 1.356925 1.359850 1.360800 1.358675 1.352425 1.352625
[25] 1.345900 1.345575 1.347575 1.342175 1.331250 1.340775 1.345000 1.354025
[33] 1.361875 1.374725 1.378775 1.380000 1.368475 1.370975 1.374175 1.353075
[41] 1.348275 1.362775 1.364450 1.363325 1.385025 1.389475 1.407450 1.417950
[49] 1.415125 1.411775 1.418725 1.425175 1.432850 1.443900 1.460575 1.464225
[57] 1.468225 1.482425 1.478000 1.467025 1.467850 1.450025 1.437800
> arr.sd
 [1] 0.009162014 0.001587451 0.002778489 0.005278889 0.004845875 0.003371943
 [7] 0.003087070 0.005186762 0.001090489 0.002280899 0.003057641 0.003158586
[13] 0.002965356 0.004282133 0.003777124 0.001564981 0.000917878 0.005212405
[19] 0.002482438 0.003916206 0.002434475 0.003150000 0.002707243 0.004061506
[25] 0.002107131 0.003802959 0.004203471 0.008246767 0.002436528 0.001320038
[31] 0.001219289 0.006475273 0.001802544 0.005628129 0.002192981 0.003994997
[37] 0.003504640 0.007358612 0.006841235 0.011090649 0.002320022 0.004195533
[43] 0.004158926 0.005801939 0.004876730 0.005395291 0.004125126 0.004286801
[49] 0.003706189 0.007011122 0.003293807 0.006029027 0.007542988 0.003132624
[55] 0.010717703 0.005842017 0.006880589 0.001623525 0.006337718 0.008369140
[61] 0.002878078 0.013175071 0.002076857
> arr.range
 [1] 0.0186 0.0034 0.0060 0.0118 0.0104 0.0066 0.0066 0.0126 0.0026 0.0054
[11] 0.0067 0.0069 0.0063 0.0099 0.0086 0.0034 0.0021 0.0114 0.0052 0.0092
[21] 0.0055 0.0057 0.0063 0.0097 0.0046 0.0089 0.0096 0.0164 0.0058 0.0030
[31] 0.0023 0.0134 0.0044 0.0122 0.0049 0.0090 0.0071 0.0155 0.0144 0.0246
[41] 0.0054 0.0090 0.0095 0.0116 0.0102 0.0115 0.0083 0.0105 0.0086 0.0162
[51] 0.0076 0.0133 0.0161 0.0072 0.0234 0.0121 0.0146 0.0036 0.0136 0.0187
[61] 0.0069 0.0290 0.0049
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
   -0.01870      0.01688  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
       -7.144          4.980  

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

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

Coefficients:
(Intercept)     arr.mean  
   -0.04171      0.03751  

> postscript(file="/var/www/html/freestat/rcomp/tmp/1fjzy1210539307.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/freestat/rcomp/tmp/2y9fn1210539307.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/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/3k4821210539307.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/42bbc1210539307.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/5ejlz1210539307.tab") 
> 
> system("convert tmp/1fjzy1210539307.ps tmp/1fjzy1210539307.png")
> system("convert tmp/2y9fn1210539307.ps tmp/2y9fn1210539307.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.096   0.401   1.167