我迫切需要帮助。我有一个连续的Response变量，一个固定的Treatment效果（具有三个级别的因子变量），一个固定的Time效应（具有9个固定时间点的因子变量）和一个Subject效应，它是随机效应（我有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）