使用稀疏矩阵时，使用矩阵响应得到最小二乘拟合的最佳方法是什么？我的回归矩阵和我的回答矩阵都很稀疏。在我的特定情况下，它们是多值因子，这些因素与计数变量相互作用。
来自slm.fit包的 SparseM是一个选项，但似乎没有简单的方法来构建SparseM模型矩阵。所有示例都使用密集矩阵创建它们，然后转换或从文件中读取它。这意味着我可以使用Matrix构建模型矩阵，然后将其写入文件，然后使用SparseM读取它。这不是一个很好的解决方案。
speedlm，sparse.lm.fit和glm4似乎都无法处理Matrix响应。我可以简单地为响应矩阵的每一列运行speedlm一次，但效率很低。