sesame oracle: compare without bead-count masking on both sides

This commit is contained in:
kemal committed 2026-09-29 00:29:25 +01:00
1 parent 1aa9459525
commit 3a816a1c44
2 files changed
+5 -2

No files matched your search

+3 -1
View File
@@ -24,7 +24,9 @@ inband <- function(s) {
list(m = m, u = u)
}
for (prefix in args[-(1:3)]) {
read <- function() readIDATpair(prefix, platform = "EPIC", manifest = mft, controls = controls, min_beads = 3)
# no bead-count masking here (sesame's min_beads also masks every Infinium II probe, whose M address is NA);
# the comparison script runs simurg with the same setting
read <- function() readIDATpair(prefix, platform = "EPIC", manifest = mft, controls = controls)
sdf <- tryCatch(read(), error = function(e) {
message("first read failed (", conditionMessage(e), "); caching idatSignature and retrying")
sesameDataCache("idatSignature")