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. > 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/1zmlj1228171827.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 1 2 1 0 -1 0 1 2 1 0 -1 0 1 [19] 0 -1 -2 -1 0 1 0 1 2 3 2 1 0 1 2 3 2 3 [37] 2 3 2 1 0 -1 -2 -3 -2 -1 0 -1 0 -1 -2 -3 -4 -3 [55] -2 -3 -2 -1 -2 -1 -2 -1 0 -1 0 -1 -2 -3 -2 -3 -2 -1 [73] -2 -1 -2 -3 -2 -1 0 1 2 3 2 3 4 5 6 7 6 7 [91] 6 5 6 7 6 7 6 7 6 5 6 7 6 5 4 3 4 3 [109] 2 1 0 1 2 1 2 1 0 -1 -2 -3 -4 -3 -4 -3 -4 -5 [127] -6 -5 -6 -7 -6 -7 -8 -9 -10 -9 -10 -11 -10 -11 -10 -9 -10 -9 [145] -8 -9 -8 -7 -8 -7 -8 -7 -6 -5 -4 -3 -4 -3 -2 -1 0 -1 [163] 0 1 0 -1 0 1 0 1 2 3 2 1 0 -1 0 1 2 1 [181] 0 1 2 1 0 1 0 -1 -2 -3 -2 -3 -4 -3 -2 -3 -2 -1 [199] 0 1 2 1 2 1 2 3 2 1 0 -1 -2 -1 0 -1 -2 -3 [217] -2 -1 -2 -1 -2 -1 0 1 2 3 2 3 4 5 6 7 8 9 [235] 8 9 8 7 6 7 8 7 6 7 8 9 10 9 10 11 10 11 [253] 10 11 12 11 10 11 12 13 12 11 10 9 10 11 10 11 12 11 [271] 10 11 10 11 12 11 10 11 10 9 8 9 8 7 8 7 8 7 [289] 6 7 6 5 4 5 4 3 2 3 2 1 0 -1 -2 -1 0 -1 [307] -2 -1 -2 -3 -4 -5 -4 -5 -6 -5 -4 -5 -4 -5 -4 -3 -2 -3 [325] -4 -3 -4 -5 -6 -7 -8 -7 -6 -7 -6 -7 -8 -7 -8 -7 -6 -7 [343] -6 -5 -4 -3 -4 -3 -2 -3 -2 -1 -2 -3 -2 -3 -4 -3 -2 -3 [361] -2 -1 -2 -1 -2 -1 -2 -3 -2 -1 -2 -1 -2 -1 -2 -3 -2 -1 [379] -2 -1 -2 -1 0 1 2 3 2 3 2 1 2 3 2 1 2 3 [397] 2 1 0 1 0 -1 0 1 0 -1 0 -1 -2 -1 -2 -1 -2 -1 [415] 0 1 0 1 0 -1 0 -1 0 1 2 1 2 1 2 3 4 3 [433] 2 3 4 5 6 5 6 7 8 7 6 7 8 9 10 9 10 9 [451] 10 11 10 11 12 11 12 11 10 11 10 9 10 9 8 7 6 5 [469] 6 7 8 7 6 7 6 5 6 5 4 5 4 5 4 3 2 1 [487] 0 1 0 1 2 3 4 3 2 3 4 5 4 5 > x <- b > postscript(file="/var/www/html/rcomp/tmp/2kzim1228171827.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] 2.761642e+01 2.115292e+03 2.000308e+03 5.552342e+01 6.171810e+02 [6] 6.860622e+02 5.969108e+02 5.463956e+01 6.976862e+01 3.754420e+01 [11] 9.401190e+00 1.724380e+01 3.454240e+01 3.078627e+01 2.662771e+01 [16] 9.858917e+00 2.997120e+01 2.866516e+01 3.316593e-01 1.053664e+01 [21] 3.263562e+01 1.564865e+01 3.838527e+00 7.418115e+00 2.966200e+00 [26] 6.271725e-02 9.119089e+00 2.553303e+00 5.246393e+00 2.283219e+00 [31] 8.044064e-01 1.056688e-01 7.351573e+00 4.354106e+00 3.113655e+00 [36] 9.885385e-01 5.141412e+00 3.468335e+00 1.574184e+00 5.031523e+00 [41] 9.778024e-01 5.246428e-01 4.854153e-01 4.322933e-02 5.500682e+00 [46] 4.778280e+00 1.445651e+00 7.354677e-01 9.340740e-01 2.144047e+00 [51] 1.040292e-01 2.859984e+00 6.766968e-01 2.222771e+00 2.702351e+00 [56] 2.423373e+00 3.257565e-01 4.410813e-01 3.369768e+00 2.264248e+00 [61] 1.835662e+00 1.007233e+00 2.290353e-01 1.538178e+00 2.640939e+00 [66] 1.204960e+00 2.410106e+00 8.198817e-01 3.311731e-01 4.019885e-02 [71] 8.371792e-01 1.920494e+00 1.833292e+00 1.828618e-01 1.007725e+00 [76] 2.345523e+00 2.364457e-01 7.105816e-01 1.475869e-02 5.296359e-01 [81] 7.689057e-01 4.181903e-02 2.281736e-01 9.624433e-01 1.227673e+00 [86] 1.586517e+00 3.256197e+00 2.950083e-01 1.801548e-01 6.879361e-01 [91] 3.047184e-01 5.813914e-01 3.165122e-01 1.291013e-01 5.573561e-01 [96] 1.235139e-01 1.032775e+00 1.285886e-01 1.283938e+00 7.062889e-02 [101] 2.534758e-01 1.465061e+00 1.039515e+00 1.712631e-01 1.770364e+00 [106] 3.023747e-02 5.153351e-01 9.332463e-02 4.317371e-01 2.649131e+00 [111] 4.818010e-01 2.945968e-01 6.205217e-02 1.780192e+00 1.426791e-01 [116] 3.858472e-01 2.130645e-01 1.391741e-01 1.996071e-01 8.701842e-01 [121] 1.055609e-01 6.681005e-02 1.134996e+00 1.975934e-02 8.140318e-02 [126] 1.921153e-01 1.458151e+00 2.624318e-01 1.518140e+00 7.543460e-01 [131] 1.330713e+00 8.744531e-01 2.972414e-01 9.870550e-01 1.460863e-01 [136] 1.291847e+00 2.434558e-01 8.128509e-01 7.366159e-01 2.086208e-01 [141] 8.145372e-02 2.433064e-01 5.544960e-02 3.832047e-01 4.541334e-02 [146] 1.790844e-02 2.368931e-01 1.808003e-01 3.368767e-01 6.026089e-01 [151] 1.560284e-01 9.091840e-01 1.163741e-01 4.677704e-01 4.189389e-01 [156] 2.303634e-01 1.488869e+00 1.000742e-01 3.779279e-01 8.584180e-02 [161] 1.920540e-02 9.008371e-02 4.753008e-01 5.443069e-02 1.034558e-01 [166] 3.836402e-01 1.226236e+00 5.359109e-01 1.200972e-01 1.807807e-01 [171] 7.023027e-01 1.127732e+00 3.297103e-01 3.729826e-01 4.869657e-01 [176] 7.093582e-01 6.880011e-01 7.182733e-01 3.325704e-01 1.474229e+00 [181] 1.173356e+00 4.981019e-01 5.692586e-01 7.342672e-01 1.172900e+00 [186] 4.264669e-01 9.606950e-01 7.708661e-01 4.729529e-01 8.182980e-02 [191] 1.940827e-01 4.089284e-01 3.465486e-02 2.553500e-01 8.135021e-01 [196] 6.382652e-02 1.838759e-01 1.047180e-01 1.481773e-01 1.180351e-01 [201] 1.116913e-01 7.771952e-03 4.130225e-01 3.016310e-02 2.552211e-02 [206] 1.101680e-03 1.984631e-01 3.934719e-01 4.621070e-01 2.892628e-01 [211] 1.412139e-01 3.916779e-02 6.310927e-02 2.350471e-01 1.921786e-01 [216] 2.337777e-01 1.237136e-01 1.012419e-01 2.545543e-01 3.507830e-02 [221] 3.619462e-01 5.467670e-01 2.168948e-02 7.592450e-02 1.663919e-01 [226] 1.480352e-01 2.789365e-01 3.410554e-01 3.371396e-01 1.291557e-02 [231] 8.646787e-02 1.080660e-03 4.112233e-01 2.692005e-01 3.174356e-02 [236] 4.616174e-01 2.853601e-01 1.135874e-01 1.037633e-01 1.716005e-01 [241] 3.747460e-02 2.947333e-01 2.789400e-01 7.299584e-02 1.490374e-01 [246] 9.041317e-01 1.759253e-01 1.242132e-01 1.239223e-01 7.074310e-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/3ltfc1228171827.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/1zmlj1228171827.ps tmp/1zmlj1228171827.png") > system("convert tmp/2kzim1228171827.ps tmp/2kzim1228171827.png") > system("convert tmp/3ltfc1228171827.ps tmp/3ltfc1228171827.png") > > > proc.time() user system elapsed 0.818 0.606 1.064