Entering edit mode
tanyabioinfo
▴
20
@tanyabioinfo-14091
Last seen 5.8 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
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