我有一个 lmer 模型,想用点数据绘制效果趋势,但是函数 plot_model() 虽然很擅长绘制我想要的东西,但在定制方面似乎很不灵活。
我想更改以下图中点的颜色和 alpha:
model1<-lmer(NDVI~ fall_ai+spring_ai+SI+(1|colony),data=clm,REML=T)
plot_model(model1, type = "pred", show.data = T)
一种解决方法是使用 lattice 来做到这一点:
plot(Effect(c("fall_ai"),model1),
colors="black",lwd=1.5,ylab="Spring NDVI",xlab="Autumn Aridity Index",
ci.style="lines",main="",ylim=c(0.20,0.55),xlim=c(0,45))
trellis.focus("panel", 1, 1, highlight = FALSE)
print(panel.points(clm$fall_ai,pch = 16, cex = 0.7, col = "black"))
trellis.unfocus()