R/llik.R

Defines functions llik

Documented in llik

llik = function(pars,tree){
  b = c(pars[1],(pars[1]-pars[2])/pars[3],pars[2])
  ct = sum(tree$wt)
  t = tree$wt[1:(length(tree$wt)-1)]
  lastn = tree$n[length(tree$n)]
  n = tree$n[1:(length(tree$n)-1)]
  E = tree$E
  sigma = n*(b[1]-b[2]*n + b[3]) #n-dimentional
  rho = pmax(b[1]*E-b[2]*n*E+b[3]*(1-E),0)
  #if(sum)
  l = -sum(-sigma*t+log(rho))
  if(min(b)<0){l = Inf}
  return(l)
}
franciscorichter/dmea02 documentation built on March 21, 2020, 3:46 a.m.