@@ -41,7 +41,7 @@ library(flowCore)
4141library(flowWorkspace)
4242library(CytoML)
4343library(openCyto)
44- library(ggcyto)
44+ # library(ggcyto)
4545library(data.table)
4646library(dplyr)
4747library(purrr)
@@ -79,10 +79,10 @@ Coereba primarily relies on the infrastructure provided `r Biocpkg("flowWorkspac
7979
8080# CRAN packages: install.packages("ThePackageName")
8181
82- # library(dplyr)
83- # library(purrr)
84- # library(stringr)
85- # library(ggplot2)
82+ library(dplyr)
83+ library(purrr)
84+ library(stringr)
85+ library(ggplot2)
8686# library(gt)
8787# library(plotly)
8888# library(htmltools)
@@ -91,11 +91,13 @@ Coereba primarily relies on the infrastructure provided `r Biocpkg("flowWorkspac
9191# Bioconductor packages: BiocManager::install("ThePackageName")
9292
9393library(Coereba)
94- # library(flowCore)
95- # library(flowWorkspace)
96- # library(openCyto)
97- # library(ggcyto)
94+ library(flowCore)
95+ library(flowWorkspace)
96+ library(openCyto)
97+ library(ggcyto)
9898
99+ # GitHub packages: remotes::install_github("DavidRach/Luciernaga")
100+ library(Luciernaga)
99101```
100102
101103## Locate your files
@@ -177,9 +179,6 @@ plot(UnmixedGatingSet)
177179We can additionally verify that the gating of the cell populations of interest was correct, by visualizing the gates using ` r Biocpkg("ggcyto") ` or the ` Luciernaga::Utility_GatingPlots ` function from the ` Luciernaga ` package.
178180
179181``` {r}
180- #| label: "Checking the openCyto gate"
181- library(Luciernaga)
182-
183182MyPlot <- Utility_GatingPlots(x=UnmixedGatingSet[4], sample.name=c("GROUPNAME", "TUBENAME"),
184183 removestrings=".fcs", subset="root", gtFile=UnmixedGates,
185184 DesiredGates=NULL, returnType="plots", outpath = getwd(),
@@ -879,7 +878,7 @@ Specimens <- InitialMetadata |> pull(specimen) |> unique()
879878# Adding Adult Normalization Controls Not Present in Study Metadata
880879Specimens <- c(Specimens, "NY068_02", "NY068_03", "NY068_03", "NY068_4",
881880 "NY068_5", "NY068_6", "NY068_7", "NY068_8")
882- StudyMetadata <- StudyMetadata |> filter(bid %in% Specimens)
881+ StudyMetadata <- StudyMetadata |> dplyr:: filter(bid %in% Specimens)
883882StudyMetadata <- StudyMetadata |> rename(specimen=bid)
884883CoerebaMetadata <- left_join(InitialMetadata, StudyMetadata, by="specimen")
885884```
@@ -969,7 +968,7 @@ We can start by viewing all markers:
969968#| label: MarkerPlots
970969#| warning: FALSE
971970
972- CordOnly <- Data |> filter(ptype %in% c("HU", "HEU-lo", "HEU-hi"))
971+ CordOnly <- Data |> dplyr:: filter(ptype %in% c("HU", "HEU-lo", "HEU-hi"))
973972
974973
975974ThePlot <- Utility_MarkerPlots(data=CordOnly, myfactor="ptype", shape_palette = shape_ptype, fill_palette = fill_ptype, panel=ThePanel, XAxisLevels=c("Vd2", "CD16", "CXCR5", "HLA-DR"), cex=3, size =3)
0 commit comments