Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
12.56147528 12.56584977 12.55335058 12.53064218 12.53507059 12.53223477 12.53214456 12.52466091 12.50673264 12.49552814 12.48745869 12.56813953 12.58171752 12.58553595 12.57136138 12.54780602 12.55198589 12.54322322 12.54400096 12.53118055 12.52531152 12.52216723 12.51102104 12.57695321 12.58013652 12.58553595 12.53719488 12.50553259 12.49072889 12.50259304 12.48437083 12.44979499 12.44253905 12.4133185 12.36892252 12.46286758 12.47224887 12.44769638 12.43178489 12.40255184 12.41756863 12.4233631 12.40917266 12.37615292 12.37617401 12.32631488 12.33199579 12.42091786 12.42252587 12.4170181 12.40941341 12.41415535 12.45201433 12.48491199 12.5000362 12.51739556 12.52074834 12.49661313 12.51353161 12.58693764 12.59771263 12.58894785
Data Y:
12.69284735 12.6946864 12.70162827 12.6681701 12.66835592 12.65631033 12.65164269 12.64236113 12.62132305 12.61259053 12.62489811 12.71732618 12.72990728 12.72466404 12.68828581 12.65691575 12.65869166 12.65465471 12.64324312 12.62136596 12.60816206 12.60626722 12.64636099 12.70613556 12.72213354 12.69959747 12.6414816 12.59783091 12.58249678 12.58293936 12.559114 12.52987568 12.51563607 12.4874134 12.48101137 12.60083401 12.62345974 12.56721659 12.52970205 12.4928443 12.49557681 12.50016289 12.47809017 12.45114206 12.44176908 12.4041158 12.43217581 12.54413291 12.56114581 12.50978554 12.47256683 12.45138074 12.46988202 12.48333689 12.48480605 12.48095819 12.47828061 12.45358485 12.48454115 12.59085874 12.6139116 12.56801768
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