Rmatsolve               package:adlift               R Documentation

_R_m_a_t_s_o_l_v_e

_D_e_s_c_r_i_p_t_i_o_n:

     This function calculates matrix inverses for symmetric matrices.

_U_s_a_g_e:

     Rmatsolve(m)

_A_r_g_u_m_e_n_t_s:

       m: a (symmetric) matrix.

_D_e_t_a_i_l_s:

     This function uses the eigenvalue decomposition of a matrix 'm' to
     work out its inverse.  The function is used here since standard
     matrix inverse algorithms do not cope well with matrices which are
     near singular (this often happens in the regression stages of the
     forward transforms.

_V_a_l_u_e:

     inv: the matrix inverse of 'm'.

_A_u_t_h_o_r(_s):

     Matt Nunes (matt.nunes@bristol.ac.uk), Marina Popa
     (Marina.Popa@bristol.ac.uk)

_E_x_a_m_p_l_e_s:

     #
     #create a 4x4 matrix
     m<-matrix(runif(16),4,4)

     temp<-crossprod(m)

     #i.e. temp is t(m)

     inv<-Rmatsolve(temp)

