如果您运行三个不同的回归,则无法正式测试系数的差异。它可能会提供更多信息:
- 标准化分数(标准正常可能是一个好的开始,这里有更多信息https://datascience.stackexchange.com/questions/1240/methods-for-standardizing-normalizing-different-rank-scales)
- 堆叠您的数据,以便每个观察结果显示 3 次。分数将在同一个变量中。
- 创建 3 个新变量,您可以在其中将分数与分数机构的虚拟指标进行交互。您的数据将如下所示:
rbind(
data.frame(iso3 = "USA", year = 2001, gdp = 13, score_sp = 1, score_moody = 0, score_sfse = 0 ),
data.frame(iso3 = "USA", year = 2001, gdp = 13, score_sp = 0, score_moody = 2, score_sfse = 0 ),
data.frame(iso3 = "USA", year = 2001, gdp = 13, score_sp = 0, score_moody = 0, score_sfse = 3 )
)
iso3 year gdp score_sp score_moody score_sfse
1 USA 2001 13 1 0 0
2 USA 2001 13 0 2 0
3 USA 2001 13 0 0 3
- 估计你的模型:
plm(CDS ~ GDP+Inflation+...+ score_sp + score_moody + score_ftse, data, model="within")
现在您可以简单地将系数与 t 检验进行比较。
如果您想查看使用不同分数时“控制”变量的系数如何变化,那么您的数据将如下所示:
iso3 year score_sp score_moody score_sfse gdp_sp gdp_moody gdp_sfse
1 USA 2001 1 0 0 13 0 0
2 USA 2001 0 2 0 0 13 0
3 USA 2001 0 0 3 0 0 13
现在您可以使用 t-test 来检查使用某个特定分数时 gdp 的系数是否更大(如果您对此感兴趣)。