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/rcomp/tmp/1yzwh1228068352.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  2  3  4  5  4  5  6  7  8  7  8  7  6  5  6  5  6  7  8  7  6
 [26]  5  6  5  6  7  8  9 10  9 10  9 10 11 10 11 12 13 14 15 14 13 12 13 14 15
 [51] 14 15 16 17 18 19 20 19 18 19 18 17 16 15 16 17 18 17 16 15 16 15 14 13 14
 [76] 13 14 13 12 11 12 11 10 11 10  9 10  9  8  9 10  9 10 11 12 11 12 13 12 13
[101] 14 13 14 15 14 15 14 15 16 17 16 15 14 15 16 15 16 17 16 17 16 15 14 15 14
[126] 15 14 15 14 13 14 15 14 13 12 11 12 11 10 11 10 11 10  9 10 11 12 13 12 11
[151] 10  9  8  9  8  7  6  5  4  5  6  5  6  7  6  7  6  7  8  7  8  9  8  9 10
[176]  9  8  7  6  7  6  5  6  7  8  9  8  7  6  5  6  7  8  9  8  7  6  7  8  9
[201] 10  9  8  9 10  9  8  7  8  7  6  7  8  9  8  9 10 11 12 11 12 11 10 11 10
[226] 11 12 13 12 13 14 15 14 13 14 15 14 15 16 17 18 17 18 17 18 19 20 19 20 19
[251] 20 19 20 21 22 23 22 23 22 21 20 19 18 19 20 21 22 23 22 23 22 21 22 21 22
[276] 21 22 23 22 23 22 23 24 23 22 23 24 25 26 27 26 27 28 29 30 29 28 29 30 31
[301] 32 31 32 33 32 31 30 29 30 29 28 27 28 29 28 29 28 27 28 29 28 29 30 31 32
[326] 31 32 33 32 31 30 29 28 27 26 25 26 25 26 27 28 29 28 29 28 29 28 27 28 27
[351] 26 27 26 25 24 25 24 23 22 23 22 21 20 19 18 17 16 17 16 15 14 13 12 11 12
[376] 13 14 15 16 15 14 15 16 17 16 15 16 17 18 17 16 15 16 15 14 15 16 17 18 17
[401] 18 17 16 15 14 13 14 15 14 13 14 13 14 15 14 15 16 15 16 17 18 19 20 19 18
[426] 19 18 19 18 17 18 19 20 19 20 19 20 21 20 21 22 23 24 23 22 21 22 23 24 25
[451] 26 25 24 23 22 21 20 19 20 19 20 21 20 19 18 19 18 19 20 19 20 19 18 17 18
[476] 17 18 17 18 17 18 17 16 17 16 15 16 17 18 17 18 17 16 17 18 17 16 15 16 15
> x <- b
> postscript(file="/var/www/html/rcomp/tmp/2mlqs1228068352.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] 8.806820e+02 5.161867e+03 1.130264e+03 2.265546e+02 3.382427e+02
  [6] 7.578729e+00 9.062888e+01 2.046491e+02 1.651805e+02 1.367056e+02
 [11] 4.562630e+01 3.109003e+01 2.555574e+00 1.214911e+01 6.885324e+00
 [16] 4.347561e+01 2.165976e+01 2.514840e+00 1.092126e+01 5.142066e+01
 [21] 3.219539e+00 1.323399e+01 1.082108e+01 5.282766e+00 4.142310e+00
 [26] 1.782528e+01 2.681001e+00 3.188187e+01 1.544397e+00 2.861142e+00
 [31] 5.844176e-01 1.007531e+00 5.017071e+00 6.875761e-01 3.234134e+00
 [36] 7.939753e+00 3.694494e+00 1.045840e+01 2.132733e-01 5.428369e+00
 [41] 1.099991e+01 1.030528e+01 2.673773e+00 8.311095e-03 1.232368e+00
 [46] 5.217778e+00 6.175971e+00 2.769536e+00 7.916322e+00 1.197038e+00
 [51] 1.320547e+00 2.331414e-01 8.315779e-01 9.738062e+00 3.263767e+00
 [56] 6.970387e+00 5.727997e-01 2.606867e-01 3.155378e-01 2.489806e+00
 [61] 1.298607e-01 8.202403e-02 1.397387e+00 7.640134e-01 1.818385e+00
 [66] 2.337682e+00 1.512798e+00 3.388960e+00 3.061132e-01 7.233638e-01
 [71] 1.590458e-01 3.080493e-01 5.121338e-01 4.077750e+00 1.597307e+00
 [76] 1.386852e+00 9.563246e-01 9.594877e-01 1.312961e+00 2.197303e-01
 [81] 7.931128e-01 6.997907e-01 4.351374e-01 1.298073e+00 4.100765e-01
 [86] 7.699564e-01 1.425278e-02 1.898260e-01 1.312931e-01 3.084991e+00
 [91] 1.167440e+00 6.985039e-01 1.564231e+00 5.703086e-01 7.196779e-01
 [96] 1.314843e+00 5.779426e-01 1.490730e-01 1.210801e-01 1.783894e-01
[101] 3.461447e+00 5.684927e-01 4.352803e-01 2.824963e-01 2.578743e+00
[106] 1.095152e+00 9.147498e-01 5.313252e-01 1.774085e+00 6.631826e-01
[111] 2.883394e-01 3.660423e-01 2.123097e-01 2.153678e-02 5.441377e-01
[116] 2.750574e-01 3.274793e-01 5.871858e-01 1.593125e-02 6.720616e-03
[121] 6.290360e-02 8.249880e-01 1.851963e-01 7.607826e-02 4.828863e-01
[126] 1.629009e-01 4.573205e-01 5.936089e-01 7.192614e-02 2.617854e-01
[131] 5.097632e-01 6.629841e-01 4.952452e-01 2.048102e-01 1.128712e+00
[136] 8.299411e-01 1.303661e-01 5.011730e-01 7.157827e-02 4.490205e-02
[141] 7.315331e-02 3.205726e-02 1.315748e-02 2.193388e-01 9.799585e-02
[146] 2.366356e-01 3.450666e-01 7.817057e-01 5.890019e-01 5.242927e-01
[151] 7.897009e-01 5.983205e-02 3.418591e-01 1.994646e-02 2.902741e-01
[156] 1.545208e-01 9.275931e-03 3.543532e-01 1.641994e-01 1.571902e-01
[161] 5.191006e-01 1.773498e-02 6.685367e-02 2.049725e+00 4.549669e-01
[166] 5.348243e-02 4.882457e-01 3.308955e-01 8.294134e-01 3.861756e-01
[171] 1.912861e-01 3.155716e-01 3.786774e-01 1.070055e+00 5.409173e-01
[176] 3.861975e-01 6.842535e-01 5.259165e-01 1.003918e-01 3.886998e-02
[181] 2.235132e-01 1.633606e-01 1.856434e+00 2.879741e-01 8.668176e-02
[186] 1.987106e-01 2.467501e-02 4.120318e-01 3.549258e-01 5.459094e-01
[191] 9.551229e-02 3.346961e-02 3.151987e-01 1.059167e-01 1.800211e-01
[196] 2.113885e-01 1.738585e-01 1.916213e-01 1.090408e-01 1.162906e-01
[201] 3.085304e-02 6.667080e-01 5.781712e-01 3.744162e-01 5.390528e-02
[206] 5.774408e-01 9.998521e-02 5.154489e-01 3.706899e-01 4.514334e-01
[211] 3.102056e-01 4.007153e-01 3.201785e-01 5.259641e-02 8.012669e-01
[216] 1.083706e-01 1.874756e-01 1.105412e-02 3.836814e-02 5.043802e-02
[221] 4.906132e-02 1.340998e-01 2.699564e-02 1.526890e-01 4.432685e-01
[226] 1.859032e-01 2.173571e-01 1.715346e-01 1.934871e-03 5.561935e-02
[231] 1.803028e-01 1.958185e-01 3.694445e-01 1.735909e-01 3.580368e-01
[236] 7.567931e-02 1.712670e-01 1.069976e-01 2.666171e-01 1.016592e-01
[241] 2.466916e-01 4.160432e-01 2.239608e-01 3.796379e-01 1.293805e-01
[246] 1.111248e-01 2.979900e-01 2.135346e-01 3.985298e-01 8.298685e-04

$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/rcomp/tmp/3ai2y1228068352.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/1yzwh1228068352.ps tmp/1yzwh1228068352.png")
> system("convert tmp/2mlqs1228068352.ps tmp/2mlqs1228068352.png")
> system("convert tmp/3ai2y1228068352.ps tmp/3ai2y1228068352.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.839   0.611   0.982