我已经使用因子变量拟合了 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"
有谁知道为什么这会有所作为?