2017-06-13 7 views
0

Ich muss die Umkehrung einer gegebenen lognormalen Verteilung finden. Da es keine eingebaute Funktion in R für inverse Lognormal gibt, muss ich meine eigene entwerfen.Invers der logarithmischen Normalverteilung

Ich habe diese Lognormalverteilung für eine Zufallsvariable 'x'

f_lambda <- function(x,mu,sig) {dlnorm(x, meanlog = mu, sdlog = sig,log=FALSE)} 

auf Wikipedia sagt

G(y) = 1- F(1/y) 

wo G (Y) n, die inverse Verteilung F (X) und X = 1/Y.

Aber ich bin verwirrt, wie man F (1/y) in r kodiert und was man benutzt, um diese Verteilung zu definieren - mu oder 1/mu.

Ich habe Schätzungen von mu und Sigma für F (x).

Vielen Dank im Voraus.

+2

Von [eine alte Bulletin Board] (https verwenden: // stat.ethz.ch/pipermail/r-help/2002-Januar/017625.html): "Inverse cdfs werden auch Quantilfunktionen genannt, und in R heißen sie qxxxx; cdfs sind pxxxx, Zufallszahlen sind rxxxx, Dichte ist dxxxx, wobei xxxx benennt die Verteilerfamilie. ZB qnorm, pnorm, rnorm, dnorm ". Also sollte' qlnorm' wahrscheinlich die Aufgabe erfüllen. – p0bs

Antwort

2

Im Allgemeinen ist die Quantilverteilung die Umkehrung einer kumulativen Verteilung. Das ist wirklich bedeutet:

F(x) = p

x = F^-1(p)

was bedeutet, dass die Umkehrung der Lognormalverteilung finden Sie

qlnorm()