适合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的帮助。