DNAml_rates_1_0
Gary J. Olsen
August 14, 1992
The DNAml_rates program takes a set of sequences and corresponding
phylogenetic tree and produces and maximum likelihood estimate of the
rate of nucleotides substitution at each sequence position.
Input is read from standard input. The format is very much like that
of the fastDNAml program. The first line of the input file gives the
number of sequences and the number of bases per sequence. Also on
this line are the requested program option letters. Any auxiliary
data required by the options follow on subsequent lines. Either the
user must specify the empirical base frequencies (F) option, or
immediately preceding the data matrix there must be a line of data
with the frequencies of A, C, G and T. Next, the program expects a
data matrix. The first 10 characters of the first line of data for a
given sequence in interpreted as the name (blanks are counted).
Elsewhere in the data matrix, blanks and numbers are ignored. The
default data matrix format is interleaved. If all the data for a
sequence are on one input line, then interleaved and noninterleaved
are equivalent. Following the data matrix there must be a line with
the number of user-specified trees for which rates are to be estimated
(as with the U option is fastDNAml). The rest of the input file is
one or more user-specified trees with branch lengths (as with the U
and L options in fastDNAml).
The program writes to standard output. The output lists the estimated
rate of change at every site in the sequence, or "Undefined" if there
are not sufficient unambiguous data at the site.
If the C option is specified, the program also categorizes the rates
into the requested number of categories. The current categorization
algorithm is rather crude, but is probably adequate if the number of
categories is large enough. A weighting mask is also created in which
sites with Undefined rates are assigned a weight of zero.
If the Y option is specified, the program writes the weights and
categories data to a file in a format appropriate for use by
fastDNAml.
Options summary:
1 - print data. Toggles print data option (default = noprint).
C - write categories. Requires auxiliary line with a C and the desired
number of categories.
F - empirical base frequencies. Calculates base frequencies from data matrix,
rather than expecting a base frequency input line.
I - interleave. Toggles the data interleave option (default = interleave).
L - userlengths. This is implicit in the program, so the option is ignored.
M - minimum informative sequences. Requires an auxiliary data line with an
M and the minimum number of sequences in which a sequence position
(alignment column) must have unambiguous information in order for the rate
at the site to be defined (default = 4).
T - transitions/transversion ratio. Requires auxiliary line with a T and
the ration of observed transitions to transversions (default = 2.0).
U - user trees. This is implicit in the program, so the option is ignored.
W - user weights. Requires weights auxiliary data.
Y - categories file. Writes the weights and categories to a file.
The option scripts usertree, weights, n_categories and categories_file are
useful for adding the appropriate options to the input data matrix.
The option script weights_categories is useful for adding the resulting
outfile to a fastDNAml input file.
|