Question: inferential replicates in Kallisto when summarising values via tximport
0
gravatar for Mozart
29 days ago by
Mozart10
Mozart10 wrote:

hello there, I am not fully sure what I am doing here. Basically, I have always known that when using tsv files, there should not be any inferential replicates considered. But now, in the latest version of tximport, I guess something has changed?

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


 Note: importing `abundance.h5` is typically faster than `abundance.tsv`
 reading in files with read.delim (install 'readr' package for speed up)
 1 2 3 4 5 6 
 transcripts missing from tx2gene: 22002
 summarizing abundance
 summarizing counts
 summarizing length
 summarizing inferential replicates

I would like to ask if, in this way, inferential replicates are automatically taken into account for the downstream analysis or it is just a new warning message that can also be ignored, because:

> rownames(sampleTable) <- colnames(txi.kallisto.tsv$counts)
> dds <- DESeqDataSetFromTximport(txi.kallisto.tsv, sampleTable, ~condition)
using counts and average transcript lengths from tximport
> dds <- DESeq(dds)
estimating size factors
using 'avgTxLength' from assays(dds), correcting for library size
estimating dispersions
gene-wise dispersion estimates
mean-dispersion relationship
final dispersion estimates
fitting model and testing

thanks

tximport • 40 views
ADD COMMENTlink modified 28 days ago by Michael Love24k • written 29 days ago by Mozart10
Answer: inferential replicates in Kallisto when summarising values via tximport
0
gravatar for Michael Love
28 days ago by
Michael Love24k
United States
Michael Love24k wrote:

This is a few versions back I think.

tximport imports inferential replicates if they are present, and you can prevent this with dropInfReps=TRUE. By the way, those are messages above, not warnings.

The inferential replicates are not used by DESeq2.

ADD COMMENTlink written 28 days ago by Michael Love24k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 190 users visited in the last hour