From mathieu.ribatet at epfl.ch  Thu Aug  7 11:52:28 2008
From: mathieu.ribatet at epfl.ch (Mathieu Ribatet)
Date: Thu, 07 Aug 2008 11:52:28 +0200
Subject: [Rd] Numerical optimisation and "non-feasible" regions
In-Reply-To: <loom.20080807T004607-309@post.gmane.org>
References: <4899D089.6020208@epfl.ch>
	<loom.20080807T004607-309@post.gmane.org>
Message-ID: <489AC5DC.10508@epfl.ch>

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 !!!
Best,

Mathieu

Ben Bolker a ?crit :
> 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
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>   

-- 
Institute of Mathematics
Ecole Polytechnique F?d?rale de Lausanne
STAT-IMA-FSB-EPFL, Station 8
CH-1015 Lausanne   Switzerland
http://stat.epfl.ch/
Tel: + 41 (0)21 693 7907



[Rd] Numerical optimisation and "non-feasible" regions
Mathieu Ribatet