y <- c(1200,916,878,841,824,819,823,825,773,836,862,886,1010,846,911,856,881,830,830,827,773,797,826,947,1110,896,917,873,845,807,841,829,781,861,831,969,991,891,945,911,847,823,838,862,822,864,862,1044,1035,858,889,832,810,792,812,783,773,840,820,945)
x <- c(1770,2203,2836,1976,2837,2150,2180,2631,1781,2327,2260,2051,2250,2102,2957,2485,2871,2447,2570,2622,1840,2682,2369,2119,2531,2214,3206,2709,2734,2348,2702,2642,2064,2647,2534,2297,2718,2321,3112,2664,2808,2668,2934,2616,2228,2463,2416,2407,2582,2101,3305,2818,2401,3019,2507,2948,2210,2467,2596,2451)
par1 = '0'
par1 <- as.numeric(par1)
library(lattice)
z <- as.data.frame(cbind(x,y))
m <- lm(y~x)
summary(m)
postscript(file="/var/wessaorg/rcomp/tmp/1d5121324412088.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()
postscript(file="/var/wessaorg/rcomp/tmp/28ewm1324412088.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
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()
postscript(file="/var/wessaorg/rcomp/tmp/3dkxh1324412088.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
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()
postscript(file="/var/wessaorg/rcomp/tmp/4x3kf1324412088.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
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()
postscript(file="/var/wessaorg/rcomp/tmp/5iq5b1324412088.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)
abline(m5)
grid()
dev.off()
postscript(file="/var/wessaorg/rcomp/tmp/67xb31324412088.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()
postscript(file="/var/wessaorg/rcomp/tmp/7lsib1324412088.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()
postscript(file="/var/wessaorg/rcomp/tmp/87gpk1324412088.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
acf(m$residuals,main='Residual Autocorrelation Function')
dev.off()
postscript(file="/var/wessaorg/rcomp/tmp/9lrnt1324412088.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
qqnorm(x)
qqline(x)
grid()
dev.off()

#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/106kfj1324412088.tab") 

try(system("convert tmp/1d5121324412088.ps tmp/1d5121324412088.png",intern=TRUE))
try(system("convert tmp/28ewm1324412088.ps tmp/28ewm1324412088.png",intern=TRUE))
try(system("convert tmp/3dkxh1324412088.ps tmp/3dkxh1324412088.png",intern=TRUE))
try(system("convert tmp/4x3kf1324412088.ps tmp/4x3kf1324412088.png",intern=TRUE))
try(system("convert tmp/5iq5b1324412088.ps tmp/5iq5b1324412088.png",intern=TRUE))
try(system("convert tmp/67xb31324412088.ps tmp/67xb31324412088.png",intern=TRUE))
try(system("convert tmp/7lsib1324412088.ps tmp/7lsib1324412088.png",intern=TRUE))
try(system("convert tmp/87gpk1324412088.ps tmp/87gpk1324412088.png",intern=TRUE))
try(system("convert tmp/9lrnt1324412088.ps tmp/9lrnt1324412088.png",intern=TRUE))