Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
100.30 101.90 102.10 103.20 103.70 106.20 107.70 109.90 111.70 114.90 116.00 118.30 120.40 126.00 128.10 130.10 130.80 133.60 134.20 135.50 136.20 139.10 139.00 139.60 138.70 140.90 141.30 141.80 142.00 144.50 144.60 145.50 146.80 149.50 149.90 150.10 150.90 152.80 153.10 154.00 154.90 156.90 158.40 159.70 160.20 163.20 163.70 164.40 163.70 165.50 165.60 166.80 167.50 170.60 170.90 172.00 171.80 173.90 174.00 173.80 173.90 176.00 176.60 178.20 179.20 181.30 181.80 182.90 183.80 186.30 187.40 189.20 189.70 191.90 192.60 193.70 194.20 197.60 199.30 201.40 203.00 206.30 207.10 209.80 211.10 215.30 217.40 215.50 210.90 212.60
Data Y:
100.00 102.83 109.50 115.91 107.94 110.86 118.89 123.38 113.33 116.38 122.04 125.47 115.62 117.91 122.40 125.05 114.18 114.74 120.63 123.68 112.84 115.64 122.32 124.59 116.33 117.45 125.64 128.38 119.87 121.22 128.98 131.35 121.35 123.72 131.06 134.55 125.93 128.90 136.19 140.34 130.48 134.68 141.05 145.44 136.21 139.85 147.13 151.44 143.62 148.55 153.54 159.79 152.55 155.84 160.38 164.22 156.40 160.05 165.60 171.15 161.90 167.21 171.34 176.83 166.27 172.30 176.71 182.99 172.07 178.17 182.20 188.49 176.88 182.13 185.32 192.86 180.27 184.92 187.82 194.94 184.36 188.80 193.42 199.76 188.78 191.49 194.87 198.28 183.24 204.87
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation