我们项目的目标是在与栅栏互动时评估鹿的2种行为反应(穿越成功和穿越决定)。我们在我们的研究中使用了改进的BACI设计,我们在预修改期间监控围栏部分,然后在我们称之为修改后期间进行了某些修改。
我们使用逻辑回归来评估在这些围栏部分的骡鹿和白尾鹿的交叉成功和交叉决定。我们使用具有logit链接的广义线性模型(即逻辑回归)模拟某些人口统计因素和围栏特征对我们的2个响应变量的影响。
现在我们已经确定了两种分析的顶级模型,我们正在尝试创建一个毛虫图,以图形方式显示修改前期和修改后期间的交叉成功的变化,我们标题为& #34; Type_Period&#34 ;.
然而,我们无法弄清楚如何改变我们的情节(我已经附上)的规模。我们使用GGally,broom和ggplot2来创建Caterpillar图。
我们已经查看了多个GGally和ggplot2帮助页面,但我还没有找到任何关于如何操纵比例的提及。我们已经尝试过"缩放"在图上但ggplot2无法读取命令。我们尝试更改xlim,但收到一条错误消息,表示"忽略未知参数xlim。"
这是我们仅用于毛虫图的脚本。如有必要,我将提供完整的脚本和示例数据库。
library(GGally)
library(broom)
library(ggplot2)
log.reg <- glm(Decision~ Type_Period, data=analysis, weights = analysis$Freq, x=TRUE,family=binomial(link="logit"))
ggcoef(log.reg, xlim=c(-1, 2))