Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
11.57467886 11.60253762 11.66544052 11.72235136 11.65106 11.67780411 11.74772878 11.78474489 11.6997867 11.72634994 11.77385023 11.80157956 11.71983394 11.73940744 11.77681933 11.7982247 11.70729087 11.71213691 11.76222787 11.7871879 11.695472 11.71994776 11.77614321 11.79451904 11.72595403 11.73551686 11.80294979 11.82449739 11.75594223 11.76714124 11.82912939 11.84736707 11.76821089 11.78755269 11.84519603 11.87147365 11.80520694 11.82858256 11.88355065 11.91355892 11.84076051 11.87240248 11.91861721 11.94929475 11.88371631 11.91010299 11.96081129 11.98970608 11.93670189 11.97041993 12.0034545 12.04337135 11.99701073 12.0183264 12.04707691 12.07071673 12.02193344 12.04496449 12.07909118 12.11206961 12.0565164 12.08872951 12.11315111 12.14467585 12.0831307 12.11873212 12.14403233 12.17895835 12.11743813 12.15222873 12.17460778 12.20856454 12.14495228 12.17423104 12.1915984 12.23144795 12.16395293 12.189435 12.20499206 12.24221169 12.18638511 12.21018418 12.2343831 12.2666455 12.21011447 12.22432014 12.24185476 12.25918678 12.18031885 12.2918833
Data Y:
4.608165695 4.62399194 4.625952725 4.636668853 4.641502115 4.665324109 4.679349584 4.699570861 4.715816706 4.744062185 4.753590191 4.773223771 4.790819533 4.836281907 4.852811209 4.868303386 4.873669439 4.894850261 4.899331225 4.90897164 4.914124394 4.935193099 4.934473933 4.93878119 4.932313327 4.948050419 4.95088529 4.954417614 4.955827058 4.973279508 4.97397131 4.980176087 4.989071116 5.007296393 5.009968405 5.011301739 5.016617366 5.029129877 5.031091303 5.036952602 5.042779747 5.05560866 5.065123479 5.073297055 5.076423035 5.094976443 5.098035484 5.102302483 5.098035484 5.108971195 5.109575242 5.11679549 5.120983351 5.139321635 5.14107859 5.147494477 5.14633101 5.158480421 5.159055299 5.157905213 5.158480421 5.170483995 5.173887288 5.182906515 5.188502501 5.200153118 5.202907182 5.208939555 5.21384821 5.227358278 5.23324537 5.242804657 5.245443877 5.256974403 5.260615499 5.26631057 5.268888556 5.286244785 5.294811227 5.30529298 5.313205979 5.32933142 5.333201768 5.346154696 5.352331955 5.372032404 5.381738975 5.37296091 5.351384087 5.359412466
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