Adding Multiple quadrantGates of Same Channels to workFlow (flowCore)
1
0
Entering edit mode
Aric Gregson ▴ 270
@aric-gregson-3058
Last seen 9.6 years ago
Hello, What I would like to do seems like a common application in analysis and I am probably missing a very easy way to accomplish this. I would like to apply a 'memory gate' generated by quadrantGate to both CD4 and CD8 lymphocyte populations. However, there appears to be no way to give the second gate a unique name for the workFlow. I use a quadrantGate to create filters for channels PE.Cy7.A and PerCP.Cy5.5.A. The quadrantGates are slightly different for the CD4 and CD8 populations (though that does not matter here). The 'action_' name in the workFlow is different for each one, but the workFlow node name appears to be the same thus there is no easy way to call the population for a plot. After adding the gate I must use an arbitrary gateview ID for plotting, which is not at all convenient especially with Sweave. Is there a way around this? Below is a short example. Thanks, Aric ----------------- quadGate.memory.CD4 <- quadrantGate(Data(wf[['CD4+CD8-']]), c("PerCP.Cy5.5.A","PE.Cy7.A"), plot=TRUE, alpha=c(0.5, 1), filterId="CD4_MEMORY" ) add(wf, quadGate.memory.CD4, parent="CD4+CD8-") quadGate.memory.CD8.base <- quadrantGate(Data(wf[['CD4-CD8+']])[[5]], c("PerCP.Cy5.5.A","PE.Cy7.A"), plot=TRUE, alpha=c(0.5, 0.9), filterId="CD8_MEMORY" ) add(wf, quadGate.memory.CD8.base, parent="CD4-CD8+") Warning: 'CD45ra+CCR7+' is not a unique view name in this workFlow object. You will only be able to access the view using its UID 'gateViewRef_OkTFWgIFTF'
• 939 views
ADD COMMENT
0
Entering edit mode
Aric Gregson ▴ 270
@aric-gregson-3058
Last seen 9.6 years ago
On Mon, 14 Nov 2011 15:58:09 -0800 Aric Gregson <a.gregson at="" ucla.edu=""> wrote: > Is there a way around this? OK. I have come up with a way around this, although considering a more global answer for working with the workFlow and perhaps a different naming scheme for nodes/gateViews may be better in the long run. I can plot the gate using the filter results without adding it to the workFlow: > memory.cd8.quadgate.base.result <- filter(Data(wf[['CD4-CD8+']]), quadGate.memory.CD8.base) > xyplot(`PE.Cy7.A` ~ `PerCP.Cy5.5.A` | SampleID, Data(wf[['CD4-CD8+']]), filter=memory.cd8.quadgate.base.result) I can extract the 'boundary' results from the quadGate result and create multiple rectangleGates to add to the workFlow, for example: > memory.cd8.naive <- rectangleGate("PerCP.Cy5.5.A"=c(quadGate.memory.CD8.base at boundary[[1]], Inf), "PE.Cy7.A"=c(quadGate.memory.CD8.base at boundary[[2]], Inf), filterId="CD8_NAIVE")
ADD COMMENT

Login before adding your answer.

Traffic: 812 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6