带vcovCL的POLR ordered logit上的聚集标准错误给出；“不一致的论点”；错误

我试图使用MASS中的POLR函数获得有序logit回归的聚集标准误差。回归有大约2900个固定效应（美国各州和年份）

当我在没有固定效果的情况下运行回归时，我能够毫无问题地使用sandwich的vcovCL函数，（）但是，当我包括固定效果时，我得到以下错误：

Error in xmat %*% x$coefficients : non-conformable arguments

我会发布数据，但是数据集非常大，回归需要大约2天的时间来运行所有的固定效应

是否有人遇到过类似的问题，您是否提出了解决方案

编辑：根据@ben bolker在下面的评论，我尝试了1000行随机样本的回归。某些固定效果似乎会产生一个方形矩阵，产生以下错误：

   1: glm.fit: fitted probabilities numerically 0 or 1 occurred 
   2: design appears to be rank-deficient, so dropping some coefs

当我尝试使用vcovCL函数对标准错误进行聚类时，我得到了相同的“非一致参数”错误。所以我怀疑，由于一些系数是NA，矩阵vcovCL与POLR对象的大小不同，这就是导致错误的原因

---

有没有办法去除这些缺失的固定效果？无论如何，我不会在回归输出中显示它们。

你能用更小的数据集重现问题吗，即一些固定效应，但只有少数（不是很多观察值）？好主意。我采样了1000行，发现了同样的问题。如果我发布被截断的数据会有帮助吗？