tximport samples with separate tx2gene maps
Entering edit mode
Last seen 22 months ago
United States

I have have tx2gene files that are output from Orthofinder. Kallisto abundances were generated from for two different assemblies from different tissues. Ultimately I am hoping to use edgeR TMM means to get a normalization factor across the tissues.

I was able to upload the abundance files according to the vignette https://bioconductor.org/packages/release/bioc/vignettes/tximport/inst/doc/tximport.html#Import_transcript-level_estimates

I used cat to merge the two tx2gene files so that the transcript IDs from both tissues were included in the map, however when I ran the following code I got the error because tximport requires that the tx ids are in the same order for the abundance files and the tx2gene map. Is there a way to deal with this issue where I have separate tx2gene maps for each sample using tximport?

Thank you so much!!

txi.kallisto.tsv <- tximport(files, type = "kallisto", tx2gene = tx2gene)

#tximport is throwing the error 
1 2 Error in tximport(files, type = "kallisto", tx2gene = tx2gene) :
  all(txId == raw[[txIdCol]]) is not TRUE
In addition: Warning message:
In txId == raw[[txIdCol]] :
  longer object length is not a multiple of shorter object length

Code should be placed in three backticks as shown below

# include your problematic code here with any corresponding output 
# please also include the results of running the following in an R session 

sessionInfo( )
tximport • 947 views
Entering edit mode
Last seen 2 days ago
United States

Missed this but followed up on the later post:

tximport error: all(txId == raw[[txIdCol]]) is not TRUE


Login before adding your answer.

Traffic: 551 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6