model_limma: Model limma

In latlio/tidyde: Tidy Differential Expression

Description

'model_limma()' is a wrapper function for 'lmFit'. Fits a linear model for each gene.

Usage

model_limma(.data, .f = limma::lmFit, ...)

Arguments

.f	limma::lmFit
...	additional arguments to .f
voom	a voom object

Details

Please refer to lmFit for more information.

Value

a 'list' ('MArrayLM') object containing the results of the fit

Examples

counts <- readr::read_delim("data/GSE60450_Lactation-GenewiseCounts.txt", delim = "\t")
meta <- readr::read_delim("data/SampleInfo_Corrected.txt", delim = "\t") %>%
  mutate(FileName = stringr::str_replace(FileName, "\\.", "-"))

id <- as.character(counts$EntrezGeneID)
check_sample_names(counts, c(1,2), meta, FileName) %>%
  purrr::pluck("mod_count") %>%
  filter_genes(., id, "edgeR") %>%
  make_voom(., my_design) %>%
  model_limma()

latlio/tidyde documentation built on Dec. 21, 2021, 9:40 a.m.