Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
-7.2267351 -6.290090218 -5.490850323 -5.390850323 -5.503065284 -5.440470271 -5.403065284 -5.253445336 -4.978635363 -5.21604035 -6.427495205 -6.602305179 -6.564900192 -5.840470271 -5.390850323 -5.32825531 -5.265660297 -5.065660297 -4.790850323 -4.51604035 -4.71604035 -5.01604035 -6.415280244 -6.890090218 -6.615280244 -5.92825531 -5.41604035 -5.21604035 -5.253445336 -5.253445336 -5.353445336 -5.490850323 -5.241230376 -4.954205442 -5.091610428 -4.879395468 -4.52977552 -4.604585494 -4.604585494 -4.541990481 -4.204585494 -3.892370533 -3.505345599 -3.293130639 -3.180915678 -3.530535625 -4.716800455 -5.303825389 -5.054205442 -4.904585494 -4.454965547 -4.467180507 -4.554205442 -4.391610428 -4.116800455 -3.72977552 -3.342750586 -3.480155573 -4.279395468 -4.629015415
Data Y:
19.41535803 17.7617828 16.28082704 15.58082704 14.89868545 14.57130492 14.49868545 14.70820757 14.16296863 13.8355881 17.13440227 17.47964121 17.50702174 17.07130492 16.18082704 15.75344651 14.72606598 14.42606598 13.68082704 13.5355881 13.1355881 13.1355881 16.71654386 17.3617828 17.11654386 15.65344651 14.2355881 13.1355881 12.90820757 12.70820757 12.90820757 14.18082704 13.59034916 12.92725181 15.59987128 15.08201288 13.791535 12.73677394 11.63677394 11.70939341 12.03677394 11.51891553 10.15581818 9.237959766 8.420101357 10.21057924 16.65463234 18.11772969 16.72725181 14.13677394 11.34629606 11.96415447 14.32725181 15.69987128 16.05463234 15.291535 13.92843765 14.40105712 18.08201288 18.47249075
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