R中享乐价格分析中的邻域固定效应

Question

我正在使用R进行享乐价格分析，目的是使用差异法来估计当地震荡对县的影响。

我的因变量是房屋的价格，我想将未观察到的邻域因素作为我估计的参数的一部分。我认为我应该在邮政编码级别包含固定效果，并具有进行此操作的数据。

我当前的型号规格如下：

reg_1 = lm(logAppr18 ~ Dist_CBD_Miles + SqFeet + Bedrooms + BathsTotal + YearBuilt.f + year_dum + Bio + year_dum:Bio, data = reg_1)

我如何将邮政编码级别的固定效果合并到R中的此模型中？

Answer 1

1）提供可复制的示例是获得快速而准确的响应的最佳选择：[How to make a great R reproducible example?

2）如果您询问如何编码要在固定效果模型中使用的虚拟变量，则可以使用model.matrix。下面是一个如何使用它的示例。

dat <- data.frame("value" = c(4,3,2,3), "zip"=c("20002","20021","20021","20202"))

dummy <- model.matrix(~zip,dat)

newdat <- cbind(dat,dummy)

>newdat
  value   zip (Intercept) zip20021 zip20202
1     4 20002           1        0        0
2     3 20021           1        1        0
3     2 20021           1        1        0
4     3 20202           1        0        1

然后您可以从模型中排除“ zip”功能。