※上記の広告は60日以上更新のないWIKIに表示されています。更新することで広告が下部へ移動します。

NB <- 1000
scpue.boot <- matrix(0,length(scpue),NB)

for(i in 1:NB){
sres3.boot <- glm.nb(number~factor(year)  +
                depth2*factor(Region), data=sdata[sample(1:nrow(sdata),nrow(sdata),replace=T),])
#dataからsampleしているsample(サンプリングの範囲(この場合は行番号),個数,繰り返しの有無)

pres3.boot <- predict(sres3.boot,newdata=ndata,se=TRUE,
                       type="response")

scpue.boot[,i] <- tapply(pres3.boot$fit,ndata$year,mean,na.rm=T)
}


par(mfrow=c(1,1))

#図の書き方もちょっと勉強
matplot(t(apply(scpue.boot,1,quantile,probs=c(0.05,0.5,0.95))),type=c("l","b","l"),
col=1,lty=c(2,1,2),pch=1,xlab="CPUE")
|