Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
81.71 87.703 90.09 100.639 83.042 89.956 89.561 105.38 86.554 93.131 92.812 102.195 88.925 94.184 94.196 108.932 91.134 97.149 96.415 112.432 92.47 98.61410515 97.80117197 111.8560178 95.63981455 104.1120262 104.0148224 118.1743476 102.033431 109.3138852 108.1523649 121.30381 103.8725146 112.7185207 109.0381253 122.4434864 106.6325686 113.8153852 111.1071252 130.039536 109.6121057 116.8592117 113.8982545 128.9375926 111.8120023 119.9689463 117.018539 132.4743387 116.3369106 124.6405636 121.025249 137.2054829 120.0187687 127.0443429 124.349043 143.6114438
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
4
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