具有聚类误差的回归模型中固定效应的联合显着性检验

Question

我正在尝试测试在城市级别聚集的回归模型中的固定效果（在下面的代码中用“ lg”表示）是否显着，F检验或卡方检验是否适合尝试。 R是我使用的语言。

问题是我尝试过的所有可能的测试代码都无法在具有“聚集错误”的模型下进行F或卡方检验。

我尝试了以下测试：

<td onclick="displayComment($row['num'])">Column 3</td>

model <- miceadds::lm.cluster(y ~ x1 + x2 + factor(year) + factor(lg) + factor(cityid) , 
                              cluster="cityid", data = data)
linearHypothesis(model, c("factor(lg)=0"))

或

model1 <- miceadds::lm.cluster(y ~ x1 + x2 + factor(year) + factor(lg) + factor(cityid), 
                               cluster="cityid", data = data)
model2 <- miceadds::lm.cluster(y ~ x1 + x2 + factor(year) + factor(cityid) , 
                               cluster="cityid", data = data)
anova(model1, model2, test="Chisq")

我收到以下错误：

  

没有适用于“ lm.cluster”类对象的'anova'方法

我想知道是否有任何方法可以让我测试lg固定效果是否共同显着？也许我可能不得不自己计算卡方统计量，但是我不熟悉如何使用自己的代码进行统计。