2017-06-13 193 views
0

我需要找到給定的對數正態分佈的逆。 由於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有估計值。

在此先感謝。

+2

從[老公告板(https://開頭stat.ethz.ch/pipermail/r-help/2002-January/017625.html):「逆cdfs也稱爲分位數函數,在R中它們被命名爲 qxxxx; cdfs是pxxxx,隨機數字是rxxxx,密度是dxxxx,其中 xxxx命名分配系列。例如qnorm,pnorm,rnorm,dnorm「,所以'qlnorm'應該可以完成這項工作。 – p0bs

回答

2

一般來說,分位數分佈是累積分佈的倒數。這實際上意味着:

F(x) = p

x = F^-1(p)

這意味着找對數正態分佈的逆可以用

qlnorm()