Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
571000 584000 599000 582000 530000 528000 536000 546000 559000 562000 541000 539000 548000 563000 581000 572000 519000 521000 531000 540000 548000 556000 551000 549000 564000 586000 604000 601000 545000 537000 552000 563000 575000 580000 575000 558000 564000 581000 597000 587000 536000 524000 537000 536000 533000 528000 516000 502000 506000 518000 534000 528000 478000 469000 490000 493000 508000 517000 514000 510000 527000 542000 565000 555000 499000 511000 526000 532000 549000 561000 557000 566000 588000 620000 626000 620000 573000 573000 574000 580000 590000 593000 597000 595000 612000 628000 629000 621000 569000 567000 573000 584000 589000 591000 595000 594000 611000 613000 611000 594000 543000 537000 544000 555000 561000 562000 555000 547000 565000 578000 580000 569000 507000 501000 509000 510000 517000 519000 512000 509000 519000 523000 525000 517000 456000 455000 461000 470000 475000 476000 471000 471000 503000 513000 510000 484000 431000 436000 443000 448000 460000 467000 460000 464000 485000 501000 521000 488000 439000 442000 457000 462000 481000 493000
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,var(myx)) a<-table.element(a,'Range',header=TRUE) a<-table.element(a,max(myx)-min(myx)) a<-table.element(a,'Trim Var.',header=TRUE) smyx <- sort(myx) sn <- length(smyx) a<-table.element(a,var(smyx[smyx>quantile(smyx,0.05) & smyx<quantile(smyx,0.95)])) 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