R/likeli.asynorm.R

In dirttee: Distributional Regression for Time to Event Data

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