r - 因子変数の問題のレベルの名前

Question

因子変数を使用して lme モデルを適合させましたstudytreat。これはモデルの 2 つの場所で、固定効果として表示されます。また、残差分散の不均一性をモデル化して、各レベルが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"

なぜこれが違いを生むべきなのか、誰にも分かりますか？