Dayhoff (Original) Calculation

Calling Sequences:
CreateOrigDayMatrix(mutations, counts, PAM)
CreateOrigDayMatrix(mutations, counts, 1..UpperPam)
Parameters:

mutations	:	array(real, real)
counts	:	array(real)
PAM, UpperPam	:	real

Returns: {DayMatrix, list(DayMatrix)}

Synposis: This function computes a Dayhoff matrix (structured type DayMatrix) computed by the method first given by Dayhoff et. all [] given an observed mutation matrix mutations, a frequency vector counts and a PAM distance PAM (or range of PAM distances beginning at 1).

Examples:

> OrigTot := [87, 41, 40, 47, 33, 38, 50, 89, 34, 37, 
>             85, 81, 15, 40, 51, 70, 58, 10, 30, 65];
> OrigFreq := OrigTot/sum(OrigTot);
> OrigDM := CreateOrigDayMatrix(Mutations1978, OrigFreq, 250);
> print(OrigDM);

See also: CreateDayMatrix