-
Notifications
You must be signed in to change notification settings - Fork 7
Description
Hi Boyle-Lab,
I'm able to run the tutorial data successfully; however, when I apply the commands to my actual data I'm getting an error about an -inf gamma. Could you provide advice on what this means and practically how I could resolve it?
For context, the error is in BaumWelch.c, ComputeGamma_P() line 588. The error is only occurring after the reestimation of the variance or covariance matrix, on line 224. I'll attach my toy example input files here. For this exact test data, the error is occurring on the 12th iteration of delta.
`
./TRACE toy_seq.txt toy_count.txt toy_slope_2.txt --initial-model ctcf_init_model.txt --final-model ctcf_final_model.txt --peak-file toy_peaks.bed --thread 1
#ERROR: -inf gamma t:0
`
ctcf_init_model.txt
toy_slope_2.txt
toy_seq.txt
toy_slope_1.txt
toy_count.txt
toy_peaks.txt
Thank you!
Amanda