Skip to content

Commit 22b6818

Browse files
committed
fix FA function
1 parent 142091f commit 22b6818

File tree

2 files changed

+7
-7
lines changed
  • inst/templates/rnaseq

2 files changed

+7
-7
lines changed

inst/templates/rnaseq/00_libs/FA.R

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -40,15 +40,15 @@ get_databases=function(sps="human"){
4040
}
4141

4242
run_fora_v2=function(input, uni, all_in_life){
43-
total_deg=length(unique(input$ENTREZID))/length(unique(uni$ENTREZID))
43+
total_deg=length(unique(input$entrez))/length(unique(uni$entrez))
4444
pathways_ora_all = lapply(names(all_in_life), function(database){
4545
p = all_in_life[[database]]
4646
#browser()
4747
pathway = split(x = p$entrez_gene, f = p$gs_name)
4848
db_name = database
4949
respath <- fora(pathways = pathway,
50-
genes = unique(input$ENTREZID),
51-
universe = unique(uni$ENTREZID),
50+
genes = unique(input$entrez),
51+
universe = unique(uni$entrez),
5252
minSize = 15,
5353
maxSize = 500)
5454
respath %>%
@@ -57,11 +57,11 @@ run_fora_v2=function(input, uni, all_in_life){
5757
mutate(analysis="ORA")
5858
ora_tb = pathways_ora_all %>% unnest(overlapGenes) %>%
5959
group_by(pathway) %>%
60-
left_join(uni, by =c("overlapGenes"="ENTREZID")) %>%
61-
dplyr::select(pathway, padj, NES, SYMBOL, analysis,
60+
left_join(uni, by =c("overlapGenes"="entrez")) %>%
61+
dplyr::select(pathway, padj, NES, gene_name, analysis,
6262
database) %>%
6363
group_by(pathway,padj,NES,database,analysis) %>%
64-
summarise(genes=paste(SYMBOL,collapse = ","))
64+
summarise(genes=paste(gene_name,collapse = ","))
6565
ora_tb
6666

6767
}

inst/templates/rnaseq/02_differential_expression/DEG.Rmd

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -660,7 +660,7 @@ fa_list=lapply(de_list,function(contrast){
660660
#change to the right species
661661
input_entrezid <- rdata %>% filter(gene_id %in% ora_input, !is.na(entrez))
662662
# AnnotationDbi::select(ann.org, ora_input, 'ENSEMBL', columns = c('ENTREZID', 'SYMBOL'))
663-
all=run_fora_v2(input_entrezid, universe_mapping,all_in_life)
663+
all=run_fora_v2(input_entrezid, universe_mapping, all_in_life)
664664
665665
ora_input = res %>% filter(!is.na(padj), padj<0.01, lfc> 0.3) %>% pull(gene_id)
666666
#change to the right species

0 commit comments

Comments
 (0)