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

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

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] [] [2008-12-02 14:20:28] [b5935c41f1031f8c061510fc5ad27e97] [Current]
Feedback Forum
2008-12-07 17:05:58 [Samira Zeroual] [reply] 
Op de eerste grafiek kunnen we een langzaam dalende trend waarnemen. Op de x-as op wordt de frequentie weergegeven, terwijl de y-as het spectrum voorstelt. Op deze grafiek zijn duidelijk lage frequenties en grote periodes vast te stellen. Ook zien we duidelijk dat de varianties steeds verschillen bij de verscheidene frequenties.

Op het cumulatief periodogram kunnen we ook zien dat de cumulatieve frequentiecurve buiten het betrouwbaarheidsinterval ligt.
2008-12-08 18:05:12 [Tinneke De Bock] [reply] 
De uitleg bij de Raw Periodogram ontbreekt. De dalende trend die we hierbij kunnen waarnemen duidt op een lange termijn. Een lage frequentie komt overeen met een lange periode, dit verschijnsel kunnen we vergelijken met een sinus die langzaam op en neer gaat. De lengte van 1 golf (= de periode) is hierbij lang en het aantal golven (=de frequentie) klein. Golfbewegingen met een lage frequentie moeten vaak voorkomen aangezien het hier gaat over een lange termijn.

De definitie die de student geeft aan het cumulatieve periodogram is juist, maar hieruit kunnen we nog een en andere conclusies trekken. We kunnen deze grafiek ook interpreteren als r kwadraat.
Indien we een groot aantal waarnemingen (bijvoorbeeld 80%) van de tijdreeks willen verklaren, dan is de overeenstemmende frequentie laag. Aangezien een lage frequentie overeenkomt met een lange periode, wijst een stijl stijgend verloop zoals dit inderdaad op een lange termijn trend.
Wat je ook nog kon vermelden is dat de 2 stippellijnen het 95%-betrouwbaarheidsinterval weergeven. Het verloop van onze grafiek bevindt zich heel duidelijk buiten dit betrouwbaarheidsinterval, waardoor het geen toeval kan zijn.
2008-12-08 19:04:36 [Marlies Polfliet] [reply] 
Het dalende patroon op de raw periodogram en de zeer stijgende lijn (links) op de cumulative periodogram wijzen op een lange termijn trend. De student heeft voor het cumulative periodogram een goede definitie gegeven, maar we kunnen de cumulative periodogram ook beschouwen als R kwadraat, hieruit kunnen we concluderen dat bijvoorbeeld 80% (y-as) kan verklaard worden door een zeer lage frequentie (een lage frequentie betekent een lange periode/lange termijn tijdreeks). Op de Y-as kunnen we de intensiteit aflezen en op de X-as de frequentie.
Bij het cumulatieve periodogram kon er nog iets gezegd worden over de seizoenaliteit. Die is hier afwezig omdat er geen duidelijke trapfunctie te zien is.
Ook liggen de waarden/grafiek buiten het betrouwbaarheidsinterval van 95% (de twee stippelijnen), hierdoor kunnen we opmerken dat er geen sprake is van toeval.

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 time3 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 3 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=27840&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]3 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=27840&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=27840&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 time3 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24


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