Question: EBSeqHMM error with many replicates per time point
0
gravatar for Akos
11 months ago by
Akos20
Akos20 wrote:

Hi All,

Anyone using EBSeqHMM out there?

I would like to use EBSeqHMM on a time - series RNA-seq data, where I have 5 time points and more than 10 replicates for each time point. However, the time points have different number of replicates.

First, when I tested EBSeqHMM on a dummy dataset with only triplicates per time-point it was working. However, with my real data It gives error message at the following step:

EBSeqHMMGeneOut <- EBSeqHMMTest(Data=matrix2, sizeFactors=Sizes, Conditions=Conditions, UpdateRd=5)

Error in if (max(Pi.up) > PIBound) { : 
    Missing value where TRUE/FALSE needed

Anyone knows what might be the problem?

Or is there a way to use EBSeqHMM with many and uneven replicates per time points? Thanks in advance.

rna-seq time-course ebseqhmm • 275 views
ADD COMMENTlink modified 9 months ago by Biostar ♦♦ 20 • written 11 months ago by Akos20
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 2.3.0
Traffic: 833 users visited in the last hour