Skip to content

The hmr command fails when methylation levels are at extremes #251

@andrewdavidsmith

Description

@andrewdavidsmith

Describe the bug
When methylation levels are either 0 or 1 at all sites, there seems to be a bug.

To Reproduce
This was found by simulating a methylome with those extremes: a set of intervals used as HMRs was used and sites within the HMRs had methylation level of 0 while those outside had 1.0. Then reads were simulated with high coverage (e.g., 200).

Expected behavior
If the estimation is not possible, or does not make sense, then it should be reported.

Desktop (please complete the following information):
Ubuntu 22.04

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions