tximport samples with separate tx2gene maps
0
0
Entering edit mode
@nicolettesipperly-16835
Last seen 1 day 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 • 25 views

Login before adding your answer.

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