我需要找到给定对数正态分布的倒数。 由于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。
提前致谢。
答案 0 :(得分:2)
通常,分位数分布是累积分布的倒数。这实际上意味着:
这意味着要找到对数正态分布的倒数,您可以使用
qlnorm()