2017-01-15 7 views
1

Ich habe ein Problem Programmierung dieser Wahrscheinlichkeitsfunktion in R:Ich kann nicht meine likehood Funktion erhalten in R arbeiten

Likelihood function

ich es geschafft haben, die Funktionen h und lambda zu codieren. Es ist die Summe, die ich brauche, um herauszufinden, so weit mein Versuch, die folgende war:

log.lik<-function(x){ 
    sigma<-x[1]; tau<-x[2] 
    for(j in (1:6)){ 
    t<-rep(0,7) 
    t[j]<- -0.5((I-j)*log((1/p[j])*sigma)-log(1+lambda1(j,sigma,tau)*sum1[j])-h_j.1(j,sigma,tau) 
} 
return(-0.5*sum(t)) 
} 

Aber das scheint nicht zu funktionieren, da das Ergebnis zu seltsam ist.

+0

Welche Art von Objekt ist 'x'? – yoland

+0

x ist ein Vektor von 2 Parametern: x <- c (x_1, x_2) –

+0

Schließen Sie die Abstimmung, weil keine Definition von 'h_j.1' oder' lambda' oder einem der anderen Vektoren vorliegt. –

Antwort

1

Ohne etwas über Ihre eigentliche Problem zu wissen, vermute ich, dass Sie die Linie

t<-rep(0,7) 

aus der for-Schleife verschieben möchten.

+0

hmm, das hat tatsächlich ein besseres Ergebnis ergeben, zumindest meine hessische Matrix hat jetzt Werte. Nur nicht sicher, ob meine Parameter stimmen, aber danke! –

Verwandte Themen