如何将 R 用于密度 f 为混合物的对数似然函数？答案

【问题标题】：How to use R for a log-likelihood function with density f being a mixture?如何将 R 用于密度 f 为混合物的对数似然函数？
【发布时间】：2020-05-01 15:25:08
【问题描述】：

考虑从混合物 f(x) = πf1(x) + (1 − π)f2(x) 得到的密度 f，其中 f1 ∼ N(0, 1) 和 f2 ∼N(μ,1 )。为参数 (π, μ) 的对数似然函数编写一个 R 函数，该函数将来自 f 的观测值 x 向量作为输入

我知道要将 f(x) 合并到 log 函数中，您可以使用 dnorm 函数来获取 f1 和 f2。

f1 = dnorm(x, 0, 1)
f2 = dnorm(x, u, 1)

然后我们会这样做

LL=sum(log(f(x))

但总的来说，我不确定如何实际执行这个问题。任何帮助将不胜感激。谢谢！

【问题讨论】：

也许你可以从这里找到有用的信息stephens999.github.io/fiveMinuteStats/intro_to_em.html

标签： r statistics

【解决方案1】：

通过链接https://stephens999.github.io/fiveMinuteStats/intro_to_em.html，我猜你可能已经知道如何构造对数似然函数了。

以下是可能对您的目标有用的示例

f <- function(X,p,u) sum(log(p*dnorm(X,0,1) + (1-p)*dnorm(X,u,1)))

【讨论】：