如何在R中更改此图的颜色和子标题（冬季，夏季）？

Question

glmtom  <- glm(Organic.Matter ~ Season*Bay, data = nuttom)
tefftom <- allEffects(glmtom)
plot(tefftom, ylab = "Mean Organic Matter [%]", main = "GLM Organic Matter")