
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:05 02/24/2011

Input data from file "/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_17858.phylip" (19 taxa, 589 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


(Output sorted by AICc. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AICc		wAICc	cwAICc
HKYG	40	-1390.8315	2867.6484	0.4447	0.4447
HKYI	40	-1391.0162	2868.0179	0.3697	0.8143
HKYIG	41	-1390.7090	2869.7141	0.1583	0.9727
GTRG	44	-1389.7507	2874.7808	0.0126	0.9852
GTRI	44	-1389.9289	2875.1373	0.0105	0.9957
GTRIG	45	-1389.6590	2876.9424	0.0043	1.0000
HKY	39	-1405.0190	2893.7211	0.0000	1.0000
F81G	39	-1408.0556	2899.7942	0.0000	1.0000
F81I	39	-1408.2503	2900.1837	0.0000	1.0000
K2PG	37	-1410.7245	2900.5525	0.0000	1.0000
GTR	43	-1403.8928	2900.7287	0.0000	1.0000
K2PI	37	-1410.8971	2900.8977	0.0000	1.0000
F81IG	40	-1407.9373	2901.8601	0.0000	1.0000
K2PIG	38	-1410.6417	2902.6725	0.0000	1.0000
SYMG	41	-1409.3081	2906.9124	0.0000	1.0000
SYMI	41	-1409.4561	2907.2083	0.0000	1.0000
SYMIG	42	-1409.2560	2909.1273	0.0000	1.0000
F81	38	-1421.6416	2924.6724	0.0000	1.0000
K2P	36	-1424.0823	2924.9907	0.0000	1.0000
JC69G	36	-1426.0942	2929.0145	0.0000	1.0000
JC69I	36	-1426.2725	2929.3710	0.0000	1.0000
SYM	40	-1422.3334	2930.6521	0.0000	1.0000
JC69IG	37	-1426.0072	2931.1178	0.0000	1.0000
JC69	35	-1439.2546	2953.0662	0.0000	1.0000



(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1390.8315	2861.6630	0.4223	0.4223
HKYI	40	-1391.0162	2862.0325	0.3511	0.7733
HKYIG	41	-1390.7090	2863.4179	0.1756	0.9489
GTRG	44	-1389.7507	2867.5014	0.0228	0.9717
GTRI	44	-1389.9289	2867.8579	0.0191	0.9908
GTRIG	45	-1389.6590	2869.3181	0.0092	1.0000
HKY	39	-1405.0190	2888.0381	0.0000	1.0000
GTR	43	-1403.8928	2893.7856	0.0000	1.0000
F81G	39	-1408.0556	2894.1112	0.0000	1.0000
F81I	39	-1408.2503	2894.5007	0.0000	1.0000
K2PG	37	-1410.7245	2895.4491	0.0000	1.0000
K2PI	37	-1410.8971	2895.7942	0.0000	1.0000
F81IG	40	-1407.9373	2895.8747	0.0000	1.0000
K2PIG	38	-1410.6417	2897.2834	0.0000	1.0000
SYMG	41	-1409.3081	2900.6162	0.0000	1.0000
SYMI	41	-1409.4561	2900.9121	0.0000	1.0000
SYMIG	42	-1409.2560	2902.5120	0.0000	1.0000
F81	38	-1421.6416	2919.2833	0.0000	1.0000
K2P	36	-1424.0823	2920.1646	0.0000	1.0000
JC69G	36	-1426.0942	2924.1884	0.0000	1.0000
JC69I	36	-1426.2725	2924.5449	0.0000	1.0000
SYM	40	-1422.3334	2924.6667	0.0000	1.0000
JC69IG	37	-1426.0072	2926.0143	0.0000	1.0000
JC69	35	-1439.2546	2948.5092	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1390.8315	3036.8000	0.5325	0.5325
HKYI	40	-1391.0162	3037.1695	0.4427	0.9751
HKYIG	41	-1390.7090	3042.9334	0.0248	0.9999
K2PG	37	-1410.7245	3057.4508	0.0000	1.0000
K2PI	37	-1410.8971	3057.7960	0.0000	1.0000
HKY	39	-1405.0190	3058.7967	0.0000	1.0000
GTRG	44	-1389.7507	3060.1521	0.0000	1.0000
GTRI	44	-1389.9289	3060.5086	0.0000	1.0000
K2PIG	38	-1410.6417	3063.6636	0.0000	1.0000
F81G	39	-1408.0556	3064.8698	0.0000	1.0000
F81I	39	-1408.2503	3065.2593	0.0000	1.0000
GTRIG	45	-1389.6590	3066.3472	0.0000	1.0000
F81IG	40	-1407.9373	3071.0117	0.0000	1.0000
K2P	36	-1424.0823	3077.7880	0.0000	1.0000
SYMG	41	-1409.3081	3080.1317	0.0000	1.0000
SYMI	41	-1409.4561	3080.4276	0.0000	1.0000
JC69G	36	-1426.0942	3081.8117	0.0000	1.0000
GTR	43	-1403.8928	3082.0579	0.0000	1.0000
JC69I	36	-1426.2725	3082.1683	0.0000	1.0000
F81	38	-1421.6416	3085.6635	0.0000	1.0000
SYMIG	42	-1409.2560	3086.4059	0.0000	1.0000
JC69IG	37	-1426.0072	3088.0161	0.0000	1.0000
SYM	40	-1422.3334	3099.8037	0.0000	1.0000
JC69	35	-1439.2546	3101.7541	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_17858.phylip = 1 - 589;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_17858.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_17858.phylip = 1 - 589;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_17858.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_17858.phylip = 1 - 589;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_17858.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

