タラバガニの成熟時のサイズを予測するために使用しているロジスティック回帰モデルがありますが、ブート パッケージを使用したブートストラップ用のコードの設定に問題があります。これは私が持っているものです:
#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")
カニの 50% が成熟するサイズに関心があるため、統計 CL=98.97 をブートストラップしたいと思いますが、関数を設定してその統計を指定する方法がわかりません。一般的なブートストラップ関数は言うまでもありません。私の 95% CI を取得するには、どんな助けでも大歓迎です! ありがとう！