[R] differential gene expression

[Duncan Murdoch]

On 21/08/2021 6:55 p.m., Anas Jamshed wrote:
> 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

I think you're asking in the wrong place.  You need statistical help, 
not R help.  Once you know what you want to do, ask here about how to do 
it; in the meantime, look for statistical help.  Some people suggest 
stats.stackexchange.com for that; I haven't read it enough to know if it 
is useful or not.  In many cases, talking to a local expert is the best 
solution.

Duncan Murdoch