Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
258881.9875 255059.9883 247697.9873 244778.9879 240900.9888 239932.9876 270246.9882 283892.9885 282347.986 273569.9859 254755.9873 254353.9863 255842.9886 254489.9882 251994.9864 246338.9879 244018.9884 245952.9859 279805.9879 283110.9877 281096.9855 275963.9864 270693.9865 271900.9855 274411.9877 272432.9881 268360.9865 268585.9875 264767.9875 269973.9865 304743.989 309364.9879 308346.9857 298426.9868 289230.9865 291974.9869 294911.9875 293487.9873 290554.9853 284735.9877 281817.9866 287853.986 316262.9885 325411.9873 326010.9856 328281.9854 317479.9858 317538.9865 313736.987 312275.9872 309390.9854 302949.9874 300315.9868 304034.9857 333475.9874 337697.9871 335931.9864 323930.9852 313926.9859 314484.9864 313217.9867 309663.9868 302962.9864 298988.9856 298422.9873 301630.986 329764.987 335082.9882 327615.9853 309118.9855 295915.9872 291412.9859 291458.4609 284594.461 276391.4593 272482.4605 264897.461 263206.4597 296722.4612
Data Y:
208069.33 204924.3308 200204.3298 198002.3304 195100.3314 191309.3301 213682.3307 225694.331 230493.3285 229175.3283 216142.3298 216627.3287 220120.3312 220029.3308 218358.3289 214826.3304 210619.331 209587.3283 236955.3304 241995.3303 241792.3279 242535.3288 238459.329 240251.3279 244666.3303 244490.3306 241486.3289 240455.33 236003.33 236911.3289 264493.3316 270263.3304 269559.3281 266951.3293 258027.3289 262727.3294 267327.33 267280.3298 264691.3277 258777.3302 254758.3291 254782.3285 277181.3311 285265.3298 286516.328 282956.3278 276601.3282 277829.329 277042.3294 277017.3297 274951.3278 270064.3299 267054.3293 264907.3281 287173.3299 291100.3296 292214.3288 288100.3275 281391.3283 282570.3289 280104.3292 280322.3292 276750.3289 275130.328 274266.3298 271225.3284 289716.3295 290640.3307 292214.3276 278420.3278 269740.3297 265775.3283 268948.3304 264090.3305 255112.3287 253267.33 245971.3306 235286.3291 258470.3308
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