Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
64.42592896 47.69064634 61.95277662 69.48479846 66.46888521 72.11049407 75.08901594 51.9287114 74.0142332 70.04903748 61.37521144 54.3296681 52.12640731 44.72803771 43.39112469 55.16936356 56.31634195 60.38853759 46.72545061 53.26821152 60.80579823 51.78432011 41.53551598 37.18268971 43.88575516 44.40186375 50.43666803 42.19860296 52.12382022 60.64808078 42.37779853 39.76284198 67.37262435 52.22823305 46.60512454 54.40742862 37.18527681 42.77616814 56.51960273 52.31634195 40.54625504 48.2149069 52.60771163 42.56773317 44.51960273 39.83503763 64.54366795 55.07568979 35.35373332 45.55440701 41.06495073 41.3296681 59.84060251 47.51960273 52.35373332 52.12382022 39.80838532 48.20158074 43.60512454 40.0168203
Data Y:
0.169087676 -0.973912324 -0.195912324 3.597087676 1.985087676 -1.740912324 2.351087676 -5.182912324 4.311087676 -0.247912324 4.116087676 3.221087676 -3.128912324 -2.453912324 6.728087676 1.687087676 0.688087676 -2.350912324 -0.532912324 2.714087676 1.803087676 -5.104912324 0.650087676 -4.779912324 1.001087676 0.182087676 2.623087676 -3.167912324 3.792087676 7.348087676 -1.804912324 -2.012912324 3.037087676 -4.884912324 -1.467912324 -5.521912324 -0.700912324 5.520087676 1.038087676 -3.311912324 -4.895912324 -1.417912324 -0.388912324 -0.987912324 2.038087676 0.948087676 1.025087676 -2.181912324 0.208087676 -3.520912324 -0.635912324 -0.778912324 4.244087676 -2.961912324 2.208087676 -2.207912324 -2.117912324 0.049087676 -2.467912324 2.390087676
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