Type II功能响应曲线随h、a参数变化绘图结果异常求助

核心错误点 #

• 公式运算优先级错误：你编写的a*x / 1 + a*h*x会按照R的运算顺序先算除法再算加法，实际等价于(a*x/1) + a*h*x，完全不符合Type II功能响应的定义。正确公式需要给分母整体加括号：a*x / (1 + a*h*x)
• 变量取值范围不符合生物学意义：功能响应中x代表猎物密度，不可能为负数，你设置的x范围从-5开始没有实际意义，还会导致曲线出现不符合预期的负数值部分。
• y轴范围设置不合理：修正后的Type II功能响应曲线值域在0到1/h之间，你设置的-10到10的范围包含大量无效区间，不利于观察曲线形态。

修正后的完整代码 #

# 定义修正后的Type II功能响应函数
fr <- function(x, a = 1, h = 1){
  a*x / (1 + a*h*x)
}

# 绘制不同处理时间h对应的曲线
plot(NA, xlim=c(0, 10), ylim=c(0, 1.2), 
     xlab = "猎物密度", ylab = "捕食者摄食率",
     main = "不同处理时间下的Type II功能响应曲线")
for (h in 1:5){
  curve(expr = fr(x, a=1, h=h), from = 0, to = 10, n = 100, add= TRUE, col = h, lwd=2)
}
legend("topright", legend = paste0("h=",1:5), col=1:5, lwd=2)

# 绘制不同攻击率a对应的曲线
plot(NA, xlim=c(0, 10), ylim=c(0, 1.2), 
     xlab = "猎物密度", ylab = "捕食者摄食率",
     main = "不同攻击率下的Type II功能响应曲线")
for (a in 1:5){
  curve(expr = fr(x, a=a, h=1), from = 0, to = 10, n = 100, add= TRUE, col = a, lwd=2)
}
legend("bottomright", legend = paste0("a=",1:5), col=1:5, lwd=2)

修正后效果说明 #

修正后曲线会呈现典型的Type II特征：随着猎物密度升高，摄食率先快速上升，之后增速放缓，最终趋近于最大值1/h，完全符合功能响应的生物学定义。

内容的提问来源于stack exchange，提问作者Alison Meeth