tximport samples with separate tx2gene maps
1
0
Entering edit mode
@nicolettesipperly-16835
Last seen 9 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 • 507 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 1 day ago
United States

Missed this but followed up on the later post:

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

ADD COMMENT

Login before adding your answer.

Traffic: 165 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