Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
1431 869 1536 2216 933 463 3202 372 1584 1440 1837 1549 1431 2274 1556 4120 1766 2251 1061 1955 2080 2294 1220 1290 2557 2148 2381 1639 1247 1897 1768 1592 1786 807 2600 1941 2733 1502 917 2620 2802 1413 1408 1363 2190 4398 918 1831 3558 1792 1438 496 2259 744 1161 2385 2207 2955 2064 1863 3336 2664 2130 2033 1227 1441 1499 2280 1834 1286 1944 893 1763 1403 1455 1930 1900 1546 1466 1431 1465 3176 1317 1761 870 1654 1054 3107 2105 2564 1939 1684 1499 2484 2016 1382 1445 1659 3096 2896 2071 1904 1475 602 1799 1617 2162 2279 874 1282 1446 2049 392 1605 530 1386 2395 387 1743 449 2724 1613 1271 1253 568 1512 2252 1123 2990 2023 2933 1002 956 2306 3744 1107 3351 1588 1474 1760 1374 3494 1695 2035 2468 1477 2211 1962 2 207 5 8 0 0 1678 2544 0 4 151 474 141 873 29 1370
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
1
0
1
2
Degree (D) of seasonal differencing
2
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