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(235.1
+ ,280.7
+ ,264.6
+ ,240.7
+ ,201.4
+ ,240.8
+ ,241.1
+ ,223.8
+ ,206.1
+ ,174.7
+ ,203.3
+ ,220.5
+ ,299.5
+ ,347.4
+ ,338.3
+ ,327.7
+ ,351.6
+ ,396.6
+ ,438.8
+ ,395.6
+ ,363.5
+ ,378.8
+ ,357
+ ,369
+ ,464.8
+ ,479.1
+ ,431.3
+ ,366.5
+ ,326.3
+ ,355.1
+ ,331.6
+ ,261.3
+ ,249
+ ,205.5
+ ,235.6
+ ,240.9
+ ,264.9
+ ,253.8
+ ,232.3
+ ,193.8
+ ,177
+ ,213.2
+ ,207.2
+ ,180.6
+ ,188.6
+ ,175.4
+ ,199
+ ,179.6
+ ,225.8
+ ,234
+ ,200.2
+ ,183.6
+ ,178.2
+ ,203.2
+ ,208.5
+ ,191.8
+ ,172.8
+ ,148
+ ,159.4
+ ,154.5
+ ,213.2
+ ,196.4
+ ,182.8
+ ,176.4
+ ,153.6
+ ,173.2
+ ,171
+ ,151.2
+ ,161.9
+ ,157.2
+ ,201.7
+ ,236.4
+ ,356.1
+ ,398.3
+ ,403.7
+ ,384.6
+ ,365.8
+ ,368.1
+ ,367.9
+ ,347
+ ,343.3
+ ,292.9
+ ,311.5
+ ,300.9
+ ,366.9
+ ,356.9
+ ,329.7
+ ,316.2
+ ,269
+ ,289.3
+ ,266.2
+ ,253.6
+ ,233.8
+ ,228.4
+ ,253.6
+ ,260.1
+ ,306.6
+ ,309.2
+ ,309.5
+ ,271
+ ,279.9
+ ,317.9
+ ,298.4
+ ,246.7
+ ,227.3
+ ,209.1
+ ,259.9
+ ,266
+ ,320.6
+ ,308.5
+ ,282.2
+ ,262.7
+ ,263.5
+ ,313.1
+ ,284.3
+ ,252.6
+ ,250.3
+ ,246.5
+ ,312.7
+ ,333.2
+ ,446.4
+ ,511.6
+ ,515.5
+ ,506.4
+ ,483.2
+ ,522.3
+ ,509.8
+ ,460.7
+ ,405.8
+ ,375
+ ,378.5
+ ,406.8
+ ,467.8
+ ,469.8
+ ,429.8
+ ,355.8
+ ,332.7
+ ,378
+ ,360.5
+ ,334.7
+ ,319.5
+ ,323.1
+ ,363.6
+ ,352.1
+ ,411.9
+ ,388.6
+ ,416.4
+ ,360.7
+ ,338
+ ,417.2
+ ,388.4
+ ,371.1
+ ,331.5
+ ,353.7
+ ,396.7
+ ,447
+ ,533.5
+ ,565.4
+ ,542.3
+ ,488.7
+ ,467.1
+ ,531.3
+ ,496.1
+ ,444
+ ,403.4
+ ,386.3
+ ,394.1
+ ,404.1
+ ,462.1
+ ,448.1
+ ,432.3
+ ,386.3
+ ,395.2
+ ,421.9
+ ,382.9
+ ,384.2
+ ,345.5
+ ,323.4
+ ,372.6
+ ,376
+ ,462.7
+ ,487
+ ,444.2
+ ,399.3
+ ,394.9
+ ,455.4
+ ,414
+ ,375.5
+ ,347
+ ,339.4
+ ,385.8
+ ,378.8
+ ,451.8
+ ,446.1
+ ,422.5
+ ,383.1
+ ,352.8
+ ,445.3
+ ,367.5
+ ,355.1
+ ,326.2
+ ,319.8
+ ,331.8
+ ,340.9
+ ,394.1
+ ,417.2
+ ,369.9
+ ,349.2
+ ,321.4
+ ,405.7
+ ,342.9
+ ,316.5
+ ,284.2
+ ,270.9
+ ,288.8
+ ,278.8
+ ,324.4
+ ,310.9
+ ,299
+ ,273
+ ,279.3
+ ,359.2
+ ,305
+ ,282.1
+ ,250.3
+ ,246.5
+ ,257.9
+ ,266.5
+ ,315.9
+ ,318.4
+ ,295.4
+ ,266.4
+ ,245.8
+ ,362.8
+ ,324.9
+ ,294.2
+ ,289.5
+ ,295.2
+ ,290.3
+ ,272
+ ,307.4
+ ,328.7
+ ,292.9
+ ,249.1
+ ,230.4
+ ,361.5
+ ,321.7
+ ,277.2
+ ,260.7
+ ,251
+ ,257.6
+ ,241.8
+ ,287.5
+ ,292.3
+ ,274.7
+ ,254.2
+ ,230
+ ,339
+ ,318.2
+ ,287
+ ,295.8
+ ,284
+ ,271
+ ,262.7
+ ,340.6
+ ,379.4
+ ,373.3
+ ,355.2
+ ,338.4
+ ,466.9
+ ,451
+ ,422
+ ,429.2
+ ,425.9
+ ,460.7
+ ,463.6
+ ,541.4
+ ,544.2
+ ,517.5
+ ,469.4
+ ,439.4
+ ,549
+ ,533
+ ,506.1
+ ,484
+ ,457
+ ,481.5
+ ,469.5
+ ,544.7
+ ,541.2
+ ,521.5
+ ,469.7
+ ,434.4
+ ,542.6
+ ,517.3
+ ,485.7
+ ,465.8
+ ,447
+ ,426.6
+ ,411.6
+ ,467.5
+ ,484.5
+ ,451.2
+ ,417.4
+ ,379.9
+ ,484.7
+ ,455
+ ,420.8
+ ,416.5
+ ,376.3
+ ,405.6
+ ,405.8
+ ,500.8
+ ,514
+ ,475.5
+ ,430.1
+ ,414.4
+ ,538
+ ,526
+ ,488.5
+ ,520.2
+ ,504.4
+ ,568.5
+ ,610.6
+ ,818
+ ,830.9
+ ,835.9
+ ,782
+ ,762.3
+ ,856.9
+ ,820.9
+ ,769.6
+ ,752.2
+ ,724.4
+ ,723.1
+ ,719.5
+ ,817.4
+ ,803.3
+ ,752.5
+ ,689
+ ,630.4
+ ,765.5
+ ,757.7
+ ,732.2
+ ,702.6
+ ,683.3
+ ,709.5
+ ,702.2
+ ,784.8
+ ,810.9
+ ,755.6
+ ,656.8
+ ,615.1
+ ,745.3
+ ,694.1
+ ,675.7
+ ,643.7
+ ,622.1
+ ,634.6
+ ,588
+ ,689.7
+ ,673.9
+ ,647.9
+ ,568.8
+ ,545.7
+ ,632.6
+ ,643.8
+ ,593.1
+ ,579.7
+ ,546
+ ,562.9
+ ,572.5)
> 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] 372
> (np <- floor(n / par1))
[1] 93
> 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,] 235.1 201.4 206.1 299.5 351.6 363.5 464.8 326.3 249.0 264.9 177.0 188.6
[2,] 280.7 240.8 174.7 347.4 396.6 378.8 479.1 355.1 205.5 253.8 213.2 175.4
[3,] 264.6 241.1 203.3 338.3 438.8 357.0 431.3 331.6 235.6 232.3 207.2 199.0
[4,] 240.7 223.8 220.5 327.7 395.6 369.0 366.5 261.3 240.9 193.8 180.6 179.6
     [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
[1,] 225.8 178.2 172.8 213.2 153.6 161.9 356.1 365.8 343.3 366.9 269.0 233.8
[2,] 234.0 203.2 148.0 196.4 173.2 157.2 398.3 368.1 292.9 356.9 289.3 228.4
[3,] 200.2 208.5 159.4 182.8 171.0 201.7 403.7 367.9 311.5 329.7 266.2 253.6
[4,] 183.6 191.8 154.5 176.4 151.2 236.4 384.6 347.0 300.9 316.2 253.6 260.1
     [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36]
[1,] 306.6 279.9 227.3 320.6 263.5 250.3 446.4 483.2 405.8 467.8 332.7 319.5
[2,] 309.2 317.9 209.1 308.5 313.1 246.5 511.6 522.3 375.0 469.8 378.0 323.1
[3,] 309.5 298.4 259.9 282.2 284.3 312.7 515.5 509.8 378.5 429.8 360.5 363.6
[4,] 271.0 246.7 266.0 262.7 252.6 333.2 506.4 460.7 406.8 355.8 334.7 352.1
     [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48]
[1,] 411.9 338.0 331.5 533.5 467.1 403.4 462.1 395.2 345.5 462.7 394.9 347.0
[2,] 388.6 417.2 353.7 565.4 531.3 386.3 448.1 421.9 323.4 487.0 455.4 339.4
[3,] 416.4 388.4 396.7 542.3 496.1 394.1 432.3 382.9 372.6 444.2 414.0 385.8
[4,] 360.7 371.1 447.0 488.7 444.0 404.1 386.3 384.2 376.0 399.3 375.5 378.8
     [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60]
[1,] 451.8 352.8 326.2 394.1 321.4 284.2 324.4 279.3 250.3 315.9 245.8 289.5
[2,] 446.1 445.3 319.8 417.2 405.7 270.9 310.9 359.2 246.5 318.4 362.8 295.2
[3,] 422.5 367.5 331.8 369.9 342.9 288.8 299.0 305.0 257.9 295.4 324.9 290.3
[4,] 383.1 355.1 340.9 349.2 316.5 278.8 273.0 282.1 266.5 266.4 294.2 272.0
     [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72]
[1,] 307.4 230.4 260.7 287.5 230.0 295.8 340.6 338.4 429.2 541.4 439.4 484.0
[2,] 328.7 361.5 251.0 292.3 339.0 284.0 379.4 466.9 425.9 544.2 549.0 457.0
[3,] 292.9 321.7 257.6 274.7 318.2 271.0 373.3 451.0 460.7 517.5 533.0 481.5
[4,] 249.1 277.2 241.8 254.2 287.0 262.7 355.2 422.0 463.6 469.4 506.1 469.5
     [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84]
[1,] 544.7 434.4 465.8 467.5 379.9 416.5 500.8 414.4 520.2 818.0 762.3 752.2
[2,] 541.2 542.6 447.0 484.5 484.7 376.3 514.0 538.0 504.4 830.9 856.9 724.4
[3,] 521.5 517.3 426.6 451.2 455.0 405.6 475.5 526.0 568.5 835.9 820.9 723.1
[4,] 469.7 485.7 411.6 417.4 420.8 405.8 430.1 488.5 610.6 782.0 769.6 719.5
     [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93]
[1,] 817.4 630.4 702.6 784.8 615.1 643.7 689.7 545.7 579.7
[2,] 803.3 765.5 683.3 810.9 745.3 622.1 673.9 632.6 546.0
[3,] 752.5 757.7 709.5 755.6 694.1 634.6 647.9 643.8 562.9
[4,] 689.0 732.2 702.2 656.8 675.7 588.0 568.8 593.1 572.5
> 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] 255.275 226.775 201.150 328.225 395.650 367.075 435.425 318.575 232.750
[10] 236.200 194.500 185.650 210.900 195.425 158.675 192.200 162.250 189.300
[19] 385.675 362.200 312.150 342.425 269.525 243.975 299.075 285.725 240.575
[28] 293.500 278.375 285.675 494.975 494.000 391.525 430.800 351.475 339.575
[37] 394.400 378.675 382.225 532.475 484.625 396.975 432.200 396.050 354.375
[46] 448.300 409.950 362.750 425.875 380.175 329.675 382.600 346.625 280.675
[55] 301.825 306.400 255.300 299.025 306.925 286.750 294.525 297.700 252.775
[64] 277.175 293.550 278.375 362.125 419.575 444.850 518.125 506.875 473.000
[73] 519.275 495.000 437.750 455.150 435.100 401.050 480.100 491.725 550.925
[82] 816.700 802.425 729.800 765.550 721.450 699.400 752.025 682.550 622.100
[91] 645.075 603.800 565.275
> arr.sd
 [1] 21.235564 18.749644 19.175940 20.773280 35.605383  9.227992 50.126731
 [8] 40.182200 18.983941 31.339113 18.352475 10.465658 23.204885 13.431896
[15] 10.510748 16.295603 11.451201 37.208153 21.292624 10.186592 22.119901
[22] 23.510476 14.787241 15.261799 18.761907 30.291734 26.998688 26.049568
[29] 26.624221 43.875002 32.597175 27.544267 17.125298 53.279139 21.748008
[36] 21.668006 25.557647 33.118915 50.964130 32.134548 37.717138  8.451183
[43] 32.932456 18.095027 24.752155 37.071373 34.134098 22.965118 31.210295
[50] 43.894068  8.946275 29.474622 41.018482  7.692150 21.839166 37.034668
[57]  8.843830 24.067180 49.474935 10.151026 33.663865 56.556697  8.360373
[64] 17.022999 47.451133 14.552978 17.650189 57.220181 20.056670 34.622476
[71] 48.340830 12.402957 34.591750 46.625815 23.667911 28.604254 45.120949
[78] 17.267793 36.963315 55.694726 48.230652 24.331461 44.707745 15.076472
[85] 58.148345 62.343484 11.244258 67.382113 53.750628 24.396858 53.690867
[92] 44.420491 14.576779
> arr.range
 [1]  45.6  39.7  45.8  47.9  87.2  21.8 112.6  93.8  43.5  71.1  36.2  23.6
[13]  50.4  30.3  24.8  36.8  22.0  79.2  47.6  21.1  50.4  50.7  35.7  31.7
[25]  38.5  71.2  56.9  57.9  60.5  86.7  69.1  61.6  31.8 114.0  45.3  44.1
[37]  55.7  79.2 115.5  76.7  87.3  17.8  75.8  39.0  52.6  87.7  79.9  46.4
[49]  68.7  92.5  21.1  68.0  89.2  17.9  51.4  79.9  20.0  52.0 117.0  23.2
[61]  79.6 131.1  18.9  38.1 109.0  33.1  38.8 128.5  37.7  74.8 109.6  27.0
[73]  75.0 108.2  54.2  67.1 104.8  40.2  83.9 123.6 106.2  53.9  94.6  32.7
[85] 128.4 135.1  26.2 154.1 130.2  55.7 120.9  98.1  33.7
> (lm1 <- lm(arr.sd~arr.mean))

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

Coefficients:
(Intercept)     arr.mean  
   11.39881      0.04515  

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

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

Coefficients:
  (Intercept)  log(arr.mean)  
      -0.6361         0.6552  

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

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

Coefficients:
(Intercept)     arr.mean  
    24.8897       0.1023  

> postscript(file="/var/www/html/rcomp/tmp/1gysi1228583423.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/2lko51228583423.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/3y9ad1228583423.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/46dgq1228583423.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/5i2ov1228583423.tab") 
> 
> system("convert tmp/1gysi1228583423.ps tmp/1gysi1228583423.png")
> system("convert tmp/2lko51228583423.ps tmp/2lko51228583423.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.672   0.312   0.809