Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
4880 5170 4430 4430 3940 4160 4530 5580 7200 6550 5870 5570 5180 4630 4330 4180 4030 4630 4550 5270 6990 6510 5550 5040 5350 4540 3950 3930 3710 3840 4050 4940 5830 5250 5400 4870 4320 4240 4320 4200 3630 3740 3960 5110 4230 4790 5210 4660 4800 3610 4170 3500 3490 4390 4230 3730 4140 4660 5160 3890 4400 3610 3410 3640 3520 3680 4120 4420 4320 5360 4800 4070 4000 3840 3580 3400 3340 3790 4080 3710 4080 4330 3850 3700 3620 3710 3290 3470 3230 3840 3910 4190 4140 4490 3990 3690 3940 3680 3460 3240 3280 3410 3560 3550 3870 4450 4260 4050
Seasonal Period
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