Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
16.1387531550 14.1499707810 12.8051002773 12.2826650254 11.5938826513 11.1163179032 11.0163179032 11.1163179032 10.7826650254 10.4826650254 14.4377945217 14.4826650254 14.4826650254 13.7938826513 12.7051002773 12.3051002773 11.1275355291 10.6499707810 9.7611884069 9.5611884069 9.1611884069 9.4275355291 13.7602297736 14.8153592699 14.7817063921 13.0929240181 11.4153592699 10.2265768958 9.9377945217 9.6490121477 9.9377945217 11.4153592699 10.6265768958 10.0041416440 13.0592711403 12.1929240181 10.8817063921 10.0592711403 9.0480535144 9.2368358884 9.2704887662 8.4041416440 6.9041416440 6.2592711403 5.6256182625 7.7919653848 14.3256182625 14.9490121477 13.3602297736 11.0153592699 8.8256182625 9.7919653848 11.9368358884 12.5377945217 12.3051002773 11.2275355291 9.9051002773 10.6714473995 14.4826650254 14.5387531550 13.6499707810
Data Y:
5.2015218082 4.9781165096 4.1717377038 3.4185483008 3.3951430023 3.6483324052 3.7483324052 3.6483324052 3.3185483008 3.1185483008 3.3121694949 4.3185483008 4.5185483008 4.3951430023 3.7717377038 3.3717377038 3.4249271067 3.4781165096 3.4547112111 3.2547112111 3.0547112111 3.1249271067 3.2653588979 4.1589800920 4.4291959876 4.2057906890 3.6589800920 3.3355747935 3.3121694949 3.3887641964 3.4121694949 3.4589800920 3.5355747935 3.1823853905 2.7760065846 2.8057906890 2.5291959876 2.1760065846 2.3994118831 2.4228171817 2.4526012861 2.2823853905 1.9823853905 1.5760065846 1.3462224802 1.1164383758 1.4462224802 2.4887641964 2.9653588979 2.6589800920 2.4462224802 2.1164383758 2.3228171817 2.6121694949 2.5717377038 2.3249271067 1.8717377038 1.4419535993 1.5185483008 2.2015218082 2.5781165096
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