R/rf_llik.R

Defines functions rf_estot rf_lesk rf_lesk_int

rf_estot <- function(pars,qq)
{
   ph = pars[1];
   ii = pars[2];
   estot <- -ii * ph * log(1 - qq)/(ph - ii * log(1 - qq));
   return(estot);
}

rf_lesk <- function(pars,qq,k)
{
   logfun <- function(x) rf_lesk_int(x,pars,qq,k);
   lesk <- integral_peak(logfun);
   return(lesk)
}

rf_lesk_int <- function(x,pars,qq,k)
{
   ph = pars[1];
   ii = pars[2];
   y <- rep(0,length(x));
   for(cnt in 1:length(x))
   {
      y[cnt] <- sum(log(ii * x[cnt] + (0:(k - 1)))) + ii * x[cnt] * log(1 - qq) - ph * x[cnt];
   }
   y <- y + k * log(qq) - lgamma(k + 1) + 2 * log(ph);
   return(y);
}

Try the SADISA package in your browser

Any scripts or data that you put into this service are public.

SADISA documentation built on Oct. 30, 2019, 10:25 a.m.