## ----------------------------------------------------------------------------------------------------------------------- # install.packages("BiocManager") # BiocManager::install("MEB") ## -------------------------------------------------------------------------------------------------------------------库（MEB）## ----------------------------------------------------------------------------------------------------------------------------------------------------------- 数据（sim_data_sp）sim_data_sp ## ----------------------------------------------------------------------------------------------------------------------------------------- data(real_data_sp) real_data_sp ## -------------------------------------------------------------------------------------------------------------------数据（SIM_DATA_DSP）SIM_DATA_DSP ## --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------数据（real_data_dsp）real_data_dsp ## ----消息= false，警告= false -------------------------------------------------------------------------------------------库（摘要折磨）## ------------------------------------------------------------------------------------------------------------------------------------------------------------------- data(sim_data_sp) gamma<-seq（1E-06,5E-05,1E-06）sim_model_sp <-nimeb（countstable = assay（sim_data_sp），train_id = 1：1000，伽马，nu = 0.01，reffec_rate = 0.05，ds = false）## ------------------------------------------------------------------------------------------------------------------------------------------------------------------------- data(real_data_sp) gamma<-seq（1E-06,5E-05,1E-06）real_model_sp <-nimeb（countstable = assay（real_data_sp），train_id = 1：530，gamma，nu = 0.01，nu = 0.01，reffect_rate = 0.1--------------------------------------------------------------------------------------------------------------- data(sim_data_dsp) gamma <- seq(1e-07,2e-05,1e-06) sim_model_dsp <- NIMEB(countsTable=assay(sim_data_dsp), train_id=1:1000, gamma, nu = 0.01, reject_rate = 0.1, ds = TRUE) ## ----------------------------------------------------------------------------- data(real_data_dsp) gamma <- seq(5e-08,5e-07,1e-08) real_model_dsp <- NIMEB(countsTable=assay(real_data_dsp), train_id=1:143, gamma, nu = 0.01, reject_rate = 0.1, ds = TRUE) ## ----------------------------------------------------------------------------- sim_model_sp_pred <- predict(sim_model_sp$model, assay(sim_data_sp)) summary(sim_model_sp_pred)