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, 12 Dec 2008 05:50:04 -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/12/t122908624712u9d4shl6qtzng.htm/, Retrieved Fri, 17 May 2024 02:39:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=32645, Retrieved Fri, 17 May 2024 02:39:27 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact183
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Spectral Analysis] [specturm vrouwen ...] [2008-12-12 12:36:12] [e43247bc0ab243a5af99ac7f55ba0b41]
-   P   [Spectral Analysis] [vrouwen met aanpa...] [2008-12-12 12:41:08] [e43247bc0ab243a5af99ac7f55ba0b41]
-   P       [Spectral Analysis] [vrouwen met juist...] [2008-12-12 12:50:04] [f24298b2e4c2a19d76cf4460ec5d2246] [Current]
-   P         [Spectral Analysis] [spectrum met juis...] [2008-12-17 17:58:18] [e43247bc0ab243a5af99ac7f55ba0b41]
Feedback Forum

Post a new message
Dataseries X:
9.0
9.1
8.7
8.2
7.9
7.9
9.1
9.4
9.5
9.1
9.0
9.3
9.9
9.8
9.4
8.3
8.0
8.5
10.4
11.1
10.9
9.9
9.2
9.2
9.5
9.6
9.5
9.1
8.9
9.0
10.1
10.3
10.2
9.6
9.2
9.3
9.4
9.4
9.2
9.0
9.0
9.0
9.8
10.0
9.9
9.3
9.0
9.0
9.1
9.1
9.1
9.2
8.8
8.3
8.4
8.1
7.8
7.9
7.9
8.0
7.9
7.5
7.2
6.9
6.6
6.7
7.3
7.5




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=32645&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=32645&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=32645&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)0
Degree of non-seasonal differencing (d)2
Degree of seasonal differencing (D)1
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0185 (54)0.000108
0.037 (27)0.000116
0.0556 (18)9.4e-05
0.0741 (13.5)0.000308
0.0926 (10.8)0.003157
0.1111 (9)0.000159
0.1296 (7.7143)0.010001
0.1481 (6.75)0.003846
0.1667 (6)0.003073
0.1852 (5.4)0.001632
0.2037 (4.9091)0.001507
0.2222 (4.5)0.000485
0.2407 (4.1538)0.000382
0.2593 (3.8571)0.000431
0.2778 (3.6)0.001068
0.2963 (3.375)0.000968
0.3148 (3.1765)4.4e-05
0.3333 (3)6.7e-05
0.3519 (2.8421)0.001368
0.3704 (2.7)0.001557
0.3889 (2.5714)4e-06
0.4074 (2.4545)0.000157
0.4259 (2.3478)0.000867
0.4444 (2.25)0.001088
0.463 (2.16)0.00157
0.4815 (2.0769)8.2e-05
0.5 (2)0.003076

\begin{tabular}{lllllllll}
\hline
Raw Periodogram \tabularnewline
Parameter & Value \tabularnewline
Box-Cox transformation parameter (lambda) & 0 \tabularnewline
Degree of non-seasonal differencing (d) & 2 \tabularnewline
Degree of seasonal differencing (D) & 1 \tabularnewline
Seasonal Period (s) & 12 \tabularnewline
Frequency (Period) & Spectrum \tabularnewline
0.0185 (54) & 0.000108 \tabularnewline
0.037 (27) & 0.000116 \tabularnewline
0.0556 (18) & 9.4e-05 \tabularnewline
0.0741 (13.5) & 0.000308 \tabularnewline
0.0926 (10.8) & 0.003157 \tabularnewline
0.1111 (9) & 0.000159 \tabularnewline
0.1296 (7.7143) & 0.010001 \tabularnewline
0.1481 (6.75) & 0.003846 \tabularnewline
0.1667 (6) & 0.003073 \tabularnewline
0.1852 (5.4) & 0.001632 \tabularnewline
0.2037 (4.9091) & 0.001507 \tabularnewline
0.2222 (4.5) & 0.000485 \tabularnewline
0.2407 (4.1538) & 0.000382 \tabularnewline
0.2593 (3.8571) & 0.000431 \tabularnewline
0.2778 (3.6) & 0.001068 \tabularnewline
0.2963 (3.375) & 0.000968 \tabularnewline
0.3148 (3.1765) & 4.4e-05 \tabularnewline
0.3333 (3) & 6.7e-05 \tabularnewline
0.3519 (2.8421) & 0.001368 \tabularnewline
0.3704 (2.7) & 0.001557 \tabularnewline
0.3889 (2.5714) & 4e-06 \tabularnewline
0.4074 (2.4545) & 0.000157 \tabularnewline
0.4259 (2.3478) & 0.000867 \tabularnewline
0.4444 (2.25) & 0.001088 \tabularnewline
0.463 (2.16) & 0.00157 \tabularnewline
0.4815 (2.0769) & 8.2e-05 \tabularnewline
0.5 (2) & 0.003076 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=32645&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]0[/C][/ROW]
[ROW][C]Degree of non-seasonal differencing (d)[/C][C]2[/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.0185 (54)[/C][C]0.000108[/C][/ROW]
[ROW][C]0.037 (27)[/C][C]0.000116[/C][/ROW]
[ROW][C]0.0556 (18)[/C][C]9.4e-05[/C][/ROW]
[ROW][C]0.0741 (13.5)[/C][C]0.000308[/C][/ROW]
[ROW][C]0.0926 (10.8)[/C][C]0.003157[/C][/ROW]
[ROW][C]0.1111 (9)[/C][C]0.000159[/C][/ROW]
[ROW][C]0.1296 (7.7143)[/C][C]0.010001[/C][/ROW]
[ROW][C]0.1481 (6.75)[/C][C]0.003846[/C][/ROW]
[ROW][C]0.1667 (6)[/C][C]0.003073[/C][/ROW]
[ROW][C]0.1852 (5.4)[/C][C]0.001632[/C][/ROW]
[ROW][C]0.2037 (4.9091)[/C][C]0.001507[/C][/ROW]
[ROW][C]0.2222 (4.5)[/C][C]0.000485[/C][/ROW]
[ROW][C]0.2407 (4.1538)[/C][C]0.000382[/C][/ROW]
[ROW][C]0.2593 (3.8571)[/C][C]0.000431[/C][/ROW]
[ROW][C]0.2778 (3.6)[/C][C]0.001068[/C][/ROW]
[ROW][C]0.2963 (3.375)[/C][C]0.000968[/C][/ROW]
[ROW][C]0.3148 (3.1765)[/C][C]4.4e-05[/C][/ROW]
[ROW][C]0.3333 (3)[/C][C]6.7e-05[/C][/ROW]
[ROW][C]0.3519 (2.8421)[/C][C]0.001368[/C][/ROW]
[ROW][C]0.3704 (2.7)[/C][C]0.001557[/C][/ROW]
[ROW][C]0.3889 (2.5714)[/C][C]4e-06[/C][/ROW]
[ROW][C]0.4074 (2.4545)[/C][C]0.000157[/C][/ROW]
[ROW][C]0.4259 (2.3478)[/C][C]0.000867[/C][/ROW]
[ROW][C]0.4444 (2.25)[/C][C]0.001088[/C][/ROW]
[ROW][C]0.463 (2.16)[/C][C]0.00157[/C][/ROW]
[ROW][C]0.4815 (2.0769)[/C][C]8.2e-05[/C][/ROW]
[ROW][C]0.5 (2)[/C][C]0.003076[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=32645&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=32645&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)0
Degree of non-seasonal differencing (d)2
Degree of seasonal differencing (D)1
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0185 (54)0.000108
0.037 (27)0.000116
0.0556 (18)9.4e-05
0.0741 (13.5)0.000308
0.0926 (10.8)0.003157
0.1111 (9)0.000159
0.1296 (7.7143)0.010001
0.1481 (6.75)0.003846
0.1667 (6)0.003073
0.1852 (5.4)0.001632
0.2037 (4.9091)0.001507
0.2222 (4.5)0.000485
0.2407 (4.1538)0.000382
0.2593 (3.8571)0.000431
0.2778 (3.6)0.001068
0.2963 (3.375)0.000968
0.3148 (3.1765)4.4e-05
0.3333 (3)6.7e-05
0.3519 (2.8421)0.001368
0.3704 (2.7)0.001557
0.3889 (2.5714)4e-06
0.4074 (2.4545)0.000157
0.4259 (2.3478)0.000867
0.4444 (2.25)0.001088
0.463 (2.16)0.00157
0.4815 (2.0769)8.2e-05
0.5 (2)0.003076



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