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(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)
> x <- 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)
> 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 
-5.372 -2.838 -1.083  1.411 25.889 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)  5.37185    0.51940   10.34   <2e-16 ***
x            0.24966    0.00636   39.26   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 4.205 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/1f5e61354457650.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/2cfgt1354457650.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 
-3.288e-14 -9.760e-16  1.040e-16  7.730e-16  1.090e-14 

Coefficients:
             Estimate Std. Error   t value Pr(>|t|)    
(Intercept) 5.372e+00  3.067e-16 1.752e+16   <2e-16 ***
x           2.497e-01  3.755e-18 6.649e+16   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 2.483e-15 on 287 degrees of freedom
Multiple R-squared:     1,	Adjusted R-squared:     1 
F-statistic: 4.421e+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/37iqv1354457650.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 
-5.372 -2.838 -1.083  1.411 25.889 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)
(Intercept) 2.381e-17  5.194e-01       0        1
x           3.863e-19  6.360e-03       0        1

Residual standard error: 4.205 on 287 degrees of freedom
Multiple R-squared: 1.307e-34,	Adjusted R-squared: -0.003484 
F-statistic: 3.752e-32 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/4qmha1354457650.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 
-17.929  -7.194  -1.451   7.537  31.255 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)    
(Intercept)  2.330e+01  5.732e-01   40.65   <2e-16 ***
m$residuals -1.380e-17  1.368e-01    0.00        1    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 9.745 on 287 degrees of freedom
Multiple R-squared: 1.389e-32,	Adjusted R-squared: -0.003484 
F-statistic: 3.987e-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/5sqpc1354457650.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 
-5.402 -2.854 -1.032  1.407 25.919 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.004057   0.248191  -0.016    0.987
Residuals   -0.009944   0.059125  -0.168    0.867

Residual standard error: 4.212 on 286 degrees of freedom
  (2 observations deleted due to missingness)
Multiple R-squared: 9.89e-05,	Adjusted R-squared: -0.003397 
F-statistic: 0.02829 on 1 and 286 DF,  p-value: 0.8666 

> abline(m5)
> grid()
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6kmy11354457650.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/7m6p91354457650.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/8xfkx1354457650.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/9li5a1354457650.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/10a2pa1354457650.tab") 
> 
> try(system("convert tmp/1f5e61354457650.ps tmp/1f5e61354457650.png",intern=TRUE))
character(0)
> try(system("convert tmp/2cfgt1354457650.ps tmp/2cfgt1354457650.png",intern=TRUE))
character(0)
> try(system("convert tmp/37iqv1354457650.ps tmp/37iqv1354457650.png",intern=TRUE))
character(0)
> try(system("convert tmp/4qmha1354457650.ps tmp/4qmha1354457650.png",intern=TRUE))
character(0)
> try(system("convert tmp/5sqpc1354457650.ps tmp/5sqpc1354457650.png",intern=TRUE))
character(0)
> try(system("convert tmp/6kmy11354457650.ps tmp/6kmy11354457650.png",intern=TRUE))
character(0)
> try(system("convert tmp/7m6p91354457650.ps tmp/7m6p91354457650.png",intern=TRUE))
character(0)
> try(system("convert tmp/8xfkx1354457650.ps tmp/8xfkx1354457650.png",intern=TRUE))
character(0)
> try(system("convert tmp/9li5a1354457650.ps tmp/9li5a1354457650.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  5.380   0.850   6.225