如何拟合四参数逻辑回归(泊松)R

时间:2021-02-18 18:53:52

标签: r logistic-regression poisson

我正在尝试拟合四参数逻辑回归来模拟鸟类物种丰富度 (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)

有什么建议吗?谢谢

0 个答案:

没有答案