RMark vs. MARK model likelihood error

posts related to the RMark library, which may not be of general interest to users of 'classic' MARK

RMark vs. MARK model likelihood error

Postby jwarren » Fri Jun 05, 2009 12:20 am

When I append model results from RMark to MARK for running GOF tests I receive the following error message: "Have you changed the data, because the saturated model's likelihood has changed!". I also noticed that the model deviance in MARK is considerably lower than the deviance in RMark for the same model. Any ideas/assistance would be greatly appreciated.
jwarren
 
Posts: 11
Joined: Tue Nov 25, 2008 11:33 am
Location: Lakeview, Montana

Appending RMark models to MARK results browser

Postby dhewitt » Fri Jun 05, 2009 1:56 am

I just dealt with this same issue, I think. Are you using both export.chdata() and export.model() to get both the input file and the model results for appending in MARK? Chances are you need to because the original INP you imported into RMark may not be appropriate for building the MARK project to go with the RMark model. If you're using the current R and RMark versions, export.chdata() won't work right now (which led me to just use my original INP to build the MARK project, but that doesn't work). Jeff Laake just patched the function for me, so if this is your problem you might contact him for it (or I can pass it along).

The reason is that the ordering of groups is set by different defaults in RMark than MARK, so you need to export the processed data from RMark to pull into MARK so that the INP matches the model structure.

If this isn't your problem, I wasted a lot of keystrokes.
dhewitt
 
Posts: 150
Joined: Tue Nov 06, 2007 12:35 pm
Location: Fairhope, AL 36532

Postby jwarren » Fri Jun 05, 2009 9:35 pm

You're exactly right in how I'm using the export.chdata() and export.model() functions. Thanks for the help, I'll check in with Jeff Laake.
jwarren
 
Posts: 11
Joined: Tue Nov 25, 2008 11:33 am
Location: Lakeview, Montana


Return to RMark

Who is online

Users browsing this forum: No registered users and 1 guest

cron