Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data X:
104 71 55 87 40 32 85 5 43 49 73 58 34 100 55 95 91 53 99 101 73 74 23 92 90 64 87 77 39 57 60 66 73 14 107 56 70 58 34 69 61 64 45 34 36 47 21 68 94 63 80 7 86 21 31 85 62 127 111 25 87 116 44 51 49 54 71 59 83 50 58 21 52 71 94 87 76 44 66 56 36 81 28 93 59 27 66 53 41 90 41 105 43 60 72 58 95 108 74 142 45 97 22 49 56 92 64 86 41 49 62 74 11 70 30 46 93 4 72 39 70 46 13 46 16 48 86 82 87 105 124 62 25 72 111 6 70 41 76 57 19 69 103 16 71 34 74 131 0 5 0 0 0 0 80 93 0 0 6 13 3 18 0 33
Data Y:
78973 46146 46492 60656 21898 36555 74680 22807 61282 37981 41553 45081 38557 51641 30658 52924 79256 53462 68950 53639 67819 48333 28001 51665 39019 46221 65792 39858 19574 41829 78688 36781 44314 24874 56911 37048 48426 33388 26998 46502 41507 40001 33144 29501 43059 43249 29272 49821 98341 44372 42448 5950 64839 32551 30767 62046 71930 67328 67253 35373 85544 88087 30621 50580 49670 25456 69245 43787 53638 35683 38008 18801 44324 51408 53880 55708 63858 183643 35660 41664 29883 62047 33321 46553 56622 15430 49379 58215 38253 77786 21331 55292 30105 37651 59370 46216 73122 93927 55935 93308 74344 78094 25625 43750 28995 47336 57582 60875 165877 32984 61638 36367 1168 40530 21427 15024 39088 855 80455 14116 43915 76705 40112 41821 8773 52045 51491 53470 53211 63091 131634 41745 23656 51442 54574 35708 66627 39585 50029 25266 34860 62759 62307 37238 42452 59820 75075 97567 0 6023 0 0 0 0 42420 31116 0 0 1644 6179 3926 23238 0 38818
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