> lm.M<-lm(y~x1+x2+x4+x5+x6,data=ms)
> summary(lm.M)

Call:
lm(formula = y ~ x1 + x2 + x4 + x5 + x6, data = ms)

Residuals:
        1         2         3         4         5         6         7
-0.006811 -3.602484  0.138659  0.652129  0.776480  0.657713  1.384316

Coefficients:
              Estimate Std. Error t value Pr(>|t|)  
(Intercept)  4.065e+02  5.229e+01   7.773   0.0815 .
x1          -4.763e-01  7.405e-02  -6.431   0.0982 .
x2          -5.759e-01  2.261e-01  -2.547   0.2382  
x4          -3.843e+00  1.303e+00  -2.948   0.2082  
x5          -3.109e+02  6.736e+01  -4.615   0.1358  
x6          -3.635e+04  3.649e+04  -0.996   0.5012  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 4.047 on 1 degrees of freedom
Multiple R-squared:  0.9979,    Adjusted R-squared:  0.9876
F-statistic: 96.68 on 5 and 1 DF,  p-value: 0.07705