Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
12.28411643 12.27179952 12.25231194 12.22609584 12.20725694 12.22423668 12.39452537 12.46961698 12.40214917 12.40680698 12.3612585 12.37673061 12.3683743 12.34737254 12.33351554 12.30898116 12.2883068 12.29888881 12.4546734 12.49066871 12.42369699 12.39302513 12.34254658 12.35365178 12.29809105 12.29825521 12.26721478 12.27732762 12.26028617 12.26177695 12.38859418 12.3910407 12.34835315 12.24910819 12.19562829 12.17868089 12.17860381 12.13181235 12.09013416 12.05877841 12.02755584 12.03203468 12.21750702 12.21794193 12.11459312 12.06379454 12.0224142 12.05115417 12.04179334 12.00873085 11.99482061 11.98050115 11.92504174 11.97107157 12.13634736 12.13923281 12.06870533 12.03979371 12.02471263 12.08974667 12.13571456 12.15069735 12.16297202 12.14900482 12.13352883 12.1823382 12.31730986 12.33226026 12.27431991
Data Y:
12.46443275 12.46321948 12.45411952 12.43120622 12.42071214 12.42742259 12.47848625 12.4852789 12.51740655 12.51517736 12.49711407 12.49107112 12.48346572 12.48659792 12.48137592 12.46452542 12.45813013 12.45868567 12.51166585 12.52266282 12.54708456 12.51603222 12.50508034 12.49958892 12.52976355 12.51446502 12.50066578 12.48145185 12.48226016 12.49089438 12.5485305 12.56867548 12.57564456 12.54697079 12.52043803 12.50254844 12.51587084 12.50723964 12.47015466 12.4770459 12.44037768 12.38280779 12.4624111 12.49395095 12.48096959 12.4609629 12.42451717 12.42232853 12.44262983 12.43657307 12.38628364 12.38718097 12.34057294 12.34374981 12.4292122 12.44515254 12.43391713 12.40055331 12.39453365 12.40486203 12.41530387 12.41845045 12.41956578 12.42514793 12.393904 12.40068091 12.48097718 12.49928316 12.48350739
Sample Range:
(leave blank to include all observations)
From:
To:
bandwidth of density plot
(?)
# lags (autocorrelation function)
(?)
36
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
Chart options
Label y-axis:
Label x-axis:
R Code
par1 <- as.numeric(par1) par2 <- as.numeric(par2) x <- as.ts(x) y <- as.ts(y) mylm <- lm(y~x) cbind(mylm$resid) library(lattice) bitmap(file='pic1.png') plot(y,type='l',main='Run Sequence Plot of Y[t]',xlab='time or index',ylab='value') grid() dev.off() bitmap(file='pic1a.png') plot(x,type='l',main='Run Sequence Plot of X[t]',xlab='time or index',ylab='value') grid() dev.off() bitmap(file='pic1b.png') plot(x,y,main='Scatter Plot',xlab='X[t]',ylab='Y[t]') grid() dev.off() bitmap(file='pic1c.png') plot(mylm$resid,type='l',main='Run Sequence Plot of e[t]',xlab='time or index',ylab='value') grid() dev.off() bitmap(file='pic2.png') hist(mylm$resid,main='Histogram of e[t]') dev.off() bitmap(file='pic3.png') if (par1 > 0) { densityplot(~mylm$resid,col='black',main=paste('Density Plot of e[t] bw = ',par1),bw=par1) } else { densityplot(~mylm$resid,col='black',main='Density Plot of e[t]') } dev.off() bitmap(file='pic4.png') qqnorm(mylm$resid,main='QQ plot of e[t]') qqline(mylm$resid) grid() dev.off() if (par2 > 0) { bitmap(file='pic5.png') acf(mylm$resid,lag.max=par2,main='Residual Autocorrelation Function') grid() dev.off() } summary(x) load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Model: Y[t] = c + b X[t] + e[t]',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'c',1,TRUE) a<-table.element(a,mylm$coeff[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'b',1,TRUE) a<-table.element(a,mylm$coeff[[2]]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Descriptive Statistics about e[t]',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'# observations',header=TRUE) a<-table.element(a,length(mylm$resid)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'minimum',header=TRUE) a<-table.element(a,min(mylm$resid)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Q1',header=TRUE) a<-table.element(a,quantile(mylm$resid,0.25)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'median',header=TRUE) a<-table.element(a,median(mylm$resid)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mean',header=TRUE) a<-table.element(a,mean(mylm$resid)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Q3',header=TRUE) a<-table.element(a,quantile(mylm$resid,0.75)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'maximum',header=TRUE) a<-table.element(a,max(mylm$resid)) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation