covariate adjustment for metabolomics data

Hi,

I haven't used Python before and am completely new to its world. Can you please help me with the covariate adjustment for my features? I have 70 targets and ~500 features which I want to adjust for 3 covariates.

I tried by seeing the examples for resAdj but was not able to replicate it on my data. I am confused about how should I tell my code which is needed to be adjusted for what.

I would be very grateful for any help provided.

Thank you very much!!