Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
3280 3444 3855 3811 3785 4075 3547 3863 4064 4176 4191 4307 4179 4622 4798 4673 4635 4875 4097 4262 4135 4238 3891 3573 3963 4192 4306 4316 4249 4408 3731 4096 4102 3962 3845 3734 3933 4176 4150 4137 4016 4113 3611 3474 3654 3712 3394 3348 3476 3908 4009 4102 4253 4532 4080 4402 4597 4844 4877 4735 4768 5251 5553 5548 5519 5798 4918 5271 5492 5547 5244 5149 5453 5584 5773 5811 5687 5647 4892 5235 5311 5378 4994 4559 4895 5104 5477 5302 5360 5540 4877 5241 5233 5561 5049 4482 4846 4636 4431 4702 4775 4834 4344 4800 4981 5069 4655 4254 4753 4888 5048 4991 4962 5150 4444 4815
Seasonal Period
12
12
1
2
3
4
6
12
Chart options
R Code
par1 <- as.numeric(par1) n.orig <- length(x) x <- na.omit(x) 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) if(n.orig!=n) { a<-table.row.start(a) a<-table.element(a,'Warning: NAs were removed from the time series! The results shown below will only be correct if the NAs are all located at the start and/or end of the time series.',6,F) 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,signif(var(myx), digits=6)) a<-table.element(a,'Range',header=TRUE) a<-table.element(a,signif(max(myx)-min(myx), digits=6)) a<-table.element(a,'Trim Var.',header=TRUE) smyx <- sort(myx) sn <- length(smyx) a<-table.element(a,signif(var(smyx[smyx>quantile(smyx,0.05) & smyx<quantile(smyx,0.95)]), digits=6)) 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