Breast cancer classification using limma (Bioconductor)

knitr::opts_chunk$set(fig.align="center", out.width="80%", echo = TRUE)

install.packages("tidyverse", quietly=T)
install.packages("gplots", quietly=T)
install.packages("readxl", quietly=T)
install.packages("gplots", quietly=T)

install.packages("BiocManager", quietly=T)
BiocManager::install("limma", quietly=T, force=TRUE)
BiocManager::install("BiocGenerics", quietly=T, force=TRUE)
BiocManager::install("Biobase", quietly=T, force=TRUE)

library("tidyverse", quietly=T)
library("parallel", quietly=T)
library("limma", quietly=T)
library("readxl", quietly=T)
library("BiocGenerics", quietly=T)
library("gplots", quietly=T)
library("Biobase", quietly=T)

url_file <- "https://github.com/friedrfe/breast_cancer_protein_abundance/raw/main/20211212_SubtypeExample.xlsx"
download.file(url_file, "file.xlsx")

e <- read_xlsx("file.xlsx")

head(e)

glimpse(e)

summary(e)

sapply(e, function(x) sum(is.na(x)))