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(10
+ ,10
+ ,10
+ ,10
+ ,10
+ ,9.94
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.06
+ ,9.94
+ ,9.94
+ ,9.94
+ ,9.94
+ ,9.94
+ ,9.94
+ ,10.06
+ ,10.06
+ ,9.94
+ ,10.06
+ ,10.06
+ ,10.06
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.28
+ ,10.18
+ ,10.18
+ ,10.28
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.18
+ ,10.28
+ ,10.28
+ ,10.28
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.34
+ ,10.34
+ ,10.34
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.34
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.34
+ ,10.34
+ ,10.34
+ ,10.34
+ ,10.42
+ ,10.42
+ ,10.34
+ ,10.34
+ ,10.34
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.34
+ ,10.34
+ ,10.34
+ ,10.34
+ ,10.34
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.34
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.34
+ ,10.34
+ ,10.34
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.42
+ ,10.34
+ ,10.34
+ ,10.42
+ ,10.42
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.49
+ ,10.49
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.49
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.49
+ ,10.49
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.49
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.49
+ ,10.49
+ ,10.49
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.49
+ ,10.55
+ ,10.55
+ ,10.55
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.57
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.57
+ ,10.57
+ ,10.57
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.57
+ ,10.63
+ ,10.63
+ ,10.57
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.63
+ ,10.57
+ ,10.63
+ ,10.6
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.6
+ ,10.6
+ ,10.6
+ ,10.6
+ ,10.7
+ ,10.7
+ ,10.6
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.8
+ ,10.8
+ ,10.8
+ ,10.8
+ ,10.7
+ ,10.7
+ ,10.7
+ ,10.85
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.75
+ ,10.75
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.85
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75
+ ,10.75)
> par1 = '6'
> #'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] 350
> (np <- floor(n / par1))
[1] 58
> 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] [,11] [,12]
[1,] 10.00 10.06 10.06 10.06  9.94 10.06 10.28 10.18 10.18 10.28 10.18 10.28
[2,] 10.00 10.06 10.06 10.06  9.94 10.06 10.28 10.18 10.18 10.28 10.18 10.28
[3,] 10.00 10.06 10.06  9.94 10.06 10.18 10.28 10.18 10.18 10.18 10.28 10.28
[4,] 10.00 10.06 10.06  9.94 10.06 10.28 10.18 10.28 10.18 10.18 10.18 10.28
[5,] 10.00 10.06 10.06  9.94  9.94 10.28 10.28 10.28 10.18 10.18 10.18 10.28
[6,]  9.94 10.06 10.06  9.94 10.06 10.18 10.28 10.18 10.28 10.28 10.18 10.28
     [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
[1,] 10.18 10.18 10.28 10.28 10.28 10.28 10.28 10.42 10.42 10.42 10.42 10.42
[2,] 10.18 10.18 10.28 10.28 10.28 10.28 10.28 10.42 10.34 10.42 10.42 10.42
[3,] 10.18 10.18 10.28 10.28 10.18 10.28 10.28 10.42 10.34 10.34 10.42 10.42
[4,] 10.18 10.28 10.28 10.28 10.18 10.28 10.42 10.42 10.34 10.42 10.42 10.42
[5,] 10.18 10.28 10.28 10.18 10.28 10.28 10.42 10.42 10.42 10.42 10.42 10.42
[6,] 10.18 10.28 10.18 10.28 10.28 10.18 10.42 10.42 10.42 10.42 10.42 10.42
     [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36]
[1,] 10.42 10.42 10.42 10.34 10.42 10.42 10.34 10.42 10.42 10.34 10.55 10.55
[2,] 10.42 10.42 10.42 10.42 10.42 10.42 10.34 10.42 10.42 10.42 10.55 10.49
[3,] 10.34 10.34 10.34 10.42 10.42 10.42 10.34 10.42 10.42 10.42 10.55 10.49
[4,] 10.34 10.34 10.34 10.42 10.42 10.42 10.42 10.42 10.42 10.55 10.55 10.55
[5,] 10.34 10.34 10.34 10.34 10.42 10.42 10.42 10.42 10.42 10.55 10.55 10.55
[6,] 10.34 10.42 10.34 10.42 10.42 10.42 10.42 10.42 10.34 10.55 10.55 10.55
     [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48]
[1,] 10.55 10.55 10.55 10.55 10.55 10.55 10.63 10.63 10.63 10.63 10.57  10.7
[2,] 10.55 10.55 10.55 10.55 10.49 10.49 10.63 10.63 10.63 10.57 10.63  10.7
[3,] 10.55 10.55 10.55 10.55 10.49 10.55 10.63 10.63 10.63 10.63 10.60  10.7
[4,] 10.49 10.55 10.49 10.55 10.49 10.55 10.63 10.57 10.63 10.63 10.70  10.7
[5,] 10.55 10.49 10.55 10.55 10.55 10.55 10.57 10.57 10.57 10.63 10.70  10.7
[6,] 10.55 10.49 10.55 10.55 10.55 10.63 10.63 10.57 10.63 10.63 10.70  10.7
     [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58]
[1,]  10.7  10.7  10.7  10.7 10.80 10.75 10.75 10.85 10.85 10.85
[2,]  10.6  10.6  10.7  10.7 10.80 10.75 10.75 10.85 10.85 10.85
[3,]  10.6  10.7  10.7  10.7 10.70 10.75 10.75 10.85 10.85 10.85
[4,]  10.6  10.7  10.7  10.7 10.70 10.75 10.75 10.85 10.75 10.75
[5,]  10.6  10.7  10.7  10.8 10.70 10.75 10.75 10.85 10.75 10.75
[6,]  10.7  10.7  10.7  10.8 10.85 10.75 10.85 10.85 10.85 10.75
> 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]  9.99000 10.06000 10.06000  9.98000 10.00000 10.17333 10.26333 10.21333
 [9] 10.19667 10.23000 10.19667 10.28000 10.18000 10.23000 10.26333 10.26333
[17] 10.24667 10.26333 10.35000 10.42000 10.38000 10.40667 10.42000 10.42000
[25] 10.36667 10.38000 10.36667 10.39333 10.42000 10.42000 10.38000 10.42000
[33] 10.40667 10.47167 10.55000 10.53000 10.54000 10.53000 10.54000 10.55000
[41] 10.52000 10.55333 10.62000 10.60000 10.62000 10.62000 10.65000 10.70000
[49] 10.63333 10.68333 10.70000 10.73333 10.75833 10.75000 10.76667 10.85000
[57] 10.81667 10.80000
> arr.sd
 [1] 0.02449490 0.00000000 0.00000000 0.06196773 0.06572671 0.09852242
 [7] 0.04082483 0.05163978 0.04082483 0.05477226 0.04082483 0.00000000
[13] 0.00000000 0.05477226 0.04082483 0.04082483 0.05163978 0.04082483
[19] 0.07668116 0.00000000 0.04381780 0.03265986 0.00000000 0.00000000
[25] 0.04131182 0.04381780 0.04131182 0.04131182 0.00000000 0.00000000
[31] 0.04381780 0.00000000 0.03265986 0.09064583 0.00000000 0.03098387
[37] 0.02449490 0.03098387 0.02449490 0.00000000 0.03286335 0.04457204
[43] 0.02449490 0.03286335 0.02449490 0.02449490 0.05796551 0.00000000
[49] 0.05163978 0.04082483 0.00000000 0.05163978 0.06645801 0.00000000
[55] 0.04082483 0.00000000 0.05163978 0.05477226
> arr.range
 [1] 0.06 0.00 0.00 0.12 0.12 0.22 0.10 0.10 0.10 0.10 0.10 0.00 0.00 0.10 0.10
[16] 0.10 0.10 0.10 0.14 0.00 0.08 0.08 0.00 0.00 0.08 0.08 0.08 0.08 0.00 0.00
[31] 0.08 0.00 0.08 0.21 0.00 0.06 0.06 0.06 0.06 0.00 0.06 0.14 0.06 0.06 0.06
[46] 0.06 0.13 0.00 0.10 0.10 0.00 0.10 0.15 0.00 0.10 0.00 0.10 0.10
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
    0.17970     -0.01407  

> (lnlm1 <- lm(log(arr.sd)~log(arr.mean)))
Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 
  NA/NaN/Inf in foreign function call (arg 4)
Calls: lm -> lm.fit -> .Fortran
Execution halted