Support RleList genomic tracks for compartment eigenvector phasing#18
Conversation
Agent-Logs-Url: https://github.com/js2264/HiContacts/sessions/b522298c-8b9e-4b0a-9cc0-70727bf49c6b Co-authored-by: js2264 <19528818+js2264@users.noreply.github.com>
Agent-Logs-Url: https://github.com/js2264/HiContacts/sessions/b522298c-8b9e-4b0a-9cc0-70727bf49c6b Co-authored-by: js2264 <19528818+js2264@users.noreply.github.com>
RleList genomic tracks for compartment eigenvector phasing
|
In the test for phasing with coverage track, use a coverage track of GC content (calculated from the yeast genomic sequence, and smoothed along chromosomes with a bon width of 100) |
| test_that("compartments can be phased with an RleList track", { | ||
| full_contacts_yeast <- HiCExperiment::contacts_yeast(full = TRUE) | ||
| full_contacts_yeast_VI <- full_contacts_yeast["VI"] | ||
| VI_regions <- HiCExperiment::regions(full_contacts_yeast_VI) |
There was a problem hiding this comment.
@copilot In the test for phasing with coverage track, use a coverage track of GC content (calculated from the yeast genomic sequence, and smoothed along chromosomes with a bon width of 100). Also, IRanges::Rle() does not exist
Agent-Logs-Url: https://github.com/js2264/HiContacts/sessions/f02bcc18-320e-4ba4-9316-8b5381c7eedc Co-authored-by: js2264 <19528818+js2264@users.noreply.github.com>
Agent-Logs-Url: https://github.com/js2264/HiContacts/sessions/f02bcc18-320e-4ba4-9316-8b5381c7eedc Co-authored-by: js2264 <19528818+js2264@users.noreply.github.com>
Updated in commit cfe7482 (plus follow-up cleanup 06a8417): the test now builds an |
Agent-Logs-Url: https://github.com/js2264/HiContacts/sessions/f02bcc18-320e-4ba4-9316-8b5381c7eedc Co-authored-by: js2264 <19528818+js2264@users.noreply.github.com>
Agent-Logs-Url: https://github.com/js2264/HiContacts/sessions/f02bcc18-320e-4ba4-9316-8b5381c7eedc Co-authored-by: js2264 <19528818+js2264@users.noreply.github.com>
Agent-Logs-Url: https://github.com/js2264/HiContacts/sessions/f02bcc18-320e-4ba4-9316-8b5381c7eedc Co-authored-by: js2264 <19528818+js2264@users.noreply.github.com>
getCompartments()already supported phasing with sequence/gene-derived inputs, but not with coverage-style genomic tracks. This PR adds support for BigWig-imported tracks provided asRleList, so compartment sign phasing can use arbitrary per-base signal tracks (e.g. GC coverage tracks).Core fix:
RleListphasing path.eigChr(), theRleListbranch now correctly uses the provided track as the coverage source before computing per-bin means.API/docs alignment
getCompartments()parameter docs to explicitly list supportedgenomeinput types:BSgenome,DNAStringSet,TxDb, andRleList.Focused test coverage
RleListtrack and runs:getCompartments(full_contacts_yeast, genome = cov_track, chromosomes = "VI")RleListphasing route executes without error for a chromosome-scoped run.