Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
-6.399530 1.957360 2.478650 1.589720 0.184681 0.918505 -1.995250 1.586390 -0.140506 1.356850 1.727480 5.067280 -2.558530 -0.947617 1.134260 0.264652 2.593240 0.743810 0.743969 -1.112320 1.219600 1.306220 0.553112 0.291992 2.656690 -0.855357 1.923590 1.744140 0.735418 1.048160 0.426579 3.418350 -0.513488 0.833442 -1.064710 -0.237221 -0.759509 1.165110 -7.109370 1.726660 1.589090 -1.344360 1.481350 -0.236383 1.388380 1.018800 -1.946000 -2.255320 1.808120 2.125650 -1.159330 7.908360 0.926642 0.590183 2.058870 0.760719 -0.455908 -2.371740 1.112290 2.440940 -2.108970 -1.545510 -1.471520 3.059330 -0.878558 2.184800 -1.854340 -2.484340 -1.072920 -2.315740 2.722060 -0.029346 -3.150900 0.404606 1.887430 -2.627560 -3.574700 0.754203 -2.192110 1.409100 -3.193550 2.848590 -0.385986 1.743400 -2.023100 0.413897 1.122490 0.182360 -0.302073 0.411861 -2.430770 -0.489293 -0.852392 1.413970 0.056108 1.001380 3.783360 3.091870 -0.146159 -0.833078 -1.827790 0.067725 -0.429261 -0.310928 1.809600 0.624170 -0.700127 -0.723910 1.868270 -3.616350 -1.970360 1.675560 1.876060 3.632780 1.012410 2.040960 3.995540 2.484770 1.588610 2.972830 -0.499202 -3.911090 -2.829110 -6.120530 0.339289 -0.944523 0.035212 -1.042880 -1.323190 -1.133120 0.183183 1.967530 -2.511840 -0.287708 -0.369008 0.685670 0.758509 -1.400430 -3.497010 -0.776372 0.760543 2.310800 -2.673360 0.700963 -4.503390 -5.186020 -1.880700 -7.006000 -0.421422 -2.080020 -3.165650 0.737049 1.532660 0.811355 -0.429889 0.602124 1.069580 -1.168290 -0.541346 -0.852089 1.349240 1.763480 0.356035 -0.821847 0.901744 -1.230250
Seasonal Period
12
12
1
2
3
4
6
12
Chart options
R Code
par1 <- as.numeric(par1) 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) 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