Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
4870 4240 3800 3990 3290 4710 4210 4440 5040 5070 4900 4790 3890 3450 4080 3280 3130 3310 3860 4570 5110 4820 4250 4210 3610 3710 2760 2710 2710 3290 2670 3620 4440 3910 4610 3760 3460 3020 3360 2610 2670 2480 2610 3320 2800 3030 3740 3060 3040 2620 3190 2750 2630 3290 2430 2730 3690 2980 2590 3360 2370 2200 2330 2370 2200 2430 2400 2840 2870 3320 3090 2680 2420 2550 2420 2430 2330 2520 2630 2570 2800 2680 2430 2790 2420 2750 2350 2330 2290 2330 2490 2480 2760 2590 2950 2570 2960 2540 2400 2470 2390 2310 2470 2490 2510 2690 3060 2690 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
Seasonal Period
4
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