bupdate | R Documentation |
A function to update the value of bstart
, the starting value for the
b fixed effect, to minimise the correlation between the random effects b and
c.
bupdate(x)
x |
a |
Returns an updated value of the b fixed effect, based on the random effect covariance matrix.
Tim Cole tim.cole@ucl.ac.uk
## fit sitar model with b fixed effect starting value defaulting to 'mean'
m1 <- sitar(x=age, y=height, id=id, data=heights, df=5)
print(fixef(m1)['b'])
## refit with starting value chosen to minimise b-c correlation and df increased
m2 <- update(m1, bstart=bupdate(m1), df=6)
print(fixef(m2)['b'])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.