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.

> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = '0.5'
> par1 = '500'
> ylab = ''
> xlab = ''
> main = ''
> n <- as.numeric(par1)
> p <- as.numeric(par2)
> heads=rbinom(n-1,1,p)
> a=2*(heads)-1
> b=diffinv(a,xi=0)
> c=1:n
> pheads=(diffinv(heads,xi=.5))/c
> postscript(file="/var/www/html/freestat/rcomp/tmp/1v17m1228230535.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> op=par(mfrow=c(2,1))
> plot(c,b,type='n',main='Law of Averages',xlab='Toss Number',ylab='Excess of Heads',lwd=2,cex.lab=1.5,cex.main=2)
> lines(c,b,col='red')
> lines(c,rep(0,n),col='black')
> plot(c,pheads,type='n',xlab='Toss Number',ylab='Proportion of Heads',lwd=2,cex.lab=1.5)
> lines(c,pheads,col='blue')
> lines(c,rep(.5,n),col='black')
> par(op)
> dev.off()
null device 
          1 
> b
  [1]  0 -1  0 -1  0  1  2  3  4  5  6  7  6  5  6  5  6  7  8  7  8  7  6  7  6
 [26]  5  6  5  4  5  6  7  6  5  4  3  2  3  4  5  6  5  6  5  6  5  4  5  4  3
 [51]  2  1  2  3  2  1  0 -1  0  1  0  1  2  3  2  3  2  3  2  3  4  5  4  5  4
 [76]  5  6  7  8  9 10 11 12 13 14 15 16 17 16 17 16 17 16 17 18 19 18 17 16 17
[101] 16 17 18 19 18 19 20 19 18 19 20 21 20 19 20 21 20 21 20 19 20 21 20 19 20
[126] 21 22 23 24 23 24 25 26 27 28 29 30 29 30 31 30 31 32 33 34 33 32 33 34 33
[151] 32 31 32 31 30 31 30 29 30 29 30 29 30 29 28 27 26 27 28 27 28 27 28 29 30
[176] 31 30 31 32 31 32 33 32 33 32 31 30 31 32 33 32 33 34 35 34 33 32 33 34 35
[201] 34 33 34 33 34 35 36 37 36 37 38 39 40 39 38 39 38 37 38 37 38 39 38 39 40
[226] 39 38 37 36 35 34 35 36 37 38 39 40 41 42 41 40 41 40 39 38 37 36 35 36 37
[251] 38 37 36 35 34 35 34 33 32 31 32 31 30 29 30 29 28 27 26 27 28 29 30 31 32
[276] 31 30 31 32 31 30 29 28 27 28 27 26 27 28 27 28 27 28 29 28 27 26 25 26 27
[301] 28 27 26 25 24 25 26 27 28 27 26 25 26 27 26 25 24 25 26 27 26 25 24 25 24
[326] 23 24 25 26 25 24 25 24 23 22 23 24 25 24 25 24 23 24 25 26 25 24 25 24 25
[351] 24 23 22 21 20 21 22 21 22 21 22 23 22 21 22 21 20 19 20 21 20 21 20 19 18
[376] 17 16 15 14 15 16 17 16 17 18 17 18 19 18 19 18 17 16 17 18 17 18 19 20 21
[401] 20 19 18 17 16 17 18 17 18 17 18 19 20 19 20 19 18 17 18 17 18 17 16 17 18
[426] 19 18 17 18 19 20 21 20 21 22 23 22 23 22 23 24 25 24 25 24 25 26 27 26 27
[451] 28 29 30 31 30 31 32 33 34 35 34 33 32 31 30 31 32 33 32 31 30 31 32 33 32
[476] 33 34 33 34 33 32 33 32 33 32 31 30 29 28 29 28 27 28 27 28 27 26 25 24 23
> x <- b
> postscript(file="/var/www/html/freestat/rcomp/tmp/2ft141228230535.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> r <- spectrum(x,main='Raw Periodogram')
> dev.off()
null device 
          1 
> r
$freq
  [1] 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 0.018 0.020 0.022 0.024
 [13] 0.026 0.028 0.030 0.032 0.034 0.036 0.038 0.040 0.042 0.044 0.046 0.048
 [25] 0.050 0.052 0.054 0.056 0.058 0.060 0.062 0.064 0.066 0.068 0.070 0.072
 [37] 0.074 0.076 0.078 0.080 0.082 0.084 0.086 0.088 0.090 0.092 0.094 0.096
 [49] 0.098 0.100 0.102 0.104 0.106 0.108 0.110 0.112 0.114 0.116 0.118 0.120
 [61] 0.122 0.124 0.126 0.128 0.130 0.132 0.134 0.136 0.138 0.140 0.142 0.144
 [73] 0.146 0.148 0.150 0.152 0.154 0.156 0.158 0.160 0.162 0.164 0.166 0.168
 [85] 0.170 0.172 0.174 0.176 0.178 0.180 0.182 0.184 0.186 0.188 0.190 0.192
 [97] 0.194 0.196 0.198 0.200 0.202 0.204 0.206 0.208 0.210 0.212 0.214 0.216
[109] 0.218 0.220 0.222 0.224 0.226 0.228 0.230 0.232 0.234 0.236 0.238 0.240
[121] 0.242 0.244 0.246 0.248 0.250 0.252 0.254 0.256 0.258 0.260 0.262 0.264
[133] 0.266 0.268 0.270 0.272 0.274 0.276 0.278 0.280 0.282 0.284 0.286 0.288
[145] 0.290 0.292 0.294 0.296 0.298 0.300 0.302 0.304 0.306 0.308 0.310 0.312
[157] 0.314 0.316 0.318 0.320 0.322 0.324 0.326 0.328 0.330 0.332 0.334 0.336
[169] 0.338 0.340 0.342 0.344 0.346 0.348 0.350 0.352 0.354 0.356 0.358 0.360
[181] 0.362 0.364 0.366 0.368 0.370 0.372 0.374 0.376 0.378 0.380 0.382 0.384
[193] 0.386 0.388 0.390 0.392 0.394 0.396 0.398 0.400 0.402 0.404 0.406 0.408
[205] 0.410 0.412 0.414 0.416 0.418 0.420 0.422 0.424 0.426 0.428 0.430 0.432
[217] 0.434 0.436 0.438 0.440 0.442 0.444 0.446 0.448 0.450 0.452 0.454 0.456
[229] 0.458 0.460 0.462 0.464 0.466 0.468 0.470 0.472 0.474 0.476 0.478 0.480
[241] 0.482 0.484 0.486 0.488 0.490 0.492 0.494 0.496 0.498 0.500

$spec
  [1] 1.306305e+04 4.810343e+03 1.513047e+03 1.343228e+03 1.174777e+02
  [6] 2.573327e+00 1.651135e+02 2.285844e+02 8.281104e+01 1.356647e+02
 [11] 6.639729e+01 5.034508e-01 6.415991e+00 2.315442e+01 4.543163e+01
 [16] 4.277996e+01 3.645494e+00 6.815629e-01 3.347709e+00 9.639760e+00
 [21] 1.450842e+00 2.042361e+01 3.579887e+01 3.436282e+00 4.305061e+00
 [26] 6.983125e+00 1.429471e+01 1.216054e+01 3.867869e+00 5.561016e+00
 [31] 7.923944e+00 8.288675e+00 1.770867e+00 1.267452e+01 2.293895e-01
 [36] 7.582911e+00 1.251495e+01 9.938744e+00 1.435247e+00 5.023072e+00
 [41] 9.846493e+00 1.899864e-01 1.631651e+00 1.010125e+00 4.616225e-01
 [46] 9.473782e-01 7.324928e+00 4.016964e-02 2.056596e+00 3.040882e+00
 [51] 4.604716e+00 1.122274e+00 8.293863e-02 3.095029e+00 8.113798e-01
 [56] 1.668254e+00 2.574913e+00 4.400807e+00 3.634341e-01 6.559078e-02
 [61] 1.626423e+00 4.630315e+00 3.067828e+00 7.757485e-01 6.300650e-02
 [66] 2.504860e-01 1.541908e+00 6.050008e-01 7.138099e+00 8.537311e-01
 [71] 1.851003e+00 1.799059e+00 9.872915e-01 5.090623e-01 1.561091e+00
 [76] 4.842649e-01 1.016089e+00 1.105842e-01 3.831155e+00 3.370208e+00
 [81] 1.948094e+00 1.868926e+00 1.349561e+00 4.031773e-02 1.196815e+00
 [86] 1.592737e+00 4.508441e-01 1.429584e-02 2.947127e+00 9.446612e-01
 [91] 4.770207e-01 2.953624e-01 5.328869e-01 2.050514e+00 3.615283e-01
 [96] 6.689593e-01 2.805342e+00 2.790071e-01 1.376581e+00 4.613162e-01
[101] 7.485103e-02 6.239506e-02 1.207897e+00 4.411123e+00 1.278101e-01
[106] 7.279272e-01 6.626837e-02 1.522844e-01 4.498600e-01 8.329963e-01
[111] 1.417139e+00 1.713531e+00 7.282015e-01 3.800646e-03 7.656023e-02
[116] 2.647660e-03 7.444772e-01 3.445266e-01 3.940563e-02 1.076736e+00
[121] 2.951826e-01 9.768562e-01 2.799649e-01 6.427093e-01 2.687866e-01
[126] 4.852926e-01 9.643706e-02 4.038779e-01 8.475176e-02 8.849412e-02
[131] 2.904444e-01 1.598535e-01 5.056590e-01 4.685906e-01 2.347546e-01
[136] 3.611410e-01 4.435141e-01 6.799689e-02 9.425555e-01 5.678849e-01
[141] 4.165628e-01 1.833079e-01 2.527495e-01 1.723677e-01 9.474319e-02
[146] 1.386987e-01 4.623145e-02 1.170225e-01 1.448180e-01 2.366948e-01
[151] 6.677000e-01 1.838763e-01 3.000206e-02 5.882798e-01 1.980038e-02
[156] 2.376375e-01 2.180589e-01 7.410489e-01 2.391509e-01 2.021854e+00
[161] 6.898193e-02 2.271014e-01 1.165690e-02 1.067495e-01 2.440619e-01
[166] 1.921550e-01 3.592735e-03 6.733364e-02 8.234278e-01 6.466831e-01
[171] 9.273520e-02 1.556312e-02 4.257988e-01 1.019617e-01 6.034959e-01
[176] 4.675261e-01 3.425962e-02 1.790151e-01 1.324494e-01 4.352137e-01
[181] 2.309539e-01 6.193558e-01 1.147119e-01 6.898924e-02 4.644315e-01
[186] 1.087227e-01 3.504701e-01 2.340363e-01 4.490405e-01 1.140642e+00
[191] 3.169013e-01 6.144619e-01 1.934648e-01 7.878421e-02 2.535419e-02
[196] 6.249311e-01 8.413134e-02 3.907525e-02 2.879325e-01 9.118811e-02
[201] 1.608861e-01 2.840973e-01 3.681614e-01 5.313029e-01 1.210893e-01
[206] 2.986105e-01 5.816545e-02 9.016424e-02 9.960744e-02 6.850491e-03
[211] 2.449830e-01 2.895230e-01 3.497907e-01 1.362503e-01 2.080958e-02
[216] 5.941308e-02 2.189899e-01 1.173585e-01 1.362624e-01 2.608459e-01
[221] 2.255694e-02 2.664910e-01 3.730897e-02 1.611842e-02 1.303482e-01
[226] 9.244090e-02 6.717221e-02 1.499812e-01 9.536950e-02 8.133883e-02
[231] 1.157740e-01 1.036074e-01 2.581993e-02 3.501270e-01 4.827454e-02
[236] 3.570899e-02 8.703207e-02 6.021805e-01 1.050149e-01 3.967910e-01
[241] 1.510503e-01 7.712900e-01 1.777987e-01 9.769201e-02 4.625119e-02
[246] 1.896373e+00 1.315213e-01 8.823116e-02 3.125809e-01 5.414023e-01

$coh
NULL

$phase
NULL

$kernel
NULL

$df
[1] 1.791590

$bandwidth
[1] 0.0005773503

$n.used
[1] 500

$orig.n
[1] 500

$series
[1] "x"

$snames
NULL

$method
[1] "Raw Periodogram"

$taper
[1] 0.1

$pad
[1] 0

$detrend
[1] TRUE

$demean
[1] FALSE

attr(,"class")
[1] "spec"
> postscript(file="/var/www/html/freestat/rcomp/tmp/3focb1228230535.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> cpgram(x,main='Cumulative Periodogram')
> dev.off()
null device 
          1 
> 
> system("convert tmp/1v17m1228230535.ps tmp/1v17m1228230535.png")
> system("convert tmp/2ft141228230535.ps tmp/2ft141228230535.png")
> system("convert tmp/3focb1228230535.ps tmp/3focb1228230535.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.050   0.738   1.162