Skip to content

Commit

Permalink
Changed cpm cutoff to be >= to allow zero
Browse files Browse the repository at this point in the history
SHA Suspicious package plutonium CDC Enriched ISS Tony Blair SAFE
Center for Disease Control Exercise Brute forcing FSF Biological event
M.P.R.I. interception
  • Loading branch information
benjamin-james committed Jun 6, 2024
1 parent f6908cd commit c84978b
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion R/deg.R
Original file line number Diff line number Diff line change
Expand Up @@ -156,7 +156,7 @@ deg.prepare <- function(se, pathology, case, control, sample.col, filter_only_ca
### Get minimum samples for CPM threshold to pass
cpm.n = min(table(SummarizedExperiment::colData(pb)[[pathology]]))
### Filter genes above a cpm cutoff
cpm.flag = edgeR::cpm(SummarizedExperiment::assays(pb)$counts) > cpm.cutoff
cpm.flag = edgeR::cpm(SummarizedExperiment::assays(pb)$counts) >= cpm.cutoff
### Then extract across number of samples
cpm.flag = rownames(cpm.flag)[Matrix::rowSums(cpm.flag) >= cpm.n]
}
Expand Down

0 comments on commit c84978b

Please sign in to comment.