# #设置,包括= FALSE --------------------------------------------------- knitr: opts_chunk美元集(崩溃= TRUE,发表评论 = "#>" ) ## ---- getPackage eval = FALSE --------------------------------------------------- # 如果(!requireNamespace(“BiocManager”,悄悄地= TRUE)) # install.packages (BiocManager) # BiocManager::安装(“PAIRADISE ") ## ---- getDevel eval = FALSE ----------------------------------------------------- # BiocManager::安装(“hubentu / PAIRADISE ") ## ---- 加载消息= FALSE ------------------------------------------------------ 库(PAIRADISE ) ## ----------------------------------------------------------------------------- 库(abind) icount < -矩阵(1:4,1)scount < -矩阵(8,1) account <- abind(icount, scount, along = 3) account design <- data.frame(sample = rep(c("s1"， "s2")， 2)， group = rep(c("T"， "N")， each =2)) lens <- data.frame(sLen=1L, iLen=2L) PDseDataSet(account, design，镜头 ) ## ----------------------------------------------------------------------------- sample_dataset数据(“sample_dataset”) ## ----------------------------------------------------------------------------- 数据(“sample_dataset_CEU”)(“sample_dataset_LUSC数据 ") ## ----------------------------------------------------------------------------- pdat < - PDseDataSetFromMat pdat (sample_dataset) ## ----------------------------------------------------------------------------- pairadise_output < -pairadise (pdat numCluster = 2 ) ## ----------------------------------------------------------------------------- res < -结果(pairadise_output p.adj =“黑洞”,sig.level = 0.01) res ## ----------------------------------------------------------------------------- res < -结果(pairadise_output细节= TRUE) colnames (res) res潜[3美元 ,] ## ----------------------------------------------------------------------------- sessionInfo ()