Re: [Rd] Numerical optimisation and "non-feasible" regions

From: Mathieu Ribatet <>
Date: Thu, 07 Aug 2008 11:52:28 +0200

Thanks Ben for your tips.
I'm not sure it'll be so easy to do (as the non-feasible regions depend on the model parameters), but I'm sure it's worth giving a try. Thanks !!!


Ben Bolker a écrit :
> Mathieu Ribatet <mathieu.ribatet <at>> writes:
>> Dear list,
>> I'm currently writing a C code to compute the (composite) likelihood -
>> well this is done but not really robust. The C code is wrapped in an R
>> one which call the optimizer routine - optim or nlm. However, the
>> fitting procedure is far from being robust as the parameter space
>> depends on the parameter - I have a covariance matrix that should be a
>> valid one for example.
> One reasonably straightforward hack to deal with this is
> to add a penalty that is (e.g.) a quadratic function of the
> distance from the feasible region, if that is reasonably
> straightforward to compute -- that way your function will
> get gently pushed back toward the feasible region.
> Ben Bolker
> ______________________________________________
> mailing list

Institute of Mathematics
Ecole Polytechnique Fédérale de Lausanne
CH-1015 Lausanne   Switzerland
Tel: + 41 (0)21 693 7907

______________________________________________ mailing list
Received on Thu 07 Aug 2008 - 10:03:14 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 07 Aug 2008 - 12:36:10 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive