gcell.cell.celltype.GETHydraCellType.get_powerlaw_at_distance#
- GETHydraCellType.get_powerlaw_at_distance(distances, min_distance=5000)[source]#
Get the powerlaw trained on k562 Hi-C data at a given distance.
- Parameters:
- distances array-like
The distances to the TSS in base pairs.
- min_distance int, default=5000
Minimum distance threshold in base pairs. Distances below this value will be clipped to this value.
- Returns:
numpy.ndarray Powerlaw contact frequency values corresponding to the input distances.
Notes
The powerlaw parameters are derived from K562 Hi-C data with: - gamma = 1.024238616787792 - scale = 5.9594510043736655
For distances below min_distance (5kb), values are clipped as the model doesn’t accurately represent contact frequencies at short distances.