我正在分析一些面板数据。对于结果，我对不同国家/地区的效果感兴趣。由于实施效果的时间不同和治疗者的年龄不同，每个国家的假人都观察到了这种效果。

例如，我有两个国家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结合使用的变量？