1

私は lme4 を学ぼうとしていますが、20 年以上前に、SAS (ノースカロライナ州立大学) の本拠地で、統計学の教育訓練が伝統的な ANOVA 設計であったことを前もって認めます。私は現在、R で約 4 ～ 5 年の経験がありますが、ほとんどの場合、すべての因子が固定効果として扱われる従来の一般線形モデルを使用しています。

ここで、淡水ムール貝を使った実験の von Bertalanffy 成長モデルのパラメーターを推定したいと思います。既知個体(ID)を対象に年1回、4年間殻長を測定したため、繰り返し測定設計です。最終的には、主な治療因子と分割プロット因子を追加したいと考えています。しかし、現時点では単純化されたモデルを実行できません。私が使用しているRコードは次のとおりです。

VBGF <- function(Linf, k, age, age.0) {
  Linf * (1-exp(-k*((age-age.0))))
}
VBGF <- deriv(L ~ VBGF, namevec=c("Lin", "k", "age.0"), function.arg=VBGF)
mod1 <- nlmer(L ~ VBGF(Linf,k,AGE,age.0) ~ (Linf+k+age.0|ID), 
              mussels,
              start=c(50,0.1,-11),
              verbose=TRUE)

次のエラーが表示されますが、それが何を意味するのかについての参照を見つけることができませんでした。

Error: s > 0 is not TRUE

見落としているのはおそらく単純なことなので、事前に無知であることをお詫びします。

4

1 に答える 1

2

関数でエラーが発生する行nlmerは

stopifnot(length(start$fixef) > 0, s > 0, inherits(data, 
    "data.frame"), nrow(data) > 1)

sコードを見ると、リストの名前の数が多い印象を受けますstart。あなたの場合、vector名前なしで提供したため、おそらくエラーが発生しました。

ヘルプページには

start: a named list of starting values for the parameters in the
      model. (..)

したがって、彼らは名前が必要であるという事実について非常に具体的です。例では、パラメータにも名前が付けられていることがわかります。

## nonlinear mixed models
(nm1 <- nlmer(circumference ~ SSlogis(age, Asym, xmid, scal) ~ Asym|Tree,
             Orange, start = c(Asym = 200, xmid = 725, scal = 350)))
于 2012-06-02T18:13:12.317 に答える