Home » date » 2010 » Dec » 04 »

Cumulatieve periodigram: uitvoer

*The author of this computation has been verified*

R Software Module: /rwasp_spectrum.wasp (opens new window with default values)

Title produced by software: Spectral Analysis

Date of computation: Sat, 04 Dec 2010 08:52:41 +0000

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Dec/04/t1291452659ns0arl7jhurf905.htm/, Retrieved Sun, 19 May 2013 14:50:46 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

System-generated keywords (parent):

t1291022416gbwwep1jkc6r0zm (pk = 102761)

Estimated Impact

Dataseries X:

» Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability «

Output produced by software:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	2 seconds
R Server	'Sir Ronald Aylmer Fisher' @ 193.190.124.24

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)	16815217.946585
0.0333 (30)	34336342.256614
0.05 (20)	9567080.01611
0.0667 (15)	4061096.912144
0.0833 (12)	1076476.051033
0.1 (10)	1572034.643123
0.1167 (8.5714)	36570.782799
0.1333 (7.5)	320236.060064
0.15 (6.6667)	152265.821227
0.1667 (6)	8231247.258462
0.1833 (5.4545)	1094046.844298
0.2 (5)	176435.490157
0.2167 (4.6154)	2082268.290889
0.2333 (4.2857)	1621912.832365
0.25 (4)	13144711.869345
0.2667 (3.75)	1354708.266804
0.2833 (3.5294)	78867.210145
0.3 (3.3333)	310836.444992
0.3167 (3.1579)	455104.465438
0.3333 (3)	6592283.310466
0.35 (2.8571)	3778623.043561
0.3667 (2.7273)	942616.859875
0.3833 (2.6087)	313114.929052
0.4 (2.5)	1689079.330869
0.4167 (2.4)	64649.06712
0.4333 (2.3077)	421607.721836
0.45 (2.2222)	368902.871189
0.4667 (2.1429)	466694.143279
0.4833 (2.069)	449967.45734
0.5 (2)	4257788.495537

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/Dec/04/t1291452659ns0arl7jhurf905/1psqq1291452758.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Dec/04/t1291452659ns0arl7jhurf905/1psqq1291452758.ps (opens in new window)

Click here to open pdf file.

http://www.freestatistics.org/blog/date/2010/Dec/04/t1291452659ns0arl7jhurf905/2ijqt1291452758.png (opens in new window)

http://www.freestatistics.org/blog/date/2010/Dec/04/t1291452659ns0arl7jhurf905/2ijqt1291452758.ps (opens in new window)

Click here to open pdf file.

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