RNA seq tximport
1
0
Entering edit mode
tanyabioinfo ▴ 20
@tanyabioinfo-14091
Last seen 5.1 years ago

Hi

I am using tximport to generate the count matrix for RNA seq data:

txi <- tximport(files, type="salmon", tx2gene=tx2gene, ignoreTxVersion=TRUE,dropInfReps=TRUE)

I am using following to generate tx2gene

 

library(EnsDb.Mmusculus.v79)

txdf <- transcripts(EnsDb.Mmusculus.v79, return.type="DataFrame")

 

tx2gene <- as.data.frame(txdf[,c("tx_id", "gene_id")])

 

Is there a way so that I can get gene name instead of gene ids.

 

Regards

 

Tanya

 

 

 

tximport rnaseq • 864 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 1 day ago
United States
> txdf <- transcripts(EnsDb.Mmusculus.v79, return.type = "DataFrame")
> txdf$symbol <- mapIds(EnsDb.Mmusculus.v79, txdf$gene_id, "GENENAME", "GENEID")
> tx2gene <- as.data.frame(txdf[,c("tx_id","symbol")])
> head(tx2gene)
               tx_id  symbol
1 ENSMUST00000082387   mt-Tf
2 ENSMUST00000179436    Cd99
3 ENSMUST00000082388 mt-Rnr1
4 ENSMUST00000177695  Ccl21c
5 ENSMUST00000082389   mt-Tv
6 ENSMUST00000082390 mt-Rnr2
ADD COMMENT
0
Entering edit mode

Hi James

Thanks for your help. I did this and got the tx2gene. However when I run tximport

 txi <- tximport(files, type="salmon", tx2gene=tx2gene ,ignoreTxVersion=TRUE,dropInfReps=TRUE)

it gives me an extra row at the start:

   

    0 wt    0 wt    0 wt    0 wt    6 wt    6 wt
    16    172.3771    2548.5088    2603.6048    30.00003    127.9545
0610007P14Rik    0    10.12589    669.19047    578.66497    6    22.00608
0610009B22Rik    0    6    223.9998    207    0    11

Can you help me to rectify this.

 

Tanya

0610007P14Rik    0    10.12589    669.19047    578.66497    6    22.00608    544.64281    502.03345    12.12808    409.58344    85.20258    115.00295    15.16641    11    0    403.9625    19.01683    183.04368    56.00142    24.00381    27.0143    152.0066    244.04903    36.6826

 

 

ADD REPLY

Login before adding your answer.

Traffic: 486 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6