＃＃ - - 在里面 - - - - - - - - - - - - - - - - - - - - - --------------------------- library(profileScoreDist) ## ----inr--------------------------------------------------------------------------------------------------------数据（INR）INR ## ----正规化------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- inr.reg <- regularizeMatrix(INR) inr.reg ## ----DIST参数----------------------------------------------------------------------------------------------------------------------------------------------------------------------粒度<-0.05 gcgran <-0.01 gcmin <-0.01 gcmax <-0.99 ## gc分数考虑gcs <-seq（gcgran*round（gcmin/gcgran），gcmin/gcgran），gcgran*round(gcmax/gcgran), gcgran) ## ----dist--------------------------------------------------------------------- ## compute probability distributions distlist <- lapply(gcs, function（x）computeScoredist（inr.reg，x，粒度））## ----图，图。cap =“ Rahmann等人在文章中对图1的复制。情节派（distlist [[50]]）## ----截止------------------------------------------------------------------------------------------------------------------------------------------------------------ AB5 <-ScoreistCutoffs（diStlist [[50]，500，1，C = 1，0.05）## 5％FDR AB5 $ CUTOFFA ## 5％FNR AB5 $ catoffb ## fdr = fnr ab5 $ cutoffopt