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

您可以使用 uniroot :

f <- function(x) dnorm(x, m=14, sd=0.9) * .6 - dnorm(x, m=16, sd=0.6) * .4

uniroot(f, interval=c(12, 16))

$root
[1] 15.19999

$f.root
[1] 2.557858e-06

$iter
[1] 5

$estim.prec
[1] 6.103516e-05

ETA一些说明:
uniroot是一个单变量根查找器，即给定一个函数 f一个变量 x ，它找到 x 的值求解方程 f(x) = 0 .

要使用它，您需要提供函数 f ，以及假设解值所在的区间。在这种情况下， f只是两个密度之间的差异；它们相交的点将是 f为零。在这个例子中，我通过绘制一个图并看到它们在 x=15 附近相交得到了间隔 (12, 16)。