|
28 | 28 | #' you use a list of gene names use "genes", if you use a target bed use "target". |
29 | 29 | #' @export |
30 | 30 | #' @import rlist |
31 | | -#' @import OrganismDbi |
32 | 31 | #' @import Rsamtools |
33 | | -#' @import org.Hs.eg.db |
34 | 32 | #' @import TxDb.Hsapiens.UCSC.hg19.knownGene |
35 | 33 | #' @import TxDb.Hsapiens.UCSC.hg38.knownGene |
36 | 34 | #' @importFrom utils write.table |
37 | 35 | #' @import GenomicRanges |
| 36 | +#' @import OrganismDbi |
| 37 | +#' @import org.Hs.eg.db |
38 | 38 | #' @importFrom S4Vectors queryHits |
39 | 39 | #' @importFrom S4Vectors subjectHits |
40 | 40 | #' @examples |
@@ -80,8 +80,9 @@ buildInput<- function(geneList,genome,type_bam,bamList,outDir,type_input, |
80 | 80 | ###check control: stop when gene names in list are incorrect |
81 | 81 | if (type_input == "target"){ |
82 | 82 | return=NULL}else{ |
83 | | - my_gene_name<- OrganismDbi::select(org.Hs.eg.db, key= gene.List, columns="ENTREZID", |
| 83 | + my_gene_name_df<- OrganismDbi::select(org.Hs.eg.db, key= gene.List, columns="ENTREZID", |
84 | 84 | keytype="ALIAS") |
| 85 | + my_gene_name<-my_gene_name_df[!duplicated(my_gene_name_df$ALIAS),] |
85 | 86 | for (i in my_gene_name$ENTREZID){ |
86 | 87 | if (is.na(i)){ |
87 | 88 | s<-print(subset(my_gene_name$ALIAS,is.na(my_gene_name$ENTREZID))) |
@@ -246,7 +247,8 @@ buildInput<- function(geneList,genome,type_bam,bamList,outDir,type_input, |
246 | 247 | start.1, end.1, width.1, strand.1)) |
247 | 248 |
|
248 | 249 | colnames(statistiche)[1:3]<- c("chromosome","start","end") |
249 | | - merge_g<- dplyr::full_join(for_bed,statistiche, by="SYMBOL", all=TRUE) |
| 250 | + #merge_g<- dplyr::full_join(for_bed,statistiche, by="SYMBOL", all=TRUE) |
| 251 | + merge_g<- dplyr::full_join(for_bed,statistiche, by="SYMBOL") |
250 | 252 | col_name<- colnames(merge_g) |
251 | 253 |
|
252 | 254 | col.sub<- col_name[grepl("sample_", col_name)] |
|
0 commit comments