flowViz xyplot smooth=F gate outline ONLY
1
0
Entering edit mode
Aric Gregson ▴ 270
@aric-gregson-3058
Last seen 7.7 years ago
Hello, I am trying to use code from a presentation of Deepayan that I found on the web. The code does a nice job a presenting density information in the flowViz xyplot utilizing a more traditional appearance as found in other flow software while still ploting all points. I would like to use this as an occasional alternative to relying upon the alpha setting to depict density. With this code, I would like to: 1) ONLY outline the gated cells rather than color and outline the cells. I have not found a way to do this. The best I have managed is in the example below, which overdraws with a light rose color. Setting fill = "transparent" appears to have no effect. 2) Apply this type of plot to a flowSet. I have not figured this out yet either. Replacing panel.xyplot.flowframe with panel.xyplot.flowset fails. Any hints would be greatly appreciated. xyplot(`Indo.1..Violet..A` ~ `APC.H7.A`, Data(wf[["logicle"]]) [[4]], smooth=F, pch=21, cex=0.3, aspect="iso", panel = function(x, y, col, ...) { cols <- densCols(x, y, nbin=128, colramp = colori) panel.xyplot(x, y, col = cols, ...) panel.xyplot.flowframe(x, y, filter=livefilter, outline=TRUE, names=TRUE, col="transparent", ...) }, par.settings=list(gate=list(col= cols, fill= cols, alpha=.9), gate.text=list(col="black", alpha=0.7, cex=0.6))) Thanks in advance. Aric ============================ R version 2.8.0 (2008-10-20) i386-pc-solaris2.11 locale: C attached base packages: [1] tools stats graphics grDevices utils datasets methods [8] base other attached packages: [1] geneplotter_1.20.0 annotate_1.20.1 xtable_1.5-4 [4] AnnotationDbi_1.4.1 colorspace_1.0-0 flowViz_1.7.10 [7] lattice_0.17-20 flowCore_1.9.18 rrcov_0.4-07 [10] robustbase_0.4-3 Biobase_2.2.1 Hmisc_3.4-4 loaded via a namespace (and not attached): [1] DBI_0.2-4 KernSmooth_2.22-22 MASS_7.2-44 RColorBrewer_1.0-2 [5] RSQLite_0.7-1 cluster_1.11.11 feature_1.1-14 graph_1.20.0 [9] grid_2.8.0 latticeExtra_0.5-4 stats4_2.8.0 tcltk_2.8.0
flowViz flowViz • 855 views
ADD COMMENT
0
Entering edit mode
@deepayan-sarkar-1436
Last seen 7.7 years ago
On Wed, Jun 24, 2009 at 4:15 PM, Aric Gregson<a.gregson at="" ucla.edu=""> wrote: > Hello, > > I am trying to use code from a presentation of Deepayan that I found > on the web. The code does a nice job a presenting density information in > the flowViz xyplot utilizing a more traditional appearance as found in > other flow software while still ploting all points. I would like to > use this as an occasional alternative to relying upon the alpha setting > to depict density. With this code, I would like to: > > 1) ONLY outline the gated cells rather than color and outline the > cells. I have not found a way to do this. The best I have managed is in > the example below, which overdraws with a light rose color. Setting > fill = "transparent" appears to have no effect. > > 2) Apply this type of plot to a flowSet. I have not figured this out > yet either. Replacing panel.xyplot.flowframe with panel.xyplot.flowset > fails. Looks like you are calling panel.xyplot.flowframe() only to draw the gate. That's not necessary; you can call glpolygon() directly. For example, library(flowViz) data(GvHD) my.panel <- function(x, y, ..., frame, my.filter) { cols <- densCols(x, y, nbin=128, colramp = colorRampPalette(rainbow(7))) panel.xyplot(x, y, col = cols, ...) glpolygon(filter(frame, my.filter)) } xyplot(`FSC-H` ~ `SSC-H`, GvHD[[1]], my.filter = norm2Filter("SSC-H", "FSC-H"), pch = 21, cex = 0.3, panel = my.panel) The details may need to be tweaked depending on your filter. This is harder to generalize for flowSets so that all kinds of filters are handled, but specific types are simple: xyplot(`FSC-H` ~ `SSC-H`, GvHD[1:4], my.filter = norm2Filter("SSC-H", "FSC-H"), pch = 21, cex = 0.3, panel = function(x, frames, channel.x, channel.y, ...) { nm <- as.character(x) x <- flowViz:::evalInFlowFrame(channel.x, frames[[nm]]) y <- flowViz:::evalInFlowFrame(channel.y, frames[[nm]]) my.panel(x, y, frame = frames[[nm]], ...) }) -Deepayan > > Any hints would be greatly appreciated. > > xyplot(`Indo.1..Violet..A` ~ `APC.H7.A`, Data(wf[["logicle"]]) > [[4]], > ? ? ? smooth=F, > ? ? ? pch=21, > ? ? ? cex=0.3, > ? ? ? aspect="iso", > ? ? ? panel = function(x, y, col, ...) { > ? ? ? ? cols <- > ? ? ? ? ? densCols(x, y, > ? ? ? ? ? ? ? ? nbin=128, > ? ? ? ? ? ? ? ? colramp = colori) > ? ? ? ? panel.xyplot(x, y, col = cols, ...) > ? ? ? ? panel.xyplot.flowframe(x, y, filter=livefilter, outline=TRUE, > ? ? ? ? ? ? ? ?names=TRUE, col="transparent", ...) }, > ? ? ? par.settings=list(gate=list(col= cols, fill= cols, alpha=.9), > ? ? ? ? gate.text=list(col="black", alpha=0.7, cex=0.6))) > > Thanks in advance. > > Aric > ============================ > > R version 2.8.0 (2008-10-20) > i386-pc-solaris2.11 > > locale: > C > > attached base packages: > [1] tools ? ? stats ? ? graphics ?grDevices utils ? ? datasets > methods [8] base > > other attached packages: > ?[1] geneplotter_1.20.0 ?annotate_1.20.1 ? ? xtable_1.5-4 > ?[4] AnnotationDbi_1.4.1 colorspace_1.0-0 ? ?flowViz_1.7.10 > ?[7] lattice_0.17-20 ? ? flowCore_1.9.18 ? ? rrcov_0.4-07 > [10] robustbase_0.4-3 ? ?Biobase_2.2.1 ? ? ? Hmisc_3.4-4 > > loaded via a namespace (and not attached): > ?[1] DBI_0.2-4 ? ? ? ? ?KernSmooth_2.22-22 MASS_7.2-44 > RColorBrewer_1.0-2 [5] RSQLite_0.7-1 ? ? ?cluster_1.11.11 > feature_1.1-14 ? ? graph_1.20.0 [9] grid_2.8.0 > latticeExtra_0.5-4 stats4_2.8.0 ? ? ? tcltk_2.8.0 > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor >
ADD COMMENT

Login before adding your answer.

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