R version 2.12.0 (2010-10-15)
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(315.42
+ ,316.32
+ ,316.49
+ ,317.56
+ ,318.13
+ ,318.00
+ ,316.39
+ ,314.66
+ ,313.68
+ ,313.18
+ ,314.66
+ ,315.43
+ ,316.27
+ ,316.81
+ ,317.42
+ ,318.87
+ ,319.87
+ ,319.43
+ ,318.01
+ ,315.75
+ ,314.00
+ ,313.68
+ ,314.84
+ ,316.03
+ ,316.73
+ ,317.54
+ ,318.38
+ ,319.31
+ ,320.42
+ ,319.61
+ ,318.42
+ ,316.64
+ ,314.83
+ ,315.15
+ ,315.95
+ ,316.85
+ ,317.78
+ ,318.40
+ ,319.53
+ ,320.41
+ ,320.85
+ ,320.45
+ ,319.44
+ ,317.25
+ ,316.12
+ ,315.27
+ ,316.53
+ ,317.53
+ ,318.58
+ ,318.92
+ ,319.70
+ ,321.22
+ ,322.08
+ ,321.31
+ ,319.58
+ ,317.61
+ ,316.05
+ ,315.83
+ ,316.91
+ ,318.20
+ ,319.41
+ ,320.07
+ ,320.74
+ ,321.40
+ ,322.06
+ ,321.73
+ ,320.27
+ ,318.54
+ ,316.54
+ ,316.71
+ ,317.53
+ ,318.55
+ ,319.27
+ ,320.28
+ ,320.73
+ ,321.97
+ ,322.00
+ ,321.71
+ ,321.05
+ ,318.71
+ ,317.65
+ ,317.14
+ ,318.71
+ ,319.25
+ ,320.46
+ ,321.43
+ ,322.22
+ ,323.54
+ ,323.91
+ ,323.59
+ ,322.26
+ ,320.21
+ ,318.48
+ ,317.94
+ ,319.63
+ ,320.87
+ ,322.17
+ ,322.34
+ ,322.88
+ ,324.25
+ ,324.83
+ ,323.93
+ ,322.39
+ ,320.76
+ ,319.10
+ ,319.23
+ ,320.56
+ ,321.80
+ ,322.40
+ ,322.99
+ ,323.73
+ ,324.86
+ ,325.41
+ ,325.19
+ ,323.97
+ ,321.92
+ ,320.10
+ ,319.96
+ ,320.97
+ ,322.48
+ ,323.52
+ ,323.89
+ ,325.04
+ ,326.01
+ ,326.67
+ ,325.96
+ ,325.13
+ ,322.90
+ ,321.61
+ ,321.01
+ ,322.08
+ ,323.37
+ ,324.34
+ ,325.30
+ ,326.29
+ ,327.54
+ ,327.54
+ ,327.21
+ ,325.98
+ ,324.42
+ ,322.91
+ ,322.90
+ ,323.85
+ ,324.96
+ ,326.01
+ ,326.51
+ ,327.01
+ ,327.62
+ ,328.76
+ ,328.40
+ ,327.20
+ ,325.28
+ ,323.20
+ ,323.40
+ ,324.64
+ ,325.85
+ ,326.60
+ ,327.47
+ ,327.58
+ ,329.56
+ ,329.90
+ ,328.92
+ ,327.89
+ ,326.17
+ ,324.68
+ ,325.04
+ ,326.34
+ ,327.39
+ ,328.37
+ ,329.40
+ ,330.14
+ ,331.33
+ ,332.31
+ ,331.90
+ ,330.70
+ ,329.15
+ ,327.34
+ ,327.02
+ ,327.99
+ ,328.48
+ ,329.18
+ ,330.55
+ ,331.32
+ ,332.48
+ ,332.92
+ ,332.08
+ ,331.02
+ ,329.24
+ ,327.28
+ ,327.21
+ ,328.29
+ ,329.41
+ ,330.23
+ ,331.24
+ ,331.87
+ ,333.14
+ ,333.80
+ ,333.42
+ ,331.73
+ ,329.90
+ ,328.40
+ ,328.17
+ ,329.32
+ ,330.59
+ ,331.58
+ ,332.39
+ ,333.33
+ ,334.41
+ ,334.71
+ ,334.17
+ ,332.88
+ ,330.77
+ ,329.14
+ ,328.77
+ ,330.14
+ ,331.52
+ ,332.75
+ ,333.25
+ ,334.53
+ ,335.90
+ ,336.57
+ ,336.10
+ ,334.76
+ ,332.59
+ ,331.41
+ ,330.98
+ ,332.24
+ ,333.68
+ ,334.80
+ ,335.22
+ ,336.47
+ ,337.59
+ ,337.84
+ ,337.72
+ ,336.37
+ ,334.51
+ ,332.60
+ ,332.37
+ ,333.75
+ ,334.79
+ ,336.05
+ ,336.59
+ ,337.79
+ ,338.71
+ ,339.30
+ ,339.12
+ ,337.56
+ ,335.92
+ ,333.74
+ ,333.70
+ ,335.13
+ ,336.56
+ ,337.84
+ ,338.19
+ ,339.90
+ ,340.60
+ ,341.29
+ ,341.00
+ ,339.39
+ ,337.43
+ ,335.72
+ ,335.84
+ ,336.93
+ ,338.04
+ ,339.06
+ ,340.30
+ ,341.21
+ ,342.33
+ ,342.74
+ ,342.07
+ ,340.32
+ ,338.27
+ ,336.52
+ ,336.68
+ ,338.19
+ ,339.44
+ ,340.57
+ ,341.44
+ ,342.53
+ ,343.39
+ ,343.96
+ ,343.18
+ ,341.88
+ ,339.65
+ ,337.80
+ ,337.69
+ ,339.09
+ ,340.32
+ ,341.20
+ ,342.35
+ ,342.93
+ ,344.77
+ ,345.58
+ ,345.14
+ ,343.81
+ ,342.22
+ ,339.69
+ ,339.82
+ ,340.98
+ ,342.82
+ ,343.52
+ ,344.33
+ ,345.11
+ ,346.88
+ ,347.25
+ ,346.61
+ ,345.22
+ ,343.11
+ ,340.90
+ ,341.17
+ ,342.80
+ ,344.04
+ ,344.79
+ ,345.82
+ ,347.25
+ ,348.17
+ ,348.75
+ ,348.07
+ ,346.38
+ ,344.52
+ ,342.92
+ ,342.63
+ ,344.06
+ ,345.38
+ ,346.12
+ ,346.79
+ ,347.69
+ ,349.38
+ ,350.04
+ ,349.38
+ ,347.78
+ ,345.75
+ ,344.70
+ ,344.01
+ ,345.50
+ ,346.75
+ ,347.86
+ ,348.32
+ ,349.26
+ ,350.84
+ ,351.70
+ ,351.11
+ ,349.37
+ ,347.97
+ ,346.31
+ ,346.22
+ ,347.68
+ ,348.82
+ ,350.29
+ ,351.58
+ ,352.08
+ ,353.45
+ ,354.08
+ ,353.66
+ ,352.25
+ ,350.30
+ ,348.58
+ ,348.74
+ ,349.93
+ ,351.21
+ ,352.62
+ ,352.93
+ ,353.54
+ ,355.27
+ ,355.52
+ ,354.97
+ ,353.74
+ ,351.51
+ ,349.63
+ ,349.82
+ ,351.12
+ ,352.35
+ ,353.47
+ ,354.51
+ ,355.18
+ ,355.98
+ ,356.94
+ ,355.99
+ ,354.58
+ ,352.68
+ ,350.72
+ ,350.92
+ ,352.55
+ ,353.91)
> 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] 384
> (np <- floor(n / par1))
[1] 32
> 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,] 315.42 316.27 316.73 317.78 318.58 319.41 319.27 320.46 322.17 322.40
 [2,] 316.32 316.81 317.54 318.40 318.92 320.07 320.28 321.43 322.34 322.99
 [3,] 316.49 317.42 318.38 319.53 319.70 320.74 320.73 322.22 322.88 323.73
 [4,] 317.56 318.87 319.31 320.41 321.22 321.40 321.97 323.54 324.25 324.86
 [5,] 318.13 319.87 320.42 320.85 322.08 322.06 322.00 323.91 324.83 325.41
 [6,] 318.00 319.43 319.61 320.45 321.31 321.73 321.71 323.59 323.93 325.19
 [7,] 316.39 318.01 318.42 319.44 319.58 320.27 321.05 322.26 322.39 323.97
 [8,] 314.66 315.75 316.64 317.25 317.61 318.54 318.71 320.21 320.76 321.92
 [9,] 313.68 314.00 314.83 316.12 316.05 316.54 317.65 318.48 319.10 320.10
[10,] 313.18 313.68 315.15 315.27 315.83 316.71 317.14 317.94 319.23 319.96
[11,] 314.66 314.84 315.95 316.53 316.91 317.53 318.71 319.63 320.56 320.97
[12,] 315.43 316.03 316.85 317.53 318.20 318.55 319.25 320.87 321.80 322.48
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,] 323.52 324.34 326.01 326.60 328.37 329.18 330.23 331.58 332.75 334.80
 [2,] 323.89 325.30 326.51 327.47 329.40 330.55 331.24 332.39 333.25 335.22
 [3,] 325.04 326.29 327.01 327.58 330.14 331.32 331.87 333.33 334.53 336.47
 [4,] 326.01 327.54 327.62 329.56 331.33 332.48 333.14 334.41 335.90 337.59
 [5,] 326.67 327.54 328.76 329.90 332.31 332.92 333.80 334.71 336.57 337.84
 [6,] 325.96 327.21 328.40 328.92 331.90 332.08 333.42 334.17 336.10 337.72
 [7,] 325.13 325.98 327.20 327.89 330.70 331.02 331.73 332.88 334.76 336.37
 [8,] 322.90 324.42 325.28 326.17 329.15 329.24 329.90 330.77 332.59 334.51
 [9,] 321.61 322.91 323.20 324.68 327.34 327.28 328.40 329.14 331.41 332.60
[10,] 321.01 322.90 323.40 325.04 327.02 327.21 328.17 328.77 330.98 332.37
[11,] 322.08 323.85 324.64 326.34 327.99 328.29 329.32 330.14 332.24 333.75
[12,] 323.37 324.96 325.85 327.39 328.48 329.41 330.59 331.52 333.68 334.79
       [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27]  [,28]  [,29]  [,30]
 [1,] 336.05 337.84 339.06 340.57 341.20 343.52 344.79 346.12 347.86 350.29
 [2,] 336.59 338.19 340.30 341.44 342.35 344.33 345.82 346.79 348.32 351.58
 [3,] 337.79 339.90 341.21 342.53 342.93 345.11 347.25 347.69 349.26 352.08
 [4,] 338.71 340.60 342.33 343.39 344.77 346.88 348.17 349.38 350.84 353.45
 [5,] 339.30 341.29 342.74 343.96 345.58 347.25 348.75 350.04 351.70 354.08
 [6,] 339.12 341.00 342.07 343.18 345.14 346.61 348.07 349.38 351.11 353.66
 [7,] 337.56 339.39 340.32 341.88 343.81 345.22 346.38 347.78 349.37 352.25
 [8,] 335.92 337.43 338.27 339.65 342.22 343.11 344.52 345.75 347.97 350.30
 [9,] 333.74 335.72 336.52 337.80 339.69 340.90 342.92 344.70 346.31 348.58
[10,] 333.70 335.84 336.68 337.69 339.82 341.17 342.63 344.01 346.22 348.74
[11,] 335.13 336.93 338.19 339.09 340.98 342.80 344.06 345.50 347.68 349.93
[12,] 336.56 338.04 339.44 340.32 342.82 344.04 345.38 346.75 348.82 351.21
       [,31]  [,32]
 [1,] 352.62 353.47
 [2,] 352.93 354.51
 [3,] 353.54 355.18
 [4,] 355.27 355.98
 [5,] 355.52 356.94
 [6,] 354.97 355.99
 [7,] 353.74 354.58
 [8,] 351.51 352.68
 [9,] 349.63 350.72
[10,] 349.82 350.92
[11,] 351.12 352.55
[12,] 352.35 353.91
> 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] 315.8267 316.7483 317.4858 318.2967 318.8325 319.4625 319.8725 321.2117
 [9] 322.0200 322.8317 323.9325 325.2700 326.1567 327.2950 329.5108 330.0817
[17] 330.9842 331.9842 333.7300 335.3358 336.6808 338.5142 339.7608 340.9583
[25] 342.6092 344.2450 345.7283 346.9908 348.7883 351.3458 352.7517 353.9525
> arr.sd
 [1] 1.615039 2.041291 1.775681 1.844684 2.047505 1.902960 1.663770 1.974256
 [9] 1.848749 1.877304 1.846432 1.665184 1.794392 1.640834 1.768386 1.947333
[17] 1.887343 2.012355 1.852522 1.894213 1.903797 1.915271 2.106062 2.128341
[25] 1.969058 2.085907 2.046189 1.920641 1.767298 1.842471 1.983545 1.974027
> arr.range
 [1] 4.95 6.19 5.59 5.58 6.25 5.52 4.86 5.97 5.73 5.45 5.66 4.64 5.56 5.22 5.29
[16] 5.71 5.63 5.94 5.59 5.47 5.60 5.57 6.22 6.27 5.89 6.35 6.12 6.03 5.48 5.50
[31] 5.89 6.22
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
   0.430656     0.004399  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
      -3.9902         0.7968  

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

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

Coefficients:
(Intercept)     arr.mean  
    1.12293      0.01374  

> postscript(file="/var/www/rcomp/tmp/181091293661033.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/2isiu1293661033.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/34byi1293661033.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/4pbf61293661033.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/5sudc1293661033.tab") 
> 
> try(system("convert tmp/181091293661033.ps tmp/181091293661033.png",intern=TRUE))
character(0)
> try(system("convert tmp/2isiu1293661033.ps tmp/2isiu1293661033.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.680   0.400   1.051