Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
9256885 10333752 10720612 11907834 9690590 10609448 11977254 13244600 11240162 11576825 11807926 13588068 11691588 12562273 12438117 12928814 11548838 12776845 11997156 13330088 11435231 12572825 11345956 13470822 10775234 12519177 11934407 13812730 11072269 12875557 12839205 14102759 10644606 13551116 12790714 13436225 10543361 13588113 13910100 13688184 9460957 12983955 12072813 13410011 11595556 12711511 13192150 11446149 11595556 12188025 13224051 13065877 9609338 12607393 12501842 12689869 10650116 12591209 13880292 13514738 9758680 11408706 11300823 11987461 9644061 11129888 11516053 10674049 8831349 10896072 10658420 11588748 9034816 10147881 11040127 11968843 9204166 9828988 10599002 10383324 9344956 10236895 9834064 10900904 9573817 10054837 9099351 10830753 8890756 10240704 9407973 10746225 8685878 10120696 9299008 8928589 7499163 8814962 9215177 10158324 8373334 9485738 9333946 10910330 8816168 8923734
Data Y:
420535353 488878453 498041239 564072495 457511553 534305471 580509081 648770437 585105127 630821361 635048475 684750205 624579278 708146235 679748791 727949555 711704976 823818767 751357794 852585004 769520139 832010662 718771367 872040833 743812598 818570610 886246705 976236126 819407290 1026133596 953889829 1060970548 831565599 1082691438 1037403716 1100450194 877855090 1190802579 1145401598 1124835873 828092945 1113564364 1079819494 1220707994 1027966674 1255047042 1229675334 1039095435 1027966674 1243925959 1291053814 1369959521 1032152788 1312286401 1353955493 1423890029 1198277525 1459407892 1612018595 1528714149 1169666679 1403102002 1420269219 1510080957 1275608999 1471568598 1483636755 1482944107 1274746307 1636335358 1618764870 1787598977 1436519436 1673395701 1798072680 1944544177 1530417042 1726170719 1781938866 1819651590 2023539182 2310956246 2296138765 2489571272 2342282823 2622730494 2443617370 2758409400 2534556428 2871825467 2811245838 2962815256 2623492171 3127575470 2817693911 2647347460 2229602701 2625109728 2803549524 2950848868 2646481372 2928423819 2993317766 3259160572 2879505487 2938057843
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