R/mu.R

In touchard: Touchard Model and Regression

 mu0 <- function(lambda, delta, N=NULL, eps=sqrt(.Machine$double.eps)){
        #stopifnot (lambda >= 0)
        if(lambda < 0) stop("touchard:::mu cannot compute with negative lambda")
        t0        <- tau(lambda,delta,N,eps)
        m0        <- tau(lambda,delta+1,N,eps)/t0
        return(m0 - 1)
}
   #val <- sapply(lambda, mu0, delta=delta, N=N, eps=eps) 
  
mu <- Vectorize(mu0, vectorize.args = c("lambda", "delta"))