R version 2.8.0 (2008-10-20)
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(358.59
+ ,362.96
+ ,362.42
+ ,364.97
+ ,364.04
+ ,361.06
+ ,358.48
+ ,352.96
+ ,359.59
+ ,360.39
+ ,357.40
+ ,362.93
+ ,364.55
+ ,365.73
+ ,364.70
+ ,364.65
+ ,359.43
+ ,362.14
+ ,356.97
+ ,354.82
+ ,353.17
+ ,357.06
+ ,356.18
+ ,355.01
+ ,355.65
+ ,357.31
+ ,357.07
+ ,357.91
+ ,358.48
+ ,358.97
+ ,351.77
+ ,352.16
+ ,359.08
+ ,360.35
+ ,359.53
+ ,359.30
+ ,358.41
+ ,359.68
+ ,355.31
+ ,357.08
+ ,349.71
+ ,354.13
+ ,345.49
+ ,341.69
+ ,344.25
+ ,340.17
+ ,342.47
+ ,344.43
+ ,333.23
+ ,339.72
+ ,342.61
+ ,346.36
+ ,339.09
+ ,339.73
+ ,341.12
+ ,335.94
+ ,333.46
+ ,335.66
+ ,341.12
+ ,342.21
+ ,342.62
+ ,346.06
+ ,344.43
+ ,346.65
+ ,343.74
+ ,335.67
+ ,342.75
+ ,341.77
+ ,345.84
+ ,346.52
+ ,350.79
+ ,345.44
+ ,345.87
+ ,338.48
+ ,337.21
+ ,340.81
+ ,339.86
+ ,342.86
+ ,343.33
+ ,341.73
+ ,351.38
+ ,351.13
+ ,345.99
+ ,347.55
+ ,346.02
+ ,345.29
+ ,347.03
+ ,348.01
+ ,345.48
+ ,349.40
+ ,351.05
+ ,349.70
+ ,350.86
+ ,354.45
+ ,355.30
+ ,357.48
+ ,355.24
+ ,351.79
+ ,355.22
+ ,351.02
+ ,350.28
+ ,350.17
+ ,348.16
+ ,340.30
+ ,343.75
+ ,344.71
+ ,344.13
+ ,342.14
+ ,345.04
+ ,346.02
+ ,346.43
+ ,347.07
+ ,339.33
+ ,339.10
+ ,337.19
+ ,339.58
+ ,327.85
+ ,326.81
+ ,321.73
+ ,320.45
+ ,327.69
+ ,323.95
+ ,320.47
+ ,322.13
+ ,316.34
+ ,314.78
+ ,308.90
+ ,308.62
+ ,314.41
+ ,306.88
+ ,310.60
+ ,321.60
+ ,321.50
+ ,325.68
+ ,324.35
+ ,320.01
+ ,326.88
+ ,332.39
+ ,331.48
+ ,332.62
+ ,324.79
+ ,327.12
+ ,328.91
+ ,328.37
+ ,324.83
+ ,325.90
+ ,326.18
+ ,328.94
+ ,333.78
+ ,328.06
+ ,325.87
+ ,325.41
+ ,318.86
+ ,319.13
+ ,310.16
+ ,311.73
+ ,306.54
+ ,311.16
+ ,311.98
+ ,306.72
+ ,308.05
+ ,300.76
+ ,301.90
+ ,293.09
+ ,292.76
+ ,294.58
+ ,289.90
+ ,296.69
+ ,297.21
+ ,293.31
+ ,296.25
+ ,298.60
+ ,296.87
+ ,301.02
+ ,304.73
+ ,301.92
+ ,295.72
+ ,293.18
+ ,298.35
+ ,297.99
+ ,299.85
+ ,299.85
+ ,304.45
+ ,299.45
+ ,298.14
+ ,298.78
+ ,297.02
+ ,301.33
+ ,294.96
+ ,296.69
+ ,300.73
+ ,301.96
+ ,297.38
+ ,293.87
+ ,285.96
+ ,285.41
+ ,283.70
+ ,284.76
+ ,277.11
+ ,274.73
+ ,274.73
+ ,274.73
+ ,274.73
+ ,274.69
+ ,275.42
+ ,264.15
+ ,276.24
+ ,268.88
+ ,277.97
+ ,280.49
+ ,281.09
+ ,276.16
+ ,272.58
+ ,270.94
+ ,284.31
+ ,283.94
+ ,284.18
+ ,282.83
+ ,283.84
+ ,282.71
+ ,279.29
+ ,280.70
+ ,274.47
+ ,273.44
+ ,275.49
+ ,279.46
+ ,280.19
+ ,288.21
+ ,284.80
+ ,281.41
+ ,283.39
+ ,287.97
+ ,290.77
+ ,290.60
+ ,289.67
+ ,289.84
+ ,298.55
+ ,296.07
+ ,297.14
+ ,295.34
+ ,296.25
+ ,294.30
+ ,296.15
+ ,296.49
+ ,298.05
+ ,301.03
+ ,300.52
+ ,301.50
+ ,296.93
+ ,289.84
+ ,291.44
+ ,286.88
+ ,286.74
+ ,288.93
+ ,292.19
+ ,295.39
+ ,295.86
+ ,293.36
+ ,292.86
+ ,292.73
+ ,296.73
+ ,285.02
+ ,285.24
+ ,288.62
+ ,283.36
+ ,285.84
+ ,291.48
+ ,291.41
+ ,287.77
+ ,284.97
+ ,286.05
+ ,278.19
+ ,281.21
+ ,277.92
+ ,280.08
+ ,269.24
+ ,268.48
+ ,268.83
+ ,269.54
+ ,262.37
+ ,265.12
+ ,265.34
+ ,263.32
+ ,267.18
+ ,260.75
+ ,261.78
+ ,257.27
+ ,255.63
+ ,251.39
+ ,259.49
+ ,261.18
+ ,261.65
+ ,262.01
+ ,265.23
+ ,268.10
+ ,262.27
+ ,263.59
+ ,257.85
+ ,265.69
+ ,271.15
+ ,266.69
+ ,265.77
+ ,262.32
+ ,270.48
+ ,273.03
+ ,269.13
+ ,280.65
+ ,282.75
+ ,281.44
+ ,281.99
+ ,282.86
+ ,287.21
+ ,283.11
+ ,280.66
+ ,282.39
+ ,280.83
+ ,284.71
+ ,279.99
+ ,283.50
+ ,284.88
+ ,288.60
+ ,284.80
+ ,287.20
+ ,286.22
+ ,286.54
+ ,279.58
+ ,283.08
+ ,288.88
+ ,280.18
+ ,284.16
+ ,290.57
+ ,286.82
+ ,273.00
+ ,278.69
+ ,264.54
+ ,271.92
+ ,283.60
+ ,269.25
+ ,263.58
+ ,264.16
+ ,268.85
+ ,269.67
+ ,249.41
+ ,268.99
+ ,268.65
+ ,260.16
+ ,256.55
+ ,251.47
+ ,234.93
+ ,232.96
+ ,215.49
+ ,213.68
+ ,236.07
+ ,235.41
+ ,214.77
+ ,225.85
+ ,224.64
+ ,238.26
+ ,232.44
+ ,222.50
+ ,225.28
+ ,220.49
+ ,216.86
+ ,234.70
+ ,230.06
+ ,238.27
+ ,238.56
+ ,242.70
+ ,249.14
+ ,234.89
+ ,227.78
+ ,234.04
+ ,230.70
+ ,230.17
+ ,218.23
+ ,232.20
+ ,220.76
+ ,215.60
+ ,217.69
+ ,204.35
+ ,191.44
+ ,203.84
+ ,211.86
+ ,210.57
+ ,219.57
+ ,219.98
+ ,226.01
+ ,207.04
+ ,212.52
+ ,217.92
+ ,210.45
+ ,218.53
+ ,223.32
+ ,218.76
+ ,217.63)
> 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] 395
> (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,] 358.59 364.55 355.65 358.41 333.23 342.62 345.87 346.02 355.24 345.04
 [2,] 362.96 365.73 357.31 359.68 339.72 346.06 338.48 345.29 351.79 346.02
 [3,] 362.42 364.70 357.07 355.31 342.61 344.43 337.21 347.03 355.22 346.43
 [4,] 364.97 364.65 357.91 357.08 346.36 346.65 340.81 348.01 351.02 347.07
 [5,] 364.04 359.43 358.48 349.71 339.09 343.74 339.86 345.48 350.28 339.33
 [6,] 361.06 362.14 358.97 354.13 339.73 335.67 342.86 349.40 350.17 339.10
 [7,] 358.48 356.97 351.77 345.49 341.12 342.75 343.33 351.05 348.16 337.19
 [8,] 352.96 354.82 352.16 341.69 335.94 341.77 341.73 349.70 340.30 339.58
 [9,] 359.59 353.17 359.08 344.25 333.46 345.84 351.38 350.86 343.75 327.85
[10,] 360.39 357.06 360.35 340.17 335.66 346.52 351.13 354.45 344.71 326.81
[11,] 357.40 356.18 359.53 342.47 341.12 350.79 345.99 355.30 344.13 321.73
[12,] 362.93 355.01 359.30 344.43 342.21 345.44 347.55 357.48 342.14 320.45
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,] 327.69 321.50 324.83 306.54 297.21 299.85 297.38 275.42 284.18 284.80
 [2,] 323.95 325.68 325.90 311.16 293.31 299.85 293.87 264.15 282.83 281.41
 [3,] 320.47 324.35 326.18 311.98 296.25 304.45 285.96 276.24 283.84 283.39
 [4,] 322.13 320.01 328.94 306.72 298.60 299.45 285.41 268.88 282.71 287.97
 [5,] 316.34 326.88 333.78 308.05 296.87 298.14 283.70 277.97 279.29 290.77
 [6,] 314.78 332.39 328.06 300.76 301.02 298.78 284.76 280.49 280.70 290.60
 [7,] 308.90 331.48 325.87 301.90 304.73 297.02 277.11 281.09 274.47 289.67
 [8,] 308.62 332.62 325.41 293.09 301.92 301.33 274.73 276.16 273.44 289.84
 [9,] 314.41 324.79 318.86 292.76 295.72 294.96 274.73 272.58 275.49 298.55
[10,] 306.88 327.12 319.13 294.58 293.18 296.69 274.73 270.94 279.46 296.07
[11,] 310.60 328.91 310.16 289.90 298.35 300.73 274.73 284.31 280.19 297.14
[12,] 321.60 328.37 311.73 296.69 297.99 301.96 274.69 283.94 288.21 295.34
       [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27]  [,28]  [,29]  [,30]
 [1,] 296.25 286.74 283.36 268.48 251.39 266.69 283.11 286.54 283.60 234.93
 [2,] 294.30 288.93 285.84 268.83 259.49 265.77 280.66 279.58 269.25 232.96
 [3,] 296.15 292.19 291.48 269.54 261.18 262.32 282.39 283.08 263.58 215.49
 [4,] 296.49 295.39 291.41 262.37 261.65 270.48 280.83 288.88 264.16 213.68
 [5,] 298.05 295.86 287.77 265.12 262.01 273.03 284.71 280.18 268.85 236.07
 [6,] 301.03 293.36 284.97 265.34 265.23 269.13 279.99 284.16 269.67 235.41
 [7,] 300.52 292.86 286.05 263.32 268.10 280.65 283.50 290.57 249.41 214.77
 [8,] 301.50 292.73 278.19 267.18 262.27 282.75 284.88 286.82 268.99 225.85
 [9,] 296.93 296.73 281.21 260.75 263.59 281.44 288.60 273.00 268.65 224.64
[10,] 289.84 285.02 277.92 261.78 257.85 281.99 284.80 278.69 260.16 238.26
[11,] 291.44 285.24 280.08 257.27 265.69 282.86 287.20 264.54 256.55 232.44
[12,] 286.88 288.62 269.24 255.63 271.15 287.21 286.22 271.92 251.47 222.50
       [,31]  [,32]
 [1,] 225.28 230.70
 [2,] 220.49 230.17
 [3,] 216.86 218.23
 [4,] 234.70 232.20
 [5,] 230.06 220.76
 [6,] 238.27 215.60
 [7,] 238.56 217.69
 [8,] 242.70 204.35
 [9,] 249.14 191.44
[10,] 234.89 203.84
[11,] 227.78 211.86
[12,] 234.04 210.57
> 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] 360.4825 359.5342 357.2983 349.4017 339.1875 344.3567 343.8500 350.0058
 [9] 348.0758 336.3833 316.3642 327.0083 323.2375 301.1775 297.9292 299.4342
[17] 281.8167 276.0142 280.4008 290.4625 295.7817 291.1392 283.1267 263.8008
[25] 262.4667 275.3600 283.9075 280.6633 264.5283 227.2500 232.7308 215.6175
> arr.sd
 [1]  3.350653  4.576091  2.792411  7.150406  3.950096  3.639861  4.630751
 [8]  4.015544  5.003521  9.737154  6.799237  4.034161  6.982641  7.692122
[15]  3.378979  2.569710  7.997813  6.111709  4.351839  5.531761  4.528750
[22]  4.111683  6.299986  4.479623  5.061969  8.341718  2.680960  7.725978
[29]  9.297536  9.013365  9.206917 12.216293
> arr.range
 [1] 12.01 12.56  8.58 19.51 13.13 15.12 14.17 12.19 14.94 26.62 20.81 12.61
[13] 23.62 22.08 11.55  9.49 22.69 20.16 14.77 17.14 14.62 11.71 22.24 13.91
[25] 19.76 24.89  8.61 26.03 34.19 24.58 32.28 40.76
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
   15.75697     -0.03282  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
       10.527         -1.551  

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

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

Coefficients:
(Intercept)     arr.mean  
    54.0357      -0.1172  

> postscript(file="/var/www/html/rcomp/tmp/1ge4u1229208378.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/2b7n61229208378.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/3k36r1229208378.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/4wj6f1229208378.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/5qfac1229208378.tab") 
> 
> system("convert tmp/1ge4u1229208378.ps tmp/1ge4u1229208378.png")
> system("convert tmp/2b7n61229208378.ps tmp/2b7n61229208378.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.578   0.293   0.656