我有一个逻辑回归模型,我用它来预测帝王蟹成熟时的大小,但是我在使用引导包设置引导代码时遇到了麻烦。这就是我所拥有的:
#FEMALE GKC SAM#
LowerChatham<-read.table(file=file.choose(),header=TRUE)
#LOGISTIC REGRESSION FIT#
glm.out<-glm(Mature~CL,family=binomial(link=logit),data=LowerChatham)
plot(Mature~CL,data=LowerChatham)
lines(LowerChatham$CL,glm.out$fitted,col="red")
title(main="Lower Chatham")
summary(glm.out)
segments(98.9,0,98.9,0.5,col=1,lty=3,lwd=3)
SAM<-data.frame(CL=98.97)
predict(glm.out,SAM,type="response")
我想引导统计 CL=98.97 因为我对 50% 的螃蟹成熟的大小感兴趣,但我不知道如何设置我的函数来指定该统计数据,更不用说引导函数了获得我的 95% CI 任何帮助将不胜感激!谢谢!