Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
9492 8767 1423 15387 11936 22386 7703 7694 14513 12552 10893 14584 5469 8334 8651 13657 22452 6023 10602 12476 783 16597 5506 5852 6630 4065 12579 10432 50116 8628 3482 7195 11938 15129 5050 0 7537 8665 3710 11128 7980 6851 15966 22214 7713 6154 2928 10805 17412 9243 6656 3585 11016 14667 11708 7619 11853 9540 10081 2574 3772 3886 8984 5956 10878 9852 13736 8033 9776 1536 6948 5846 7318 4898 1358 0 9580 17612 2941 6924 6387 3690 7153 8765 5761 8809 10686 8055 8352 6030 14070 17349 5645 4120 13983 3853 19150 3895 4408 7498 2325 11429 17265 6265 3087 3979 8931 3721 0 6185 14254 7500 5890 0 0 0 10411 28040 103 14377 2781 8759 2229 2423 13656 2790 10255 4845 5264 0 0 0 0 0 0 0 0 7711 797 0 0 4010 5959 2338
Data Y:
2845 2234 603 6113 4648 4632 2404 2159 5660 3571 5285 5774 2597 2264 1233 2292 7207 1 4131 1824 495 5345 594 881 2398 1818 3241 4053 4373 3887 1 2662 3044 31 1658 0 2233 3223 783 3011 2196 2591 1647 1890 1353 1883 1095 2079 4174 2808 1783 1780 1917 2865 2290 2489 2584 2098 3715 76 2356 967 1091 1588 1881 2636 5431 1833 3009 516 1574 416 1377 2019 437 0 2881 1997 1706 1651 40 1730 1825 3202 302 2316 4208 1584 2204 2391 2571 2458 1191 471 2280 1978 2540 786 2282 1806 593 1235 2242 2395 953 460 3009 983 0 2178 349 1830 1992 0 0 0 4975 1768 52 2245 4 1335 624 84 358 1186 3070 1748 1815 0 0 0 0 0 0 0 0 1769 361 0 0 1239 1993 753
Sample Range:
(leave blank to include all observations)
From:
To:
bandwidth of density plot
(?)
Chart options
Label y-axis:
Label x-axis:
R Code
par1 <- as.numeric(par1) library(lattice) z <- as.data.frame(cbind(x,y)) m <- lm(y~x) summary(m) bitmap(file='test1.png') plot(z,main='Scatterplot, lowess, and regression line') lines(lowess(z),col='red') abline(m) grid() dev.off() bitmap(file='test2.png') m2 <- lm(m$fitted.values ~ x) summary(m2) z2 <- as.data.frame(cbind(x,m$fitted.values)) names(z2) <- list('x','Fitted') plot(z2,main='Scatterplot, lowess, and regression line') lines(lowess(z2),col='red') abline(m2) grid() dev.off() bitmap(file='test3.png') m3 <- lm(m$residuals ~ x) summary(m3) z3 <- as.data.frame(cbind(x,m$residuals)) names(z3) <- list('x','Residuals') plot(z3,main='Scatterplot, lowess, and regression line') lines(lowess(z3),col='red') abline(m3) grid() dev.off() bitmap(file='test4.png') m4 <- lm(m$fitted.values ~ m$residuals) summary(m4) z4 <- as.data.frame(cbind(m$residuals,m$fitted.values)) names(z4) <- list('Residuals','Fitted') plot(z4,main='Scatterplot, lowess, and regression line') lines(lowess(z4),col='red') abline(m4) grid() dev.off() bitmap(file='test5.png') myr <- as.ts(m$residuals) z5 <- as.data.frame(cbind(lag(myr,1),myr)) names(z5) <- list('Lagged Residuals','Residuals') plot(z5,main='Lag plot') m5 <- lm(z5) summary(m5) abline(m5) grid() dev.off() bitmap(file='test6.png') hist(m$residuals,main='Residual Histogram',xlab='Residuals') dev.off() bitmap(file='test7.png') if (par1 > 0) { densityplot(~m$residuals,col='black',main=paste('Density Plot bw = ',par1),bw=par1) } else { densityplot(~m$residuals,col='black',main='Density Plot') } dev.off() bitmap(file='test8.png') acf(m$residuals,main='Residual Autocorrelation Function') dev.off() bitmap(file='test9.png') qqnorm(x) qqline(x) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Simple Linear Regression',5,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Statistics',1,TRUE) a<-table.element(a,'Estimate',1,TRUE) a<-table.element(a,'S.D.',1,TRUE) a<-table.element(a,'T-STAT (H0: coeff=0)',1,TRUE) a<-table.element(a,'P-value (two-sided)',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'constant term',header=TRUE) a<-table.element(a,m$coefficients[[1]]) sd <- sqrt(vcov(m)[1,1]) a<-table.element(a,sd) tstat <- m$coefficients[[1]]/sd a<-table.element(a,tstat) pval <- 2*(1-pt(abs(tstat),length(x)-2)) a<-table.element(a,pval) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'slope',header=TRUE) a<-table.element(a,m$coefficients[[2]]) sd <- sqrt(vcov(m)[2,2]) a<-table.element(a,sd) tstat <- m$coefficients[[2]]/sd a<-table.element(a,tstat) pval <- 2*(1-pt(abs(tstat),length(x)-2)) a<-table.element(a,pval) 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