[R] Help: Maximum likelihood estimation

dave fournier davef at otter-rsch.com
Fri Oct 22 20:56:17 CEST 2010


Actually it is not that difficult to parameterize the covariance matrix 
so that the
optimization is unconstrained. first parameterize the correlation matrix 
and the
standard deviations separately. the std devs can be parameterized as

       sigma_i=exp(x_i)     1<=i<=n

For the correlation matrix parameterize it in terms of its choleski 
decomposition.

this is a lower triangular matrix

    1      0       0
    a_1  a_2    0
    b _1 b_2 b_3
      ....


such that  the norm  of each row is 1

to ensure this  form start with

        1
       y_1  1
       y_2  y_3  1
         ...

  and normalize each row to have norm 1. There are n*(n-1)/2  of these y's.
  together with the n x_i's you have n*(n+1)/2 parameters as you should.



More information about the R-help mailing list