R/likeli.ald.R

In dirttee: Distributional Regression for Time to Event Data

likeli.ald <-
function(par,tau,ycens,delta,X,K)
  { -sum(delta * log(dald(ycens,X%*%par[-1],abs(par[1])+0.0001,tau)) + (1-delta) * log(1-pald(ycens,X%*%par[-1],abs(par[1])+0.0001,tau))) + t(par[-1])%*%K%*%par[-1]}