@@ -956,29 +956,29 @@ color_mirtop <- if (params$input_type == 'smrnaseq') "#f8c383" else "#f0f0f0"
956956# Generate DOT string with conditional color
957957dot_code <- sprintf('
958958digraph flowchart {
959- graph [layout = dot, rankdir = LR]
960-
961- node [shape = box, style = filled, fontname = Helvetica, fontcolor = black]
962-
963- A [label = "Primary Bioinformatics\\nnf-core/rnaseq", fillcolor = "%s"]
964- B [label = "Primary Bioinformatics\\nnf-core/smrnaseq", fillcolor = "%s"]
965- C [label = "Quantification", shape = diamond, fillcolor = "#f0f0f0", fontcolor = "#333333" width = 1.5, height = 1.5 ]
966- D [label = "featureCounts", fillcolor = "%s", fontcolor = "#333333"]
967- E [label = "STAR-Salmon", fillcolor = "%s", fontcolor = "#333333"]
968- F [label = "RSEM", fillcolor = "%s", fontcolor = "#333333"]
969- G [label = "miRTop", fillcolor = "%s", fontcolor = "#333333"]
970- H [label = "Secondary Bioinformatics\\nqbic-pipelines/rnadeseq", fillcolor = "#2b64b9", fontcolor=white]
971-
972- A -> C
973- B -> C
974- C -> D
975- C -> E
976- C -> F
977- C -> G
978- D -> H
979- E -> H
980- F -> H
981- G -> H
959+ graph [layout = dot, rankdir = LR]
960+
961+ node [shape = box, style = filled, fontname = Helvetica, fontcolor = black]
962+
963+ A [label = "Primary Bioinformatics\\nnf-core/rnaseq", fillcolor = "%s"]
964+ B [label = "Primary Bioinformatics\\nnf-core/smrnaseq", fillcolor = "%s"]
965+ C [label = "Quantification", shape = diamond, fillcolor = "#f0f0f0", fontcolor = "#333333" width = 1.5, height = 1.5 ]
966+ D [label = "featureCounts", fillcolor = "%s", fontcolor = "#333333"]
967+ E [label = "STAR-Salmon", fillcolor = "%s", fontcolor = "#333333"]
968+ F [label = "RSEM", fillcolor = "%s", fontcolor = "#333333"]
969+ G [label = "miRTop", fillcolor = "%s", fontcolor = "#333333"]
970+ H [label = "Secondary Bioinformatics\\nqbic-pipelines/rnadeseq", fillcolor = "#2b64b9", fontcolor=white]
971+
972+ A -> C
973+ B -> C
974+ C -> D
975+ C -> E
976+ C -> F
977+ C -> G
978+ D -> H
979+ E -> H
980+ F -> H
981+ G -> H
982982}
983983', color_nfcore_rnaseq, color_nfcore, color_featurecounts, color_salmon, color_rsem, color_mirtop)
984984
@@ -3013,18 +3013,18 @@ and if their log2FC was greater than `r logFC_text`. We compared the following c
30133013``` {r summary, warning=FALSE, message=FALSE, results='asis'}
30143014for (file in allgenes_files) {
30153015 # Reading DE genes list
3016- fname <- tools::file_path_sans_ext(basename(file))
3017- DE_genes <- read.csv(file = paste0(wd, "/differential_gene_expression/allgenes/", file), sep="\t", header = T)
3018- DE_genes$contrast <- rep(fname, nrow(DE_genes))
3019- num_all <- DE_genes %>% nrow()
3020- numDE <- DE_genes %>%
3021- filter(padj <= params$adj_pval_threshold, abs(log2FoldChange) >= params$logFC_threshold) %>%
3022- nrow()
3023-
3024- contrast_len <- sapply(strsplit(fname, "_"), length)
3025- current_condition <- sapply(strsplit(contrast_names[i], '_'), function(x) paste(x[2], collapse=' '))
3026- current_contrast_name <- sapply(strsplit(fname, '_'), function(x) paste(x[(contrast_len-2):contrast_len], collapse=' '))
3027- cat(paste("* Comparison:", current_condition, ":", current_contrast_name, "has", numDE, "DE genes out of", num_all, "tested genes.\n"))
3016+ fname <- tools::file_path_sans_ext(basename(file))
3017+ DE_genes <- read.csv(file = paste0(wd, "/differential_gene_expression/allgenes/", file), sep="\t", header = T)
3018+ DE_genes$contrast <- rep(fname, nrow(DE_genes))
3019+ num_all <- DE_genes %>% nrow()
3020+ numDE <- DE_genes %>%
3021+ filter(padj <= params$adj_pval_threshold, abs(log2FoldChange) >= params$logFC_threshold) %>%
3022+ nrow()
3023+
3024+ contrast_len <- sapply(strsplit(fname, "_"), length)
3025+ current_condition <- sapply(strsplit(contrast_names[i], '_'), function(x) paste(x[2], collapse=' '))
3026+ current_contrast_name <- sapply(strsplit(fname, '_'), function(x) paste(x[(contrast_len-2):contrast_len], collapse=' '))
3027+ cat(paste("* Comparison:", current_condition, ":", current_contrast_name, "has", numDE, "DE genes out of", num_all, "tested genes.\n"))
30283028}
30293029```
30303030
@@ -3180,7 +3180,7 @@ cat(paste0(
31803180
31813181The secondary analysis of differential gene expression was performed by the using the [rnadeseq pipeline](https://github.com/qbic-pipelines/rnadeseq), revision ",
31823182as.character(params$revision),"
3183- , which was written using the nf-core template [@ewels2020nf]. The read quantification data resulting from `",
3183+ , which was written using the nf-core template [@ewels2020nf]. The read quantification data resulting from `",
31843184quant_tool, "` was processed with the R package [DESeq2 v",
31853185packageVersion("DESeq2"),
31863186"](https://bioconductor.org/packages/release/bioc/html/DESeq2.html) [@love2014differential].
0 commit comments