Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
60635408600.00 64463389370.00 66750859382.00 69032584255.95 72841303906.58 72838685607.53 75887977453.76 78936158549.66 81986612544.38 85038924238.38 87328862789.67 85794984680.72 88082838620.17 91885510657.40 96355573323.80 101321342608.80 105791395836.80 113241494103.00 117215849652.80 120940898880.60 125658810316.50 130873880200.80 139583536623.40 148226520849.00 153789461350.00 161871513310.40 171781295610.40 178996585758.20 176620960002.00 186604680395.00 187772917033.20 193109744878.60 197630528464.00 206483683756.40 203906590621.60 206783719069.34 206759082201.76 211878483671.40 214009149959.00 217203709135.40 222331811922.60 232825724880.00 241180274038.70 248494123822.10 253049202253.08 256922518700.16 254451243638.75 262662316800.94 268926171539.67 272037080259.13 281118338865.04 286389613939.39 296190694318.91 307294826961.69 309697986635.60 314369521231.48 317533640477.42 327005697520.69 332458473876.00 339794119726.27
Data Y:
79717000 82855667 79346667 76270333 89532667 83541333 82456000 85818333 92106667 96873333 96580000 87769000 88396000 90995667 92785000 98112667 105779667 103660333 105435000 105204000 107466333 118554333 123614333 125620000 121524333 130830333 138871333 135215667 122096333 130005333 126566000 134808667 138761333 133778333 122598667 116446000 100745333 104419333 103477000 102054333 102446667 99528000 106549667 107587333 117091333 114029667 109721333 114697000 115063667 121234667 118393000 121308000 117579000 117414000 117740333 109417000 116442333 113311000 109391333 107198667
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