Skip to content

Commit

Permalink
Merge pull request #21 from yerkes-gencore/voom-edgeR
Browse files Browse the repository at this point in the history
`bulk$contrasts` -> `bulk$fit.contrasts` see ca6bf39 of #19
  • Loading branch information
derrik-gratz authored Feb 28, 2024
2 parents 4421b93 + 4d8e8f1 commit e896ee7
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions R/fitModels.R
Original file line number Diff line number Diff line change
Expand Up @@ -57,10 +57,10 @@ fitVoomLm <- function(bulkObj, contr.matrix, block = NULL, sample.weights = TRUE
var.design = var.design,
var.group = var.group,
plot = plotVoom)
bulkObj$contrasts <- limma::contrasts.fit(bulkObj$fit, contrasts = contr.matrix)
bulkObj$contrasts <- limma::eBayes(bulkObj$contrasts, robust=TRUE)
limma::plotSA(bulkObj$contrasts, main="Final model: Mean-variance trend")
bulkObj$de <- limma::decideTests(bulkObj$contrasts, lfc = 0)
bulkObj$fit.contr <- limma::contrasts.fit(bulkObj$fit, contrasts = contr.matrix)
bulkObj$fit.contr <- limma::eBayes(bulkObj$fit.contr, robust=TRUE)
limma::plotSA(bulkObj$fit.contr, main="Final model: Mean-variance trend")
bulkObj$de <- limma::decideTests(bulkObj$fit.contr, lfc = 0)
return(bulkObj)
}

Expand Down

0 comments on commit e896ee7

Please sign in to comment.