Re: difficulty with FOCE
From: Justin Wilkins justin.wilkins@farmbio.uu.se
Subject: Re: [NMusers] difficulty with FOCE
Date: Thu, 13 Oct 2005 17:48:30 +0200
Hi Jing,
Without seeing your model it's hard to narrow down the specific cause of your problem, but
a typical rule of thumb to follow when this sort of thing happens is to make sure that your
model is not overparameterized. FOCE is much more sensitive to unneeded model parameters than FO.
I would test each of your parameters - starting with ETAs, and the ones with the largest RSE -
to make sure it belongs in the model. If the 95% CI of a parameter (calculated from the standard error)
overlaps zero, it's probably not needed.
Check your FOCE output file for zero gradients (or very low gradient values) in any of the parameters
as another clue as to which to start investigating first. If any of the FOCE parameter estimates
look odd or very low, remove them and try again. If you're using ADVAN6, you might try using ADVAN8
instead, especially if your system is stiff (if there are large differences in the quantities used
in your differential equations).
It would be very helpful if you could post your model code so we can get a better idea of
where the fault lies!
Use of FO is associated with inflated rates of type I error (see publications on this topic by Whlby
et al). Despite FOCE's relatively poor stability and greatly-increased run times, if you can get you
model to converge using it - preferably with the INTER option - I would be much more confident in your
parameter estimates. FO is useful for 'quick and dirty' model-based analysis but it is very definitely
not as good as FOCE. I would feel satisfied with FO only where the time needed for FOCE was excessive.
Justin
_______________________________________________________