I am working with RNA-seq count data. when I am using EBSeqHMM EBSeqHMMTest function, I got an error when UpdateRd is larger than a number. Do you know how to fix it or why it generates an error like this? Thank you! This is the function: EBSeqHMMTest (Data=GeneExampleData, sizeFactors=Sizes, Conditions=Conditions, UpdateRd=5) Here is the error: Error in if (max(PI.up) > PIBound) { : missing value where TRUE/FALSE needed