make design data problem with RDHFHet

[darryl]

Hi All,
I've finally launched into the great unknown of RMark, and wonder why I didn't start a long time ago! Thanks for this Jeff. However, when I try to create the design data for the above model type I get the error back:

> snails.proc<-process.data(snails,model="RDHFHet",time.intervals=time.intervals,groups="site")
> snails.ddl<-make.design.data(snails.proc)

Error in rbind(design.data, add.design.data) :
number of columns of matrices must match (see arg 2)

If I change model type (eg RDFullHet or RDHHet) everything goes fine, it seems specific to this model type. Am I missing a newbie step or is there a problem with this particular model type?

Cheers
Darryl

[darryl]

Answering my own question here for prosperity. Jeff contacted me off list advising there was a slight error in the parameter definition file that RMark uses for this model (and provided me with the correction). I'm sure Jeff will implement this with the next update of RMark.