Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
4419 4336 4214 4294 4650 4608 4650 4625 4739 5010 4808 4474 4527 4652 4677 4904 4851 4956 4819 4940 5217 5305 5265 5256 5671 5617 5811 5728 5629 5490 5605 4944 5555 5956 5872 5795 6033 6337 6396 6244 6200 6082 5866 5917 6134 6428 6187 6228 6269 6586 6223 6724 6294 6445 6163 6207 6816 6850 6439 6401 6913 6969 7064 6987 6882 6683 6530 6748 6773 7375 7208 6676 7167 7146 7193 7162 7145 6819 6702 6702 6782 7307 6818 6966 7012 7754 7462 7183 7165 7299 7103 6950 7506 7708 7693 7495 7955 8316 9230 8654 8307 7940 7509 7752 8310 8616 8358 8150 8664 8817 8927 8537 8497 8270 7658 8049 8365 8971 8854 8540 8878 9184
Seasonal Period
Default
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