パッケージで使用gnls
して、ロジスティック成長モデルをRに適合させようとしています。nlme
私は以前にモデルをうまく適合させました:
mod1 <- gnls(Weight ~ I(A/(1+exp(b + v0*Age + v1*Sum.T))),
data = df,
start = c(A= 13.157132, b= 3, v0= 0.16, v1= -0.0059),
na.action=na.omit)
ただし、モデルに適合しないように制約したいb
ので、2 番目のモデルを適合させてみました。
mod2 <- gnls(Weight ~ I(A/(1+exp(log((A/1.022)-1) + v0*Age + v1*Sum.T))),
data = df,
start = c(A= 13.157132, v0= 0.16, v1= -0.0059),
na.action=na.omit)
このモデルは次のエラーを返しました:
Error in gnls(Weight ~ A/(1 + exp(log((A/1.022) - 1) + v0 * Age + :
approximate covariance matrix for parameter estimates not of full rank
Warning messages:
1: In log((A/1.022) - 1) : NaNs produced
エラーを検索すると、問題がモデルの対称性によって引き起こされていることが示唆され、特定の問題の解決策には、さまざまなパラメーターを使用して式を適応させることが含まれます。残念ながら、私の統計的知識は、a) 問題を完全に理解するか、b) 式を自分で適応させるのに十分ではありません。
警告メッセージ (全部で 15 ありましたが、すべて同じです) については、モデルのこのセクションが単独で機能するため (例の数字を使用)、なぜそれらが発生するのかわかりません。
これらのクエリのいずれかに関するヘルプをいただければ幸いです。