如何解决为什么在进行参数自举拟合优度检验时我不能做超过 nsim=99 的事情？
这是我的适合度测试：
(fm <- distsamp(~1 ~ndvi_avg,keyfun="hazard",umf))
# Function returning three fit-statistics.
fitstats <- function(fm) {
observed <- getY(fm@data)
expected <- fitted(fm)
resids <- residuals(fm)
sse <- sum(resids^2)
chisq <- sum((observed - expected)^2 / expected)
freeTuke <- sum((sqrt(observed) - sqrt(expected))^2)
out <- c(SSE=sse,Chisq=chisq,freemanTukey=freeTuke)
return(out)
}
(pb <- parboot(fm,fitstats,nsim=100,report=1))
我最多可以做 nsim=99，然后我尝试 nsim=100，结果如下：
> (pb <- parboot(fm,report=1))
t0 = 103.385 540.1714 70.38591
Running in parallel on 3 cores. Bootstrapped statistics not reported.
解决方法
在 parallel = FALSE 中尝试 parboot()。
文档：https://www.rdocumentation.org/packages/unmarked/versions/1.1.0/topics/parboot
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。
