在R中拟合重复测量[关闭]

Question

在 R 中拟合重复测量，收敛问题。我有以下拟合，它是许多数据集之一，它不会收敛。我做其他的套装。这个数据集和模型在 SAS 中工作......我可以为在 R 中进行这项工作找到一些方向吗？要看的东西（矩阵、选项设置、关于这个主题的 r/splus 参考...）。

fit.gls <- gls(resp~grpnum+grpnum/day,data=long, corr=cormat,na.action=na.omit)
Error in glsEstimate(object, control = control) : 
  computed "gls" fit is singular, rank 62

我已阅读以下内容，但仍在努力解决... Converting Repeated Measures mixed model formula from SAS to R

Answer 1

问题在于数据。 gls 需要反转矩阵才能工作（请参阅Wikipedia 了解估计协变量的公式）。对于您的特定数据集，该矩阵是不可逆的。

您可以使用 control 参数允许奇异值：

fit.gls <- gls(resp~grpnum+grpnum/day,data=long, corr=cormat,na.action=na.omit, control = list(singular.ok = TRUE))

小心这个，因为你可能会得到不好的结果！之后总是检查模型的拟合度。

查看gls 和glsConrol 的帮助以了解有关选项的更多详细信息。