Problem fixing gamma parameters in RMark

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

Problem fixing gamma parameters in RMark

Postby mks » Thu Dec 13, 2012 3:00 pm

Hi,

I'm a newcomer to RMark. I'm using the Robust design and I've been trying to set up a "no emigration" model (gammas=0) in RMark to compare with Markovian and random emigration models. I've tried fixing gamma double prime and gamma prime to 0, following examples in the RMark chapter and in a similar thread from the forum (viewtopic.php?f=21&t=2265&p=6997&hilit=fixed+parameters#p6994); I used the following formulae for the gammas:

GammaPrime.zero=list(formula=~1,fixed=0)
GammaDoublePrime.zero=list(formula=~1,fixed=0)

When I run my code with these formulae the Beta value for gamma double prime (intercept) equals 0, as expected, but the output for the real parameters for the gammas is blank and does not contain 0's. Here is a snippet of the output for gamma double prime:

Real Parameter GammaDoublePrime
Group:initialage0.sexFemale
2008 2009
2008
2009

Group:initialage1.sexFemale
2008 2009
2008
2009

Group:initialage2.sexFemale
2008 2009
2008
2009

I've tried debugging but I'm at a standstill. Any suggestions as to why the output is blank would be greatly appreciated.

-MKS
mks
 
Posts: 5
Joined: Thu Dec 13, 2012 12:49 am

Re: Problem fixing gamma parameters in RMark

Postby jlaake » Thu Dec 13, 2012 3:17 pm

The default for the display of fixed values is to show a blank. See the argument show.fixed in the help for summary.mark. If you set it to TRUE then the fixed values will be shown.

--jeff
jlaake
 
Posts: 1479
Joined: Fri May 12, 2006 12:50 pm
Location: Escondido, CA

Re: Problem fixing gamma parameters in RMark

Postby mks » Thu Dec 13, 2012 5:19 pm

Got it! Thanks, Jeff.
mks
 
Posts: 5
Joined: Thu Dec 13, 2012 12:49 am


Return to RMark

Who is online

Users browsing this forum: No registered users and 1 guest

cron