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 computationFri, 11 Dec 2009 07:34:45 -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/2009/Dec/11/t1260542405ca2whjd9atyfrs1.htm/, Retrieved Mon, 29 Apr 2024 02:11:47 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=66256, Retrieved Mon, 29 Apr 2024 02:11:47 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact120
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Variance Reduction Matrix] [VRM - vaste rente...] [2008-12-12 10:51:08] [c5a66f1c8528a963efc2b82a8519f117]
- RM D  [Standard Deviation-Mean Plot] [SDMP inschrijving...] [2008-12-12 11:03:14] [c5a66f1c8528a963efc2b82a8519f117]
- RM      [Variance Reduction Matrix] [VRM - inschrijvin...] [2008-12-12 11:08:27] [c5a66f1c8528a963efc2b82a8519f117]
- RMP       [(Partial) Autocorrelation Function] [ACF - inschrijvin...] [2008-12-12 11:14:36] [c5a66f1c8528a963efc2b82a8519f117]
- RMP         [Spectral Analysis] [SA - inschrijving...] [2008-12-12 11:17:12] [c5a66f1c8528a963efc2b82a8519f117]
-  MPD            [Spectral Analysis] [spectraalanalyse ...] [2009-12-11 14:34:45] [557d56ec4b06cd0135c259898de8ce95] [Current]
Feedback Forum

Post a new message
Dataseries X:
17.8
17.9
17.4
16.7
16
16.6
19.1
17.8
17.2
18.6
16.3
15.1
19.2
17.7
19.1
18
17.5
17.8
21.1
17.2
19.4
19.8
17.6
16.2
19.5
19.9
20
17.3
18.9
18.6
21.4
18.6
19.8
20.8
19.6
17.7
19.8
22.2
20.7
17.9
20.9
21.2
21.4
23
21.3
23.9
22.4
18.3
22.8
22.3
17.8
16.4
16
16.4
17.7
16.6
16.2
18.3
17.6
15.1




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

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

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







Raw Periodogram
ParameterValue
Box-Cox transformation parameter (lambda)-1.1
Degree of non-seasonal differencing (d)1
Degree of seasonal differencing (D)1
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0208 (48)2e-06
0.0417 (24)5e-06
0.0625 (16)4e-06
0.0833 (12)5e-06
0.1042 (9.6)0
0.125 (8)2e-06
0.1458 (6.8571)2e-06
0.1667 (6)1e-06
0.1875 (5.3333)1e-06
0.2083 (4.8)7e-06
0.2292 (4.3636)1e-06
0.25 (4)3e-06
0.2708 (3.6923)6e-06
0.2917 (3.4286)5e-06
0.3125 (3.2)1.1e-05
0.3333 (3)8e-06
0.3542 (2.8235)0.000104
0.375 (2.6667)1e-05
0.3958 (2.5263)6e-05
0.4167 (2.4)3e-06
0.4375 (2.2857)8e-06
0.4583 (2.1818)0
0.4792 (2.087)3.7e-05
0.5 (2)2e-06

\begin{tabular}{lllllllll}
\hline
Raw Periodogram \tabularnewline
Parameter & Value \tabularnewline
Box-Cox transformation parameter (lambda) & -1.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) & 2e-06 \tabularnewline
0.0417 (24) & 5e-06 \tabularnewline
0.0625 (16) & 4e-06 \tabularnewline
0.0833 (12) & 5e-06 \tabularnewline
0.1042 (9.6) & 0 \tabularnewline
0.125 (8) & 2e-06 \tabularnewline
0.1458 (6.8571) & 2e-06 \tabularnewline
0.1667 (6) & 1e-06 \tabularnewline
0.1875 (5.3333) & 1e-06 \tabularnewline
0.2083 (4.8) & 7e-06 \tabularnewline
0.2292 (4.3636) & 1e-06 \tabularnewline
0.25 (4) & 3e-06 \tabularnewline
0.2708 (3.6923) & 6e-06 \tabularnewline
0.2917 (3.4286) & 5e-06 \tabularnewline
0.3125 (3.2) & 1.1e-05 \tabularnewline
0.3333 (3) & 8e-06 \tabularnewline
0.3542 (2.8235) & 0.000104 \tabularnewline
0.375 (2.6667) & 1e-05 \tabularnewline
0.3958 (2.5263) & 6e-05 \tabularnewline
0.4167 (2.4) & 3e-06 \tabularnewline
0.4375 (2.2857) & 8e-06 \tabularnewline
0.4583 (2.1818) & 0 \tabularnewline
0.4792 (2.087) & 3.7e-05 \tabularnewline
0.5 (2) & 2e-06 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=66256&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.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]2e-06[/C][/ROW]
[ROW][C]0.0417 (24)[/C][C]5e-06[/C][/ROW]
[ROW][C]0.0625 (16)[/C][C]4e-06[/C][/ROW]
[ROW][C]0.0833 (12)[/C][C]5e-06[/C][/ROW]
[ROW][C]0.1042 (9.6)[/C][C]0[/C][/ROW]
[ROW][C]0.125 (8)[/C][C]2e-06[/C][/ROW]
[ROW][C]0.1458 (6.8571)[/C][C]2e-06[/C][/ROW]
[ROW][C]0.1667 (6)[/C][C]1e-06[/C][/ROW]
[ROW][C]0.1875 (5.3333)[/C][C]1e-06[/C][/ROW]
[ROW][C]0.2083 (4.8)[/C][C]7e-06[/C][/ROW]
[ROW][C]0.2292 (4.3636)[/C][C]1e-06[/C][/ROW]
[ROW][C]0.25 (4)[/C][C]3e-06[/C][/ROW]
[ROW][C]0.2708 (3.6923)[/C][C]6e-06[/C][/ROW]
[ROW][C]0.2917 (3.4286)[/C][C]5e-06[/C][/ROW]
[ROW][C]0.3125 (3.2)[/C][C]1.1e-05[/C][/ROW]
[ROW][C]0.3333 (3)[/C][C]8e-06[/C][/ROW]
[ROW][C]0.3542 (2.8235)[/C][C]0.000104[/C][/ROW]
[ROW][C]0.375 (2.6667)[/C][C]1e-05[/C][/ROW]
[ROW][C]0.3958 (2.5263)[/C][C]6e-05[/C][/ROW]
[ROW][C]0.4167 (2.4)[/C][C]3e-06[/C][/ROW]
[ROW][C]0.4375 (2.2857)[/C][C]8e-06[/C][/ROW]
[ROW][C]0.4583 (2.1818)[/C][C]0[/C][/ROW]
[ROW][C]0.4792 (2.087)[/C][C]3.7e-05[/C][/ROW]
[ROW][C]0.5 (2)[/C][C]2e-06[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=66256&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=66256&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.1
Degree of non-seasonal differencing (d)1
Degree of seasonal differencing (D)1
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0208 (48)2e-06
0.0417 (24)5e-06
0.0625 (16)4e-06
0.0833 (12)5e-06
0.1042 (9.6)0
0.125 (8)2e-06
0.1458 (6.8571)2e-06
0.1667 (6)1e-06
0.1875 (5.3333)1e-06
0.2083 (4.8)7e-06
0.2292 (4.3636)1e-06
0.25 (4)3e-06
0.2708 (3.6923)6e-06
0.2917 (3.4286)5e-06
0.3125 (3.2)1.1e-05
0.3333 (3)8e-06
0.3542 (2.8235)0.000104
0.375 (2.6667)1e-05
0.3958 (2.5263)6e-05
0.4167 (2.4)3e-06
0.4375 (2.2857)8e-06
0.4583 (2.1818)0
0.4792 (2.087)3.7e-05
0.5 (2)2e-06



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