MIXTURE - Dirichlet Mixture Simulation

MIXTURE simulates a process governed by a Dirichlet Mixture probability density function. Each realization of the process involves choosing one of a specified set of Dirichlet PDF's, with given weight, and then sampling that PDF.

The current code is pretty rough, and simply does a simulation. Future work will involve using the simulation to test procedures for estimating the number of Dirichlet densities in the mixture, and their parameters, using a Bayesian approach.

Files you may copy include:

• mixture.f90, the source code;
• mixture.dat, a mixture definition file;
• mixture.out, the sample output.

The list of routines includes:

• MAIN is the main program for the mixture simulation.
• WEIGHT_INIT initializes the estimated weights.
• WEIGHT_UPDATE updates the estimated weights.
• AMINO_PRINT prints the amino acid parameters.
• 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_AMINO_NAME converts a character to an amino acid name.
• C_TO_DIGIT returns the integer value of a base 10 digit.
• COMP_PARAM_PRINT prints the parameters for the mixture components.
• COMP_STATS_PRINT prints the mean and variance for the mixture components.
• FAVOR_RATIO_COMPUTE computes the ratio by which a component density favors an amino acid.
• FAVOR_RATIO_PRINT prints the favor ratios.
• I_NEXT "reads" integers from a string, one at a time.
• MIXTURE_PRINT prints the Dirichlet mixture parameters.
• MIXTURE_READ reads the Dirichlet mixture parameters from a file.
• R_NEXT "reads" real numbers from a string, one at a time.
• R_RANDOM returns a random real in a given range.
• RROW_MEAN returns the means of rows of a real array.
• RROW_VARIANCE returns the variances of the rows of a real array.
• RVEC_SUM sums the entries of a real vector.
• S_BEGIN is TRUE if one string matches the beginning of the other.
• S_EQI is a case insensitive comparison of two strings for equality.
• S_TO_I reads an integer value from a string.
• S_TO_R reads a real number from a string.
• SAMPLE_ANALYZE analyzes the samples from the Dirichlet mixture PDF.
• SAMPLE_PROJECT computes the projection of a sample onto the mixture components.
• SAMPLE_PROJECT computes the projection of a sample onto the mixture components.

Return to the biomedical software page.