From bolker at ufl.edu  Thu Aug  7 02:47:11 2008
From: bolker at ufl.edu (Ben Bolker)
Date: Thu, 7 Aug 2008 00:47:11 +0000 (UTC)
Subject: [Rd] Numerical optimisation and "non-feasible" regions
References: <4899D089.6020208@epfl.ch>
Message-ID: <loom.20080807T004607-309@post.gmane.org>

Mathieu Ribatet <mathieu.ribatet <at> epfl.ch> 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



[Rd] Numerical optimisation and "non-feasible" regions
Ben Bolker