library(linearModel)
y=c(30,40,60,70,100);y
x=c(2,3,5,6,8);x
z=c(6,7,8,10,15);z
t=lm(y~x+z)
anova(t)
plot(t)