## ---------------------------------------------------------------------------------------------------------------------------＃if（！sireseenamespace（“ biocmanager”，悄悄= true））＃install.packages（“ biocmanager”）＃biocmanager :: install（“”gscreend”）## ----消息= false ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------库（GSCREEND）库（摘要Exexperiment）## --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- raw_counts < - read.table（system.file（“ extdata”，“ simulated_counts.txt”，package =“ gscreend”），header = true）## ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- cbind（raw_counts $ library0，raw_counts $ r0_0，raw_counts $ r1_0）rowdata <-data.frame（sgrna_id = raw_counts $ sgrna_id $ sgrna_id，gene gene gene gene gene= raw_counts $ gene）coldata < - data.frame（samplename = c（“库”，“ r1”，“ r2”），＃timepoint命名命名约定：＃t0-> reference，＃t1->增殖时间点= c（“ t0”，“ t1”，“ t1”））se < - summarizedExperiment（assays = list = counts = counts_matrix），rowdata = rowdata，coldata = coldata）## ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- pse <- createPoolScreenExp(se) ## --------------------------------------------------------------------------------------------------------------- pse_an <- RunGscreend(pse) ## ---- fig.width=4, fig.height=4.5--------------------------------------------- plotReplicateCorrelation(pse_an) ## --------------------------------------------------------------------------------------------------------------- plotModelParameters(pse_an) ## --------------------------------------------------------------------------------------------------------------- res < - 结果稳定（pse_an，direction =“否定”）头（res）## ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- sessioninfo（）