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 computationTue, 07 Dec 2010 09:00:05 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2010/Dec/07/t12917123560vy8bltg0wgvfpw.htm/, Retrieved Sat, 04 May 2024 00:06:42 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=106036, Retrieved Sat, 04 May 2024 00:06:42 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact128
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [data set] [2008-12-01 19:54:57] [b98453cac15ba1066b407e146608df68]
- RMP   [Spectral Analysis] [Unemployment] [2010-11-29 09:21:38] [b98453cac15ba1066b407e146608df68]
-    D      [Spectral Analysis] [WS9 Spectrum] [2010-12-07 09:00:05] [67e3c2d70de1dbb070b545ca6c893d5e] [Current]
-   PD        [Spectral Analysis] [WS9 Spectrum D=d=1] [2010-12-07 09:24:10] [07a238a5afc23eb944f8545182f29d5a]
Feedback Forum

Post a new message
Dataseries X:
562.325 
560.854 
555.332 
543.599 
536.662 
542.722 
593.530 
610.763 
612.613 
611.324 
594.167 
595.454 
590.865 
589.379 
584.428 
573.100 
567.456 
569.028 
620.735 
628.884 
628.232 
612.117 
595.404 
597.141 
593.408 
590.072 
579.799 
574.205 
572.775 
572.942 
619.567 
625.809 
619.916 
587.625 
565.742 
557.274 
560.576 
548.854 
531.673 
525.919 
511.038 
498.662 
555.362 
564.591 
541.657 
527.070 
509.846 
514.258 
516.922 
507.561 
492.622 
490.243 
469.357 
477.580 
528.379 
533.590 
517.945 
506.174 
501.866 
516.141 
528.222 
532.638 
536.322 
536.535 
523.597 
536.214 
586.570 
596.594 
580.523 
564.478 
557.560 
575.093 




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk

\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 & 'RServer@AstonUniversity' @ vre.aston.ac.uk \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=106036&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]'RServer@AstonUniversity' @ vre.aston.ac.uk[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=106036&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=106036&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 time1 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk







Raw Periodogram
ParameterValue
Box-Cox transformation parameter (lambda)1
Degree of non-seasonal differencing (d)0
Degree of seasonal differencing (D)0
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0139 (72)12900.281427
0.0278 (36)6917.386817
0.0417 (24)938.096366
0.0556 (18)1518.185475
0.0694 (14.4)207.130335
0.0833 (12)7135.534685
0.0972 (10.2857)482.223604
0.1111 (9)167.725283
0.125 (8)27.17485
0.1389 (7.2)51.813418
0.1528 (6.5455)40.177759
0.1667 (6)5530.202698
0.1806 (5.5385)123.023946
0.1944 (5.1429)21.204353
0.2083 (4.8)20.050946
0.2222 (4.5)1.41595
0.2361 (4.2353)5.585053
0.25 (4)787.553517
0.2639 (3.7895)110.857821
0.2778 (3.6)4.984687
0.2917 (3.4286)10.979276
0.3056 (3.2727)4.087461
0.3194 (3.1304)9.210075
0.3333 (3)457.698943
0.3472 (2.88)34.046295
0.3611 (2.7692)27.5949
0.375 (2.6667)1.300353
0.3889 (2.5714)4.068703
0.4028 (2.4828)2.830611
0.4167 (2.4)227.091892
0.4306 (2.3226)21.964992
0.4444 (2.25)13.022081
0.4583 (2.1818)2.200669
0.4722 (2.1176)9.885945
0.4861 (2.0571)14.534196
0.5 (2)64.393594

\begin{tabular}{lllllllll}
\hline
Raw Periodogram \tabularnewline
Parameter & Value \tabularnewline
Box-Cox transformation parameter (lambda) & 1 \tabularnewline
Degree of non-seasonal differencing (d) & 0 \tabularnewline
Degree of seasonal differencing (D) & 0 \tabularnewline
Seasonal Period (s) & 12 \tabularnewline
Frequency (Period) & Spectrum \tabularnewline
0.0139 (72) & 12900.281427 \tabularnewline
0.0278 (36) & 6917.386817 \tabularnewline
0.0417 (24) & 938.096366 \tabularnewline
0.0556 (18) & 1518.185475 \tabularnewline
0.0694 (14.4) & 207.130335 \tabularnewline
0.0833 (12) & 7135.534685 \tabularnewline
0.0972 (10.2857) & 482.223604 \tabularnewline
0.1111 (9) & 167.725283 \tabularnewline
0.125 (8) & 27.17485 \tabularnewline
0.1389 (7.2) & 51.813418 \tabularnewline
0.1528 (6.5455) & 40.177759 \tabularnewline
0.1667 (6) & 5530.202698 \tabularnewline
0.1806 (5.5385) & 123.023946 \tabularnewline
0.1944 (5.1429) & 21.204353 \tabularnewline
0.2083 (4.8) & 20.050946 \tabularnewline
0.2222 (4.5) & 1.41595 \tabularnewline
0.2361 (4.2353) & 5.585053 \tabularnewline
0.25 (4) & 787.553517 \tabularnewline
0.2639 (3.7895) & 110.857821 \tabularnewline
0.2778 (3.6) & 4.984687 \tabularnewline
0.2917 (3.4286) & 10.979276 \tabularnewline
0.3056 (3.2727) & 4.087461 \tabularnewline
0.3194 (3.1304) & 9.210075 \tabularnewline
0.3333 (3) & 457.698943 \tabularnewline
0.3472 (2.88) & 34.046295 \tabularnewline
0.3611 (2.7692) & 27.5949 \tabularnewline
0.375 (2.6667) & 1.300353 \tabularnewline
0.3889 (2.5714) & 4.068703 \tabularnewline
0.4028 (2.4828) & 2.830611 \tabularnewline
0.4167 (2.4) & 227.091892 \tabularnewline
0.4306 (2.3226) & 21.964992 \tabularnewline
0.4444 (2.25) & 13.022081 \tabularnewline
0.4583 (2.1818) & 2.200669 \tabularnewline
0.4722 (2.1176) & 9.885945 \tabularnewline
0.4861 (2.0571) & 14.534196 \tabularnewline
0.5 (2) & 64.393594 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=106036&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]0[/C][/ROW]
[ROW][C]Degree of seasonal differencing (D)[/C][C]0[/C][/ROW]
[ROW][C]Seasonal Period (s)[/C][C]12[/C][/ROW]
[ROW][C]Frequency (Period)[/C][C]Spectrum[/C][/ROW]
[ROW][C]0.0139 (72)[/C][C]12900.281427[/C][/ROW]
[ROW][C]0.0278 (36)[/C][C]6917.386817[/C][/ROW]
[ROW][C]0.0417 (24)[/C][C]938.096366[/C][/ROW]
[ROW][C]0.0556 (18)[/C][C]1518.185475[/C][/ROW]
[ROW][C]0.0694 (14.4)[/C][C]207.130335[/C][/ROW]
[ROW][C]0.0833 (12)[/C][C]7135.534685[/C][/ROW]
[ROW][C]0.0972 (10.2857)[/C][C]482.223604[/C][/ROW]
[ROW][C]0.1111 (9)[/C][C]167.725283[/C][/ROW]
[ROW][C]0.125 (8)[/C][C]27.17485[/C][/ROW]
[ROW][C]0.1389 (7.2)[/C][C]51.813418[/C][/ROW]
[ROW][C]0.1528 (6.5455)[/C][C]40.177759[/C][/ROW]
[ROW][C]0.1667 (6)[/C][C]5530.202698[/C][/ROW]
[ROW][C]0.1806 (5.5385)[/C][C]123.023946[/C][/ROW]
[ROW][C]0.1944 (5.1429)[/C][C]21.204353[/C][/ROW]
[ROW][C]0.2083 (4.8)[/C][C]20.050946[/C][/ROW]
[ROW][C]0.2222 (4.5)[/C][C]1.41595[/C][/ROW]
[ROW][C]0.2361 (4.2353)[/C][C]5.585053[/C][/ROW]
[ROW][C]0.25 (4)[/C][C]787.553517[/C][/ROW]
[ROW][C]0.2639 (3.7895)[/C][C]110.857821[/C][/ROW]
[ROW][C]0.2778 (3.6)[/C][C]4.984687[/C][/ROW]
[ROW][C]0.2917 (3.4286)[/C][C]10.979276[/C][/ROW]
[ROW][C]0.3056 (3.2727)[/C][C]4.087461[/C][/ROW]
[ROW][C]0.3194 (3.1304)[/C][C]9.210075[/C][/ROW]
[ROW][C]0.3333 (3)[/C][C]457.698943[/C][/ROW]
[ROW][C]0.3472 (2.88)[/C][C]34.046295[/C][/ROW]
[ROW][C]0.3611 (2.7692)[/C][C]27.5949[/C][/ROW]
[ROW][C]0.375 (2.6667)[/C][C]1.300353[/C][/ROW]
[ROW][C]0.3889 (2.5714)[/C][C]4.068703[/C][/ROW]
[ROW][C]0.4028 (2.4828)[/C][C]2.830611[/C][/ROW]
[ROW][C]0.4167 (2.4)[/C][C]227.091892[/C][/ROW]
[ROW][C]0.4306 (2.3226)[/C][C]21.964992[/C][/ROW]
[ROW][C]0.4444 (2.25)[/C][C]13.022081[/C][/ROW]
[ROW][C]0.4583 (2.1818)[/C][C]2.200669[/C][/ROW]
[ROW][C]0.4722 (2.1176)[/C][C]9.885945[/C][/ROW]
[ROW][C]0.4861 (2.0571)[/C][C]14.534196[/C][/ROW]
[ROW][C]0.5 (2)[/C][C]64.393594[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=106036&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=106036&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)0
Degree of seasonal differencing (D)0
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0139 (72)12900.281427
0.0278 (36)6917.386817
0.0417 (24)938.096366
0.0556 (18)1518.185475
0.0694 (14.4)207.130335
0.0833 (12)7135.534685
0.0972 (10.2857)482.223604
0.1111 (9)167.725283
0.125 (8)27.17485
0.1389 (7.2)51.813418
0.1528 (6.5455)40.177759
0.1667 (6)5530.202698
0.1806 (5.5385)123.023946
0.1944 (5.1429)21.204353
0.2083 (4.8)20.050946
0.2222 (4.5)1.41595
0.2361 (4.2353)5.585053
0.25 (4)787.553517
0.2639 (3.7895)110.857821
0.2778 (3.6)4.984687
0.2917 (3.4286)10.979276
0.3056 (3.2727)4.087461
0.3194 (3.1304)9.210075
0.3333 (3)457.698943
0.3472 (2.88)34.046295
0.3611 (2.7692)27.5949
0.375 (2.6667)1.300353
0.3889 (2.5714)4.068703
0.4028 (2.4828)2.830611
0.4167 (2.4)227.091892
0.4306 (2.3226)21.964992
0.4444 (2.25)13.022081
0.4583 (2.1818)2.200669
0.4722 (2.1176)9.885945
0.4861 (2.0571)14.534196
0.5 (2)64.393594



Parameters (Session):
par1 = 1 ; par2 = 0 ; par3 = 0 ; par4 = 12 ;
Parameters (R input):
par1 = 1 ; par2 = 0 ; par3 = 0 ; 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')