library(flowWorkspace)
library(openCyto)
library(ggcyto)
load gating set
dataDir <- system.file("extdata",package="flowWorkspaceData")
gs <- load_gs(list.files(dataDir, pattern = "gs_manual",full = TRUE))
get parent CD3+
fr <- getData(gs[[1]], "CD3+")
visually check the target population locations
p <- autoplot(fr, x = "CD4", y = "CD8") + xlim(c(0,4e3))
p

gate out cd8+cd4- by setting target properly
cd8.gate <- flowClust.2d(fr, "<B710-A>", "<R780-A>", K = 2, target = c(1000, 3500))
flowClust gate out cd8-cd4+
cd4.gate <- flowClust.2d(fr, "<B710-A>", "<R780-A>", K = 2, target = c(3000, 100))
visualize the gates
p + geom_gate(cd4.gate) + geom_gate(cd8.gate)
