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/1z3ze1228158652.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  3  2  1  0  1  0  1  2  3  2  3  2  1  0 -1 -2
 [26] -3 -4 -3 -2 -1 -2 -1 -2 -1  0 -1 -2 -3 -2 -3 -2 -3 -4 -3 -4 -3 -2 -3 -4 -3
 [51] -2 -3 -4 -3 -4 -3 -2 -1  0 -1  0 -1  0  1  2  3  4  3  2  1  2  3  2  3  2
 [76]  3  4  3  4  5  6  7  6  7  8  9 10 11 10 11 10  9 10  9  8  7  8  9 10 11
[101] 10 11 10  9 10 11 12 13 12 11 10 11 12 11 10  9 10  9 10 11 12 13 14 15 16
[126] 17 16 15 14 13 12 13 12 13 14 15 14 15 14 15 16 17 16 17 18 19 18 17 18 17
[151] 18 19 20 21 22 21 22 23 22 23 24 23 24 23 24 23 22 23 24 25 26 25 26 25 26
[176] 27 28 29 30 29 28 29 30 29 28 29 28 29 28 29 28 27 28 29 28 29 28 29 30 31
[201] 32 33 34 35 36 37 36 37 38 39 38 39 38 37 38 37 36 37 38 37 38 37 36 37 38
[226] 37 36 37 36 37 36 35 36 35 34 33 32 33 34 33 34 35 36 35 34 35 34 33 34 35
[251] 36 35 36 35 34 33 32 31 30 31 30 31 32 31 32 33 34 35 36 35 34 35 34 33 32
[276] 31 32 33 34 35 34 35 36 35 36 37 38 37 36 37 38 37 38 39 40 41 42 41 42 43
[301] 44 43 42 43 44 45 46 47 46 47 46 47 46 45 44 45 44 45 44 45 46 45 46 47 46
[326] 47 46 45 44 45 44 45 44 43 44 45 44 45 44 43 44 45 46 47 48 49 48 47 48 49
[351] 50 49 50 51 50 49 48 47 46 47 46 47 46 45 46 45 44 45 44 45 46 47 48 49 50
[376] 51 50 49 50 51 52 51 50 49 50 51 52 51 50 49 50 49 48 49 50 49 50 51 52 53
[401] 52 51 50 49 50 51 52 53 54 55 56 55 56 55 56 57 58 59 60 61 60 59 58 59 58
[426] 59 60 59 60 61 62 63 64 65 64 65 66 67 68 67 68 67 68 67 68 69 68 67 68 69
[451] 70 69 70 69 68 67 66 67 66 65 64 63 62 63 64 65 66 65 64 65 66 67 68 67 66
[476] 67 66 65 66 67 68 67 66 65 66 67 66 65 64 63 62 61 60 59 60 61 62 63 62 61
> x <- b
> postscript(file="/var/www/html/rcomp/tmp/2tuig1228158652.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] 5.153956e+02 1.407998e+03 1.499061e+02 6.913451e+02 3.363180e+02
  [6] 3.164936e+01 3.789562e+02 3.802871e+00 6.412360e+00 3.626276e+00
 [11] 9.581398e+01 1.074515e+01 3.514516e+00 3.050380e+01 6.169818e+01
 [16] 1.350934e+01 5.983901e+01 6.896821e+00 1.934953e+00 1.266417e+01
 [21] 2.638504e+01 6.706013e+00 1.949209e+00 4.156491e+01 9.045524e+00
 [26] 1.080318e+01 2.649433e+01 8.327130e+00 1.067697e+01 4.111437e+00
 [31] 2.711867e+00 1.088690e+01 9.312307e-01 3.041730e+00 7.269159e+00
 [36] 3.745112e+00 1.104942e+00 1.266074e+01 2.453285e+00 1.341468e+00
 [41] 1.053153e-01 5.211240e+00 6.508732e-01 5.002232e+00 1.015179e+01
 [46] 4.893330e+00 2.412456e+00 2.424333e+00 5.135927e-01 7.428225e-01
 [51] 2.223050e+00 1.552657e+00 5.601015e+00 4.324370e+00 4.940652e+00
 [56] 3.362488e-01 2.694673e+00 5.160137e-01 2.292324e+00 5.782358e-01
 [61] 2.165535e+00 9.712420e-01 1.827577e+00 9.258200e-01 1.032596e-01
 [66] 6.969315e-01 7.373354e-02 3.517487e-01 6.814233e-01 2.095175e+00
 [71] 2.864439e+00 1.565233e+00 4.486468e-01 9.521569e-01 5.650147e+00
 [76] 3.062588e-01 1.076584e+00 3.406600e+00 1.003869e+00 1.127886e+00
 [81] 2.029247e+00 9.169380e-02 4.594661e-02 1.958218e+00 3.973620e-01
 [86] 1.169574e-01 2.572038e-01 2.030807e+00 1.239746e-02 1.867116e-02
 [91] 4.053437e-01 1.508360e-01 1.580544e-01 2.458193e+00 3.715687e-01
 [96] 1.075660e-01 4.085933e-01 3.824674e-01 1.228077e+00 1.498522e+00
[101] 3.159652e-01 1.723607e+00 1.129841e+00 4.050177e-02 9.908314e-02
[106] 7.896748e-01 3.545452e-01 6.095955e-01 1.860729e-01 4.324086e-01
[111] 1.027218e+00 7.400558e-01 1.195517e-01 3.217298e-01 9.229809e-01
[116] 4.246691e-01 1.371980e-01 5.568991e-01 4.308161e-01 9.588010e-01
[121] 1.580648e-01 1.543257e-01 5.464973e-01 3.010213e-01 8.706931e-01
[126] 3.413082e-02 3.422604e-02 4.283068e-01 3.033553e-01 1.333369e-01
[131] 1.598967e-01 2.242574e-01 1.897519e-02 5.638940e-01 2.924002e+00
[136] 1.727995e-01 2.085851e-01 2.009943e-01 2.062928e-01 2.678100e-01
[141] 8.220434e-02 1.766035e-01 1.594879e-01 5.323861e-01 4.181656e-02
[146] 2.560617e-01 1.362175e-01 5.014682e-01 7.304662e-01 1.344333e-01
[151] 8.065976e-02 3.883664e-02 4.635771e-01 4.982882e-01 1.119210e-01
[156] 1.030150e+00 3.802806e-01 7.780507e-01 1.890316e-01 2.236081e-02
[161] 1.734474e-01 2.374072e-02 7.652269e-01 1.639618e-01 3.328007e-01
[166] 1.827069e-01 8.933881e-01 3.809425e-02 1.733239e-01 6.796819e-01
[171] 1.239917e+00 4.502470e-01 1.577756e-01 3.274595e-02 5.838902e-02
[176] 1.451219e-03 8.868074e-02 1.102263e-01 2.941919e-01 1.706216e+00
[181] 6.367917e-01 2.786189e-01 3.972978e-01 5.217190e-01 5.181222e-02
[186] 4.828536e-01 2.894918e-01 1.182833e-01 5.225377e-01 7.976016e-02
[191] 1.189866e-01 5.391222e-01 4.110421e-01 1.690771e-02 2.767678e-01
[196] 1.093702e-01 4.142980e-01 5.251772e-01 5.755351e-01 3.543155e-01
[201] 3.253969e-01 1.581483e-01 1.302033e-02 6.220887e-01 3.981153e-01
[206] 3.666182e-01 2.159559e-01 2.850616e-01 1.717601e-02 5.278254e-01
[211] 5.244032e-02 4.738466e-02 3.620272e-01 2.817368e-01 4.209248e-01
[216] 2.380020e-01 5.304063e-02 1.598065e-01 1.518940e-01 1.431119e-01
[221] 3.591211e-01 2.471787e-01 1.267302e-01 1.435754e-01 8.056768e-02
[226] 9.676985e-02 1.941006e-01 2.475895e-02 1.081846e-01 3.178300e-01
[231] 2.855613e-01 4.903411e-02 7.346925e-01 1.703934e-01 1.519239e-01
[236] 1.745314e-02 5.336098e-02 9.220864e-02 6.732157e-01 1.462694e-01
[241] 1.003964e+00 3.778559e-01 1.695496e-01 2.443133e-02 3.326761e-01
[246] 3.295149e-01 3.085834e-01 2.104337e-01 9.462220e-02 1.195792e-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/rcomp/tmp/3lrvv1228158652.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/1z3ze1228158652.ps tmp/1z3ze1228158652.png")
> system("convert tmp/2tuig1228158652.ps tmp/2tuig1228158652.png")
> system("convert tmp/3lrvv1228158652.ps tmp/3lrvv1228158652.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.809   0.591   0.992