r - 置换测试,看两个glm是否相交

时间:2018-01-11 15:42:53

标签: r permutation

我有两个数据集,我已经进行了逻辑回归。当我使用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

0 个答案:

没有答案