R/pasynorm.R

In dirttee: Distributional Regression for Time to Event Data

Defines functions pasynorm

Documented in pasynorm

pasynorm <-
function(q,mu=0,sigma=1,tau=0.5)
{
  if(any(sigma <= 0)) stop("standard deviation must be strictly > 0.")
  if(any(tau <= 0 | tau >= 1)) stop("tau must be between 0 and 1.")
  mapply(function(z,mu)integrate(dasynorm,-Inf,z,mu=mu,sigma=sigma,tau=tau)$value,q,mu)
}

Try the dirttee package in your browser