library(linearModel)
y = c(30, 40, 60, 70, 100);y
## [1]  30  40  60  70 100
x = c(2, 3, 5, 6, 8);x
## [1] 2 3 5 6 8
z = c(6, 7, 8, 10, 15);z
## [1]  6  7  8 10 15
t = lm(y ~ x+z)
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)