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