uscbiostats / hilda Goto Github PK

Hello Zhi,

I am trying to run the HiLDA tutorial provided with 3 flanking bases instead of 5 bases on your simulated data and esophageal data. This is the code that I ran (for the simulated data, I aggregated the 5-base counts and features to 3-base):

# HiLDA simulated data
set.seed(123)
## Run hildaTest with 3-base (aggregated simulated data)
hildaGlobal <- hildaTest(inputG=newG, numSig=3, localTest=FALSE, refGroup=1:4, nIter=1000)

# Espohageal data
inputFile <- system.file("extdata/esophageal.mp.txt.gz", package="HiLDA")
G <- hildaReadMPFile(inputFile, numBases=3, trDir=TRUE)
hildaGlobal <- hildaTest(inputG=G, numSig=4, localTest=FALSE, refGroup=1:60, nIter=1000)

Both running these codes gives me errors while running JAGS:
Error in update.jags(object, n.iter, ...) : Error in node sCat1[4,76] Cannot normalize density

Error in update.jags(object, n.iter, ...) : Error in node sCat1[60,187] Cannot normalize density

Would you have any suggestions as to what might be going wrong here?

Best,
Ji-Eun