这是我第一次访问你的博客，如果可能的话，我需要你的帮助或建议来解决我将要描述的问题。
我和我的团队使用 R 估计了以下固定效应 OLS 估计：
iv_tsls_fe <- feols(mort_u65_pop ~ year + datacoverage + GDP + unemployment_rate + at_risk_poverty | clust2 | pm25 ~ rain_yr_mm + temp_c + wdsp, data = df2, vcov = vcov_cluster("clust2"))
我们使用“feols”函数估计固定效应工具变量（IV）模型。因变量是“mort_u65_pop”，自变量包括“year”、“datacoverage”、“GDP”、“unemployment_rate”和“at_risk_poverty”。工具变量是“rain_yr_mm”、“temp_c”和“wdsp”。我们使用“vcov_cluster”来计算固定对象的聚类 VCOV。在我们的例子中，“clust2”是数据集“data = df2”内的变量。
然后我们计算固定对象的拟合统计量如下：
fitstat(iv_tsls_fe, "ivwald")
fitstat(iv_tsls_fe, "cd")
fitstat(iv_tsls_fe, "kpr")
虽然 Wald 和 The Cragg-Donald 测试了弱仪器工作，但对于弱仪器的 Kleibergen-Paap 测试，我遇到了以下错误：
Error in vcov + (-1)^(i + 1) * vcovClust(index, bread, scores, adj = ssc$cluster.adj && :
non-conformable arrays
不幸的是，我们无法找出问题所在，因为前两个测试都有效。
有人对如何处理这个问题有什么建议吗？我将非常感激。
我希望我的演示已经很清楚了，并且可以随时提供更多详细信息，以防我遗漏了什么。
非常感谢， 洛雷娜
我认为这里的问题可能是由于集群错误与此包中的 Kleibergen-Paap 测试配合得不好。我也遇到过这个错误，并注意到它只在我使用集群错误时出现，如果它与您的情况相同，我会很感兴趣。
不幸的是，我不知道修复程序中有任何直接解决此问题的方法。我当前的解决方案是利用 ivDiag 进行 IV 诊断。虽然在我的例子中，所涉及的固定效果的数量相对于使用fixt时可能拥有的数量而言处于较低端，但ivDiag的后端是来自lfe的，因此即使对于大量的固定效果也应该没问题。最后要注意的是，ivDiag 实际上并不报告 KP 测试，但它确实报告了 Montiel Olea & Pflueger 有效 f 统计量，该统计量在刚刚确定的情况下简化为 KP 测试，并且似乎优于 kp 测试至少由本文的作者“https://scholar.harvard.edu/files/wirev_092218-_ Corrected_0.pdf”。
很抱歉没有提供更多帮助，但我希望这能有所帮助。
我认为你可以使用：
waldtest(iv_tsls_fe, vcov = vcovHC)