Ich habe ein Problem Programmierung dieser Wahrscheinlichkeitsfunktion in R:Ich kann nicht meine likehood Funktion erhalten in R arbeiten
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.
Welche Art von Objekt ist 'x'? – yoland
x ist ein Vektor von 2 Parametern: x <- c (x_1, x_2) –
Schließen Sie die Abstimmung, weil keine Definition von 'h_j.1' oder' lambda' oder einem der anderen Vektoren vorliegt. –