Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_rwalk.wasp
Title produced by softwareLaw of Averages
Date of computationSun, 30 Nov 2008 03:50:28 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Nov/30/t1228042388v3dxg13slsjcon6.htm/, Retrieved Sun, 19 May 2024 09:19:51 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=26426, Retrieved Sun, 19 May 2024 09:19:51 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact198
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Law of Averages] [Random Walk Simul...] [2008-11-25 18:05:16] [b98453cac15ba1066b407e146608df68]
F         [Law of Averages] [Non Stationary Ti...] [2008-11-30 10:50:28] [dafd615cb3e0decc017580d68ecea30a] [Current]
Feedback Forum
2008-12-07 11:52:06 [Dana Molenberghs] [reply
Het gaat hier om de autocorrelatie van de kopkant. Lag 1 = 1 X kop gegooid, lag 2 = 2de maal kop gegooid. Er werd hier berekend of er een verband tussen was.

We zien een duidelijk dalend trend in de autocorrelatie. Als ik hier een logische verklaring aan moet proberen te geven is deze het volgende: onder het begrip autocorrelatie verstaan we de correlatie tussen twee opeenvolgende waarden in dezelfde serie getallen. Hoeveel te meer we tossen wordt het duidelijk dat het verband tussen het tossen van de kop daalt. Dit komt omdat het louter toeval is wanneer kop wordt gegooid.
2008-12-08 20:51:54 [Jeroen Michel] [reply
Je analyse klopt niet helemaal. De stelling die je maakt weinig te maken met het verloop van de grafiek.

Op de onderstaande grafiek is duidelijk af te lezen dat de sterke correlatie die aanwezig is sterk afneemt (dalend verloop).

Dit aangezien autocorrelatie te maken heeft met het feit dat er een verband bestaat tussen voorgaande data en de huidge data van de tijdreeksen die worden genomen.

Post a new message




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 2 seconds \tabularnewline
R Server & 'Herman Ole Andreas Wold' @ 193.190.124.10:1001 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=26426&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Herman Ole Andreas Wold' @ 193.190.124.10:1001[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=26426&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=26426&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Herman Ole Andreas Wold' @ 193.190.124.10:1001



Parameters (Session):
par1 = 500 ; par2 = 0.5 ;
Parameters (R input):
par1 = 500 ; par2 = 0.5 ;
R code (references can be found in the software module):
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
bitmap(file='test1.png')
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()
b
bitmap(file='pic1.png')
racf <- acf(b,n/10,main='Autocorrelation',xlab='lags',ylab='ACF')
dev.off()
racf