Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
2.366550914 -2.784331965 11.17687724 -13.08690014 -3.702416453 13.53551828 9.502773534 2.81572168 0.60106142 -7.213545194 -6.130773624 2.389932679 0.66061216 -3.473042289 12.52091167 -20.54891176 -6.198884935 11.12267743 3.776128477 6.474416363 -0.522105764 -4.597921587 -3.270313182 11.43658892 -4.449553233 3.353978285 8.860826741 -21.06844013 -5.968440133 8.548788298 9.837605911 6.268638544 -4.941687785 -0.980478576 -1.983046747 10.82647717 -2.982939457 0.162753436 8.784262151 -16.11744996 -3.633875981 4.70491481 14.16269979 5.551463759 -7.174432532 1.213636315 2.506787859 -1.385400337 5.889613074 -2.333661401 4.590469131 -15.42590324 -10.26046739 8.260131619 10.20132965 -4.500489758 -8.71605972 -11.09283889 -8.636819669 1.076984534
Data Y:
-0.3160103 -6.358448377 12.27737808 -13.46500204 -12.99867146 9.566650862 6.96160958 1.918725333 -1.176092864 -2.333737351 -2.125109214 12.64884125 3.499204854 -1.09186136 14.07900638 -17.46250601 -12.86891915 9.203882529 -0.057125727 3.885171829 -3.922249573 -3.325720513 -1.493786564 17.69419692 -4.072662393 5.117089915 11.37789969 -17.20456299 -12.32456299 7.746808873 10.81456927 6.481908112 -5.881480264 -1.385653811 -4.382207478 8.629503386 -5.827860058 0.624296976 10.00992495 -10.98777749 -9.827471845 4.766701702 13.85712327 5.457709958 -7.776521546 0.226735053 -2.474074724 -2.690066305 5.981727122 -3.608777006 4.190578344 -9.771942297 -15.91903344 10.5005696 18.89977568 -0.392274179 -2.763117308 -12.47978689 -8.764241479 1.111725495
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