我想将R中的cluster.bs.glm包中的clusterSEs函数（即成对的群集自举p值配对）用于多个glm对象。
特别是，我想使用glm.nb对象来实现它，该对象是使用glm.nb包中的MASS函数估计的，但是遇到了错误。下面显示了我正在使用的代码。
cluster.bs.glm(neg_bin, dat = DATASET, cluster = ~ firm, boot.reps = 1000)
Error in `rownames<-`(`*tmp*`, value = c("(Intercept)", "treat", "post", :
length of 'dimnames' [1] not equal to array extent
我怀疑是theta的估计值导致了这一点。我最终可以使用glm来工作，但是我必须先设置theta的值（不确定是否合适）。如果需要，我将尝试获取一个可行的示例。有什么想法吗？也许删除theta，以便将置信区间打印到控制台上？