parboot and ranef work for me with the example code in the occuMS help file. Do you have missing values in your y-matrix and/or covariates? I suspect these may not be being handled properly.
On Thu, Jul 21, 2022 at 05:07:43PM -0700, Benjamin Becker wrote:
> Dear List,
> Is there any guidance for calculating *c-hat *and *goodness of fit *from an
> *occuMS()* model? I see that parboot() does not support the condbinom
> parameterization and when I run the multinom parameterization I get:
> > parboot(mod04b_stacked_fit)
> Error in sample.int
(length(x), size, replace, prob) :
> NA in probability vector
> The model object works fine with predict(), SSE(), plot(), and gives
> reasonable errors, plots, and parameter results, but not ranef().
> > ranef(mod04b_stacked_fit)
> Mean Mode 2.5% 97.5%
> [1,] NA integer,0 NA NA
> [2,] NA integer,0 NA NA
> [3,] NA integer,0 NA NA
> Thanks for any guidance.
> *** Three hierarchical modeling email lists ***
> (1) unmarked (this list): for questions specific to the R package unmarked
> (2) SCR: for design and Bayesian or non-bayesian analysis of spatial capture-recapture
> (3) HMecology: for everything else, especially material covered in the books by Royle & Dorazio (2008), Kéry & Schaub (2012), Kéry & Royle (2016, 2021) and Schaub & Kéry (2022)
> You received this message because you are subscribed to the Google Groups "unmarked" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to unmarked+u...@googlegroups.com
> To view this discussion on the web visit https://groups.google.com/d/msgid/unmarked/50e90e60-8809-43a1-831c-34f66a3435c6n%40googlegroups.com