AlphaGenes / AlphaPeel

AlphaPeel: calling, phasing, and imputing genotype and sequence data in pedigrees
MIT License
2 stars 11 forks source link

What does e=1e-06 in peelingInfo.segregationTensorXY = ProbMath.generateSegregationXYChrom(e=1e-06) mean? #123

Open gregorgorjanc opened 7 months ago

gregorgorjanc commented 7 months ago

Looking at recent PR https://github.com/AlphaGenes/AlphaPeel/pull/122/files I noticed this line of code

peelingInfo.segregationTensorXY = ProbMath.generateSegregationXYChrom(e=1e-06)

@XingerTang and @AprilYUZhang please have a look at what is the meaning of e=1e-06 as in, is this another "error/noise" rate that we would like to control via an argument or is this just a numerical convenience to start a tensor with probabilities that are not zero?