www.phidot.org

Hi all,

Is there any particular reason why Nest survival models are not implemented in RMARK?
I have been adopting the code from functions make.mark.model and run.mark model in the package, and managed to successfully run some models while using identity matrix and logit or sin links.
However, when i try to repeatedly run the same code for different samples
i sometimes get:
"
forrtl sever (161): Program exception -array bounds exceeded
"
and a whole bunch of "Mark.exe" 's and a lot of other information that seems Greek to me.

Does any one know what happens?
I suspect that this is one of the reasons why nest survival models are not (yet?) implemented in RMARK.
I can send example code and datasubset for anyone interested in testing

Any help is much appreciated
Many thanks

Caspar

OK... i figured out its not a RMARK problem, nor with the code etc. Mark fails to run this particular datasets alltogether. i.e. resulting in "Program exception - array bounds exceeded"

Anyone knows what this means? Or how to overcome?

Thanks

Caspar

If you're excited to run nest survival models in R, why not just implement a "Shaffer model" and do the analyses with GLMs? The correct link function is built in to the latest R versions, so the modeling process is about as difficult as performing any other regression. The form of the data needs to be a little different; other than that, it's a quick job. Program MARK is not needed for the GLM modeling approach. There are pros and cons to each approach; look in the recent ornithological literature for several papers addressing this issue.

With Jay Rotella's help I've coded the Nest survival model into RMark and in the process of testing and writing up some help. So far appears to work fine. It has a very different structure than most MARK models so it has taken some specialized coding.

My guess about the error you encountered would be that you have some of the recorded days that are greater than the number of occasions. Don't know for sure but MARK probably sets up an array based on number of occasions and then if you had a day value outside that limit it would give the array exception. I'll include that check into the RMark code.