Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
-1813.130876 -2101.143684 -2058.137618 -2267.706355 -2101.007565 -2368.612017 -2303.984657 -1874.249475 -2271.114711 -1868.681416 -1534.283836 -1863.492063 -1977.995434 -1826.136941 -2331.23331 -2358.156491 -2133.800693 -2364.103257 -2322.073607 -1858.604597 -2322.26296 -2088.771043 -1712.609985 -2362.32226 -2426.8728 -2251.227245 -2975.922001 -2422.926035 -2397.760266 -3206.088446 -2168.055411 -2180.274414 -2151.108645 -2596.63089 -2038.688158 -2154.002177 -2272.937488 -2554.003531 -2789.306095 -2918.99882 -3094.833728 -3074.804078 -2330.919969 -2594.956362 -2631.636955 -2755.382914 -2500.914581 -2961.141004 -3669.374169 -2764.087092 -3128.461087 -2273.937488 -2966.51971 -2613.09248 -2491.0507 -2299.564848 -2404.47051 -2541.405821 -2049.498805 -2315.724551 -2398.168622 -2314.967138 -2649.471187 -2370.594497 -2247.765654 -2880.986689 -2478.96175 -2156.025761 -1885.391659 -2651.388302 -2101.58169 -2090.634924
Data Y:
-1700.130876 -1990.143684 -1951.137618 -2164.706355 -2003.007565 -2269.612017 -2166.984657 -1727.249475 -2124.114711 -1729.681416 -1404.283836 -1735.492063 -1850.995434 -1703.136941 -2213.23331 -2245.156491 -2024.800693 -2253.103257 -2172.073607 -1699.604597 -2164.26296 -1941.771043 -1575.609985 -2225.32226 -2290.8728 -2118.227245 -2850.922001 -2302.926035 -2283.760266 -3090.088446 -2015.055411 -2018.274414 -1990.108645 -2447.63089 -1899.688158 -2020.002177 -2142.937488 -2428.003531 -2667.306095 -2801.99882 -2982.833728 -2961.804078 -2181.919969 -2437.956362 -2474.636955 -2608.382914 -2364.914581 -2829.141004 -3544.374169 -2641.087092 -3011.461087 -2159.937488 -2855.51971 -2501.09248 -2347.0507 -2149.564848 -2255.47051 -2406.405821 -1926.498805 -2200.724551 -2281.168622 -2203.967138 -2544.471187 -2268.594497 -2152.765654 -2787.986689 -2354.96175 -2025.025761 -1761.391659 -2536.388302 -1995.58169 -1985.634924
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