[Master Index] [Index for export/lyngby]

lyngby_pde_binmix

(export/lyngby/lyngby_pde_binmix.m)


Function Synopsis

[f, df] = lyngby_pde_binmix_cd33(lambda, P1, P2, N)

Help text

 lyngby_pde_binmix     - Estimation of binomial mixture model

       O = lyngby_pde_binmix(N)

       Input:  N  Observed frequencies

       Output: O  Parameters [ P1 P2 lambda ]

       Estimation of mixture of two binomials (binomial mixture)

       The first element in N should contain the observed frequency
       zero positive outcomes, the second the observed frequency of
       one positive outcomes, the third the observed frequency of two
       positive outcomes, ...
 
       The input argument 'N' can also be a matrix. The rows are then
       modeled with a common lambda.

       The parameters in the output argument is ordered as: P1, P2,
       lambda, where lambda is the mixing coefficient. The return
       argument is sorted so that P1<=P2.


       Ref: John Uebersax, MIXBIN.
            Gelfand AE, Solomon A, JASA 69:32

       Example:
         N = [ 5 0 2 2 3 ];
         lyngby_pde_binmix(N)

         % Yerushalmy radiograph data from MIXBIN
         %   P1=0.0100   P2=0.6446, lambda= 0.9808
         N = [ 13560 877 168 66 42 28 23 39 64];
         lyngby_pde_binmix(N, 'info', 5)
 
       See also LYNGBY, LYNGBY_PDE_BINMIX_LL, LYNGBY_PDE_BINMIX_DER,
                LYNGBY_OPT_CG, LYNGBY_OPT_GRADIENT.  

 $Id: lyngby_pde_binmix.m,v 1.2 2007/09/21 16:15:34 fn Exp $

Cross-Reference Information

This function calls
This function is called by

Produced by mat2html on Wed Jul 29 15:43:40 2009
Cross-Directory links are: OFF