Nothing
### 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().
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.