library(linearModel)
y=c(30,40,60,70,100)
x=c(2,3,5,6,8)
z=c(6,7,8,10,15)
t=lm(y~x+z)
cor(x,y)
## [1] 0.9941348
cor(y,z)
## [1] 0.9733995
cor(x,z)
## [1] 0.9461426
anova(t)
## Analysis of Variance Table
## 
## Response: y
##           Df  Sum Sq Mean Sq  F value    Pr(>F)    
## x          1 2964.91 2964.91 1384.386 0.0007216 ***
## z          1   30.80   30.80   14.383 0.0630239 .  
## Residuals  2    4.28    2.14                       
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
plot(t)