MLE und logit Modell < Matlab < Mathe-Software < Mathe < Vorhilfe
|
Hallo zusammen,
In meiner vorweihnachtlichen Verzweiflung wende ich mich mit folgendem Problem an euch. Ich habe einen Graphen bestehend aus 50 Werten den ich mit einem logit Modell glätten möchte. Hier das Modell:
Das logistische Modell stellt sich wie folgt da
[mm] q(x)=1-exp[-(\bruch{beta*exp(b*x)}{1+a*exp(b*x) }+c) [/mm] ]
Der Parametervektor θ wird so gewählt, dass die folgende Log-Likelihood Funktion maximal wird
ln [mm] L(Teta)=\summe_{x=65}^{120} [/mm] [T(x) * ln(q(x))+(L(x)-T(x)) * ln(1-q(x))]
[mm] T_x ~Binomial(L_x q_x,L_x q_x p_x)
[/mm]
[mm] L_x [/mm] die Lebenden des Alters x des Modellbestandes
[mm] T_x [/mm] die Zufallsvariable der im Alter x Gestorbenen im Modellbestand
q Überlebenswahrscheinlichkeit
Leider habe ich bei der Umsetzung mit dem Befehl
[...] = mle(data,'distribution',dist)
große Probleme. Laut einem Paper ergeben sich alle 4 Parameter (alpha, beta, b und c) aus der Maximum Likelihood estimation.
Habe ich das richtig verstanden?
Für dist gebe ich binomial ein.Wofür steht aber 'distribution' (anscheinend number of tails)?
Bei meinen Versuchen gebe ich die Überlebenswahrscheinlichkeit q ein und erhalte nur 2 Parameter zurück (was ich nachvollziehen kann), ich brauche aber 4?
Ich bin für jeden Kommentar dankbar.
VIELEN DANK und viele Grüße
ich
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 10:20 Di 30.12.2008 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
