Forbes<-read.table("P:/linmal/Forbes.dtr", header=TRUE)
attach(Forbes)
summary(Forbes)
T<- 1/(Temp+273.15)
m1<-lm(log(Press)~T)
names(m1)
t<-(m1$coefficients[2]+4895.36)/(sqrt(vcov(m1)[2, 2]))
p<-2*(1-pt(abs(t), m1$df))



Typpi<-read.table("P:/linmal/Stack.dtr", header=TRUE)
attach(Typpi)
summary(Typpi)
m1<-lm(Y~X2)
m2<-lm(X3~X2)
plot(m1$residuals~m2$residuals)

R_osa<-cov(m1$residuals, m2$residuals)/sqrt(var(m1$residuals)*var(m2$residuals))


m3<-lm(Y~X2+X3)
summary(m3)
t<-(m3$coefficients[3])/(sqrt(vcov(m3)[3, 3]))
p<-2*(1-pt(abs(t), m3$df))
m4<-lm(X4~X2+X3)
plot(m3$residuals~m4$residuals)
R_osa<-cov(m3$residuals, m4$residuals)/sqrt(var(m3$residuals)*var(m4$residuals))
m5<-lm(Y~X2+X3+X4)
summary(m5)
plot(m3$residuals~X3)
X5<-(X3)^2
m6<-lm(Y~X2+X3+X5)
summary(m6)