Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
6396 6509 7018 6791 6300 6737 6399 6758 6601 6891 6789 5967 6995 6771 7360 7162 7025 7188 6801 7395 7203 7534 7267 6741 7185 7295 7660 7344 7439 7264 7192 7368 7246 7625 7153 6875 7413 6981 7651 7148 7072 6943 7022 6925 6904 7666 7112 6907 7880 7261 7478 7742 7499 7563 7812 7582 7362 7779 7681 7795 8049 7541 7869 7855 7582 7723 8195 7729 7856 7984 7401 8162 8223 7743 8249 8137 7650 7689 7901 7951 7956 7908 7887 8198 8358 7837 8169 8040 8112 8348 8548 8498 8167 8331 8070 7936 8553 8108 8341 8260 8113 8042 8467 8333 8243 8517 8120 8012 8801 8115 8161 8271 8108 7992 8239 8419 8781 8992 7966 8626 7935 8166 8411 8693 8861 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
Seasonal Period
0.9TRUEFALSEFALSEFALSE0.9FALSETRUETRUE1212120.912412
12
1
2
3
4
6
12
Chart options
R Code
par1 <- '4' 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