Update the b fixed effect to minimise the b-c random effect correlation
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'])
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.