simple vs. complex models

questions concerning anlysis/theory using program DENSITY and R package secr. Focus on spatially-explicit analysis.

simple vs. complex models

Postby bear16 » Tue Aug 13, 2013 4:31 pm

I am working on determining how to handle some of my models. I found that the session with trap specific behavior model (g0~session+bk, sigma~session) ranks the highest for a number of groups.

However, when I run (g0~session, sigma~session), I get maximization errors. (Different errors for different groups.) One example is
Warning messages:
1: In secr.fit(M, model = list(g0 ~ session, sigma ~ session), :
probable maximization error: optim returned convergence 1. See ?optim
2: In secr.fit(M, model = list(g0 ~ session, sigma ~ session), :
at least one variance calculation failed

I do not find any glaring errors in the results of the more complex model (i.e. parameters appear to be estimable and within reasonable range). Is there any reason why the more complex model would work while the less complex does not work? Or should I assume that there is likely an error with the more complex model, even though it does not come up with warnings?

Is there a source someone might recommend to read about this issue further?

Thank you for your help!
bear16
 
Posts: 20
Joined: Thu Oct 14, 2010 11:18 am

Re: simple vs. complex models

Postby murray.efford » Thu Aug 22, 2013 1:49 pm

This behaviour does look screwy... Make sure you have the latest version of 'secr'. Otherwise it is a matter of poking around, trying different fitting algorithms (e.g., method = 'Nelder-Mead') to check the estimates are stable etc. If you also want session-specific density estimates there is (almost) no merit in fitting a full g0~session, sigma~session model: you are better to simply fit each session-specific model separately - it's much faster and will lead you to the session with inadequate data or model fitting problems.
Murray
murray.efford
 
Posts: 686
Joined: Mon Sep 29, 2008 7:11 pm
Location: Dunedin, New Zealand


Return to analysis help

Who is online

Users browsing this forum: No registered users and 11 guests