re-doing GLM for FP richnesss with new data
updating glm() for each response variable (FP, LM, SP, W)