Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
98101950.13 187968318.5 188999043.8 210743289 230609073.9 130468610.9 190983186.8 162537001 262991089 175244644 153536881 218451087.3 116973495.2 218177896.1 139972561 142356218 112613295.5 253547719.2 123096251.3 262749137.1 102010000 147614906.3 135590955.4 85561526.32 80654369.49 104950319.4 155190863.4 177088669.1 117497470.1 139892713.1 119376190.6 113962030.1 86439787.29 108851881.9 150342218.5 119054770.4 87131416.39 135840396.2 122766400 96828411.45 55486359.51 69933078.76 71667060.61 67583576.24 108844643.7 72887198.76 72854150.57 63959435 39707815.79 57692632.99 51846967.98 37853044.1 36775240.36 52851578.19 43275872.83 59417285.63 40974770.56 49590376.35 68830403.8 92416177.78
Data Y:
10284.5 12792 12823.61538 13845.66667 15335.63636 11188.5 13633.25 12298.46667 15353.63636 12696.15385 12213.93333 13683.72727 11214.14286 13950.23077 11179.13333 11801.875 11188.82353 16456.27273 11110.0625 16530.69231 10038.41176 11681.25 11148.88235 8631 9386.444444 9764.736842 12043.75 12948.06667 10987.125 11648.3125 10633.35294 10219.3 9037.6 10296.31579 11705.41176 10681.94444 9362.947368 11306.35294 10984.45 10062.61905 8118.583333 8867.48 8346.72 8529.307692 10697.18182 8591.84 8695.607143 8125.571429 7009.758621 7883.466667 7527.645161 6763.758621 6682.333333 7855.681818 6738.88 7895.434783 6361.884615 6935.956522 8344.454545 9107.944444
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