Page 1 of 1

Problem fixing gamma parameters in RMark

PostPosted: Thu Dec 13, 2012 3:00 pm
by mks
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

Re: Problem fixing gamma parameters in RMark

PostPosted: Thu Dec 13, 2012 3:17 pm
by jlaake
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

Re: Problem fixing gamma parameters in RMark

PostPosted: Thu Dec 13, 2012 5:19 pm
by mks
Got it! Thanks, Jeff.