#gp.st is a data frame object containing the lipid classes (columns) abundance for each individual sample (rows)
#mtdt2.cl is a data frame object containing samples' metadata (e.g: sex, maternal group, etc) organized in columns.
BiocManager::install("car")
library(car)
anv<-apply(gp.st,2,function(x){a<-cbind("lipid"= x,"sex"=mtdt2.cl$sex,"group"=mtdt2.cl$group)
a<-as.data.frame(a)
a$lipid<-as.numeric(as.character(a$lipid))
model=lm(lipid ~ group + sex + group:sex, data = a)
b<-Anova(model, type="II")
c<-b$`Pr(>F)`[1:3]
})
rownames(anv)<-c("group","sex","group:sex")
Do you have any questions about this protocol?
Post your question to gather feedback from the community. We will also invite the authors of this
article to respond.