## ----设置，包括= false --------------------------------------------------------------------------------------------------------------------------------------------------------------------------- knitr::opts_chunk$set(fig.width = 7, fig.height = 7, echo = FALSE, warning = FALSE, message = FALSE, dev = 'png'，out.extra ='style =“边框宽度：0;”）## --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- rep1_df <- idr2d:::chiapet$rep1_df rep2_df <-IDR2D ::: CHIAPET $ REP2_DF ## -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------库（DT）标题<-htmltools :: withtags（表格）（class ='display'，thead（tr（th（colspan = 3，“锚a”），th（colspan = 3，“锚b”），th（rowspan = 2，“ fdr”）），tr（lapply（rep（c（“ chr。”，“启动坐标”，“ end坐标”），2），th））））））datatable（rep1_df [seq_len（min（nrow（rep1_df），1000）），]，容器=标题，rownames = false，option------------------------------------------------------------------------------------------ dataTable（rep2_df [seq_len（min（nrow（rep2_df），1000））），]，]，容器=标题，rownames = false，options = list（searching = false））％>％格式化（“ fdr”，3）## --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- - - - - - - - - - - - - - - - - - - - - - - - 图书馆（idr2d）## ------- echo = true ---------------------------------------------------------------------------------------------------------------------------------------------- idr_results <- estimate_idr2d(rep1_df, rep2_df, value_transformation = "log_additive_inverse") rep1_idr_df <- idr_results$rep1_df ## --------------------------------------------------------------------------------------------------------------- # avoid CRAN warnings rank <- rep_rank <- value <- rep_value <- idr <- NULL header <- htmltools::withTags(table( class = 'display', thead( tr( th("rank in R1"), th("rank in R2"), th("transformed value in R1"), th("transformed value in R2"), th("IDR") ) ) )) df <- dplyr::select(rep1_idr_df, rank, rep_rank, value, rep_value, idr) datatable(df[seq_len(min(nrow(df), 1000)), ], rownames = FALSE, options = list(searching = FALSE), container = header) %>% formatRound(c("value", "rep_value", "idr"), 3) ## ---- echo = TRUE------------------------------------------------------------- summary(idr_results) ## ---- echo = TRUE------------------------------------------------------------- draw_idr_distribution_histogram(rep1_idr_df) ## ---- echo = TRUE------------------------------------------------------------- draw_rank_idr_scatterplot(rep1_idr_df) ## ---- echo = TRUE------------------------------------------------------------- draw_value_idr_scatterplot(rep1_idr_df)