*Thanks for that, David.

> boxx <- box
> f0128bi = lmer(fincr ~ icfac + (1|gap) + (1|boxx)+ (1|gap:boxx), family =
binomial)
> f0128biu = update(f0128bi, . ~ . + (1|tree))
> anova(f0128bi,f0128biu)
Data:
Models:
f0128bi: fincr ~ icfac + (1 | gap) + (1 | boxx) + (1 | gap:boxx)
f0128biu: fincr ~ icfac + (1 | gap) + (1 | boxx) + (1 | gap:boxx) + (1 |
f0128biu:     tree)
         Df    AIC    BIC  logLik Chisq Chi Df Pr(>Chisq)
f0128bi   7 615.73 655.48 -300.87
f0128biu  8 687.01 732.43 -335.51     0      1          1

So update() does not help. I guess I will need to fool with starting values.

Murray
*
On 17 May 2012 14:34, David Duffy <David.Duffy@qimr.edu.au> wrote:

> On Thu, 17 May 2012, Murray Jorgensen wrote:
>
>  *I tried using 'update' which should (I believe) initialise at the simpler
>>
>> model:
>>
>>  f0128biu = update(f0128bi, . ~ . + (1|tree))
>>>
>> Error in model.frame.default(formula = fincr ~ icfac + (1 + gap) + (1 +  :
>>  invalid type (closure) for variable 'box'
>>
>
> It's finding the box() function.
>
> --
> | David Duffy (MBBS PhD)                                         ,-_|\
> | email: davidD@qimr.edu.au  ph: INT+61+7+3362-0217 fax: -0101  /     *
> | Epidemiology Unit, Queensland Institute of Medical Research   \_,-._/
> | 300 Herston Rd, Brisbane, Queensland 4029, Australia  GPG 4D0B994A v
>



-- 
Dr Murray Jorgensen      http://www.stats.waikato.ac.nz/Staff/maj.html
Department of Statistics, University of Waikato, Hamilton, New Zealand
Email: maj@waikato.ac.nz    majorgensen@ihug.co.nz      Fax 7 838 4155
Phone  +64 7 838 4773 wk    Home +64 7 825 0441   Mobile 021 0200 8350

	[[alternative HTML version deleted]]

