vignettes/test_data/new_workflow.R

library(geometric2)

setwd("more_testing/")
data("biosamps")
gse <- "GSE6346"
get_input(gse)

exprs <- get_eset(gse)

samples <- fread("sample_review.txt")


my_design <- get_design()

fit <- lmFit(exprs, my_design)

fit <- fitting()

contrasts <- get_contrast.matrix()
fit2 <- contrasts.fit(fit, contrasts)
fit3 <- eBayes(fit2)
axelmuller/geometric2 documentation built on May 25, 2019, 6:26 p.m.