R help - Apr 2008 - fix variance parameter values for lmer estimation

Dear list-members,

a model is fit with lmer, but I want to force the variance parameter values
to be as defined by me

I thought, use 'start' to specify initial values and only allow for one
iteration ?
my question is how to do that ?

to specify the 2x2 matrix of variance parameter values:
start=list(groups=array(2,-0.5,1),dim=c(2,2))

now I need to make sure the mean structure is estimated, while the variance
parameter values are kept constant ??
e.g. lmer(resp~pred1+pred2+(1 +
pred1|groups),start=list(groups=array(2,-0.5,1),dim=c(2,2)),....)


cheers,

 
Wilfried
 


Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm

On 4/25/08, Wilfried Cools <Wilfried.Cools at ped.kuleuven.be>
wrote:
> Dear list-members,
>  a model is fit with lmer, but I want to force the variance parameter
values
>  to be as defined by me
>  I thought, use 'start' to specify initial values and only allow
for one
>  iteration ?
>  my question is how to do that ?
>  to specify the 2x2 matrix of variance parameter values:
>  start=list(groups=array(2,-0.5,1),dim=c(2,2))
Well, the problem is that if you are going to decide what the
arguments to the function should be and how they should be interpreted
you are going to need to write the function yourself.  There is
nothing in the lmer documentation that says it will fit a model with
the variance parameter values as defined by you and there is nothing
to say that you can specify a list with arbitrary names as the value
of the start parameter and somehow lmer will magically infer what you
mean by these starting values and do the right thing.  Sorry but I
don't write mind-reading software.  You are going to have to write it
yourself.

>  now I need to make sure the mean structure is estimated, while the
variance
>  parameter values are kept constant ??
>  e.g. lmer(resp~pred1+pred2+(1 +
>  pred1|groups),start=list(groups=array(2,-0.5,1),dim=c(2,2)),....)
>
>
>  cheers,
>
>
>  Wilfried
>
>
>
>  Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
>
>  ______________________________________________
>  R-help at r-project.org mailing list
>  https://stat.ethz.ch/mailman/listinfo/r-help
>  PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
>  and provide commented, minimal, self-contained, reproducible code.
>