ロジスティック分布関数の不動点を見つけて、さまざまなパラメーター値で不動点がどのように変化するかを判断しようとしています。コードは次のようになります。
nfxp.reps <- 0
err <- 10
p <- seq(0, 1, by = 0.0001)
pold <- p
gamma <- 6
k <- 3
while (err > 1E-12) {
nfxp.reps <- nfxp.reps + 1
if (nfxp.reps > 999) {
stop("The number of NFXP reps needs to be increased. \n")
}
pnew <- plogis(-k + gamma * pold)
err <- max(abs(pnew - pold))
pold <- pnew
}
上記のコードは、上記のパラメーターの選択で非常にうまく機能します: gamma と k - 3 つの固定点、2 つの安定点、1 つの不安定点 (p=0.5) を見つけます。ただし、上記のパラメーターを不均等に変更すると、中央の固定点が 0.5 より上または下にある場合、次のようになります。
gamma<-7
k<-3
ループは、p=0.3225 である中央の固定点を見つけることができません (ガンマ=7 の場合、k=3)