为什么在进行参数自举拟合优度检验时我不能做超过 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