FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

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 computationTue, 02 Dec 2008 08:07:51 -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/Dec/02/t1228230593x4ighizo7u22qfb.htm/, Retrieved Sat, 25 May 2024 13:51:52 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=27923, Retrieved Sat, 25 May 2024 13:51:52 +0000
QR Codes:

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

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 18:40:39] [b98453cac15ba1066b407e146608df68]
F         [Law of Averages] [Q4] [2008-12-02 15:07:51] [3fc0b50a130253095e963177b0139835] [Current]
Feedback Forum
2008-12-04 17:18:23 [Loïque Verhasselt] [reply] 
Q4: De juiste output is getoond maar er ontbreekt een verklaring.We voeren hier dus een spectrale analyse door.We gaan dus de oorspronkelijke tijdreeks reconstrueren op basis van golfbewegingen met een bepaald aantal frequenties en perioden(geholpen met de sinus en cosinus functie).Grafiek1: wanneer we een lage frequentie zien wil dit zeggen een lange periode tussen 2 toppen van de cosinus of sinus functie. De sinus of cosinus gaat dus langzaam op en neer. Nu passen we dit toe op onze random-walk.Welke combinatie van sinus en cosinus gaat voorkomen in het periodogram.(Lange termijn trend = lange periode en lage frequentie).Je ziet dus in de raw periodogram een langzaam dalend patroon.Deze trend moeten we eruit halen door te differentiëren. Nu de tweede grafiek of de cumulatieve periodogram. Als we bijvoorbeeld een hele kleine frequentie nemen rond de 0,1 wordt 80% van de tijdreeks verklaard. Wat wijst op een Lange Termijn (LT) beweging.
2008-12-07 17:31:10 [Samira Zeroual] [reply] 
Op de eerste grafiek kunnen we ook een langzaam dalende trend waarnemen. Op de x-as op wordt de frequentie weergegeven, terwijl de y-as het spectrum voorstelt Ook zien we duidelijk dat de varianties steeds verschillen bij de verscheidene frequenties.

Op het cumulatief periodogram, zien we een stijl stijgend verloop, wat wijst op een langetermijntrend. We kunnen ook zien dat de cumulatieve frequentiecurve buiten het betrouwbaarheidsinterval ligt. Ontbreken van het getrapte patroon wijst op het ontbreken van seizoenale invloeden.
2008-12-08 13:05:12 [Anouk Greeve] [reply] 
De snelle stijging betekent dat er een trend aanwezig is, dus D mag nul blijven omdat er geen trap aanwezig is. Als er geen trap aanwezig is dan is er ook geen seizoenaliteit.

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 time2 seconds
R Server'George Udny Yule' @ 72.249.76.132

\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 & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=27923&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]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=27923&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=27923&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 time2 seconds
R Server'George Udny Yule' @ 72.249.76.132


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