...
|
...
|
@@ -441,7 +441,7 @@ genes <- getBM(attributes=c('ensembl_gene_id', 'chromosome_name', 'start_positio
|
441
|
441
|
'end_position', 'strand', 'external_gene_name',
|
442
|
442
|
'gene_biotype'),
|
443
|
443
|
mart=ensembl)
|
444
|
|
-expr <- merge(genes, expression_lv, by='ensembl_gene_name')
|
|
444
|
+expr <- merge(genes, expression_lv, by='ensembl_gene_id')
|
445
|
445
|
# Transform to GRanges
|
446
|
446
|
expression.SHR <- GRanges(seqnames=paste0('chr',expr$chromosome_name),
|
447
|
447
|
ranges=IRanges(start=expr$start, end=expr$end),
|