Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
7.4271 7.7662 7.6289 7.5281 7.3831 7.2355 7.0617 7.1237 7.4533 7.5411 7.4978 7.3525 7.3862 7.311 7.2013 7.249 7.3321 7.59 7.9082 8.2123 8.0929 8.118 8.1206 8.2883 8.4281 8.7917 8.9168 8.9446 8.9786 9.5862 9.6533 9.4125 9.2195 9.2882 9.6774 9.6857 10.1688 10.4399 10.4675 10.149 9.9163 9.9268 10.0529 10.1622 10.083 10.1134 10.3423 10.7536 11.0967 10.8588 10.7719 10.9262 10.708 10.5062 10.0683 9.8954 9.9589 9.9177 9.7189 9.5273 9.5746 9.763 9.6117 9.6581 9.8361 10.2353 10.1285 10.1347 10.2141 10.0971 9.9651 10.1286 10.3356 10.1238 10.1326 10.2467 10.44 10.3689 10.2415 10.3899 10.3162 10.4533 10.6741 10.8957 10.7404 10.6568 10.5682 10.9833 11.0237 10.8462 10.7287 10.7809 10.2609 9.8252 9.1071 8.695 9.2205 9.0496 8.7406 8.921 9.011 9.3157 9.5786 9.6246 9.7485 9.9431 10.1152 10.1827 9.9777 9.7436 9.3462 9.2623 9.1505 8.5794 8.3245 8.6538 8.752 8.8104 9.2665 9.0895
Box-Cox transformation parameter
12
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
1
2
Degree (D) of seasonal differencing
0
1
2
Seasonal Period
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