Single-sample specimen-overlap provenance pre-flight gate
Source:R/singlesample-provenance-preflight.R
os_provenance_preflight.RdAudits a candidate accession list against a manifest of known specimen- overlap pairs and returns one of three exit conditions:
NO_OVERLAPNo accession in the input list overlaps with any other accession in the input list according to the manifest, and every supplied accession appears (as either side of a pair) somewhere in the manifest.
KNOWN_OVERLAPAt least one pair of input accessions is listed in the manifest. The cross-cohort analysis must apply block-aware multiplicity correction (
singlesample_bh_fdr_correct_blocked).UNKNOWN_ACCESSIONAt least one input accession is not in the manifest. Audit required before combining with anything else. Takes precedence over
KNOWN_OVERLAP.
Intended as the first stage in the single-sample pipeline (Methods section "Provenance auditing and specimen-overlap handling").
Arguments
- accessions
Character vector of GEO/ArrayExpress accessions (e.g.,
c("GSE211692", "GSE106817")).- manifest_path
Path to a TSV file with columns
acc1, acc2, block_id(and optionallyfraction,evidence_source). IfNULL, uses the bundled manifest atsystem.file("extdata", "provenance_manifest.tsv", package = "OmicSelector").
Value
List with components:
statusOne of
"NO_OVERLAP","KNOWN_OVERLAP","UNKNOWN_ACCESSION".overlapping_pairsdata.framesubset of the manifest with rows whoseacc1andacc2are both in the input accession list.unknown_accessionsCharacter vector of input accessions not present in the manifest.
manifest_pathResolved path of the manifest used.
Examples
res <- os_provenance_preflight(c("GSE211692", "GSE106817"))
res$status
#> [1] "KNOWN_OVERLAP"