FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_rwalk.wasp
Title produced by softwareLaw of Averages
Date of computationSun, 30 Nov 2008 10:14:23 -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/t1228065323ye338eso28isysw.htm/, Retrieved Sun, 19 May 2024 10:41:37 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=26628, Retrieved Sun, 19 May 2024 10:41:37 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact154

Tree of Dependent Computations

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 17:50:19] [b98453cac15ba1066b407e146608df68]
F R       [Law of Averages] [Non-stationary ti...] [2008-11-30 17:14:23] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum
2008-12-04 10:32:05 [Steven Vercammen] [reply] 
Deze vraag werd correct beantwoord. De analogie met de beurskoers wordt ook uitgelegd. Er worden stukjes uit de theorie opgenomen. Er wordt correct geconcludeerd dat er geen sprake is van seizonaliteit of een trend. Dit is enkel schijn, omdat er volledig willekeurig een munt wordt opgegooid.
2008-12-07 09:30:22 [Käthe Vanderheggen] [reply] 
Dit is zeer correct en uitvoerig besproken. Wanneer we de random walk verschillende keren reproduceren, merken we dat de grafieken er telkens anders uitzien. De uitkomst is volledig toevallig. Het is goed dat de student de beurskoers vergelijkt met de random walk, dit klopt wel enkel als we aannemen dat de beurskoers niet voorspelbaar is.
2008-12-08 18:58:31 [Koen Van Baelen] [reply] 
Q1: De student heeft de simulatie van de random-walk slechts 1 maal gereproducet, het was de bedoeling om dit enkele keren te doen en zo te kunnen zien dat er steeds een gelijkaardig patroon wordt verkregen bij proportion of heads. Verder wordt er wel duidelijk uitgelegd dat er geen sprake is van een lange termijntrend en van seizoenaliteit. Ook het voorbeeld van de beurs is heel toepasselijk, hiermee mee kan men duidelijk bepalen dat er geen sprake is van seizoenaliteit. Dit verhaal is analoog aan wat er op de beurs gebeurt. De beurskoers wordt volgens deze redenering bepaald door volgende formule Beurskoers t = beurskoers t-1 + et waarbij et een term is die bepaald wordt door het toeval, dit kan vanalles zijn bv. de winstverwachting, slecht nieuws,… de kans dat dit geval positief is, is 50% en dus even groot als de kans dat dit getal positief is. Misschien had men nog kort kunnen verklaren waarvoor ‘Proportion of Heads’ voor staat, Proportion of Heads geeft de proportie van het aantal keren dat er kop wordt gegooid ten opzichte van het aantal keer er let wordt gegooid. Omdat de kans fifty/fifty is dat er kop of let wordt gegooid, is het logisch dat deze naar de 50% zal neigen.

Post a new message

Dataset

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\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 & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=26628&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=26628&T=0

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

As an alternative you can also use a QR Code:  
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 time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135


Figures (Output of Computation)

Input Parameters & R Code

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()