我有两个数据集,我已经进行了逻辑回归。当我使用predict()
绘制它们时,manip.model
明显低于non.manip
模型,但是它们在x = 2.3处相交。我想做一个排列测试,看看它们在任何地方交叉的频率,或者这只是一个神器,但我无法弄明白。到目前为止,我只进行了非常简单的排列测试,并且遇到了这个更复杂的问题。
manip.model <- glm(catch~traplen, family=binomial, data=manip.data)
nonmanip.model <- glm(catch~traplen, family=binomial,
data=nonmanip.data)
predict(manip.model,data.frame(traplen=x),type="resp")
predict(nonmanip.model2,data.frame(preymass=x),type="resp")
编辑:这是捕获成功由0或1表示的数据,并且处理由y或n表示。我试图确定治疗n的预测逻辑回归线在[1,2.3]范围内低于治疗y的频率,范围为(2.3,3)。
df <- data.frame(traplen=sample(1.0:3.0, 10, replace=T),
manip=sample(c('y', 'n'), 10, replace=T), catch=sample(c('0','1'), 10,
replace=T))
df