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 computationTue, 02 Dec 2008 07:46:41 -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/t1228229235r49k1anx8fdspcd.htm/, Retrieved Sun, 19 May 2024 12:20:05 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=27882, Retrieved Sun, 19 May 2024 12:20:05 +0000
QR Codes:

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

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] [W7Q4] [2008-12-02 14:46:41] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum
2008-12-08 15:40:03 [Kevin Vermeiren] [reply] 
De student geeft een correct maar ook zeer beknopt antwoord. Er werd niets gezegd over de werking van de spectraal analyse. Nuttig had geweest te vermelden dat deze methode gebruikt wordt om de trend te onderzoeken. Deze methode gaat de tijdreeks ontbinden in regelmatige golfbewegingen. Met andere woorden, de sinus en cosinusfuncties worden uit de tijdreeks gehaald. Als er sprake is van een lange periode, is er sprake van een lage frequentie. Uiteraard want de formule voor de frequentie = 1/periode. Uit het raw periodogram blijkt inderdaad dat de golfbewegingen met een lage frequentie (golfbeweging over een lange periode) vaker voorkomen dan een hogere frequentie. Er had ook nog vermeld kunnen worden dat er duidelijk sprake is van een langzaam dalende trend. Vervolgens bekijken we het cumulative periodogram. De student geeft ook hier een goede interpretatie maar er ontbreekt nuttige informatie. Het klopt dat een lage frequentie het vaakst voorkomt. Een frequentie van 0 tot 0.05 hertz vertegenwoordigd reeds 95% van alle voorkomende frequenties. Hier had nog vermeld kunnen worden dat een steil begin duidt op de aanwezigheid van een lange termijn trend. Het klopt ook dat er geen sprake is van toeval daar de gegevens niet binnen de blauwe stippellijnen (betrouwbaarheidsinterval) liggen. Er had ook nog vermeld kunnen worden hoe de grafiek geinterpreteerd moet worden namelijk: wil ik 80% (0.8 op Y-as) van de reeks verklaren is duidelijk een lage frequentie nodig (X-as).
2008-12-08 21:40:54 [] [reply] 
Je had nog wat meer kunnen vertellen. Dat er gebruik wordt gemaakt van een spectraal analyse, om de trend te onderzoeken.
De lage frequentie wijst juist op een lange periode. Op de cumulative periodogram zien we dat we 80% van de tijdsreeks kunnen verklaren. Het is dan ook een snel stijgende lijn dat we zien. Er zijn geen trapjes die zouden kunnen wijzen op seizoenaliteit. Er is dus een lange termijnbeweging en geen seizoenaliteit.
2008-12-09 11:23:40 [Yannick Van Schil] [reply] 
correct geantwoord, je had inderdaad nog kunnen vermelden dat er gebruik werd gemaakt van de spectral analysis, hierbij verwijst de lage frequentie op een lange termijn trend.

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=27882&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'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()