Displaying 1 result from an estimated 1 matches for "term_log".
Did you mean:
term_arg
2000 Feb 04
2
constrained optimiser doesn't obey constraints (PR#411)
...hird.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 minimize -log-likelihood which is
## equivalent to maximizing the log-likelihood
return(-log.term)
}
init<-c(0.001, 0.001, 0.001, 0.220)
"sim.ns" <-
c(17, 17, 15, 15, 16, 1...