因子变量中的级别名称很重要

Question

我使用因子变量studytreat拟合了lme模型。它在模型中出现在两个位置，作为固定效应，并且还在残差方差中模拟异质性，以便studytreat的每个级别都有自己的方差。这是模型：

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

我有一个问题，不能让我对lme充满信心！如果我将study中的水平标记为11,12,21,24（按此顺序），那么该功能正常。但是，如果我将数字切换为11, 21, 12, 42，那么当我调用summary(modelLME)$apVar时，我会收到错误消息：

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

有谁知道为什么这应该有所作为？