###########################################################################
## summary of the evaluation of simulation result.
## The input is an object from evalSim function.
###########################################################################
#' @name printEvalRes
#' @aliases printEvalRes
#' @title Summary table of power assessment
#' @description This function takes as input a result object from \code{\link{evaluateSim}} and prints out a table to summarize important error-rates-related quantities. The results are marginalized, meaning that they are averaged quantities over all strata and simulations. This provides a quick view of the marginal results per sample size.
#' @usage printEvalRes(evalRes)
#' @param evalRes The result object from \code{\link{evaluateSim}}.
#' @return A matrix of results per sample size considered (rows). Columns include sample size, specified nomial type I control value (for FDR or p-values), actual error rate, marginal TPR, averaged number of true and false discoveries, and false discovery costs.
#' @author Beate Vieth
#' @seealso \code{\link{simulateDE}}, \code{\link{evaluateSim}}
#' @examples
#' \dontrun{
#' ## for example evaluation result see \code{\link{evaluateSim}}
#' printEvalRes(evalRes=evalres)
#' }
#' @rdname printEvalRes
#' @export
printEvalRes <- function(evalRes) {
nreps1 <- evalRes$n1
nreps2 <- evalRes$n2
alpha.type <- evalRes$alpha.type
if(alpha.type == "raw") {
alpha.nam <- "FPR"
alpha.mar <- rowMeans(evalRes$FPR.marginal, na.rm=TRUE)
}
if(alpha.type == "adjusted") {
alpha.nam <- "FDR"
alpha.mar <- rowMeans(evalRes$FDR.marginal, na.rm=TRUE)
}
res <- matrix(0, nrow=length(nreps1), ncol=5)
colnames(res) <- c("Sample size group 1", "Sample size group 2", paste(c("Nominal", "Marginal"), alpha.nam),
"Marginal TPR")
res[,1] <- nreps1
res[,2] <- nreps2
res[,3] <- evalRes$alpha.nominal
res[,4] <- alpha.mar
res[,5] <- rowMeans(evalRes$TPR.marginal, na.rm=TRUE)
print(signif(res,2))
return(invisible(res))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.