Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
95077872409 89058674329 83654571361 85249400625 86973087744 86135206144 84422208025 81074589696 79421385124 82859925316 1.00022E+11 1.05893E+11 1.06283E+11 1.07769E+11 1.00794E+11 1.00831E+11 98430905169 97516300176 95722790881 91778702500 90189699856 92437281225 1.11206E+11 1.1404E+11 1.1285E+11 1.04931E+11 98550161329 98900815225 98105515524 95891792896 91786579369 89394422121 89056286929 90981260161 1.08745E+11 1.12281E+11 1.07332E+11 95554556161 87566279056 84921536569 84996737764 81041563684 76438425625 74292224356 70214930361 69321624100 88093801636 92171745604 82365556036 76411886329 70981747776 71370725409 72028361161 68917800484 65301713764 64088972964 59439902809 62871049081 78649398025 81371556049 73427992576 68160677776 65332893609 67795661376 69644793409 69849732681 69314252176 68944055184 65621531889 69812736841 86353699600 90428308369
Data Y:
72708426025 71308759369 66622320769 69070672969 71509712569 71484577956 70106859729 67010052769 64945464336 64957697424 76876989289 81425193201 82140706404 80112773764 76555695969 77236747225 76799928384 76786072609 75645351369 72981022500 71363779600 70221290049 82517733081 84789286596 85439290000 83051170596 79229301529 79894414336 78506436100 78628646464 76638170896 75743846656 75269019904 73609658721 83985199204 84521607076 85439290000 77565592036 72806070276 70682071321 72379293156 69788958976 65126019204 64187742609 60544047249 55399978384 66851205136 68117346049 64853243569 62821913449 59254270084 61060881025 61772628681 60044111521 56206926400 56209297225 50874606916 51455931921 61471268356 61669278889 60993686961 60073029604 60645465169 65415735225 69864533761 72010112409 74554118116 75055725369 71518804900 73980192049 85679144100 87545566161
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation