# #设置,包括= FALSE ----------------------------------------------------- knitr: opts_chunk美元集(echo = TRUE, fig.width = 8, fig.height = 6) # #——stage0,结果=“隐藏”,消息= FALSE -------------------------------- 库(flowClust)库(flowCore) # #——stage1 ------------------------------------------------------------------- 数据(利妥昔单抗)总结(利妥昔单抗)res1 < flowClust(以美罗华varNames = c(“FSC.H”、“SSC.H”),K = 1,B = 100 ) ## ---- stage2 ------------------------------------------------------------------- rituximab2 < -利妥昔单抗(美罗华%,% res1)它< flowClust (rituximab2 varNames = c(“FL1.H”、“FL3.H”),K = 1:6, B = 100 ) ## ---- stage2Result ------------------------------------------------------------- 标准(它,(BIC)总结它[[3 ]]) ## ---- stage2ChangeRule1 -------------------------------------------------------- ruleOutliers(它[[3]])< -列表(= 0.95)总结(它[[3 ]]) ## ---- stage2ChangeRule2 -------------------------------------------------------- ruleOutliers(它[[3]])< - - - - - - (z列表。截止= 0.6)总结(它[[3 ]]) ## ---- stage2Alternative -------------------------------------------------------- flowClust (rituximab2 varNames = c(“FL1.H”、“FL3.H”),K = 2, B = 100分钟= c (0, 0), max = c(400年,800年 ) ) ## ---- stage2Scatter ------------------------------------------------------------ 情节(它[[3]],数据= rituximab2 = 0.8水平,z.cutoff = 0) # #——stage2Contour ------------------------------------------------------------ res2.den < -密度(它[[3]],数据= rituximab2)情节(res2.den) # #——stage2Image -------------------------------------------------------------- 情节(res2.den type = "形象 ") ## ---- stage2Hist --------------------------------------------------------------- 嘘(它[[3]],数据= rituximab2 =“FL1.H子集 ") ## ---- stage2Hist2 -------------------------------------------------------------- 嘘(它[[3]],data = rituximab2,= 1) # #——stage2f子集 ------------------------------------------------------------------ s2filter < - tmixFilter(“s2filter”,c(“FL1.H”、“FL3.H”),K = 3, B = 100) res2f < -过滤器(rituximab2 s2filter) # #——stage2fSubsetting警告= FALSE ----------------------------------------- 子集(rituximab2 res2f)分裂(rituximab2, res2f,人口=列表(那么= 1:2,星际2 = 3 )) ## ---- stage2fRectGate ---------------------------------------------------------- rectGate < - rectangleGate (filterId =“rectRegion”、“FL1.H”= c(0, 400),“FL3.H”= c (0, 800)) MBCfilter < - tmixFilter(“MBCfilter”,c(“FL1.H”、“FL3.H”),K = 2, B = 100)过滤器(rituximab2, MBCfilter % % rectGate子集)# #——之前 -------------------------------------------------------------------- 之前set.seed(100)库(flowStats) < - flowClust2Prior(它[[2]],k = 1,Nt = 5000) prior2 <- prior prior2$Mu0[1，] <- rep(box(200, prior2$lambda)， 2) prior2$Lambda0 <- prior2$Lambda0 / 2 pfit2 <- flowClust(rituximab2, varNames = c("FL1.H"， "FL3.H")， K = 2, prior = prior2, usprior = "yes") par(mfrow = c(1,2)) plot(res2[[2]]， data = rituximab2) plot(pfit2, data = rituximab2)