R code for autompg regression

tmp<-read.table("autompg.Data", header=TRUE);

y<-tmp[,1];

n<-nrow(tmp); p<-8;

x<-matrix(0,nrow(tmp),(p-1));

for(i in 1:(p-1)) { x[,i]<-tmp[,(i+1)];}

mylm<-lm(y ~ x);

summary(mylm);

tmp<-read.table("autompg.Data", header=TRUE);

Y<-tmp[,1];

X<-matrix(0,nrow(tmp),p);

X[,1]<-1;

for(i in 2:p) { X[,i]<-tmp[,i];}

##Regression estimates by matrix

solve(t(X) %*% X) %*% t(X) %*% Y