Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
110.43 114.77 132.21 122.86 118.5 130.3 113.25 104.54 132.78 122.99 133.14 125.83 122.99 125.7 148.47 120.75 136.7 139.17 123.47 112.76 137.99 139.75 140.22 121.6 132.33 130.34 149.05 130.47 139.29 146.55 137.79 122.95 139.51 155.77 143.95 125.07 142.35 144.34 145.87 156.01 146.74 156.45 152.29 122.56 154.59 149.68 118.75 109.22 104.19 107.33 114.07 107.92 103.53 117.3 112.09 95.08 123.28 121.98 121.74 119.93 115.11
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