【发布时间】:2017-11-15 11:22:53
【问题描述】:
我需要找到给定对数正态分布的倒数。 由于 R 中没有用于逆对数正态的内置函数,因此我需要自己设计。
我有一个随机变量“x”的对数正态分布
f_lambda <- function(x,mu,sig) {dlnorm(x, meanlog = mu, sdlog = sig,log=FALSE)}
维基百科上说
G(y) = 1- F(1/y)
其中 G(Y)n 是 F(X) 和 X= 1/Y 的逆分布。
但是,我对如何在 r 中编码 F(1/y) 以及使用什么来定义该分布 - mu 或 1/mu 感到困惑。
我估计了 F(x) 的 mu 和 sigma。
提前致谢。
【问题讨论】:
-
来自an old bulletin board:“逆cdfs也称为分位数函数,在R中它们被命名为qxxxx;cdfs是pxxxx,随机数是rxxxx,密度是dxxxx,其中xxxx命名分布族。例如,qnorm、pnorm、rnorm、dnorm”。所以
qlnorm应该可以完成这项工作。
标签: r statistics distribution