One of the model fits is failing, most likely because the model is too complicated given the underlying data. But when you run all the models at once with dredge, it's impossible to tell which one specifically isn't working. You'd have to run the models one-by-one to identify the issue.
In general we recommend slowly building up the complexity of the model iteratively, so you can more easily identify the reasons for this type of error. Also, I think it's better practice to identify a candidate set of models specifically based on your hypotheses/predictions, and fit only those models, rather than using a 'dredging' approach.
Ken
On Wed, Oct 02, 2024 at 08:01:25AM -0700, Devavrat Pawar wrote:
> Hello everyone,
>
> I have the same error which Prashant has mentioned in the previous email
> while using the Dredge function. It would be great if anyone could clarify
> this.
>
> *Error in h(simpleError(msg, call)) :*
>
> *error in evaluating the argument 'x' in selecting a method for function
> 'diag': Lapack routine dgesv: system is exactly singular: U[1,1] = 0*
>
> *In addition: There were 42 warnings (use warnings() to see them)*
>
>
> Thank you
> Devavrat
>
> On Tuesday 28 February 2023 at 14:56:08 UTC+5:30 Prashant Mahajan wrote:
>
> > Hello everyone,
> >
> > I am running single-season occupancy models in R using unmarked. I have
> > 100 sites with 6 replicates each and have some missing values in them which
> > I have replaced with NA. I have multiple continuous and categorical
> > covariates. Therefore, I am using dredge function of package "MuMIn".
> > However, while running the global model with all the covariates using
> > dredge function, I am getting the following error:
> >
> > *Error in h(simpleError(msg, call)) :*
> >
> > *error in evaluating the argument 'x' in selecting a method for function
> > 'diag': Lapack routine dgesv: system is exactly singular: U[1,1] = 0*
> >
> > *In addition: There were 42 warnings (use warnings() to see them)*
> >
> > Even while running the constant model in dredge I am getting this error:
> >
> > *Warning message:*
> >
> > *In sqrt(diag(vcov(model, ...))) : NaNs produced*
> >
> > Can someone please suggest me the solution for this. Is is because of
> > missing observations or something else.
> >
> > Best regards
> >
> > Prashant
> >
> >
> >
> >
> >
> >
>
> --
> *** 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/3ea27e2e-a032-46d8-9676-cdcf32f623dfn%40googlegroups.com.