0

我已经使用因子变量拟合了 lme 模型studytreat。它出现在模型中的两个地方，作为一个固定效应，也用于对残差方差中的异质性进行建模，以便每个级别studytreat都有自己的方差。这是模型：

modelLME <- lme(logkp ~ studytreat, 
                random = ~ -1 + STUDY | ptno, 
                weights = varIdent(form = ~1 | studytreat))

我有一个问题让我对 lme 没有信心！如果我将 studytreat 中的级别标记为 11、12、21、24（按此顺序），则该功能可以正常工作。但是，如果我切换数字以获取，11, 21, 12, 42那么当我打电话时，summary(modelLME)$apVar我会收到错误消息：

[1] "Non-positive definite approximate variance-covariance"

有谁知道为什么这会有所作为？

4

0 回答 0