GETWGT - Dirichlet Mixture Estimation

GETWGT is a set of routines for Dirichlet Mixture estimation. The main requirement in the design of this library was that there be a single routine with a simple interface through which the user interacts. This routine, called getwgt, accepts a set of nucleic acid counts, updates an internal Dirichlet mixture model, and returns its current estimate for the pseudocounts. All other transactions and information are hidden from the user.

The program requires a data file, containing the parameters of the Dirichlet density distributions that make up the mixture. This file needs to be named weights.txt. Two files are provided here as possible sources of this data. These files were obtained from the UCSC computational biology page, where more information, papers, and data is available.

Files you may copy include:

• GETWGT.F90, the source code;
• GETWGT_PRB.F90, a sample calling program;
• GETWGT_PRB.OUT, the output from the sample program;
• BYST_WEIGHTS.TXT, the "current best" data defining a Dirichlet mixture.
• BLOCKS_WEIGHTS.TXT, the BLOCKS data defining a Dirichlet mixture.

The list of routines includes:

• GETWGT updates the Dirichlet mixture weights based on a set of counts.
• C_CAP capitalizes a single character.
• C_EQI is a case insensitive comparison of two characters for equality.
• C_NEXT "reads" space-separated characters from a string, one at a time.
• C_TO_DIGIT returns the integer value of a base 10 digit.
• COMP_PARAM_PRINT prints the parameters for the mixture components.
• DIRICHLET_MEAN returns the means of the Dirichlet PDF.
• DIRICHLET_MULTINOMIAL_PDF evaluates a Dirichlet Multinomial PDF.
• EVENT_PROCESS updates the mixture weight distribution parameters.
• GAMMA_LOG calculates the natural logarithm of GAMMA ( X ) for positive X.
• I_NEXT "reads" integers from a string, one at a time.
• IVEC_SUM sums the entries of an integer vector.
• MIXTURE_READ reads the Dirichlet mixture parameters from a file.
• R_NEXT "reads" real numbers from a string, one at a time.
• RVEC_COPY copies a real vector.
• RVEC_PRINT prints a real vector, with an optional title.
• RVEC_SUM sums the entries of a real vector.
• RVEC_UNIT_SUM normalizes a real vector to have unit sum.
• S_BEGIN is TRUE if one string matches the beginning of the other.
• S_TO_I reads an integer value from a string.
• S_TO_R reads a real number from a string.

Return to the biomedical software page.