Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
4.120575483 8.926511718 -4.46024537 4.727882162 4.432447952 6.950256655 -5.32280306 -14.69266743 9.590894311 23.81920504 2.589069408 -9.504994357 -4.66481116 -0.229193754 4.008248556 5.383133174 -0.630564198 4.148431753 -9.51138505 -26.78764011 7.089985394 10.64478195 -8.66297919 -18.5246209 -13.93557738 -14.03466139 7.721960062 0.405521803 -2.811832441 15.73565743 -12.38854196 -18.25109965 12.66670905 11.26625459 8.322876048 1.079497504 -5.869817274 -6.599952905 10.76077988 8.019226242 2.303242442 11.71054912 -15.01182537 -19.98990534 -5.734654326 1.602333524 10.6078153 -13.65884666 0.998690786 -4.683954969 37.83385373 8.819708971 8.615611775 48.25943771 -31.23279409 -22.29124044 -0.374802185 5.6841057 7.232050034 -26.2158943
Data Y:
9.566151798 10.1585264 8.675824146 8.49107495 6.150900994 3.328024789 2.524800673 4.548854048 2.704278807 3.633603244 12.10954987 12.70192447 8.815998101 1.57024569 -1.180777783 0.717175271 -4.162302863 -3.966704149 -2.052630515 -4.583132122 -4.220081962 -9.256162024 -4.04326556 -9.523920864 -8.538301891 -4.367933722 -2.009284847 -7.653860088 -0.568803499 -7.929806713 -6.76148549 -7.212508963 -1.135385169 -12.89756555 -10.03891668 -11.1802678 -8.313993526 -9.638046901 -12.08175237 -11.12837455 -7.61076941 -2.385846259 5.577203902 4.151973357 6.778073679 10.66486982 4.795064035 0.592402302 7.721419347 1.836362757 17.51348655 17.11783543 -0.587795425 15.36972868 6.50283962 -1.050538193 -1.805962952 1.055602646 -5.343220184 -12.14439735
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