x=read.csv(“1.csv”,header = TRUE)
fm=lm(y1~x1+x2+x3+x4+x5,data=x)
library(mvstats)
coef.sd(fm)
$coef.sd
x1 x2 x3 x4 x5
-0.01895527 0.09943056 -7.18439825 -90.72823469 -98.56855034
summary(fm)
Call:
lm(formula = y1 ~ x1 + x2 + x3 + x4 + x5, data = x)
Residuals:
Min 1Q Median 3Q Max
-285.542 -31.564 -7.255 42.318 200.318
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 156143.462 155076.350 1.007 0.3186
x1 -12.037 32.505 -0.370 0.7126
x2 9.384 4.329 2.168 0.0347 *
x3 -1685.353 1548.830 -1.088 0.2815
x4 -1536.033 1551.011 -0.990 0.3265
x5 -1560.011 1550.751 -1.006 0.3190
Signif. codes: 0 ‘’ 0.001 ‘’ 0.01 ‘’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 81.78 on 53 degrees of freedom
Multiple R-squared: 0.9415, Adjusted R-squared: 0.936
F-statistic: 170.7 on 5 and 53 DF, p-value: < 2.2e-16