Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_spectrum.wasp

Title produced by software

Spectral Analysis

Date of computation

Fri, 14 Dec 2007 03:29:11 -0700

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2007/Dec/14/t11976272867ryv6nfwi6nm215.htm/, Retrieved Thu, 02 May 2024 16:26:09 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=3809, Retrieved Thu, 02 May 2024 16:26:09 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

s0650550 s0650062

Estimated Impact

198

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-     [Standard Deviation-Mean Plot] [Inducing Stationa...] [2007-11-29 14:28:06] [68a1fecd8f1c75119cd425050381cede]
- RMPD    [Spectral Analysis] [Time series] [2007-12-14 10:29:11] [ab924f39c1cc7a5dd22761038b10db61] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

263418000000
262752000000
266433000000
267722000000
266003000000
262971000000
265521000000
264676000000
270223000000
269508000000
268457000000
265814000000
266680000000
263018000000
269285000000
269829000000
270911000000
266844000000
271244000000
269907000000
271296000000
270157000000
271322000000
267179000000
264101000000
265518000000
269419000000
268714000000
272482000000
268351000000
268175000000
270674000000
272764000000
272599000000
270333000000
270846000000
270491000000
269160000000
274027000000
273784000000
276663000000
274525000000
271344000000
271115000000
270798000000
273911000000
273985000000
271917000000
273338000000
270601000000
273547000000
275363000000
281229000000
277793000000
279913000000
282500000000
280041000000
282166000000
290304000000
283519000000

Summary of compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

\begin{tabular}{lllllllll}
\hline
Summary of compuational 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=3809&T=0

[TABLE]
[ROW][C]Summary of compuational 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=3809&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=3809&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 compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Raw Periodogram
Parameter	Value
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.0167 (60)	45401539764285145088
0.0333 (30)	1.025483352546e+18
0.05 (20)	27019663968139280384
0.0667 (15)	4072972301344728576
0.0833 (12)	27034425452962336768
0.1 (10)	4731835403865708544
0.1167 (8.5714)	2275670071695108864
0.1333 (7.5)	8422946963446760448
0.15 (6.6667)	10511475231171446784
0.1667 (6)	40039374976593371136
0.1833 (5.4545)	1236050109596857088
0.2 (5)	332739673048928064
0.2167 (4.6154)	2277307797395766016
0.2333 (4.2857)	801024781280353664
0.25 (4)	2845937906899729408
0.2667 (3.75)	1675016824189553920
0.2833 (3.5294)	478513431506132352
0.3 (3.3333)	8502532571946564608
0.3167 (3.1579)	282751513973295392
0.3333 (3)	530021280583968896
0.35 (2.8571)	3464112158545502208
0.3667 (2.7273)	1652375900394635776
0.3833 (2.6087)	106758531937505904
0.4 (2.5)	3330525795597602304
0.4167 (2.4)	4775403387517319168
0.4333 (2.3077)	885155808568176384
0.45 (2.2222)	1090112532756221184
0.4667 (2.1429)	1585208665804837376
0.4833 (2.069)	2141837200608446976
0.5 (2)	21725449193290301440

\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.0167 (60) & 45401539764285145088 \tabularnewline
0.0333 (30) & 1.025483352546e+18 \tabularnewline
0.05 (20) & 27019663968139280384 \tabularnewline
0.0667 (15) & 4072972301344728576 \tabularnewline
0.0833 (12) & 27034425452962336768 \tabularnewline
0.1 (10) & 4731835403865708544 \tabularnewline
0.1167 (8.5714) & 2275670071695108864 \tabularnewline
0.1333 (7.5) & 8422946963446760448 \tabularnewline
0.15 (6.6667) & 10511475231171446784 \tabularnewline
0.1667 (6) & 40039374976593371136 \tabularnewline
0.1833 (5.4545) & 1236050109596857088 \tabularnewline
0.2 (5) & 332739673048928064 \tabularnewline
0.2167 (4.6154) & 2277307797395766016 \tabularnewline
0.2333 (4.2857) & 801024781280353664 \tabularnewline
0.25 (4) & 2845937906899729408 \tabularnewline
0.2667 (3.75) & 1675016824189553920 \tabularnewline
0.2833 (3.5294) & 478513431506132352 \tabularnewline
0.3 (3.3333) & 8502532571946564608 \tabularnewline
0.3167 (3.1579) & 282751513973295392 \tabularnewline
0.3333 (3) & 530021280583968896 \tabularnewline
0.35 (2.8571) & 3464112158545502208 \tabularnewline
0.3667 (2.7273) & 1652375900394635776 \tabularnewline
0.3833 (2.6087) & 106758531937505904 \tabularnewline
0.4 (2.5) & 3330525795597602304 \tabularnewline
0.4167 (2.4) & 4775403387517319168 \tabularnewline
0.4333 (2.3077) & 885155808568176384 \tabularnewline
0.45 (2.2222) & 1090112532756221184 \tabularnewline
0.4667 (2.1429) & 1585208665804837376 \tabularnewline
0.4833 (2.069) & 2141837200608446976 \tabularnewline
0.5 (2) & 21725449193290301440 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=3809&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.0167 (60)[/C][C]45401539764285145088[/C][/ROW]
[ROW][C]0.0333 (30)[/C][C]1.025483352546e+18[/C][/ROW]
[ROW][C]0.05 (20)[/C][C]27019663968139280384[/C][/ROW]
[ROW][C]0.0667 (15)[/C][C]4072972301344728576[/C][/ROW]
[ROW][C]0.0833 (12)[/C][C]27034425452962336768[/C][/ROW]
[ROW][C]0.1 (10)[/C][C]4731835403865708544[/C][/ROW]
[ROW][C]0.1167 (8.5714)[/C][C]2275670071695108864[/C][/ROW]
[ROW][C]0.1333 (7.5)[/C][C]8422946963446760448[/C][/ROW]
[ROW][C]0.15 (6.6667)[/C][C]10511475231171446784[/C][/ROW]
[ROW][C]0.1667 (6)[/C][C]40039374976593371136[/C][/ROW]
[ROW][C]0.1833 (5.4545)[/C][C]1236050109596857088[/C][/ROW]
[ROW][C]0.2 (5)[/C][C]332739673048928064[/C][/ROW]
[ROW][C]0.2167 (4.6154)[/C][C]2277307797395766016[/C][/ROW]
[ROW][C]0.2333 (4.2857)[/C][C]801024781280353664[/C][/ROW]
[ROW][C]0.25 (4)[/C][C]2845937906899729408[/C][/ROW]
[ROW][C]0.2667 (3.75)[/C][C]1675016824189553920[/C][/ROW]
[ROW][C]0.2833 (3.5294)[/C][C]478513431506132352[/C][/ROW]
[ROW][C]0.3 (3.3333)[/C][C]8502532571946564608[/C][/ROW]
[ROW][C]0.3167 (3.1579)[/C][C]282751513973295392[/C][/ROW]
[ROW][C]0.3333 (3)[/C][C]530021280583968896[/C][/ROW]
[ROW][C]0.35 (2.8571)[/C][C]3464112158545502208[/C][/ROW]
[ROW][C]0.3667 (2.7273)[/C][C]1652375900394635776[/C][/ROW]
[ROW][C]0.3833 (2.6087)[/C][C]106758531937505904[/C][/ROW]
[ROW][C]0.4 (2.5)[/C][C]3330525795597602304[/C][/ROW]
[ROW][C]0.4167 (2.4)[/C][C]4775403387517319168[/C][/ROW]
[ROW][C]0.4333 (2.3077)[/C][C]885155808568176384[/C][/ROW]
[ROW][C]0.45 (2.2222)[/C][C]1090112532756221184[/C][/ROW]
[ROW][C]0.4667 (2.1429)[/C][C]1585208665804837376[/C][/ROW]
[ROW][C]0.4833 (2.069)[/C][C]2141837200608446976[/C][/ROW]
[ROW][C]0.5 (2)[/C][C]21725449193290301440[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=3809&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=3809&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
Parameter	Value
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.0167 (60)	45401539764285145088
0.0333 (30)	1.025483352546e+18
0.05 (20)	27019663968139280384
0.0667 (15)	4072972301344728576
0.0833 (12)	27034425452962336768
0.1 (10)	4731835403865708544
0.1167 (8.5714)	2275670071695108864
0.1333 (7.5)	8422946963446760448
0.15 (6.6667)	10511475231171446784
0.1667 (6)	40039374976593371136
0.1833 (5.4545)	1236050109596857088
0.2 (5)	332739673048928064
0.2167 (4.6154)	2277307797395766016
0.2333 (4.2857)	801024781280353664
0.25 (4)	2845937906899729408
0.2667 (3.75)	1675016824189553920
0.2833 (3.5294)	478513431506132352
0.3 (3.3333)	8502532571946564608
0.3167 (3.1579)	282751513973295392
0.3333 (3)	530021280583968896
0.35 (2.8571)	3464112158545502208
0.3667 (2.7273)	1652375900394635776
0.3833 (2.6087)	106758531937505904
0.4 (2.5)	3330525795597602304
0.4167 (2.4)	4775403387517319168
0.4333 (2.3077)	885155808568176384
0.45 (2.2222)	1090112532756221184
0.4667 (2.1429)	1585208665804837376
0.4833 (2.069)	2141837200608446976
0.5 (2)	21725449193290301440

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

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')

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code