用relevel()函数,如
dat$x5=relevel(dat$x5,ref="m")#将性别x5作为虚拟变量,以男性m作为参考水平,女性
f
fitlm=lm(y~x4+x5,data=dat) #建立y关于x4和x5的回归方程
summary(fitlm) #输出回归结果
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -11556.0 13143.5 -0.879 0.3856
x4 4434.4 915.5 4.844 2.92e-05 ***
x5f -16840.5 6344.9 -2.654 0.0121 *
---
得到回归方程:
男性:y=-11556+4434.4x4
女性:y=-11556+4434.4x4-16840.5