Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_spectrum.wasp
Title produced by softwareSpectral Analysis
Date of computationMon, 08 Dec 2008 14:26:21 -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/08/t12287716053h35p0tbhm7xh2f.htm/, Retrieved Thu, 16 May 2024 16:10:53 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=31062, Retrieved Thu, 16 May 2024 16:10:53 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact190
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Spectral Analysis] [Unemployment - St...] [2008-12-08 17:28:52] [57850c80fd59ccfb28f882be994e814e]
F    D    [Spectral Analysis] [] [2008-12-08 21:26:21] [6d40a467de0f28bd2350f82ac9522c51] [Current]
Feedback Forum
2008-12-15 08:26:47 [Nathalie Koulouris] [reply
De student heeft vervolgens de juiste berekeningsmethode gebruikt maar ligt zijn antwoord te weinig toe.

Post a new message
Dataseries X:
299,63
305,945
382,252
348,846
335,367
373,617
312,612
312,232
337,161
331,476
350,103
345,127
297,256
295,979
361,007
321,803
354,937
349,432
290,979
349,576
327,625
349,377
336,777
339,134
323,321
318,86
373,583
333,03
408,556
414,646
291,514
348,857
349,368
375,765
364,136
349,53
348,167
332,856
360,551
346,969
392,815
372,02
371,027
342,672
367,343
390,786
343,785
362,6
349,468
340,624
369,536
407,782
392,239
404,824
373,669
344,902
396,7
398,911
366,009
392,484




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 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 & 2 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=31062&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]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=31062&T=0

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







Raw Periodogram
ParameterValue
Box-Cox transformation parameter (lambda)1
Degree of non-seasonal differencing (d)1
Degree of seasonal differencing (D)1
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0208 (48)55.308268
0.0417 (24)58.416194
0.0625 (16)8.9813
0.0833 (12)36.073229
0.1042 (9.6)20.57184
0.125 (8)840.957569
0.1458 (6.8571)279.88906
0.1667 (6)4.941791
0.1875 (5.3333)15.796649
0.2083 (4.8)1105.840631
0.2292 (4.3636)2320.631364
0.25 (4)585.370156
0.2708 (3.6923)419.582525
0.2917 (3.4286)1123.755078
0.3125 (3.2)2058.2125
0.3333 (3)1215.204528
0.3542 (2.8235)5285.295119
0.375 (2.6667)5563.108289
0.3958 (2.5263)3549.139282
0.4167 (2.4)399.219213
0.4375 (2.2857)5079.019853
0.4583 (2.1818)2906.436783
0.4792 (2.087)672.275163
0.5 (2)100.873706

\begin{tabular}{lllllllll}
\hline
Raw Periodogram \tabularnewline
Parameter & Value \tabularnewline
Box-Cox transformation parameter (lambda) & 1 \tabularnewline
Degree of non-seasonal differencing (d) & 1 \tabularnewline
Degree of seasonal differencing (D) & 1 \tabularnewline
Seasonal Period (s) & 12 \tabularnewline
Frequency (Period) & Spectrum \tabularnewline
0.0208 (48) & 55.308268 \tabularnewline
0.0417 (24) & 58.416194 \tabularnewline
0.0625 (16) & 8.9813 \tabularnewline
0.0833 (12) & 36.073229 \tabularnewline
0.1042 (9.6) & 20.57184 \tabularnewline
0.125 (8) & 840.957569 \tabularnewline
0.1458 (6.8571) & 279.88906 \tabularnewline
0.1667 (6) & 4.941791 \tabularnewline
0.1875 (5.3333) & 15.796649 \tabularnewline
0.2083 (4.8) & 1105.840631 \tabularnewline
0.2292 (4.3636) & 2320.631364 \tabularnewline
0.25 (4) & 585.370156 \tabularnewline
0.2708 (3.6923) & 419.582525 \tabularnewline
0.2917 (3.4286) & 1123.755078 \tabularnewline
0.3125 (3.2) & 2058.2125 \tabularnewline
0.3333 (3) & 1215.204528 \tabularnewline
0.3542 (2.8235) & 5285.295119 \tabularnewline
0.375 (2.6667) & 5563.108289 \tabularnewline
0.3958 (2.5263) & 3549.139282 \tabularnewline
0.4167 (2.4) & 399.219213 \tabularnewline
0.4375 (2.2857) & 5079.019853 \tabularnewline
0.4583 (2.1818) & 2906.436783 \tabularnewline
0.4792 (2.087) & 672.275163 \tabularnewline
0.5 (2) & 100.873706 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=31062&T=1

[TABLE]
[ROW][C]Raw Periodogram[/C][/ROW]
[ROW][C]Parameter[/C][C]Value[/C][/ROW]
[ROW][C]Box-Cox transformation parameter (lambda)[/C][C]1[/C][/ROW]
[ROW][C]Degree of non-seasonal differencing (d)[/C][C]1[/C][/ROW]
[ROW][C]Degree of seasonal differencing (D)[/C][C]1[/C][/ROW]
[ROW][C]Seasonal Period (s)[/C][C]12[/C][/ROW]
[ROW][C]Frequency (Period)[/C][C]Spectrum[/C][/ROW]
[ROW][C]0.0208 (48)[/C][C]55.308268[/C][/ROW]
[ROW][C]0.0417 (24)[/C][C]58.416194[/C][/ROW]
[ROW][C]0.0625 (16)[/C][C]8.9813[/C][/ROW]
[ROW][C]0.0833 (12)[/C][C]36.073229[/C][/ROW]
[ROW][C]0.1042 (9.6)[/C][C]20.57184[/C][/ROW]
[ROW][C]0.125 (8)[/C][C]840.957569[/C][/ROW]
[ROW][C]0.1458 (6.8571)[/C][C]279.88906[/C][/ROW]
[ROW][C]0.1667 (6)[/C][C]4.941791[/C][/ROW]
[ROW][C]0.1875 (5.3333)[/C][C]15.796649[/C][/ROW]
[ROW][C]0.2083 (4.8)[/C][C]1105.840631[/C][/ROW]
[ROW][C]0.2292 (4.3636)[/C][C]2320.631364[/C][/ROW]
[ROW][C]0.25 (4)[/C][C]585.370156[/C][/ROW]
[ROW][C]0.2708 (3.6923)[/C][C]419.582525[/C][/ROW]
[ROW][C]0.2917 (3.4286)[/C][C]1123.755078[/C][/ROW]
[ROW][C]0.3125 (3.2)[/C][C]2058.2125[/C][/ROW]
[ROW][C]0.3333 (3)[/C][C]1215.204528[/C][/ROW]
[ROW][C]0.3542 (2.8235)[/C][C]5285.295119[/C][/ROW]
[ROW][C]0.375 (2.6667)[/C][C]5563.108289[/C][/ROW]
[ROW][C]0.3958 (2.5263)[/C][C]3549.139282[/C][/ROW]
[ROW][C]0.4167 (2.4)[/C][C]399.219213[/C][/ROW]
[ROW][C]0.4375 (2.2857)[/C][C]5079.019853[/C][/ROW]
[ROW][C]0.4583 (2.1818)[/C][C]2906.436783[/C][/ROW]
[ROW][C]0.4792 (2.087)[/C][C]672.275163[/C][/ROW]
[ROW][C]0.5 (2)[/C][C]100.873706[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=31062&T=1

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

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Raw Periodogram
ParameterValue
Box-Cox transformation parameter (lambda)1
Degree of non-seasonal differencing (d)1
Degree of seasonal differencing (D)1
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0208 (48)55.308268
0.0417 (24)58.416194
0.0625 (16)8.9813
0.0833 (12)36.073229
0.1042 (9.6)20.57184
0.125 (8)840.957569
0.1458 (6.8571)279.88906
0.1667 (6)4.941791
0.1875 (5.3333)15.796649
0.2083 (4.8)1105.840631
0.2292 (4.3636)2320.631364
0.25 (4)585.370156
0.2708 (3.6923)419.582525
0.2917 (3.4286)1123.755078
0.3125 (3.2)2058.2125
0.3333 (3)1215.204528
0.3542 (2.8235)5285.295119
0.375 (2.6667)5563.108289
0.3958 (2.5263)3549.139282
0.4167 (2.4)399.219213
0.4375 (2.2857)5079.019853
0.4583 (2.1818)2906.436783
0.4792 (2.087)672.275163
0.5 (2)100.873706



Parameters (Session):
par1 = 1 ; par2 = 1 ; par3 = 1 ; par4 = 12 ;
Parameters (R input):
par1 = 1 ; par2 = 1 ; par3 = 1 ; par4 = 12 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
par3 <- as.numeric(par3)
par4 <- as.numeric(par4)
if (par1 == 0) {
x <- log(x)
} else {
x <- (x ^ par1 - 1) / par1
}
if (par2 > 0) x <- diff(x,lag=1,difference=par2)
if (par3 > 0) x <- diff(x,lag=par4,difference=par3)
bitmap(file='test1.png')
r <- spectrum(x,main='Raw Periodogram')
dev.off()
bitmap(file='test2.png')
cpgram(x,main='Cumulative Periodogram')
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Raw Periodogram',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Parameter',header=TRUE)
a<-table.element(a,'Value',header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Box-Cox transformation parameter (lambda)',header=TRUE)
a<-table.element(a,par1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of non-seasonal differencing (d)',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Degree of seasonal differencing (D)',header=TRUE)
a<-table.element(a,par3)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Seasonal Period (s)',header=TRUE)
a<-table.element(a,par4)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Frequency (Period)',header=TRUE)
a<-table.element(a,'Spectrum',header=TRUE)
a<-table.row.end(a)
for (i in 1:length(r$freq)) {
a<-table.row.start(a)
mylab <- round(r$freq[i],4)
mylab <- paste(mylab,' (',sep='')
mylab <- paste(mylab,round(1/r$freq[i],4),sep='')
mylab <- paste(mylab,')',sep='')
a<-table.element(a,mylab,header=TRUE)
a<-table.element(a,round(r$spec[i],6))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')