Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
-4.213500441 -3.935935693 -3.536894327 -3.026635334 -2.926635334 -2.992982456 -3.08176483 -2.992982456 -2.93785296 -2.760288212 -3.249070586 -3.825676701 -4.003241449 -3.814459075 -3.18176483 -2.826635334 -2.715417708 -2.604200082 -2.504200082 -2.326635334 -2.149070586 -2.449070586 -3.049070586 -4.159329578 -4.636894327 -4.359329578 -3.715417708 -3.349070586 -3.249070586 -3.23785296 -3.33785296 -3.53785296 -3.526635334 -3.371505838 -3.327593967 -3.116376341 -2.950029219 -2.794899723 -3.072464471 -3.172464471 -2.961246845 -2.572464471 -2.306117349 -2.262205478 -2.295858356 -2.329511234 -2.38464073 -2.438811593 -2.782723464 -2.827593967 -3.272464471 -3.217334975 -3.083682097 -2.627593967 -2.116376341 -1.738811593 -1.494899723 -1.350987852 -1.339770226 -1.550029219 -1.805158715
Data Y:
-3.585860839 -3.329201771 -2.774235967 -2.247599696 -2.147599696 -2.332588297 -2.360917831 -2.332588297 -2.119270162 -1.962611095 -2.290940628 -3.0025655 -3.159224568 -3.030895034 -2.460917831 -2.047599696 -2.07592923 -2.004258764 -2.004258764 -1.847599696 -1.690940628 -1.890940628 -2.190940628 -3.017576899 -3.374235967 -3.117576899 -2.57592923 -2.290940628 -2.290940628 -2.319270162 -2.319270162 -2.419270162 -2.447599696 -2.234281561 -2.092633892 -1.920963426 -1.735974824 -1.522656689 -1.779315756 -1.779315756 -1.70764529 -1.479315756 -1.294327155 -0.952679486 -0.967690884 -0.882702283 -0.996020418 -1.464304358 -1.905952027 -1.792633892 -1.879315756 -1.665997621 -1.750986223 -1.592633892 -1.420963426 -1.164304358 -0.822656689 -0.58100902 -0.609338554 -0.935974824 -1.249292959
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