|
12 | 12 | #' @param filename_binom Name of the exported binomial file |
13 | 13 | #' @param filename_sigdeg Name of the exported sigDEG file |
14 | 14 | #' @importFrom stats pbinom median |
15 | | -#' @import org.Hs.eg.db |
16 | 15 | #' @rdname ClustDiffGenes |
17 | 16 | #' @return A list containing two tables. |
18 | 17 | #' @export |
@@ -103,50 +102,29 @@ setMethod( |
103 | 102 | DEGsS <- vector() |
104 | 103 | DEGsTable <- data.frame() |
105 | 104 |
|
106 | | - for (n in 1:K) { |
107 | | - if (length(cdiff[[n]][, 1]) == 0) { |
108 | | - next |
109 | | - } |
| 105 | + # Generate gene list |
| 106 | + geneList <- retrieve_geneList() |
| 107 | + GL <- c(1, "MTRNR2", "ENSG00000210082") |
| 108 | + GL1 <- c(1, "MTRNR1", "ENSG00000211459") |
| 109 | + geneList <- rbind(geneList, GL, GL1) |
| 110 | + gene_list <- geneList[, 3] |
110 | 111 |
|
| 112 | + for (n in 1:K) { |
111 | 113 | if (length(cdiff[[n]][, 1]) > 0) { |
112 | 114 | p.adj <- p.adjust(cdiff[[n]][, 4], method = "bonferroni") |
113 | 115 | out <- cbind(cdiff[[n]], p.adj) |
114 | 116 | out <- subset(out, out[, 5] < fdr) |
115 | 117 | if (length(out[, 1]) > 0) { |
116 | 118 | Regulation <- vector() |
117 | | - for (i in seq_len(length(out[, 1]))) { |
| 119 | + for (i in seq_along(out[, 1])) { |
118 | 120 | if (out[i, 1] > out[i, 2]) { |
119 | 121 | Regulation[i] <- "Down" |
120 | 122 | } else { |
121 | 123 | Regulation[i] <- "Up" |
122 | 124 | } |
123 | 125 | } |
124 | 126 | out <- cbind(out, Regulation) |
125 | | - if (quiet) { |
126 | | - suppressMessages( |
127 | | - geneList <- |
128 | | - AnnotationDbi::select( |
129 | | - org.Hs.eg.db, |
130 | | - keys = keys(org.Hs.eg.db), |
131 | | - columns = c("SYMBOL", "ENSEMBL") |
132 | | - ) |
133 | | - ) |
134 | | - GL <- c(1, "MTRNR2", "ENSG00000210082") |
135 | | - GL1 <- c(1, "MTRNR1", "ENSG00000211459") |
136 | | - geneList <- rbind(geneList, GL, GL1) |
137 | | - } else { |
138 | | - geneList <- |
139 | | - AnnotationDbi::select( |
140 | | - org.Hs.eg.db, |
141 | | - keys = keys(org.Hs.eg.db), |
142 | | - columns = c("SYMBOL", "ENSEMBL") |
143 | | - ) |
144 | | - GL <- c(1, "MTRNR2", "ENSG00000210082") |
145 | | - GL1 <- c(1, "MTRNR1", "ENSG00000211459") |
146 | | - geneList <- rbind(geneList, GL, GL1) |
147 | | - } |
148 | 127 | genes <- rownames(out) |
149 | | - gene_list <- geneList[, 3] |
150 | 128 | idx_genes <- is.element(gene_list, genes) |
151 | 129 | genes2 <- geneList[idx_genes, ] |
152 | 130 | Final <- cbind(genes, out) |
@@ -233,6 +211,7 @@ setMethod( |
233 | 211 | } |
234 | 212 | } |
235 | 213 | } |
| 214 | + |
236 | 215 | if (length(DEGsTable) > 0) { |
237 | 216 | colnames(DEGsTable) <- c( |
238 | 217 | "Target Cluster", "VS", "Gene number", "File name", |
|
0 commit comments