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 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

Traffic: 345 users visited in the last hour
FAQ
API
Stats

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