Есть ли способ указать значения в функции PLM?

Question

В настоящее время я работаю над проектом, для которого мне нужно отличить основную переменную. Я строю отдельную модель с фиксированными эффектами с помощью plm, и мне нужно «разрезать пополам» переменную. К сожалению, я не могу предоставить репекс, но часть кода.

low_troop <- plm(POLSTAB ~ US +  RDURAB + DEMOC + POP,  data = US, index = c('COUNTRY'), model = 'within', effect = "individual")
high_troop <- plm(POLSTAB ~ US+  RDURAB + DEMOC + POP,  data = US, index = c('COUNTRY'), model = 'within', effect = "individual"

По моим данным, США - это количество американских солдат за границей. Я хотел бы иметь один plm для всех данных, где значение для США> 100 и один для значения США <100. Возможно ли это без создания нового набора данных? </p>

Большое спасибо всем за помощь.

С уважением,

Лоуренс

Answer 1

Просто подмножество данных в аргументе data=; пример:

library(plm)
data(Grunfeld)

wi.lo <- plm(inv ~ value + capital, data=Grunfeld[Grunfeld$value >= 1000,], 
          index="firm", model="within", effect="individual")
wi.hi <- plm(inv ~ value + capital, data=Grunfeld[Grunfeld$value < 1000,], 
          index="firm", model="within", effect="individual")

summary(wi.lo)$coef
#          Estimate Std. Error   t-value     Pr(>|t|)
# value   0.1132322 0.02024847  5.592138 6.934243e-07
# capital 0.3455339 0.03022447 11.432257 2.868412e-16

summary(wi.hi)$coef
#           Estimate Std. Error  t-value     Pr(>|t|)
# value   0.09465974 0.01115671 8.484558 4.579648e-14
# capital 0.09808216 0.01249071 7.852408 1.431317e-12