library(GEOquery)

#### Set the GEO accession number
geo_acc_numbers <- "GSE77193"
groups_of_interest <- c("Mock", "WNVWT")

# Download and load the GEO dataset
gse <- getGEO(geo_acc_numbers, GSEMatrix = TRUE)

expression_data <- as.data.frame(exprs(gse[[1]]))
group_data <- pData(phenoData(gse[[1]]))

gse <- gse[[1]]
mAnnot <- fData(gse)
mAnnot <- mAnnot[c("ID", "GENE_SYMBOL")]