The HAPLOTYPE Procedure |
Continuing the example from the section Getting Started: HAPLOTYPE Procedure, suppose you are concerned that the likelihood reached a local and not a global maximum. You can request that PROC HAPLOTYPE use several different sets of initial haplotype frequencies to ensure that you find a global maximum of the likelihood. The following code invokes the EM algorithm with five different sets of initial values, including the set used in the "Getting Started" example:
proc haplotype data=markers prefix=SNP init=random seed=51220 nstart=5; var m1-m8; run;
The NSTART=5 option requests that the EM algorithm be run three times with randomly generated initial frequencies, including once using the seed 51220 that was previously used, once using uniform initial frequencies, and once using haplotype frequencies given by the product of the allele frequencies. The two tables in Output 8.2.1 are from the run that produced the best log likelihood:
Analysis Information | |
---|---|
Loci Used | SNP1 SNP2 SNP3 SNP4 |
Number of Individuals | 25 |
Number of Starts | 5 |
Convergence Criterion | 0.00001 |
Iterations Checked for Conv. | 1 |
Maximum Number of Iterations | 100 |
Number of Iterations Used | 19 |
Log Likelihood | -95.94742 |
Initialization Method | Random |
Random Number Seed | 499887544 |
Standard Error Method | Binomial |
Haplotype Frequency Cutoff | 0 |
Haplotype Frequencies | |||||
---|---|---|---|---|---|
Number | Haplotype | Freq | Standard Error |
95% Confidence Limits | |
1 | A-A-A-A | 0.14324 | 0.05005 | 0.04515 | 0.24133 |
2 | A-A-A-B | 0.07507 | 0.03764 | 0.00129 | 0.14885 |
3 | A-A-B-A | 0.00000 | 0.00001 | 0.00000 | 0.00001 |
4 | A-A-B-B | 0.00000 | 0.00010 | 0.00000 | 0.00019 |
5 | A-B-A-A | 0.09295 | 0.04148 | 0.01165 | 0.17425 |
6 | A-B-A-B | 0.05349 | 0.03214 | 0.00000 | 0.11649 |
7 | A-B-B-A | 0.00001 | 0.00052 | 0.00000 | 0.00103 |
8 | A-B-B-B | 0.07523 | 0.03768 | 0.00138 | 0.14909 |
9 | B-A-A-A | 0.08644 | 0.04014 | 0.00776 | 0.16512 |
10 | B-A-A-B | 0.08784 | 0.04044 | 0.00859 | 0.16710 |
11 | B-A-B-A | 0.07904 | 0.03854 | 0.00350 | 0.15459 |
12 | B-A-B-B | 0.10836 | 0.04441 | 0.02133 | 0.19540 |
13 | B-B-A-A | 0.10097 | 0.04304 | 0.01661 | 0.18533 |
14 | B-B-A-B | 0.00000 | 0.00000 | 0.00000 | 0.00000 |
15 | B-B-B-A | 0.09735 | 0.04235 | 0.01435 | 0.18035 |
16 | B-B-B-B | 0.00000 | 0.00000 | 0.00000 | 0.00000 |
Copyright © 2008 by SAS Institute Inc., Cary, NC, USA. All rights reserved.