I want to run ANCOVAs for my data. I want to look at how peptide levels are differing between groups (I have three different groups) while adjusting for ethnicity.
This is the function I am using (What is the difference pf using the * or the + sign?):
ancova.mod1 <- lm(peptide ~ group + ethnicity, data = ancova.df)
The summary provides the Adjusted R square and p_trends, but I also want the estimated marginal means for each group, is there a way to extract this information?