FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_spectrum.wasp
Title produced by softwareSpectral Analysis
Date of computationSat, 18 Dec 2010 18:30:58 +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/18/t1292697049x7d3vka9t8ms9mt.htm/, Retrieved Tue, 30 Apr 2024 01:18:11 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=112151, Retrieved Tue, 30 Apr 2024 01:18:11 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact136

Tree of Dependent Computations

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   [(Partial) Autocorrelation Function] [Workshop 8 Autoco...] [2010-11-28 16:39:23] [945bcebba5e7ac34a41d6888338a1ba9]
- RMP     [Classical Decomposition] [workshop 8 Klassi...] [2010-11-28 17:01:52] [945bcebba5e7ac34a41d6888338a1ba9]
- RMPD      [Exponential Smoothing] [Paper TSA Exponen...] [2010-12-18 18:11:45] [945bcebba5e7ac34a41d6888338a1ba9]
- RMP           [Spectral Analysis] [Paper TSA Spectra...] [2010-12-18 18:30:58] [514029464b0621595fe21c9fa38c7009] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
67
189
342
432
517
623
605
716
677
710
839
886
891
917
820
793
932
906
844
801
957
1159
1264
1097
1240
1411
1535
1862
1894
2239
2465
2423
2692
2856
3450
4162
4260
4225
4092
4160
3896
3628
3754
3749
3907
4449
5272
6197
6446
7157
7559
7674
6929
7156
6805
7095
7222
7593
7910

Tables (Output of Computation)





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

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

As an alternative you can also use a QR Code:  
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'Gwilym Jenkins' @ 72.249.127.135






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)1
Frequency (Period)Spectrum
0.0167 (60)10310840.788791
0.0333 (30)169176.871747
0.05 (20)2857425.895077
0.0667 (15)3699504.834075
0.0833 (12)1364201.571362
0.1 (10)217675.047079
0.1167 (8.5714)113778.371838
0.1333 (7.5)107934.299205
0.15 (6.6667)107751.290143
0.1667 (6)13855.370532
0.1833 (5.4545)13388.697779
0.2 (5)44429.936135
0.2167 (4.6154)17250.736887
0.2333 (4.2857)34017.464656
0.25 (4)28301.26499
0.2667 (3.75)35190.455579
0.2833 (3.5294)26426.92708
0.3 (3.3333)8567.095278
0.3167 (3.1579)879.950829
0.3333 (3)4563.722788
0.35 (2.8571)19187.872557
0.3667 (2.7273)3231.510441
0.3833 (2.6087)15125.818086
0.4 (2.5)6779.040633
0.4167 (2.4)12196.410134
0.4333 (2.3077)12166.174353
0.45 (2.2222)13705.236326
0.4667 (2.1429)10644.788941
0.4833 (2.069)25736.286139
0.5 (2)27685.968618

\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) & 1 \tabularnewline
Frequency (Period) & Spectrum \tabularnewline
0.0167 (60) & 10310840.788791 \tabularnewline
0.0333 (30) & 169176.871747 \tabularnewline
0.05 (20) & 2857425.895077 \tabularnewline
0.0667 (15) & 3699504.834075 \tabularnewline
0.0833 (12) & 1364201.571362 \tabularnewline
0.1 (10) & 217675.047079 \tabularnewline
0.1167 (8.5714) & 113778.371838 \tabularnewline
0.1333 (7.5) & 107934.299205 \tabularnewline
0.15 (6.6667) & 107751.290143 \tabularnewline
0.1667 (6) & 13855.370532 \tabularnewline
0.1833 (5.4545) & 13388.697779 \tabularnewline
0.2 (5) & 44429.936135 \tabularnewline
0.2167 (4.6154) & 17250.736887 \tabularnewline
0.2333 (4.2857) & 34017.464656 \tabularnewline
0.25 (4) & 28301.26499 \tabularnewline
0.2667 (3.75) & 35190.455579 \tabularnewline
0.2833 (3.5294) & 26426.92708 \tabularnewline
0.3 (3.3333) & 8567.095278 \tabularnewline
0.3167 (3.1579) & 879.950829 \tabularnewline
0.3333 (3) & 4563.722788 \tabularnewline
0.35 (2.8571) & 19187.872557 \tabularnewline
0.3667 (2.7273) & 3231.510441 \tabularnewline
0.3833 (2.6087) & 15125.818086 \tabularnewline
0.4 (2.5) & 6779.040633 \tabularnewline
0.4167 (2.4) & 12196.410134 \tabularnewline
0.4333 (2.3077) & 12166.174353 \tabularnewline
0.45 (2.2222) & 13705.236326 \tabularnewline
0.4667 (2.1429) & 10644.788941 \tabularnewline
0.4833 (2.069) & 25736.286139 \tabularnewline
0.5 (2) & 27685.968618 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=112151&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]1[/C][/ROW]
[ROW][C]Frequency (Period)[/C][C]Spectrum[/C][/ROW]
[ROW][C]0.0167 (60)[/C][C]10310840.788791[/C][/ROW]
[ROW][C]0.0333 (30)[/C][C]169176.871747[/C][/ROW]
[ROW][C]0.05 (20)[/C][C]2857425.895077[/C][/ROW]
[ROW][C]0.0667 (15)[/C][C]3699504.834075[/C][/ROW]
[ROW][C]0.0833 (12)[/C][C]1364201.571362[/C][/ROW]
[ROW][C]0.1 (10)[/C][C]217675.047079[/C][/ROW]
[ROW][C]0.1167 (8.5714)[/C][C]113778.371838[/C][/ROW]
[ROW][C]0.1333 (7.5)[/C][C]107934.299205[/C][/ROW]
[ROW][C]0.15 (6.6667)[/C][C]107751.290143[/C][/ROW]
[ROW][C]0.1667 (6)[/C][C]13855.370532[/C][/ROW]
[ROW][C]0.1833 (5.4545)[/C][C]13388.697779[/C][/ROW]
[ROW][C]0.2 (5)[/C][C]44429.936135[/C][/ROW]
[ROW][C]0.2167 (4.6154)[/C][C]17250.736887[/C][/ROW]
[ROW][C]0.2333 (4.2857)[/C][C]34017.464656[/C][/ROW]
[ROW][C]0.25 (4)[/C][C]28301.26499[/C][/ROW]
[ROW][C]0.2667 (3.75)[/C][C]35190.455579[/C][/ROW]
[ROW][C]0.2833 (3.5294)[/C][C]26426.92708[/C][/ROW]
[ROW][C]0.3 (3.3333)[/C][C]8567.095278[/C][/ROW]
[ROW][C]0.3167 (3.1579)[/C][C]879.950829[/C][/ROW]
[ROW][C]0.3333 (3)[/C][C]4563.722788[/C][/ROW]
[ROW][C]0.35 (2.8571)[/C][C]19187.872557[/C][/ROW]
[ROW][C]0.3667 (2.7273)[/C][C]3231.510441[/C][/ROW]
[ROW][C]0.3833 (2.6087)[/C][C]15125.818086[/C][/ROW]
[ROW][C]0.4 (2.5)[/C][C]6779.040633[/C][/ROW]
[ROW][C]0.4167 (2.4)[/C][C]12196.410134[/C][/ROW]
[ROW][C]0.4333 (2.3077)[/C][C]12166.174353[/C][/ROW]
[ROW][C]0.45 (2.2222)[/C][C]13705.236326[/C][/ROW]
[ROW][C]0.4667 (2.1429)[/C][C]10644.788941[/C][/ROW]
[ROW][C]0.4833 (2.069)[/C][C]25736.286139[/C][/ROW]
[ROW][C]0.5 (2)[/C][C]27685.968618[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=112151&T=1

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

As an alternative you can also use a QR Code:  
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)1
Frequency (Period)Spectrum
0.0167 (60)10310840.788791
0.0333 (30)169176.871747
0.05 (20)2857425.895077
0.0667 (15)3699504.834075
0.0833 (12)1364201.571362
0.1 (10)217675.047079
0.1167 (8.5714)113778.371838
0.1333 (7.5)107934.299205
0.15 (6.6667)107751.290143
0.1667 (6)13855.370532
0.1833 (5.4545)13388.697779
0.2 (5)44429.936135
0.2167 (4.6154)17250.736887
0.2333 (4.2857)34017.464656
0.25 (4)28301.26499
0.2667 (3.75)35190.455579
0.2833 (3.5294)26426.92708
0.3 (3.3333)8567.095278
0.3167 (3.1579)879.950829
0.3333 (3)4563.722788
0.35 (2.8571)19187.872557
0.3667 (2.7273)3231.510441
0.3833 (2.6087)15125.818086
0.4 (2.5)6779.040633
0.4167 (2.4)12196.410134
0.4333 (2.3077)12166.174353
0.45 (2.2222)13705.236326
0.4667 (2.1429)10644.788941
0.4833 (2.069)25736.286139
0.5 (2)27685.968618


Figures (Output of Computation)

Input Parameters & R Code

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