Generate dynamic rectangleGate object with flowCore and flowWorkspace

dyna_gates.R

compute_dynamic_gate <- function(gs, inv.bins) {
    df.stat <- data.frame()
    for (min.cutoff in inv.bins) {
        # Create dynamic FL1-A gate
        rg_gfp <- flowCore::rectangleGate(
            "FL1-A" = c(min.cutoff, Inf), "FSC-A" = c(-Inf, Inf),
            filterId = "_dynamic"
        )
        gates <- flowWorkspace::gs_get_pop_paths(gs)
        if ("/P1/_dynamic" %in% gates) {
            flowWorkspace::gs_pop_remove("_dynamic", gs = gs)
            flowWorkspace::gs_pop_add(gs, rg_gfp, parent = "P1")
            flowWorkspace::recompute(gs)
        } else {
            flowWorkspace::gs_pop_add(gs, rg_gfp, parent = "P1")
            flowWorkspace::recompute(gs)
        }
        stat <- flowWorkspace::gs_pop_get_stats(gs, "_dynamic", type = "percent") %>%
            dplyr::mutate(percent, cutoff = inv.log.trans(min.cutoff))
        df.stat <- rbind(df.stat, stat)
    }
    return(df.stat)
}