justpairs           package:hyperdirichlet           R Documentation

_H_y_p_e_r_d_i_r_i_c_h_l_e_t _d_i_s_t_r_i_b_u_t_i_o_n _w_i_t_h _n_o_n_z_e_r_o _c_o_e_f_f_i_c_i_e_n_t_s _f_o_r _j_u_s_t
_t_h_e _p_a_i_r_s

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

     Hyperdirichlet distribution with nonzero coefficients for terms
     consisting of a pair of ps.

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

     justpairs(x)

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

       x: A square matrix.  Upper triangular elements correspond to
          pairs of included ps, and lower triangular elements
          correspond to pairs of excluded ps

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

     The PDF is


          prod_{i<j}(p_i+p_j)^{M[i,j]} (1-p_i-p_j)^{M[j,i]}

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

     Returns a hyperdirichlet object (without normalizing factor) of
     dimension 'nrow(x)'.

_N_o_t_e:

     The case where 'nrow(x)==4' is treated specially.  In this case,
     the lower triangular elements are discarded, with a warning.

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

     Robin K. S. Hankin

_S_e_e _A_l_s_o:

     'paircomp'

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

     justpairs(matrix(1:25,5,5))

