Browse code

ensembl issue

ataudt authored on 03/11/2021 08:34:25
Showing 2 changed files

... ...
@@ -29,7 +29,7 @@
29 29
 #'                            'end_position', 'strand', 'external_gene_name',
30 30
 #'                            'gene_biotype'),
31 31
 #'               mart=ensembl)
32
-#'expr <- merge(genes, expression_lv, by='ensembl_gene_name')
32
+#'expr <- merge(genes, expression_lv, by='ensembl_gene_id')
33 33
 #'# Transform to GRanges
34 34
 #'expression.SHR <- GRanges(seqnames=paste0('chr',expr$chromosome_name),
35 35
 #'                          ranges=IRanges(start=expr$start, end=expr$end),
... ...
@@ -38,7 +38,7 @@ genes <- getBM(attributes=c('ensembl_gene_id', 'chromosome_name', 'start_positio
38 38
                            'end_position', 'strand', 'external_gene_name',
39 39
                            'gene_biotype'),
40 40
               mart=ensembl)
41
-expr <- merge(genes, expression_lv, by='ensembl_gene_name')
41
+expr <- merge(genes, expression_lv, by='ensembl_gene_id')
42 42
 # Transform to GRanges
43 43
 expression.SHR <- GRanges(seqnames=paste0('chr',expr$chromosome_name),
44 44
                          ranges=IRanges(start=expr$start, end=expr$end),