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 computationThu, 17 Dec 2009 02:47:20 -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/17/t12610433105qf883yj3ls8top.htm/, Retrieved Tue, 30 Apr 2024 03:14:14 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=68675, Retrieved Tue, 30 Apr 2024 03:14:14 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact147
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Explorative Data Analysis] [Run Sequence gebo...] [2008-12-12 13:32:37] [76963dc1903f0f612b6153510a3818cf]
- R  D  [Univariate Explorative Data Analysis] [Run Sequence gebo...] [2008-12-17 12:14:40] [76963dc1903f0f612b6153510a3818cf]
-         [Univariate Explorative Data Analysis] [Run Sequence Plot...] [2008-12-22 18:19:51] [1ce0d16c8f4225c977b42c8fa93bc163]
- RMP       [Spectral Analysis] [Identifying Integ...] [2009-11-22 12:38:17] [b98453cac15ba1066b407e146608df68]
-    D        [Spectral Analysis] [spectrumanalyse] [2009-11-28 09:48:21] [7773f496f69461f4a67891f0ef752622]
-    D          [Spectral Analysis] [koffie en thee] [2009-12-16 19:32:35] [7773f496f69461f4a67891f0ef752622]
-    D            [Spectral Analysis] [Thee] [2009-12-17 09:32:15] [7773f496f69461f4a67891f0ef752622]
-   P                 [Spectral Analysis] [Thee] [2009-12-17 09:47:20] [5c2088b06970f9a7d6fea063ee8d5871] [Current]
- RMP                   [Standard Deviation-Mean Plot] [Thee] [2009-12-17 10:10:23] [7773f496f69461f4a67891f0ef752622]
Feedback Forum

Post a new message
Dataseries X:
1.56
1.57
1.56
1.56
1.56
1.56
1.55
1.56
1.55
1.54
1.54
1.53
1.53
1.53
1.53
1.52
1.52
1.52
1.52
1.52
1.51
1.51
1.51
1.52
1.51
1.5
1.5
1.5
1.5
1.5
1.5
1.49
1.47
1.45
1.45
1.44




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=68675&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
Degree of non-seasonal differencing (d)1
Degree of seasonal differencing (D)0
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0278 (36)6.9e-05
0.0556 (18)2.2e-05
0.0833 (12)4.7e-05
0.1111 (9)2.2e-05
0.1389 (7.2)8.1e-05
0.1667 (6)9.8e-05
0.1944 (5.1429)3.4e-05
0.2222 (4.5)5.7e-05
0.25 (4)6.1e-05
0.2778 (3.6)7e-06
0.3056 (3.2727)6.6e-05
0.3333 (3)6.3e-05
0.3611 (2.7692)3.5e-05
0.3889 (2.5714)1.1e-05
0.4167 (2.4)3.2e-05
0.4444 (2.25)0.000125
0.4722 (2.1176)0
0.5 (2)9e-06

\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) & 0 \tabularnewline
Seasonal Period (s) & 12 \tabularnewline
Frequency (Period) & Spectrum \tabularnewline
0.0278 (36) & 6.9e-05 \tabularnewline
0.0556 (18) & 2.2e-05 \tabularnewline
0.0833 (12) & 4.7e-05 \tabularnewline
0.1111 (9) & 2.2e-05 \tabularnewline
0.1389 (7.2) & 8.1e-05 \tabularnewline
0.1667 (6) & 9.8e-05 \tabularnewline
0.1944 (5.1429) & 3.4e-05 \tabularnewline
0.2222 (4.5) & 5.7e-05 \tabularnewline
0.25 (4) & 6.1e-05 \tabularnewline
0.2778 (3.6) & 7e-06 \tabularnewline
0.3056 (3.2727) & 6.6e-05 \tabularnewline
0.3333 (3) & 6.3e-05 \tabularnewline
0.3611 (2.7692) & 3.5e-05 \tabularnewline
0.3889 (2.5714) & 1.1e-05 \tabularnewline
0.4167 (2.4) & 3.2e-05 \tabularnewline
0.4444 (2.25) & 0.000125 \tabularnewline
0.4722 (2.1176) & 0 \tabularnewline
0.5 (2) & 9e-06 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=68675&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]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.0278 (36)[/C][C]6.9e-05[/C][/ROW]
[ROW][C]0.0556 (18)[/C][C]2.2e-05[/C][/ROW]
[ROW][C]0.0833 (12)[/C][C]4.7e-05[/C][/ROW]
[ROW][C]0.1111 (9)[/C][C]2.2e-05[/C][/ROW]
[ROW][C]0.1389 (7.2)[/C][C]8.1e-05[/C][/ROW]
[ROW][C]0.1667 (6)[/C][C]9.8e-05[/C][/ROW]
[ROW][C]0.1944 (5.1429)[/C][C]3.4e-05[/C][/ROW]
[ROW][C]0.2222 (4.5)[/C][C]5.7e-05[/C][/ROW]
[ROW][C]0.25 (4)[/C][C]6.1e-05[/C][/ROW]
[ROW][C]0.2778 (3.6)[/C][C]7e-06[/C][/ROW]
[ROW][C]0.3056 (3.2727)[/C][C]6.6e-05[/C][/ROW]
[ROW][C]0.3333 (3)[/C][C]6.3e-05[/C][/ROW]
[ROW][C]0.3611 (2.7692)[/C][C]3.5e-05[/C][/ROW]
[ROW][C]0.3889 (2.5714)[/C][C]1.1e-05[/C][/ROW]
[ROW][C]0.4167 (2.4)[/C][C]3.2e-05[/C][/ROW]
[ROW][C]0.4444 (2.25)[/C][C]0.000125[/C][/ROW]
[ROW][C]0.4722 (2.1176)[/C][C]0[/C][/ROW]
[ROW][C]0.5 (2)[/C][C]9e-06[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=68675&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=68675&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)0
Seasonal Period (s)12
Frequency (Period)Spectrum
0.0278 (36)6.9e-05
0.0556 (18)2.2e-05
0.0833 (12)4.7e-05
0.1111 (9)2.2e-05
0.1389 (7.2)8.1e-05
0.1667 (6)9.8e-05
0.1944 (5.1429)3.4e-05
0.2222 (4.5)5.7e-05
0.25 (4)6.1e-05
0.2778 (3.6)7e-06
0.3056 (3.2727)6.6e-05
0.3333 (3)6.3e-05
0.3611 (2.7692)3.5e-05
0.3889 (2.5714)1.1e-05
0.4167 (2.4)3.2e-05
0.4444 (2.25)0.000125
0.4722 (2.1176)0
0.5 (2)9e-06



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