r - 2 つの法線曲線の交点

Question

これは基本的な質問だと思いますが、Rでこれを計算する方法を見つけることができないようです:

たとえば、次のパラメータを持つ2つ以上の正規分布(ヒストグラムに適合)の交点(x値が必要):

d=data.frame(mod=c(1,2),mean=c(14,16),sd=c(0.9,0.6),prop=c(0.6,0.4))

私の 2 つの曲線の平均と標準偏差を使用して、分布に対する各 mod の寄与率をプロップします。

Answer 1

申し訳ありませんが、受け入れられた答えは良くありません。参照: matlab の 2 つの曲線の交点

次のような関数を使用して、両方のルートを取得できます。

intersect <- function(m1, s1, m2, s2, prop1, prop2){

B <- (m1/s1^2 - m2/s2^2)
A <- 0.5*(1/s2^2 - 1/s1^2)
C <- 0.5*(m2^2/s2^2 - m1^2/s1^2) - log((s1/s2)*(prop2/prop1))

(-B + c(1,-1)*sqrt(B^2 - 4*A*C))/(2*A)
}

あなたの場合：

> intersect(14,0.9,16,0.6,0.6,0.4)
[1] 20.0 15.2