我有一个大型数据集，我一直在使用线性混合模型进行分析（约600000次观察）。通常，我估计最小二乘平均值，并使用 lsmeans 包中的 lsmeans 进行成对比较。我还使用自由度的估计来测试每个对比度的显着性，使用 t 统计量，因为我实际上有很多次数。但是，我无法在我的数据集上应用 lsmeans 。下面我展示了一个示例，其中我模拟了不同大小的数据集。
O(N)我还尝试了 lmerTest 包并使用了 difflsmeans 函数，它也崩溃了。我可以使用 multcomp 包中的一般线性假设检验 - glht 来获得平均估计值 - 但是此过程不会估计特定对比的自由度。您是否知道适用于大型数据集的这些程序的替代方法？
非常感谢你。