[R] differential gene expression

[Anas Jamshed]

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]]