FluteRRA {MAGeCKFlute} | R Documentation |
Integrative analysis pipeline using the gene summary table in MAGeCK RRA results
FluteRRA(gene_summary, sgrna_summary, lfcCutoff = c(-1, 1), organism = "hsa", limit = c(1, 120), pvalueCutoff = 0.25, prefix = "Test", width = 12, height = 6, outdir = ".")
gene_summary |
A file path or a data frame of gene summary data. |
sgrna_summary |
A file path or a data frame of sgRNA summary data. |
lfcCutoff |
A two-length vector (default: c(-1, 1)), specifying the logFC cutoff for negative selection and positive selection. |
organism |
"hsa" or "mmu". |
limit |
A two-length vector (default: c(1, 120)), specifying the minimal and maximal size of gene sets for enrichent analysis. |
pvalueCutoff |
A numeric, specifying pvalue cutoff of enrichment analysis, default 1. |
prefix |
A character, indicating the prefix of output file name. |
width |
The width of summary pdf in inches. |
height |
The height of summary pdf in inches. |
outdir |
Output directory on disk. |
MAGeCK RRA allows for the comparison between two experimental conditions. It can identify genes and sgRNAs are significantly selected between the two conditions. The most important output of MAGeCK RRA is the file 'gene_summary.txt'. MAGeCK RRA will output both the negative score and positive score for each gene. A smaller score indicates higher gene importance. MAGeCK RRA will also output the statistical value for the scores of each gene. Genes that are significantly positively and negatively selected can be identified based on the p-value or FDR.
The downstream analysis of this function includes identifying positive and negative selection genes, and performing biological functional category analysis and pathway enrichment analysis of these genes.
All of the pipeline results is output into the out.dir
/prefix
_Results,
which includes a pdf file and a folder named 'RRA'.
Wubing Zhang
data("rra.gene_summary") data("rra.sgrna_summary") ## Not run: # Run the FluteRRA pipeline FluteRRA(rra.gene_summary, rra.sgrna_summary, prefix="RRA", organism="hsa") ## End(Not run)