On 20/10/11 02:31, Serguei Kaniovski wrote:

Dear R-Users

I would like to estimate a constrained bivariate normal density, the
constraint being that the means are of equal magnitude but of opposite
signs. So I need to estimate four parameters:

mu    (meanvector (mu,-mu))
sigma_1 and sigma_2 (two sd deviations)
rho (correlation coefficient)

I have looked at several packages, including Gaussian mixture models in
Mclust, but I am not sure what is the best way, or the best package to use
for this task.

Greatly appreciate any suggestions!

I very much doubt that there is anything built-in that you can use.
However it shouldn't be *too* hard to get maximum likelihood estimates
using optim() to maximise the (log) likelihood.

For starting values I would try just using the ordinary covariance matrix estimate to get your sigma_1, sigma_2, and rho estimates,
and for mu use (x1.bar - x2.bar)/2 (in what I hope is an obvious
notation) for a starting value.

Good luck!

        cheers,

                Rolf Turner

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to