我正在分析一些面板数据。对于结果,我对不同国家/地区的效果感兴趣。由于实施效果的时间不同和治疗者的年龄不同,每个国家的假人都观察到了这种效果。
例如,我有两个国家A
和B
,并为它们制作了两个假人(dummyA
和dummyB
)。
就目前而言,我只使用了数据的一个子组来分析一个国家:
mydata = data[data$country == A,]
然后我将假人退回A郡:
lm(y ~ dummyA , data = mydata)
对于国家B
我做了同样的事情
mydata = data[data$country == B,]
lm(y ~ dummyB , data = mydata)
我的问题是:是否有办法将这两个回归合并为一个?也许我可以将mydata = data[data$country == A,]
变成可以与dummyA
结合使用的变量?