R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > y <- c(94,103,93,103,51,70,91,22,38,93,60,123,148,90,124,70,168,115,71,66,134,117,108,84,156,120,114,94,120,81,110,133,122,158,109,124,39,92,126,0,70,37,38,120,93,95,77,90,80,31,110,66,138,133,113,100,7,140,61,41,96,164,78,49,102,124,99,129,62,73,114,99,70,104,116,91,74,138,67,151,72,120,115,105,104,108,98,69,111,99,71,27,69,107,73,107,93,129,69,118,73,119,104,107,99,90,197,36,85,139,106,50,64,31,63,92,106,63,69,41,56,25,65,93,114,38,44,87,110,0,27,83,30,80,98,82,0,60,28,9,33,59,49,115,140,49,120,66,21,124,152,139,38,144,120,160,114,39,78,119,141,101,56,133,83,116,90,36,50,61,97,98,78,117,148,41,105,55,132,44,21,50,0,73,86,0,13,4,57,48,46,48,32,68,87,43,67,46,46,56,48,44,60,65,55,38,52,60,54,86,24,52,49,61,61,81,43,40,40,56,68,79,47,57,41,29,3,60,30,79,47,40,48,36,42,49,57,12,40,43,33,77,43,45,47,43,45,50,35,7,71,67,0,62,54,4,25,40,38,19,17,67,14,30,54,35,59,24,58,42,46,61,3,52,25,40,32,4,49,63,67,32,23,7,54,37,35,51,39) > x <- c(30,28,38,30,22,26,25,18,11,26,25,38,44,30,40,34,47,30,31,23,36,36,30,25,39,34,31,31,33,25,33,35,42,43,30,33,13,32,36,0,28,14,17,32,30,35,20,28,28,39,34,26,39,39,33,28,4,39,18,14,29,44,21,16,28,35,28,38,23,36,32,29,25,27,36,28,23,40,23,40,28,34,33,28,34,30,33,22,38,26,35,8,24,29,20,29,45,37,33,33,25,32,29,28,28,31,52,21,24,41,33,32,19,20,31,31,32,18,23,17,20,12,17,30,31,10,13,22,42,1,9,32,11,25,36,31,0,24,13,8,13,19,18,33,40,22,38,24,8,35,43,43,14,41,38,45,31,13,28,31,40,30,16,37,30,35,32,27,20,18,31,31,21,39,41,13,32,18,39,14,7,17,0,30,37,0,5,1,16,32,24,17,11,24,22,12,19,13,17,15,16,24,15,17,18,20,16,16,18,22,8,17,18,16,23,22,13,13,16,16,20,22,17,18,17,12,7,17,14,23,17,14,15,17,21,18,18,17,17,16,15,21,16,14,15,17,15,15,10,6,22,21,1,18,17,4,10,16,16,9,16,17,7,15,14,14,18,12,16,21,19,16,1,16,10,19,12,2,14,17,19,14,11,4,16,20,12,15,16) > par1 = '0' > par1 <- as.numeric(par1) > library(lattice) > z <- as.data.frame(cbind(x,y)) > m <- lm(y~x) > summary(m) Call: lm(formula = y ~ x) Residuals: Min 1Q Median 3Q Max -93.822 -7.161 2.969 10.839 31.178 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -6.86445 2.20096 -3.119 0.002 ** x 3.37657 0.08601 39.257 <2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 15.47 on 287 degrees of freedom Multiple R-squared: 0.843, Adjusted R-squared: 0.8425 F-statistic: 1541 on 1 and 287 DF, p-value: < 2.2e-16 > postscript(file="/var/wessaorg/rcomp/tmp/1wnq71354454083.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(z,main='Scatterplot, lowess, and regression line') > lines(lowess(z),col='red') > abline(m) > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/2e3gb1354454083.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > m2 <- lm(m$fitted.values ~ x) > summary(m2) Call: lm(formula = m$fitted.values ~ x) Residuals: Min 1Q Median 3Q Max -1.063e-13 -1.976e-15 1.339e-15 3.233e-15 1.698e-14 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -6.864e+00 1.214e-15 -5.653e+15 <2e-16 *** x 3.377e+00 4.745e-17 7.115e+16 <2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 8.533e-15 on 287 degrees of freedom Multiple R-squared: 1, Adjusted R-squared: 1 F-statistic: 5.063e+33 on 1 and 287 DF, p-value: < 2.2e-16 > 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() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/3tz921354454083.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > m3 <- lm(m$residuals ~ x) > summary(m3) Call: lm(formula = m$residuals ~ x) Residuals: Min 1Q Median 3Q Max -93.822 -7.161 2.969 10.839 31.178 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -7.419e-17 2.201e+00 0 1 x 4.366e-18 8.601e-02 0 1 Residual standard error: 15.47 on 287 degrees of freedom Multiple R-squared: 8.916e-38, Adjusted R-squared: -0.003484 F-statistic: 2.559e-35 on 1 and 287 DF, p-value: 1 > 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() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/4ubw21354454083.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > m4 <- lm(m$fitted.values ~ m$residuals) > summary(m4) Call: lm(formula = m$fitted.values ~ m$residuals) Residuals: Min 1Q Median 3Q Max -78.678 -24.652 -4.393 25.996 96.904 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 7.181e+01 2.108e+00 34.06 <2e-16 *** m$residuals -2.737e-17 1.368e-01 0.00 1 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 35.84 on 287 degrees of freedom Multiple R-squared: 1.479e-32, Adjusted R-squared: -0.003484 F-statistic: 4.245e-30 on 1 and 287 DF, p-value: 1 > 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() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5nnlw1354454083.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > 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) Call: lm(formula = z5) Residuals: Min 1Q Median 3Q Max -94.093 -7.241 2.971 10.588 31.402 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.002495 0.912345 0.003 0.998 Residuals -0.035031 0.059124 -0.593 0.554 Residual standard error: 15.48 on 286 degrees of freedom (2 observations deleted due to missingness) Multiple R-squared: 0.001226, Adjusted R-squared: -0.002266 F-statistic: 0.3511 on 1 and 286 DF, p-value: 0.554 > abline(m5) > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6gots1354454084.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(m$residuals,main='Residual Histogram',xlab='Residuals') > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7ek6t1354454084.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > 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() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/8i85u1354454084.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(m$residuals,main='Residual Autocorrelation Function') > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/9hjla1354454084.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(x) > qqline(x) > grid() > dev.off() null device 1 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/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="/var/wessaorg/rcomp/tmp/108zrb1354454084.tab") > > try(system("convert tmp/1wnq71354454083.ps tmp/1wnq71354454083.png",intern=TRUE)) character(0) > try(system("convert tmp/2e3gb1354454083.ps tmp/2e3gb1354454083.png",intern=TRUE)) character(0) > try(system("convert tmp/3tz921354454083.ps tmp/3tz921354454083.png",intern=TRUE)) character(0) > try(system("convert tmp/4ubw21354454083.ps tmp/4ubw21354454083.png",intern=TRUE)) character(0) > try(system("convert tmp/5nnlw1354454083.ps tmp/5nnlw1354454083.png",intern=TRUE)) character(0) > try(system("convert tmp/6gots1354454084.ps tmp/6gots1354454084.png",intern=TRUE)) character(0) > try(system("convert tmp/7ek6t1354454084.ps tmp/7ek6t1354454084.png",intern=TRUE)) character(0) > try(system("convert tmp/8i85u1354454084.ps tmp/8i85u1354454084.png",intern=TRUE)) character(0) > try(system("convert tmp/9hjla1354454084.ps tmp/9hjla1354454084.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.953 0.907 5.912