Geschlossene Form der Funktion in Bezug den natürlichen Parameter für die logarithmische Reihe Verteilung?

7

Während ich hier eine andere Frage beantwortete , erwähnte ich die logarithmische Reihenverteilung als mögliches Modell für Arten pro Gattung.

Als ich mir die PMF ansah, als ich antwortete, stellte ich fest, dass es sich um eine exponentielle Familie handelte (eine Tatsache, die mir vorher nicht bewusst war).

Die Parametrisierung des pmf auf der Wikipedia-Seite (die bis zur Änderung der Symbole die einzige Form ist, in der ich sie gesehen habe) ist von der Form

f(y)=c(p).py/y,0<p<1;y=1,2,3,...

Dabei ist die Normalisierungskonstante.c(p)=1log(11p)=1log(1p)

Geben Sie hier die Bildbeschreibung ein

Wenn wir es in "Standardform" für die Exponentialfamilie untersuchen, sehen wir den natürlichen Parameter (und auch, dass die ausreichende Statistik , so dass die MLE für eine Konstante sein sollte Funktion des Stichprobenmittelwerts, und in diesem Fall würden wir erwarten, dass MLE der Methode der Momente entspricht, was sich als der Fall herausstellt.θ=log(p)T(x)=xp

[Es sieht also so aus, als hätten wir ein glm mit einem Modell , dann sollte für ausreichen .]θ=logp=XβXyβ

Es ist einfach genug, als Funktion von abzuleiten (die Reihe ist trivial zu summieren) oder (in Bezug auf die Exponentialfamilie betrachtet) zu nehmen und so oder so in Begriffen zu schreiben von :μpb(θ)μθ

μ=c(p)p1p=1log(1eθ)eθ1eθ.

Es ist sogar durchaus möglich, die Varianzfunktion zu berechnen .V(μ)

Mir ist jedoch nicht sofort klar, wie diese mittlere Funktion invertiert werden kann, um (oder äquivalent ) als Funktion von zu schreiben (was praktisch wäre, da es sich um die Verknüpfungsfunktion handelt). Es mag Müdigkeit sein oder es könnte tatsächlich schwierig sein, aber ich konnte nicht sehen, wie ich damit irgendwohin komme ... also zur Frage:θpμ

Kann in geschlossener Form geschrieben werden?θ(μ)

( wird genauso gut tun)p(μ)

Ich möchte nicht übermäßig pedantisch sein, was in geschlossener Form enthalten ist. Es ist nützlich, diese Dinge in Form bekannter Funktionen mit bekannten Eigenschaften schreiben zu können, die durch Aufrufen einer Funktion relativ einfach ausgewertet werden können. Ausdrücke in Form von mehr oder weniger "Standard" -Funktionen, die üblicherweise in Paketen verfügbar sind, wie Gammafunktionen oder sogar Lambert-W-Funktionen (oder sogar die verschiedenen Formen von hypergeometrischen Funktionen, denke ich) könnten zählen - ich bin mehr interessiert daran, welche nützlichen Dinge wir über sagen können und wie leicht wir es bewerten können (über den offensichtlichen Schritt hinaus, die obige Gleichung für -als-Funktion-von- und nur nach zu lösen nach Bedarf).θ(μ)μθθ

Glen_b -State Monica
quelle

Antworten:

2

Sie haben die Antwort in Ihrer Frage: Lambert-W.

In R können Sie verwenden

# Load library
  library(LambertW)

# Define function to obtain p from mu
  p <- function(mu) {1 - exp((1/mu) + W(-1/(mu*exp(1/mu)),-1))}

# Show (visually) that the function provides the correct values
  mu <- 1 + 10*c(0:100)/100
  plot(p(mu), mu, xlab="p", ylab="mu", las=1)
  pp <- c(1:99)/100
  lines(pp, -pp/((1-pp)*log(1-pp)))

erhalten

R Lambert W.

In Mathematica können Sie verwenden

data = Table[{1 - Exp[(1/mu) + 
       ProductLog[-1, -1/(mu Exp[1/mu])]], mu}, {mu, 1.01, 5, 0.01}];
cp = -1/Log[1 - p];
ListPlot[{data, Table[{p, cp p/(1 - p)}, {p, 0.1, 0.93, 0.01}]}, 
 PlotStyle -> {{PointSize[0.03]}, {PointSize[0.01]}}]  

(Mathematica wurde verwendet, um die allgemeine Form der Beziehung zur Bestimmung von aus zu bestimmen .)pμ

JimB
quelle
Vielen Dank. Nachdem ich einige Zeit damit verbracht hatte, fand ich schließlich heraus, wie ich es in eine Form bringen konnte, die Wolfram Alpha lösen würde, bevor das Zeitlimit überschritten wurde (dies war vor ein paar Tagen), aber ich dachte mir, dass ich das Kopfgeld laufen lassen würde, damit jemand das bekommen könnte Vorteil. Ihre Antwort ist viel schöner als alles, was ich gepostet habe, also verdient sie das Kopfgeld.
Glen_b -State Monica