R/dlindley.R

In matheuscastro43/finiteMix: Finite Mixture Models

dlindley <- function(x, beta, log = FALSE){
  if(length(x) == 1){
    if(length(beta) == 1){
      if(beta > 0){
        aux = 0
        if(x >= 0){
          aux = (x + 1)/(beta*(beta + 1))*exp(-x/beta)
        }
        if(!log){
          return(aux)
        }else{
          return(log(aux))
        }
      }else{
        stop("The parametric space must be respected.")
      }
    }else{
      k <- function(x, beta) dlindley(x, beta, log = log)
      mapply(k, x, beta)
    }
  }else{
    h = function(x){dlindley(x, beta, log)}
    sapply(x, h)
  }
}