### 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().
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.