R version 2.7.0 (2008-04-22) 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. > par2 = '0.5' > par1 = '500' > 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/1udks1227954720.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 2 1 0 -1 0 -1 [19] 0 1 2 3 4 3 4 5 6 5 6 7 6 7 8 7 6 5 [37] 4 3 4 5 4 5 6 5 6 5 4 5 4 5 4 3 4 5 [55] 6 5 6 5 4 3 2 1 0 -1 0 -1 -2 -3 -4 -3 -2 -1 [73] 0 -1 0 1 2 1 2 1 0 -1 -2 -1 -2 -3 -4 -5 -4 -5 [91] -6 -5 -6 -7 -6 -5 -6 -5 -4 -3 -2 -1 -2 -1 0 1 2 1 [109] 2 1 0 1 2 1 0 1 0 1 2 1 2 1 0 -1 -2 -1 [127] -2 -3 -4 -5 -6 -7 -6 -5 -6 -5 -6 -7 -8 -9 -8 -7 -8 -9 [145] -10 -11 -12 -11 -10 -9 -8 -7 -8 -9 -8 -9 -10 -9 -8 -9 -8 -9 [163] -10 -11 -12 -13 -14 -15 -14 -13 -14 -15 -16 -15 -16 -15 -16 -17 -18 -17 [181] -16 -15 -14 -15 -14 -13 -12 -11 -12 -13 -14 -13 -14 -13 -14 -15 -14 -13 [199] -14 -13 -12 -11 -12 -11 -10 -9 -8 -9 -10 -11 -12 -11 -10 -9 -10 -9 [217] -8 -9 -10 -11 -10 -9 -8 -9 -10 -9 -10 -11 -12 -13 -14 -15 -16 -15 [235] -16 -17 -18 -17 -16 -17 -18 -19 -20 -19 -20 -21 -22 -23 -22 -23 -22 -21 [253] -20 -19 -18 -19 -20 -19 -18 -19 -18 -19 -18 -17 -16 -15 -16 -17 -16 -15 [271] -16 -15 -16 -17 -16 -17 -18 -19 -18 -19 -20 -21 -22 -21 -20 -21 -22 -23 [289] -22 -21 -20 -21 -20 -21 -20 -21 -22 -21 -20 -21 -20 -19 -20 -21 -20 -19 [307] -18 -19 -20 -21 -22 -23 -22 -23 -22 -23 -22 -23 -24 -23 -22 -21 -22 -23 [325] -22 -21 -22 -21 -22 -23 -22 -21 -22 -23 -24 -23 -24 -25 -26 -25 -26 -27 [343] -26 -25 -26 -27 -26 -27 -28 -29 -30 -29 -30 -29 -30 -31 -30 -31 -32 -33 [361] -32 -31 -32 -33 -34 -35 -36 -35 -34 -33 -34 -35 -36 -37 -38 -39 -40 -41 [379] -42 -41 -40 -41 -42 -43 -44 -43 -42 -43 -44 -45 -44 -43 -42 -41 -40 -41 [397] -40 -39 -38 -39 -40 -41 -42 -43 -44 -43 -44 -45 -44 -45 -46 -45 -44 -43 [415] -44 -45 -44 -45 -44 -43 -44 -43 -44 -43 -44 -45 -44 -43 -44 -43 -42 -43 [433] -44 -45 -46 -45 -44 -43 -44 -45 -46 -45 -46 -47 -48 -47 -48 -49 -48 -49 [451] -50 -51 -50 -51 -52 -51 -50 -49 -50 -51 -50 -49 -48 -49 -48 -49 -50 -51 [469] -50 -51 -50 -49 -50 -51 -52 -53 -52 -53 -52 -51 -52 -53 -54 -55 -56 -57 [487] -56 -55 -54 -55 -54 -55 -56 -55 -54 -53 -54 -55 -56 -55 > x <- b > postscript(file="/var/www/html/rcomp/tmp/2hdbm1227954720.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.392317e+03 7.686489e+02 1.750499e+02 1.792932e+02 4.843418e+02 [6] 7.399515e+01 2.447360e+02 3.020029e+01 1.439919e+02 9.455193e+01 [11] 3.645032e+01 1.614070e+01 5.855047e+01 1.219444e+02 2.212973e+01 [16] 3.370961e+01 6.359574e+00 2.289163e+01 4.143945e+01 1.367675e+01 [21] 1.597100e+01 2.623642e+01 1.198096e+01 9.431778e+00 5.200854e+00 [26] 7.497603e-01 1.454639e-01 8.169062e+00 1.128869e+01 8.659224e+00 [31] 2.042881e+01 2.559113e+00 1.096382e+00 3.254636e+00 1.687861e+00 [36] 3.817054e-01 1.019589e+00 1.512561e+01 3.883740e+00 2.968743e-01 [41] 1.041147e+00 1.938830e-01 3.942667e+00 2.949887e-01 9.490891e-01 [46] 2.160417e+00 4.038839e+00 1.418684e+00 1.378211e+00 1.185352e+00 [51] 1.412926e+00 2.270951e+00 9.235566e-02 3.351595e+00 2.983542e+00 [56] 6.695204e-01 2.409977e+00 3.558024e+00 3.751098e+00 1.014953e+00 [61] 2.812633e+00 3.284023e-01 1.155109e+00 6.735239e+00 7.471984e-01 [66] 9.804597e-01 4.408086e-01 1.694344e+00 2.131857e-01 6.062864e+00 [71] 3.836123e-01 8.682076e-01 1.941731e+00 2.119377e-01 1.040761e+00 [76] 7.800905e-01 4.540161e-01 4.138832e+00 2.502331e-01 2.271147e+00 [81] 9.471393e-01 3.749577e+00 1.435278e+00 9.820414e-02 1.340191e-01 [86] 1.975449e+00 1.382632e+00 1.593472e+00 7.283642e-01 5.586952e-01 [91] 5.513586e-01 1.267089e+00 6.735334e-02 3.755912e-01 2.643175e+00 [96] 3.158170e+00 1.738596e+00 2.820523e+00 4.879011e-01 2.540640e-01 [101] 2.507399e-01 2.578261e-01 2.166085e-01 9.692935e-01 3.272168e-01 [106] 5.467914e-01 1.156134e-01 9.890495e-01 5.155827e-01 1.004368e+00 [111] 7.010556e-02 6.078467e-01 3.023649e-02 1.999581e-01 6.707130e-01 [116] 5.991569e-02 5.829330e-02 2.103289e-01 2.340314e+00 5.476997e-02 [121] 2.090138e-01 1.021545e-02 7.702570e-01 1.713109e-01 1.297969e+00 [126] 3.627864e-02 1.382885e-01 6.957620e-02 8.510563e-01 6.776953e-03 [131] 5.999432e-01 3.023435e-01 7.670981e-01 2.853792e-01 1.940801e+00 [136] 1.570606e+00 8.991437e-01 9.436923e-01 3.606357e-01 7.180164e-01 [141] 3.241905e-01 4.939402e-01 5.522910e-01 1.028066e-01 2.733788e-01 [146] 6.727466e-02 4.617000e-01 1.457865e-01 3.565332e-03 5.390312e-01 [151] 7.790853e-01 1.496912e-01 1.007369e+00 4.137475e-02 7.952005e-02 [156] 8.402194e-02 1.728324e-01 2.168789e-01 1.247662e-01 1.906877e-01 [161] 1.241383e+00 1.050694e-01 3.039326e-01 2.829756e-01 3.214268e-01 [166] 8.454825e-03 4.697741e-01 1.724144e-02 4.394250e-01 3.144734e-01 [171] 6.978123e-01 7.593142e-01 8.159567e-02 1.578590e-02 1.110413e-01 [176] 7.355375e-01 2.247604e-01 3.688723e-01 1.762040e-01 1.374451e-01 [181] 4.585986e-03 8.840910e-01 1.860438e-01 2.040190e-02 2.401596e-01 [186] 6.382302e-01 2.668188e-02 9.287658e-02 1.125466e-01 5.272879e-02 [191] 8.409153e-01 4.315292e-01 5.192616e-02 3.320113e-01 2.012029e-02 [196] 4.096370e-01 1.687205e-02 2.960673e-01 2.472842e-01 3.007345e-01 [201] 4.927697e-02 3.568891e-01 7.953899e-01 7.870965e-01 2.494262e-01 [206] 1.295420e-01 2.484685e-01 1.422730e-01 2.314727e-02 2.844212e-02 [211] 1.955360e-01 3.927370e-01 4.264705e-01 2.335290e-01 1.704665e-01 [216] 5.993665e-02 3.573184e-02 4.893899e-01 1.808610e-01 1.526197e-01 [221] 3.372263e-01 6.063320e-01 1.113403e-01 2.689637e-01 2.722512e-02 [226] 1.389998e-01 1.091561e-01 6.052021e-02 2.875762e-02 1.238660e-01 [231] 1.608812e-01 4.438960e-03 7.198528e-02 1.031645e-01 2.822263e-01 [236] 1.055284e-01 2.662516e-02 1.584804e-03 2.428610e-01 4.086397e-01 [241] 1.042777e-01 2.536180e-02 9.229842e-02 3.860321e-01 2.847668e-01 [246] 1.671272e-01 1.877249e-01 5.832829e-01 4.463029e-02 1.335373e-02 $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/3fq6o1227954720.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/1udks1227954720.ps tmp/1udks1227954720.png") > system("convert tmp/2hdbm1227954720.ps tmp/2hdbm1227954720.png") > system("convert tmp/3fq6o1227954720.ps tmp/3fq6o1227954720.png") > > > proc.time() user system elapsed 1.908 1.010 2.029