I have this data: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE162562&fbclid=IwAR0iZQhttG8HzGhFIIMWbFgNszQrVDgiyVChYzQ_ypCx_d-1pn_tm7STjGs
and 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 : gset <- getGEO("GSE162562", GSEMatrix =TRUE, getGPL=FALSE) gset sampleinfo = pData(gset[[1]])[, 1:3] print(sampleinfo) but i need help to proceed futher [[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.