Displaying 1 result from an estimated 1 matches for "part_eta".
Did you mean:
part_end
2000 Feb 04
2
constrained optimiser doesn't obey constraints (PR#411)
...0:(no.inform[i]-loss[i]-1)) {
second.prod_second.prod*(1-pi1+r*omega1)
};} else {second.prod_1;}
third.prod_1
for (r in 0:(no.inform[i]-1)){
third.prod_third.prod*(1+r*omega1)
}
common.term<-choose(no.inform[i],loss[i])
first.part_eta*common.term*first.prod*second.prod/third.prod
whole.term_eta*common.term*first.prod*second.prod/third.prod+(1-eta)*common.term*p0^loss[i]*(1-p0)^(no.inform[i]-loss[i])
log.term_log.term+log(whole.term)
}
## log.term is the log-likelihood.
## For nlminb we need to m...