linear model for multiple response variables (1 to 30) and combination of one or two predictor variables ( i.e. different formula for linear model)

Question

I am relatively new in R programming and I want to run multiple linear models at a time along with diagnostic plots. For example:

cc is different chemical constituent