rateCvToAlphaBeta function

[xiasijian]

Why could not find rateCvToAlphaBeta function?

[xiasijian]

#' Converts mutation rate and coefficient of variation (CV) parameters
#' to equivalent alpha and beta parameters typically used for beta-binomial.
#'
#' @param rate mutation rate
#' @param cv coefficient of variation for mutation rate
#' @return Param list containing alpha and beta
rateCvToAlphaBeta <- function(rate, cv) {
ab <- rate * (1-rate) / (cv*rate)^2 - 1
my.alpha <- rate * ab
my.beta <- (1-rate)*ab
return(list(alpha=my.alpha, beta=my.beta))
}

#############################

Analyze power and false positives

when using a beta-binomial model

#############################
smgBbdFullAnalysis <- function(mu, cv, Leff, signif.level, effect.size,
desired.power, samp.sizes){

find the power and numer of samples needed for a desired power

powerResult <- smgBbdRequiredSampleSize(desired.power, mu, cv, samp.sizes,
effect.size, signif.level, Leff)
bbd.samp.size.min <- powerResult$samp.size.min
bbd.samp.size.max <- powerResult$samp.size.max
power.result.bbd <- powerResult$power

get alpha and beta parameterization

for beta-binomial

params <- rateCvToAlphaBeta(mu, cv)
}