cancerit / telomerecat

Telomerecat: The telomere computational analysis tool
GNU General Public License v3.0
20 stars 5 forks source link

Different results with same file as input. Results not reliable. #46

Open spnlab opened 1 year ago

spnlab commented 1 year ago

Hi there, This tools has given me different results with same script and same input but running twice. The script I was running: telomerecat bam2length -N 100 -v 2 -p 32 -e input.bam -r ref.fa --outbam_dir output_dir --output ouptu_telomerecat_length.csv Once I got Length = 2400.46 Other time I got Length = 2418.03

Any idea why?

AndyMenzies commented 1 year ago

There is a random number seed used as part of the algorithm, this can impact the outcome if there are many very good answers with very close likelyhoods.

There is a hidden command line option (--seed_randomness) where you can set the seed value, we use 42 for unit testing.

spnlab commented 1 year ago

Thanks for the updates. Will try (--seed_randomness).

Urja25 commented 7 months ago

There is a random number seed used as part of the algorithm, this can impact the outcome if there are many very good answers with very close likelyhoods.

There is a hidden command line option (--seed_randomness) where you can set the seed value, we use 42 for unit testing.

Hi @AndyMenzies,

I tried using the --seed_randomness 42 flag as an argument in the command but get the following error: error: unrecognized arguments: 42

Can you please help troubleshoot?

Thanks a lot!