Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
336.02 333.15 314.95 302.48 307.31 305.50 308.57 322.58 337.09 323.81 333.06 331.90 327.90 319.93 331.51 336.42 319.77 323.20 324.51 328.34 331.88 336.45 337.95 330.75 323.87 325.26 328.73 331.72 332.54 354.25 352.69 356.15 372.50 390.90 404.65 430.04 453.54 464.98 463.31 497.20 528.62 470.91 499.53 493.51 469.97 464.41 487.15 476.45 484.91 509.61 495.19 504.75 493.43 488.58 484.82 488.46 512.32 530.29 549.38 551.45 604.41 625.29 623.56 577.42 572.28 571.69 596.28 560.00 577.93 606.51 597.31 607.58 648.14 737.48 708.73 674.01 679.90 674.93 663.38 665.69 684.21 703.71 755.42 772.43
Seasonal Period
12
12
1
2
3
4
6
12
Chart options
R Code
par1 <- as.numeric(par1) n <- length(x) sx <- sort(x) load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Variance Reduction Matrix',6,TRUE) a<-table.row.end(a) for (bigd in 0:2) { for (smalld in 0:3) { mylabel <- 'V(Y[t],d=' mylabel <- paste(mylabel,as.character(smalld),sep='') mylabel <- paste(mylabel,',D=',sep='') mylabel <- paste(mylabel,as.character(bigd),sep='') mylabel <- paste(mylabel,')',sep='') a<-table.row.start(a) a<-table.element(a,mylabel,header=TRUE) myx <- x if (smalld > 0) myx <- diff(myx,lag=1,differences=smalld) if (bigd > 0) myx <- diff(myx,lag=par1,differences=bigd) a<-table.element(a,var(myx)) a<-table.element(a,'Range',header=TRUE) a<-table.element(a,max(myx)-min(myx)) a<-table.element(a,'Trim Var.',header=TRUE) smyx <- sort(myx) sn <- length(smyx) a<-table.element(a,var(smyx[smyx>quantile(smyx,0.05) & smyx<quantile(smyx,0.95)])) a<-table.row.end(a) } } a<-table.end(a) table.save(a,file='mytable.tab') bitmap(file='pic0.png') op <- par(mfrow=c(2,2)) plot(x,type='l',xlab='time',ylab='value',main='d=0, D=0') plot(diff(x,lag=1,differences=1),type='l',xlab='time',ylab='value',main='d=1, D=0') plot(diff(x,lag=par1,differences=1),type='l',xlab='time',ylab='value',main='d=0, D=1') plot(diff(diff(x,lag=1,differences=1),lag=par1,differences=1),type='l',xlab='time',ylab='value',main='d=1, D=1') par(op) dev.off()
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