我正在尝试拟合四参数逻辑回归来模拟鸟类物种丰富度 (Patch_Richness) 以响应森林覆盖 (FOREST500)。我需要将 km 作为协变量添加到模型(km = km),我正在尝试以下(代码)但不确定“km”是否正确包含在公式中,我只是在预测变量之后添加了 km( FOREST500)
logip=function(p,lambda,x){
a=p[1]
b=p[2]
c=p[3]
d=p[4]
Riq1 = d+(a/(1+exp((b-FOREST500+km)/c)))
-sum(dpois(x,lambda=Riq1, log=TRUE))
}
parnames(logip)=c("a","b","c","d")
modTR.log=mle2(minuslog=logip, start= c(a=30,b=30, c=3,d=20), data=list(x=Patch_Richness))
summary(modTR.log)
有什么建议吗?谢谢