Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
868 1022 438 594 946 797 775 978 793 813 782 892 841 991 394 536 906 734 789 1006 785 775 858 710 828 922 535 729 754 641 723 716 615 679 629 659 875 784 348 392 686 631 674 724 658 628 653 651 739 703 331 386 732 729 596 684 707 739 678 608 681
Data Y:
276444 289742 303725 298305 266795 259497 266148 271037 276239 279681 277509 271115 275902 287224 300713 293860 264221 256167 262572 263276 264291 263903 260376 255603 261076 270976 285257 280445 250741 243803 253158 255542 262522 268381 267153 266424 276427 286994 303598 296806 263290 264981 272566 276475 284678 291542 291413 295916 309119 327616 335083 329765 301631 298423 298989 302963 309664 313218 314485 313927 323931
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation