I have downloaded data from: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE162562&fbclid=IwAR0iZQhttG8HzGhFIIMWbFgNszQrVDgiyVChYzQ_ypCx_d-1pn_tm7STjGs
and now I want to compare: healthy vs Mild healthy vs Highly exposed seronegative (ishgl) Healthy vs Asymptomatic covid19 patient healthy vs Highly exposed seronegative (non ishgl) from this data. I started like : library(edgeR) library(limma) library(GEOquery) library(Biobase) Sys.setenv("VROOM_CONNECTION_SIZE" = 131072 * 2) setwd("D:\\") untar("GSE162562_RAW.tar") filelist = list.files(pattern = ".*.txt.gz") But after getting text files I don't know how to proceed further. I want to find degs from these files *Plz help me * [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.