r - 具有非结构化相关性的 gls 错误

Question

我迫切需要帮助。我有一个连续的响应变量、一个固定的治疗效果（具有三个水平的因子变量）、一个固定的时间效果（具有九个固定时间点的因子变量）和一个随机效果的受试者效果（我有 34 个受试者，基本上每个受试者在每个时间点尝试每个治疗）。我正在尝试使用以下命令拟合具有非结构化协方差的混合模型，但我不断收到以下错误。此外，R 需要 EVER 才能返回错误。问题似乎出在 corSymm 方面。我以前从未这样做过，R 对我来说仍然很难，所以任何帮助都将不胜感激。谢谢！

cs1 <- corSymm(form = ~ 1 | Subject)                                                                    
cs1 <- Initialize(cs1, data = Mydata)                                    
model <- gls(Response ~ Time+Treatment+Treatment*Time, correlation = cs1,                                                           data = Mydata) 
anova(model)

gls 中的错误（响应 ~ 时间 + 治疗 + 治疗 * 时间，相关性 = cs1，：在没有收敛的情况下达到函数评估极限 (9)

标签： rnlme