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 09:41:02 -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/t1228063295n0nla093jjascdg.htm/, Retrieved Sun, 19 May 2024 00:03:05 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=26615, Retrieved Sun, 19 May 2024 00:03:05 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordslaw of averages
Estimated Impact141
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:40:39] [b98453cac15ba1066b407e146608df68]
F         [Law of Averages] [law of averages q4] [2008-11-30 16:41:02] [f4b2017b314c03698059f43b95818e67] [Current]
Feedback Forum
2008-12-04 17:29:56 [339a57d8a4d5d113e4804fc423e4a59e] [reply
De student gebruikt de juiste software om deze vraag op te lossen. Bij het cumulatieve periodogram, worden op de y-as alle frequenties opgeteld. In dit geval zien we een steil stijgende lijn, die op waarde 1 horizontaal wordt. De lijn valt hier niet in het betrouwbaarheidsinterval.
2008-12-05 15:28:17 [Kristof Van Esbroeck] [reply
Correcte berekening en beknopte interpretatie.

We nemen eerst het Raw Perdiogram en vervolgens het Cumulative Perdiogram waar.

In het Raw Perdiogram noteren we op de y as de trend en op de x as nemen we bij interpretatie de frequentie waar. We concluderen op lange termijn een groot aantal lage frequenties.

Het Cumulative Perdiogram geeft een bijna verticale rechte weer die omstreeks 0,9 evolueert. Op de y as noteren we waarden tussen 0 en 1 bij een relatieve frequentie die we op de x as kunnen waarnemen.

Post a new message




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 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 & 4 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=26615&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]4 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=26615&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=26615&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 time4 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
x <- b
bitmap(file='test1.png')
r <- spectrum(x,main='Raw Periodogram')
dev.off()
r
bitmap(file='test2.png')
cpgram(x,main='Cumulative Periodogram')
dev.off()