-
Notifications
You must be signed in to change notification settings - Fork 0
/
R04 FLOW workflow 2.R
71 lines (56 loc) · 2.31 KB
/
R04 FLOW workflow 2.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
library(flowCore)
library(flowWorkspace)
library(openCyto)
library(ggcyto)
library(flowAI)
library(gridExtra)
library(tidyverse)
library(flowStats)
library(flowWorkspace)
library(CytoML)
library(Rtsne)
library(FlowSOM)
myfiles <- list.files(path="C:/Users/edmondsonef/Desktop/samp/15726 10Mar2022/", pattern = ".FCS", ignore.case = TRUE)
#fs <- read.flowSet(myfiles, path="C:/Users/edmondsonef/Desktop/samp/15726 10Mar2022/")#, truncate_max_range = FALSE)
wsp_file <- "C:/Users/edmondsonef/Desktop/samp/15726 10Mar2022/15726 10Mar2022 Simone.wsp"
fcs_file <- "C:/Users/edmondsonef/Desktop/samp/15726 10Mar2022/Samples_Tube_015 Animal 137 blood_015.fcs"
# Specify the cell types of interest for assigning one label per cell
cell_types <- c("/scatter/sing/hCD45+")
# Parse the FlowJo workspace
gatingResult <- GetFlowJoLabels(fcs_file, wsp_file,
cell_types = cell_types,
getData = TRUE)
# Check the number of cells assigned to each gate
colSums(gatingResult$matrix)
colnames(gatingResult$matrix)
pData(fs) %>% head(3)
colnames(fs)
colnames(fs)[colnames(fs)=="BB515-A"] <- "CD8"
colnames(fs)[colnames(fs)=="BB700-P-A"] <- "CD4"
colnames(fs)[colnames(fs)=="APC-A"] <- "CD11b"
colnames(fs)[colnames(fs)=="APC-Cy7-A"] <- "CD19"
colnames(fs)[colnames(fs)=="BV421-A"] <- "CD3"
colnames(fs)[colnames(fs)=="BV786-A"] <- "CD33"
colnames(fs)[colnames(fs)=="BUV395-A"] <- "mCD45"
colnames(fs)[colnames(fs)=="BUV805-A"] <- "huCD45"
colnames(fs)[colnames(fs)=="PE-A"] <- "CD56"
colnames(fs)[colnames(fs)=="PE-CF594-A"] <- "CD66b"
colnames(fs)[colnames(fs)=="PE-Cy7-A"] <- "CD25"
colnames(fs)
# Build a FlowSOM tree
fsom <- FlowSOM(gatingResult$flowFrame,
#compensate = TRUE,
#transform = TRUE,
#toTransform = c(7:8,11),
colsToUse = c(7:17),
nClus = 20,
seed = 1)
PlotStars(fsom)
PlotFlowSOM(fsom, equalNodeSize = F)
PlotPies(fsom, gatingResult$manual,
backgroundValues = fsom$metaclustering)
PlotManualBars(fsom, manualVector = gatingResult$manual,
manualOrder = c(cellTypes = "/scatter/sing/hCD45+/Q6: CD3+ , CD4 [PCP55]+"))
dups <- duplicated(fsom)
fsom <- fsom[!dups, ]
PlotDimRed(fsom, colsToUse = fsom$map$colsUsed, colorBy = "metaclusters", check_duplicates = FALSE)