R/dCFGeo.R

In ComRiskModel: Fitting of Complementary Risk Models

#' @export
dCFGeo<-function (x, a, b, lambda, log = FALSE)
{
G=(1-((1+(x/a)^b))^(-1))
g=((1+(x/a)^b))^(-2) *b*x^(b-1)*a^(-b)


  pdf <- x
  pdf[log == FALSE] <-(1-lambda)*(g)/(1-lambda*G)^2
  pdf[log == TRUE] <-(1-lambda)*log(g)-2*log(1-lambda*G)
  return(pdf)
}