我正在使用属于不同观察结果的动物行为数据来构建GLM,以识别具有不同协变量(例如距离和运动方向)的相机陷阱。由于我需要确保数据点之间的独立性,因此我每次观察均取十只动物的随机样本,并将模型拟合到该数据子集上。当然,这些模型的结果之间会有细微的差异,这就是为什么我想在我的结果中显示所有结果,然后将它们组合成一个图。
我使用效果包为单个模型制作图形:
library(effects)
plot(predictorEffect("Flash",D2_Fitdetected_1),
main="Effect of flash type on camera trap recognition",
axes=list(
y=list(type="response", grid=T,lab="Probability of camera trap detection"),
x=list(rotate=30),
grid=T),
lines=list(multiline=T,col="red", lty=1),
symbols=list(pch=16,cex=1.5))
是否有可能在一个坐标系中合并这些相同变量的多个图?
我也将欣赏其他想法,该想法如何以紧凑的方式基于数据的不同子集显示这些模型的效果之间的差异。