R中的虚拟变量

Question

我正在构建一个线性模型，以评估距栖息地边界的距离对昆虫数量的影响。 使用的设备存在一些差异，因此我将设备作为分类变量包括在内，以确保它不会对丰富度产生重大影响。

类别因子为3级，因此我要求r通过使用以下代码在lm中生成虚拟变量：

lm(Richness ~ Distances + factor(Equipment), data = Data) 

当我要求模型的摘要时，我可以看到两个级别及其系数。 我假设这意味着r使用一个级别作为“标准”来比较其他级别的系数。

我如何找到第三级的系数以查看其对模型的影响？

谢谢

Answer 1

您可以执行lm(y~x-1)来删除截距，在您的情况下，截距是因素之一的参考水平。 话虽这么说，但出于统计原因，使用其中一个级别作为参考。

Answer 2

要确定如何提取系数，这是一个简单的示例：

# load data 
data(mtcars)
head(mtcars)

# what are the means of wt given the factor carb?
(means <- with(mtcars, tapply(wt, factor(carb), mean)))

# run the lm
mod <- with(mtcars, lm(wt~factor(carb)))

# extract the coefficients
coef(mod)

# the intercept is the reference level (i.e., carb 1)
coef(mod)[1]
coef(mod)[2:6]
coef(mod)[1] + coef(mod)[2:6]
means

因此，您可以看到在这种简单情况下，系数被简单地添加到参考电平（即截距）。 但是，如果您有协变量，它将变得更加复杂

mod2 <- lm(wt ~ factor(carb) + disp, data=mtcars)
summary(mod2)

当disp = 0时，截距现在为carb 1。