Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
19.25 11.6 15.15 10.95 15.2 12.6 13.2 9.95 19.9 8.1 12.9 14.85 14.05 10.95 7.65 12.65 11.35 14.5 13.6 14.9 16.1 12.4 18.1 18.25 12.15 17.35 12.6 7.6 13.4 14.1 19.9 18.1 11.85 16.65 15.6 15.25 16.1 15.4 13.35 15.4 16.1 16.2 7.7 11.15 13.15 14.75 15.85 15.4 14.1 18.2 16.15 11.2 18.4 17.65 18.45 9.9 16.6 17.6 17.65 18.4 12.6 19.3 11.2 14.6 18.45 4.5 19.1 13.4 4.35 12.75 15.6 11.85 10.95 15.25 11.9 18.55 11.95 15.1 15.6 15.1 17.85 19.05 16.65 12.4 12.6 13.35 16.1 18.25 12.35 14.85 13.85 14.6 7.85 16 13.9 18.95 11.4 14.6 15.25 12.45 19.1 14.6 12.7 13.2 17.75 16.35 18.4 12.85 15.35 17.75 13.1 15.7 15.95 14.7 15.65 13.35 14.75 14.6 15.9 19.1 14.9 12.2 7.85 12.35 19.2 8.6 11.75 9.85 16.85 10.35 14.9 10.6 15.35 9.6 11.9 14.75 14.8 16.35 16.85 15.2 17.35 18.15 13.6 13.6 15 16.85 17.1 17.1 13.35 17.75 18.9 13.6 13.95 15.65 14.35 14.75 11.7 14.35 19.1 16.6 9.5 16.25 17.6 17.1 16.1 17.75 13.6 15.6 12.65 13.6 11.7
Box-Cox transformation parameter
1
1
-2.0
-1.9
-1.8
-1.7
-1.6
-1.5
-1.4
-1.3
-1.2
-1.1
-1.0
-0.9
-0.8
-0.7
-0.6
-0.5
-0.4
-0.3
-0.2
-0.1
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2.0
Degree (d) of non-seasonal differencing
0
0
1
2
Degree (D) of seasonal differencing
0
0
1
2
Seasonal Period
12
1
2
3
4
12
Chart options
R Code
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')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation