R/my.logLAll.r

### This either returns a logL vector or sum of the vector.
###
### These functions are for all genes.

### Get the specific function according to the options.
get.my.logLAll <- function(model.Phi){
  if(!any(model.Phi[1] %in% .CF.CT$model.Phi)){
    stop("model.Phi is not found.")
  }
  ret <- eval(parse(text = paste("my.logLAll.",
                                 model.Phi[1], sep = "")))
  assign("my.logLAll", ret, envir = .cubfitsEnv)
  ret
} # End of get.my.logLAll().


### Function to calculate complete logL for
### (phi, b, sigmaWsq) given y, n, and phi.Obs
my.logLAll.lognormal <- function(phi, phi.Obs, y, n, b, p.Curr, 
    reu13.df = NULL){
  ### Dispatch.
  sigmaW <- p.Curr[1]

  ### Return.
  dlnorm(phi.Obs, log(phi), sigmaW, log = TRUE) +
  .cubfitsEnv$my.logdmultinomCodAllR(b, phi, y, n, reu13.df = reu13.df)
} # End of my.logLAll.lognormal().

### No need to changed from my.logLAll.lognormal since prior does not count.
my.logLAll.logmixture <- function(phi, phi.Obs, y, n, b, p.Curr,
    reu13.df = NULL){
  my.logLAll.lognormal(phi, phi.Obs, y, n, b, p.Curr, reu13.df = reu13.df)
} # End of my.logLAll.logmixture().

Try the cubfits package in your browser

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

cubfits documentation built on May 2, 2019, 4:08 a.m.