R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(122.36
+ ,123.33
+ ,123.04
+ ,124.53
+ ,125.13
+ ,125.85
+ ,126.50
+ ,126.53
+ ,127.07
+ ,124.55
+ ,124.90
+ ,124.32
+ ,122.84
+ ,123.31
+ ,123.31
+ ,124.87
+ ,124.64
+ ,124.73
+ ,124.90
+ ,124.04
+ ,123.28
+ ,123.86
+ ,122.29
+ ,124.09
+ ,124.54
+ ,125.65
+ ,125.70
+ ,125.53
+ ,125.61
+ ,125.55
+ ,125.41
+ ,127.60
+ ,124.68
+ ,124.41
+ ,126.43
+ ,126.38
+ ,125.78
+ ,124.70
+ ,125.07
+ ,125.25
+ ,126.58
+ ,127.13
+ ,125.82
+ ,123.70
+ ,124.39
+ ,123.70
+ ,124.42
+ ,121.05
+ ,121.02
+ ,123.23
+ ,121.32
+ ,120.91
+ ,120.72
+ ,123.31
+ ,119.58
+ ,119.53
+ ,120.59
+ ,118.63
+ ,118.47
+ ,111.81
+ ,114.71
+ ,117.34
+ ,115.77
+ ,118.38
+ ,117.84
+ ,118.83
+ ,120.02
+ ,116.21
+ ,117.08
+ ,120.20
+ ,119.83
+ ,118.92
+ ,118.03
+ ,117.71
+ ,119.55
+ ,116.13
+ ,115.97
+ ,115.99
+ ,114.96
+ ,116.46
+ ,116.55
+ ,113.05
+ ,117.44
+ ,118.84
+ ,117.06
+ ,117.54
+ ,119.31
+ ,118.72
+ ,121.55
+ ,122.61
+ ,121.53
+ ,123.31
+ ,124.07
+ ,123.59
+ ,122.97
+ ,123.22
+ ,123.04
+ ,122.96
+ ,122.81
+ ,122.81
+ ,122.62
+ ,120.82
+ ,119.41
+ ,121.56
+ ,121.59
+ ,118.50
+ ,118.77
+ ,118.86
+ ,117.60
+ ,119.90
+ ,121.83
+ ,121.84
+ ,122.12
+ ,122.12
+ ,121.36
+ ,119.66
+ ,119.32
+ ,120.36
+ ,117.06
+ ,117.48
+ ,115.60
+ ,113.86
+ ,116.92
+ ,117.75
+ ,117.75
+ ,115.31
+ ,116.28
+ ,115.22
+ ,115.65
+ ,115.11
+ ,118.67
+ ,118.04
+ ,116.50
+ ,119.78
+ ,119.95
+ ,120.37
+ ,119.79
+ ,119.43
+ ,121.06
+ ,121.74
+ ,121.09
+ ,122.97
+ ,120.50
+ ,117.18
+ ,115.03
+ ,113.36
+ ,112.59
+ ,111.65
+ ,111.98
+ ,114.87
+ ,114.67
+ ,114.09
+ ,114.77
+ ,117.05
+ ,117.22
+ ,113.18
+ ,110.95
+ ,112.14
+ ,112.72
+ ,110.01
+ ,110.29
+ ,110.74
+ ,110.32
+ ,105.89
+ ,108.97
+ ,109.34
+ ,106.57
+ ,99.49
+ ,101.81
+ ,104.29
+ ,109.73
+ ,105.06
+ ,107.97
+ ,108.13
+ ,109.86
+ ,108.95
+ ,111.20
+ ,110.69
+ ,106.10
+ ,105.68
+ ,104.12
+ ,104.71
+ ,104.30
+ ,103.52
+ ,107.76
+ ,107.80
+ ,107.30
+ ,108.64
+ ,105.03
+ ,108.30
+ ,107.21
+ ,109.27
+ ,109.50
+ ,111.68
+ ,111.80
+ ,111.75
+ ,106.68
+ ,106.37
+ ,105.76
+ ,109.01
+ ,109.01
+ ,109.01
+ ,109.01
+ ,107.69
+ ,105.19
+ ,105.48
+ ,102.22
+ ,100.54
+ ,105.00
+ ,105.44
+ ,107.89
+ ,108.64
+ ,106.70
+ ,109.10
+ ,105.23
+ ,108.41
+ ,108.80
+ ,110.39
+ ,110.22
+ ,110.86
+ ,108.58
+ ,107.70
+ ,106.62
+ ,109.84
+ ,107.16
+ ,107.26
+ ,108.70
+ ,109.85
+ ,109.41
+ ,112.36
+ ,111.03
+ ,110.67
+ ,109.21
+ ,113.58
+ ,113.88
+ ,114.08
+ ,112.33
+ ,113.92
+ ,114.41
+ ,114.57
+ ,115.35
+ ,113.13
+ ,113.29
+ ,112.56
+ ,113.06
+ ,113.46
+ ,115.39
+ ,116.62
+ ,117.04
+ ,117.42
+ ,115.62
+ ,115.16
+ ,115.69
+ ,112.85
+ ,114.05
+ ,112.00
+ ,113.74
+ ,116.26
+ ,118.63
+ ,116.49
+ ,118.23
+ ,116.83
+ ,118.82
+ ,114.36
+ ,112.02
+ ,113.24
+ ,109.75
+ ,110.33
+ ,112.86
+ ,113.04
+ ,113.80
+ ,110.90
+ ,109.96
+ ,108.69
+ ,108.84
+ ,108.47
+ ,108.07
+ ,107.94
+ ,108.11
+ ,108.11
+ ,106.81
+ ,105.58
+ ,105.61
+ ,106.52
+ ,103.86
+ ,104.60
+ ,104.73
+ ,105.12
+ ,104.76
+ ,103.85
+ ,103.83
+ ,103.22
+ ,101.64
+ ,102.13
+ ,104.33
+ ,104.92
+ ,107.78
+ ,104.49
+ ,102.80
+ ,102.86
+ ,104.51
+ ,104.73
+ ,102.58
+ ,99.93
+ ,101.41
+ ,101.05
+ ,99.86
+ ,101.11
+ ,100.89
+ ,101.09
+ ,98.31
+ ,98.08
+ ,99.55
+ ,99.62
+ ,97.37
+ ,98.16
+ ,97.98
+ ,98.15
+ ,97.10
+ ,97.24
+ ,96.70
+ ,96.64
+ ,100.65
+ ,96.75
+ ,97.74
+ ,97.92
+ ,98.34
+ ,93.84
+ ,97.80
+ ,96.20
+ ,95.99
+ ,95.18
+ ,95.95
+ ,92.23
+ ,91.78
+ ,92.97
+ ,89.76
+ ,92.88
+ ,96.23
+ ,95.79
+ ,93.97
+ ,93.90
+ ,93.60
+ ,93.96
+ ,88.69
+ ,88.57
+ ,85.62
+ ,86.25
+ ,85.33
+ ,83.33
+ ,77.78
+ ,78.70
+ ,72.05
+ ,80.75
+ ,81.41
+ ,82.65
+ ,75.85
+ ,75.70
+ ,78.25
+ ,77.41
+ ,76.84
+ ,74.25
+ ,74.95
+ ,68.78
+ ,73.21
+ ,73.26
+ ,78.67
+ ,75.63
+ ,74.99
+ ,83.87
+ ,79.62
+ ,80.13
+ ,79.76
+ ,78.20
+ ,78.05
+ ,79.05
+ ,73.32
+ ,75.17
+ ,73.26
+ ,73.72
+ ,73.57
+ ,70.60
+ ,71.25
+ ,74.22
+ ,73.32
+ ,73.01
+ ,74.21
+ ,75.32
+ ,71.73
+ ,71.94
+ ,72.94
+ ,72.47
+ ,71.94
+ ,74.30
+ ,74.30)
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = '12'
> ylab = ''
> xlab = ''
> main = ''
> #'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,] 122.36 122.84 124.54 125.78 121.02 114.71 118.03 117.06 123.04 117.60
 [2,] 123.33 123.31 125.65 124.70 123.23 117.34 117.71 117.54 122.96 119.90
 [3,] 123.04 123.31 125.70 125.07 121.32 115.77 119.55 119.31 122.81 121.83
 [4,] 124.53 124.87 125.53 125.25 120.91 118.38 116.13 118.72 122.81 121.84
 [5,] 125.13 124.64 125.61 126.58 120.72 117.84 115.97 121.55 122.62 122.12
 [6,] 125.85 124.73 125.55 127.13 123.31 118.83 115.99 122.61 120.82 122.12
 [7,] 126.50 124.90 125.41 125.82 119.58 120.02 114.96 121.53 119.41 121.36
 [8,] 126.53 124.04 127.60 123.70 119.53 116.21 116.46 123.31 121.56 119.66
 [9,] 127.07 123.28 124.68 124.39 120.59 117.08 116.55 124.07 121.59 119.32
[10,] 124.55 123.86 124.41 123.70 118.63 120.20 113.05 123.59 118.50 120.36
[11,] 124.90 122.29 126.43 124.42 118.47 119.83 117.44 122.97 118.77 117.06
[12,] 124.32 124.09 126.38 121.05 111.81 118.92 118.84 123.22 118.86 117.48
       [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
 [1,] 115.60 116.50 115.03 110.95 101.81 104.12 109.50 105.19 108.80 109.41
 [2,] 113.86 119.78 113.36 112.14 104.29 104.71 111.68 105.48 110.39 112.36
 [3,] 116.92 119.95 112.59 112.72 109.73 104.30 111.80 102.22 110.22 111.03
 [4,] 117.75 120.37 111.65 110.01 105.06 103.52 111.75 100.54 110.86 110.67
 [5,] 117.75 119.79 111.98 110.29 107.97 107.76 106.68 105.00 108.58 109.21
 [6,] 115.31 119.43 114.87 110.74 108.13 107.80 106.37 105.44 107.70 113.58
 [7,] 116.28 121.06 114.67 110.32 109.86 107.30 105.76 107.89 106.62 113.88
 [8,] 115.22 121.74 114.09 105.89 108.95 108.64 109.01 108.64 109.84 114.08
 [9,] 115.65 121.09 114.77 108.97 111.20 105.03 109.01 106.70 107.16 112.33
[10,] 115.11 122.97 117.05 109.34 110.69 108.30 109.01 109.10 107.26 113.92
[11,] 118.67 120.50 117.22 106.57 106.10 107.21 109.01 105.23 108.70 114.41
[12,] 118.04 117.18 113.18  99.49 105.68 109.27 107.69 108.41 109.85 114.57
       [,21]  [,22]  [,23]  [,24]  [,25]  [,26]  [,27] [,28] [,29] [,30] [,31]
 [1,] 115.35 115.69 112.02 108.07 104.76 104.51  99.55 97.74 89.76 85.33 76.84
 [2,] 113.13 112.85 113.24 107.94 103.85 104.73  99.62 97.92 92.88 83.33 74.25
 [3,] 113.29 114.05 109.75 108.11 103.83 102.58  97.37 98.34 96.23 77.78 74.95
 [4,] 112.56 112.00 110.33 108.11 103.22  99.93  98.16 93.84 95.79 78.70 68.78
 [5,] 113.06 113.74 112.86 106.81 101.64 101.41  97.98 97.80 93.97 72.05 73.21
 [6,] 113.46 116.26 113.04 105.58 102.13 101.05  98.15 96.20 93.90 80.75 73.26
 [7,] 115.39 118.63 113.80 105.61 104.33  99.86  97.10 95.99 93.60 81.41 78.67
 [8,] 116.62 116.49 110.90 106.52 104.92 101.11  97.24 95.18 93.96 82.65 75.63
 [9,] 117.04 118.23 109.96 103.86 107.78 100.89  96.70 95.95 88.69 75.85 74.99
[10,] 117.42 116.83 108.69 104.60 104.49 101.09  96.64 92.23 88.57 75.70 83.87
[11,] 115.62 118.82 108.84 104.73 102.80  98.31 100.65 91.78 85.62 78.25 79.62
[12,] 115.16 114.36 108.47 105.12 102.86  98.08  96.75 92.97 86.25 77.41 80.13
      [,32]
 [1,] 79.76
 [2,] 78.20
 [3,] 78.05
 [4,] 79.05
 [5,] 73.32
 [6,] 75.17
 [7,] 73.26
 [8,] 73.72
 [9,] 73.57
[10,] 70.60
[11,] 71.25
[12,] 74.22
> 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] 124.84250 123.84667 125.62417 124.79917 119.92667 117.92750 116.72333
 [8] 121.29000 121.14583 120.05417 116.34667 120.03000 114.20500 108.95250
[15] 107.45583 106.49667 108.93917 105.82000 108.83167 112.45417 114.84167
[22] 115.66250 110.99167 106.25500 103.88417 101.12917  97.99250  95.49500
[29]  91.60167  79.10083  76.18333  75.01417
> arr.sd
 [1] 1.4713452 0.8526358 0.8909898 1.5862618 2.9770801 1.7631229 1.7535488
 [8] 2.4862239 1.8077935 1.8831620 1.4658806 1.7878021 1.7762908 3.5810261
[15] 2.8764197 2.0150314 2.0737294 2.5760082 1.4155682 1.9405222 1.6931940
[22] 2.2781616 1.9311788 1.5478519 1.6058499 2.0657178 1.3136979 2.3150123
[29] 3.6497244 3.7498229 3.9574035 3.0479992
> arr.range
 [1]  4.71  2.61  3.19  6.08 11.50  5.49  6.50  7.01  4.54  5.06  4.81  6.47
[13]  5.57 13.23  9.39  5.75  6.04  8.56  4.24  5.36  4.86  6.82  5.33  4.25
[25]  6.14  6.65  4.01  6.56 10.61 13.28 15.09  9.16
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
    6.40459     -0.03916  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
        8.603         -1.689  

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

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

Coefficients:
(Intercept)     arr.mean  
    21.6807      -0.1359  

> postscript(file="/var/www/html/rcomp/tmp/1qhxc1260967140.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/27z8o1260967140.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/349371260967140.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/4skoo1260967140.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/5i1pv1260967140.tab") 
> 
> try(system("convert tmp/1qhxc1260967140.ps tmp/1qhxc1260967140.png",intern=TRUE))
character(0)
> try(system("convert tmp/27z8o1260967140.ps tmp/27z8o1260967140.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.566   0.304   0.796