EBSeqHMM Error: Missing value

[sush...@medgenome.com]

Hi,

I am using EBseqHMM to identify genes following upregulated/downregulated trend across nine time points. I have 5 replicates. When I use 3 replicates, I am getting the desired output, i.e. EBSeqHMM runs properly and I can get the clusters of up-up-up-up-up... genes. But, when I use 5 replicates, the EBSeqHMM run doesnot complete, and throws the following error:

Error in if (max(PI.up)) > PIBound) { : missing value where TRUE/FALSE expected

I assumed, some NAs are being produced during the calculations, so I replaced all NAs with 0

GeneNormData[is.na(GeneNormData)] <- 0

Still the run did not complete, and threw the error:

Error:cannot allocate vector of size 584.5 Mb

Please suggest !!