我阅读了以下文章:MAYFIELD LOGISTIC RERESSION: A PRACTICAL 巢生存分析的方法,它提供了一个 SAS 代码(如下)来拟合逻辑模型:
proc logistic data = {data set};
model FAIL/OBSDAYS = MIDHT SNAGBA
VERTDENS;
run;
在此模型中,目标是估计巢穴的每日存活率。因此,响应变量 (FAIL) 是巢的命运(0=成功,1=失败),OBSDAYS 是巢的暴露时间,MIDHT、SNAGBA 和 VERTDENS 是协变量。我完全理解模型的第二部分,但我对在 R 中配置此模型中的响应变量有疑问。在 R 中设置如下是否合适?:
m1<-glm(fail~MIDHT+SNAGBA+VERTDENS, data=data set, family="binomial")