flowStats quadrantGate and Back-Gating
0
0
Entering edit mode
Aric Gregson ▴ 50
@aric-gregson-3687
Last seen 10.2 years ago
Hello, I have tried but have not been successful in using the quadrantGate function in flowStats for back-gating. I would like to turn the quadrantGate filter into a gate result as I do with other filters, but it will not plot the results into another population. For instance, I can create a rangeGate as shown: > cd45rarangegate <- rangeGate(Data(wf[["Live CD3+"]]),#[c(3,4)], stain="Qdot.605.A", filterId="CD45RA", plot=TRUE, sd=0.2) Then create a filter result: > cd45rarangegate_result <- filter(Data(wf[['CD4+CD8-']]), cd45rarangegate) Then demonstrate where these cells within the above filter result are in a different plot: > xyplot(`BD.Horizon.V450.A` ~ `AlexaFluor.647.A`, Data(wf[['CD4+CD8-']]), filter=cd45rarangegate_result, smooth=F, pch=21, alpha=0.5, par.settings=list(gate=list(col="purple", alpha=0.6)) ) However, I cannot do the same with a quadrantGate. I can show one population at a time by making a new rectangleGate filter result, but this is not exactly what I want and is a bit of a hassle. Create the filter result as above: > quadgate_cd4cd8_result <- filter(Data(wf[['Lymphs+']]), quadgate_cd4cd8) Then try to plot as above: > xyplot(`PE.A` ~ `Alexa.Fluor.700.A`, Data(wf[['Lymphs+']]), smooth=F, pch=21, cex=0.4, filter=quadgate_cd4cd8_result ) This fails with 'error using packet 1 in evaluating argument 'subset' in selecting a method for Subset'. I'm certain there is a way around this and I have been unable to figure it out thus far. Thanks in advance for suggestions. Thanks, Aric p.s. is it possible to make quadrantGate plot in a smooth=FALSE xyplot? sessionInfo() R version 2.12.2 (2011-02-25) Platform: amd64-portbld-freebsd8.2 (64-bit) locale: [1] C attached base packages: [1] grid splines tools stats graphics grDevices utils datasets methods [10] base other attached packages: [1] plateCore_1.8.0 heR.Misc_0.0.4 IDPmisc_1.1.10 Rgraphviz_1.26.0 [5] graph_1.26.0 flowQ_1.9.0 latticeExtra_0.6-14 RColorBrewer_1.0-2 [9] parody_1.6.0 bioDist_1.20.0 KernSmooth_2.23-4 outliers_0.13-3 [13] flowStats_1.8.0 cluster_1.13.3 mvoutlier_1.4 fda_2.2.3 [17] zoo_1.6-4 flowViz_1.14.0 lattice_0.19-17 flowCore_1.16.0 [21] rrcov_1.1-00 pcaPP_1.8-2 mvtnorm_0.9-92 robustbase_0.5-0-1 [25] Biobase_2.8.0 fortunes_1.4-0 loaded via a namespace (and not attached): [1] AnnotationDbi_1.10.2 DBI_0.2-5 MASS_7.3-11 RSQLite_0.9-2 [5] annotate_1.26.1 feature_1.2.5 geneplotter_1.26.0 ks_1.7.1 [9] stats4_2.12.2 tcltk_2.12.2 xtable_1.5-6
flowStats flowStats • 1.3k views
ADD COMMENT

Login before adding your answer.

Traffic: 702 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