Блог

Я разработал модель twfe с фиктивной переменной, чтобы сделать различие в различии. Чтобы ввести новые управляющие переменные, мне нужно было бы ввести еще один манекен (муниципалитет моря/озера). Если я добавлю данные в код, а затем подведу итоги, я не увижу оценок этих коэффициентов, я подозреваю, что R отбрасывает переменную, но не знаю, как этого избежать.

Ниже кода:

 ## read the data library(readxl) data lt;- read_excel("~/link")  ## load packages library(dplyr) ## for wrangling the data library(plm) ## for the regression library(lmtest) ## for robust standard errors  ## coerce MUNICIPALITY and YEAR variables to factor data lt;- data %gt;%  mutate_at(c("MUNICIPALITY", "YEAR"), as.factor)  ## run the regression model_beds lt;- plm(formula = beds ~ DID   sea, data = data, effect = "twoways",  model = "within", index = c("MUNICIPALITY", "YEAR"))  ## summarize the results of the regression summary(model_beds)  

И результат:

 Call: plm(formula = hotelbeds ~ DID   sea, data = data, effect = "twoways",   model = "within", index = c("MUNICIPALITY", "YEAR"))  Balanced Panel: n = 78, T = 18, N = 1404  Residuals:  Min. 1st Qu. Median 3rd Qu. Max.  -1229.853 -64.005 -11.272 71.209 1728.606   Coefficients:  Estimate Std. Error t-value Pr(gt;|t|)  DID 55.471 23.342 2.3765 0.01762 * --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1  Total Sum of Squares: 73147000 Residual Sum of Squares: 72833000 R-Squared: 0.0042991 Adj. R-Squared: -0.068019 F-statistic: 5.64753 on 1 and 1308 DF, p-value: 0.017623  

Как я могу решить эту проблему? Как избежать того, чтобы R отбрасывал двоичную переменную?

Заранее спасибо.