Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
16198.9 16896.2 0 16554.2 16698.00 0 19554.2 19691.6 0 15903.8 15930.7 0 18003.8 17444.6 0 18329.6 17699.4 0 16260.7 15189.8 0 14851.9 15672.7 0 18174.1 17180.8 0 18406.6 17664.9 0 18466.5 17862.9 0 16016.5 16162.3 0 17428.5 17463.6 0 17167.2 16772.1 0 19630.00 19106.9 0 17183.6 16721.3 0 18344.7 18161.3 0 19301.4 18509.9 0 18147.5 17802.7 0 16192.9 16409.9 0 18374.4 17967.7 0 20515.2 20286.6 0 18957.2 19537.3 0 16471.5 18021.9 0 18746.8 20194.3 0 19009.5 19049.6 0 19211.2 20244.7 0 20547.7 21473.3 0 19325.8 19673.6 0 20605.5 21053.2 0 20056.9 20159.5 0 16141.4 18203.6 0 20359.8 21289.5 0 19711.6 20432.3 1 15638.6 17180.4 1 14384.5 15816.8 1 13855.6 15071.8 1 14308.3 14521.1 1 15290.6 15668.8 1 14423.8 14346.9 1 13779.7 13881.00 1 15686.3 15465.9 1 14733.8 14238.2 1 12522.5 13557.7 1 16189.4 16127.6 1 16059.1 16793.9 1 16007.1 16014.00 1 15806.8 16867.9 1 15160.00 16014.6 0 15692.1 15878.6 0 18908.9 18664.9 0 16969.9 17962.5 0 16997.5 17332.7 0 19858.9 19542.1 0 17681.2 17203.6 0
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
1
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation