con_fragility: con_fragility

In LyonsZhang/ProgPerm: Progressive permutation for a dynamic representation of the robustness of microbiome discoveries

Usage

con_fragility(alloutputs, top_pm, alpha = 0.05, lgndcol = 2, yrange = NULL, pvtitle = NULL, savepvfile, pvpicdim)

Arguments

alloutputs
top_pm
alpha
lgndcol
yrange
pvtitle
savepvfile
pvpicdim

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (alloutputs, top_pm, alpha = 0.05, lgndcol = 2, yrange = NULL, 
    pvtitle = NULL, savepvfile, pvpicdim) 
{
    require(ggplot2)
    results <- alloutputs$alldist
    selname0 <- alloutputs$selname
    if (top_pm > length(selname)) {
        stop("top_pm should be smaller than the number of observations.")
    }
    df <- as.data.frame(results$pvdist)
    N <- length(unique(df$dif))
    df$dif <- factor(df$dif, levels = sort(unique(df$dif)))
    selname <- selname0[1:top_pm]
    if (top_pm < length(selname0)) {
        tt <- split(df, df$dif)
        ss <- lapply(tt, function(x) {
            return(x[1:top_pm, ])
        })
        df <- do.call("rbind", ss)
    }
    dat <- matrix(as.numeric(as.character(df$rep)), N, top_pm, 
        byrow = T)
    df <- data.frame(fragility = integer(top_pm), variable = character(top_pm), 
        stringsAsFactors = FALSE)
    df$fragility <- colSums(dat > -log10(alpha))
    df$variable <- selname
    df <- df[order(-df$fragility), ]
    df$variable <- factor(selname, levels = selname)
    df$group <- factor(df$fragility, levels = sort(unique(df$fragility), 
        decreasing = T))
    p <- ggplot(data = df, aes(x = variable, y = fragility, fill = group)) + 
        geom_bar(stat = "identity")
    if (is.null(yrange)) {
        y_range <- range(0, ceiling(max(df$fragility)))
    }
    else {
        y_range <- yrange
    }
    ggplot(data = df, aes(x = variable, y = fragility, fill = variable)) + 
        geom_bar(stat = "identity") + scale_fill_manual(aes(labels = variable), 
        values = ggplot_build(p)$data[[1]]$fill, drop = FALSE) + 
        ylab("Fragility Index") + xlab(paste("Top", top_pm, "significant features with decreasing order", 
        sep = " ")) + guides(fill = guide_legend(ncol = 2)) + 
        ggtitle(pvtitle) + ylim(y_range) + theme_bw() + theme(plot.title = element_text(size = 24, 
        face = "bold"), legend.title = element_text(size = 20), 
        legend.text = element_text(size = 12)) + theme(axis.text.x = element_blank(), 
        axis.ticks.x = element_blank(), axis.line = element_line(colour = "black"), 
        panel.border = element_blank(), panel.background = element_blank(), 
        axis.line.x = element_line(color = "black", size = 1), 
        axis.line.y = element_line(color = "black", size = 1), 
        axis.text = element_text(size = 16, face = "bold"), axis.title = element_text(size = 20, 
            face = "bold"))
    ggsave(filename = savepvfile, width = pvpicdim[1], height = pvpicdim[2])
    return(df)
  }

LyonsZhang/ProgPerm documentation built on July 16, 2020, 12:45 a.m.