0

在 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

我已阅读以下内容并仍在尝试通过它... 将重复测量混合模型公式从 SAS 转换为 R

4

1 回答 1

4

问题是数据。 gls需要反转一个矩阵才能工作（参见维基百科的公式来估计协变量）。对于您的特定数据集，该矩阵是不可逆的。

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

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

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

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

于 2013-11-01T19:34:44.250 回答