您如何计算 MLE 转换的标准误?

机器算法验证 最大似然 标准错误 增量法
2022-03-21 02:58:09

我需要推断一个正参数p. 为了适应我重新参数化的积极性p=exp(q). 使用 MLE 例程,我计算了点估计和 seq. MLE 的不变性直接给了我一个点估计p,但我不确定如何计算 sep. 提前感谢您的任何建议或参考。

2个回答

Delta 方法用于此目的在一些标准正则假设下,我们知道 MLE,θ^为了θ大约(即渐近地)分布为

θ^N(θ,I1(θ))

在哪里I1(θ)是整个样本的Fisher 信息的倒数,评估为θN(μ,σ2)表示具有均值的正态分布μ和方差σ2. MLE函数不变性表示g(θ), 在哪里g是一些已知的函数,是g(θ^)(正如你所指出的)并且具有近似分布

g(θ^)N(g(θ),I1(θ)[g(θ)]2)

您可以在其中插入未知量的一致估计器(即插入θ^在哪里θ出现在方差中)。我会假设您拥有的标准错误是基于 Fisher 信息(因为您有 MLE)。将标准误差表示为s. 那么标准误eθ^,如在您的示例中,是

s2e2θ^

我可能会向后解释你,实际上你有 MLE 的方差θ并且想要 MLE 的方差log(θ)在这种情况下,标准将是

s2/θ^2

Macro 给出了如何通过 delta 方法转换标准错误的正确答案。尽管 OP 专门要求提供标准误差,但我怀疑目标是为p. 除了计算估计的标准误p^您可以直接转换置信区间,[q1,q2], 在里面q-参数化到置信区间[exp(q1),exp(q2)]在里面p-参数化。这是完全有效的,甚至可能是一个更好的主意,具体取决于用于证明基于标准误差的置信区间的正态近似在q-参数化与p-参数化。此外,直接转换的置信区间将满足积极性约束。