Visualize flowSet with ggcyto

library(ggcyto)
data(GvHD)
fs <- GvHD[subset(pData(GvHD), Patient %in%5:7 & Visit %in% c(5:6))[["name"]]]
fr <- fs[[1]]

1d histogram/densityplot

ggcyto wrapper will construct the ggcyto object that inherits from ggplot class.

p <- ggcyto(fs, aes(x = `FSC-H`)) 
class(p)

## [1] "ggcyto_flowSet"
## attr(,"package")
## [1] "ggcyto"

is(p, "ggplot")

## [1] TRUE

Since only one dimension is specified, we can add any 1d geom layer

p1 <- p + geom_histogram() 
p1

As shown, data is facetted by samples name automatically (i.e facet_wrap(~name)).

We can overwrite the default faceting by any variables that are defined in pData

pData(fs)

##       Patient Visit Days Grade  name
## s5a05       5     5   19     3 s5a05
## s5a06       5     6   26     3 s5a06
## s6a05       6     5   19     3 s6a05
## s6a06       6     6   27     3 s6a06
## s7a05       7     5   21     3 s7a05
## s7a06       7     6   28     3 s7a06

p1 + facet_grid(Patient~Visit)

To display 1d density

p + geom_density()

Fill the same color

p + geom_density(fill = "black")