ruqianl / comapr

Other
2 stars 0 forks source link

feature request - interference function #1

Open genomestability opened 3 years ago

genomestability commented 3 years ago

Hi Ruqian, Would it be possible to add a function that detects interference? E.g. in order to make an "expected" distrubution for the null hypothesis of no interference: 1) take all double crossovers and randomise their labels 2) do this e.g. 1000 times to generate a distribution 3) calculate the p-value is = P(B ≤ b) = (b+1)/ (m+1) as per doi 10.2202/1544-6115.1585

I'm not sure if it's possible to make direct comparisons between groups.

Happy to discuss this request further. Thanks for considering it.

ruqianl commented 3 years ago

Thanks for the request. Yes, it would be possible to add a permute double crossover function to generate a NULL distribution of inter-double crossover distances and compare with the observed distribution.

Direct group comparison might be tricky using this method, we perhaps consider the gamma model which was compared here