estErrProbMethodOfMoments      package:ppiStats      R Documentation

_E_s_t_i_m_a_t_e _f_a_l_s_e _p_o_s_i_t_i_v_e _a_n_d _f_a_l_s_e _n_e_g_a_t_i_v_e _e_r_r_o_r _p_r_o_b_a_b_i_l_i_t_i_e_s
_b_y _m_e_t_h_o_d _m_o_m_e_n_t_s.

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

     Estimate false positive and false negative error probabilities by
     method moments.

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

     estErrProbMethodOfMoments(nint, nrec, nunr, ntot)

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

    nint: Integer vector. True number of interactions. Typically, the
          function is called for a range of these, returning all
          possible solutions for that range.

    nrec: Integer scalar. Observed number of reciprocated edges.

    nunr: Integer scalar. Observed number of unreciprocated edges.

    ntot: Integer scalar. Number of proteins which were tested twice
          (e.g. both as viable bait and as viable prey).

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

     The model is described in the vignette _Stochastic and systematic
     errors in PPI data, by looking at unreciprocated in- or out-edges_
     by W. Huber, T. Chiang and R. Gentleman.

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

     Matrix with 5 columns 'nint' (a copy of the input argument),
     'pfp1', 'pfn1', 'pfp2' and 'pfn2', and as many rows as the length
     of 'nint'.

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

     Wolfgang Huber <URL: http://www.ebi.ac.uk/huber>

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

     est = estErrProbMethodOfMoments(nint=seq(8000, 40000, by=100), nrec=9722, nunr=15856, ntot=2000)
     if(interactive()) {
       plot(est[, c("pfp2", "pfn2")], type="l", col="blue", lwd=2,
            xlab=expression(p[FP]), ylab=expression(p[FN]))
       abline(h=0, v=0, lty=2)
     }

