Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
3281 3397 3498.5 3538 3449.5 3673 3350.5 3604 3673.5 3747 3616 3580.5 3710 3994.5 4091 3954.5 4004 4287 3831 4046.5 4079.5 4029.5 3880 3855 3841.5 4123.5 4133 3958.5 4003 4151.5 3723 3957 3965.5 3861.5 3917.5 3704 3950 4140.5 4090 4162 4066 4358.5 4022.5 4285.5 4373.5 4284.5 4077.5 4122 4181.5 4535.5 4497 4420.5 4370 4712 4475 4578.5 4751.5 4746 4581.5 4645.5 4751 4952.5 4996.5 4998 4986.5 5348 4933 5263 5330.5 5301 5159 5258.5 5411.5 5536.5 5613 5505.5 5476 5782.5 5283 5451.5 5578 5548.5 5379.5 5117.5 5316.5 5505.5 5620.5 5383.5 5461.5 5658.5 5357.5 5622 5608 5604.5 5399 5185 5221 5379.5 5333 5214 5206.5 5630 5285.5 5512.5 5592.5 5554.5 5284.5 5198.5 5241.5 5455 5548.5 5375 5346 5730.5 5457 5603
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
par4 <- '12' par3 <- '1' par2 <- '1' par1 <- '1' 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