
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:04 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/chr8_3308.phylip" (19 taxa, 572 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	-2507.0757	5100.3285	0.6881	0.6881
HKYIG	41	-2507.0757	5102.6495	0.2156	0.9037
GTRG	44	-2504.8829	5105.2800	0.0579	0.9616
HKYI	40	-2510.6578	5107.4927	0.0191	0.9807
GTRIG	45	-2504.8829	5107.6365	0.0178	0.9986
GTRI	44	-2508.5705	5112.6552	0.0014	1.0000
K2PG	37	-2530.3509	5139.9678	0.0000	1.0000
K2PIG	38	-2530.3509	5142.2629	0.0000	1.0000
SYMG	41	-2529.2881	5147.0744	0.0000	1.0000
K2PI	37	-2534.1170	5147.4999	0.0000	1.0000
SYMIG	42	-2529.2881	5149.4042	0.0000	1.0000
HKY	39	-2535.1829	5154.2304	0.0000	1.0000
SYMI	41	-2533.1383	5154.7747	0.0000	1.0000
GTR	43	-2532.6981	5158.5628	0.0000	1.0000
K2P	36	-2559.8293	5196.6380	0.0000	1.0000
SYM	40	-2558.4785	5203.1340	0.0000	1.0000
F81G	39	-2586.5287	5256.9220	0.0000	1.0000
F81IG	40	-2586.5286	5259.2343	0.0000	1.0000
F81I	39	-2590.0735	5264.0117	0.0000	1.0000
JC69G	36	-2604.0044	5284.9882	0.0000	1.0000
JC69IG	37	-2604.0044	5287.2747	0.0000	1.0000
JC69I	36	-2607.6531	5292.2856	0.0000	1.0000
F81	38	-2613.5469	5308.6547	0.0000	1.0000
JC69	35	-2631.8459	5338.3932	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2507.0757	5094.1514	0.6156	0.6156
HKYIG	41	-2507.0757	5096.1513	0.2265	0.8421
GTRG	44	-2504.8829	5097.7658	0.1010	0.9432
GTRIG	45	-2504.8829	5099.7658	0.0372	0.9803
HKYI	40	-2510.6578	5101.3157	0.0171	0.9975
GTRI	44	-2508.5705	5105.1409	0.0025	1.0000
K2PG	37	-2530.3509	5134.7019	0.0000	1.0000
K2PIG	38	-2530.3509	5136.7019	0.0000	1.0000
SYMG	41	-2529.2881	5140.5763	0.0000	1.0000
K2PI	37	-2534.1170	5142.2340	0.0000	1.0000
SYMIG	42	-2529.2881	5142.5762	0.0000	1.0000
SYMI	41	-2533.1383	5148.2766	0.0000	1.0000
HKY	39	-2535.1829	5148.3658	0.0000	1.0000
GTR	43	-2532.6981	5151.3961	0.0000	1.0000
K2P	36	-2559.8293	5191.6586	0.0000	1.0000
SYM	40	-2558.4785	5196.9570	0.0000	1.0000
F81G	39	-2586.5287	5251.0573	0.0000	1.0000
F81IG	40	-2586.5286	5253.0573	0.0000	1.0000
F81I	39	-2590.0735	5258.1470	0.0000	1.0000
JC69G	36	-2604.0044	5280.0088	0.0000	1.0000
JC69IG	37	-2604.0044	5282.0088	0.0000	1.0000
JC69I	36	-2607.6531	5287.3062	0.0000	1.0000
F81	38	-2613.5469	5303.0938	0.0000	1.0000
JC69	35	-2631.8459	5333.6917	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2507.0757	5268.1170	0.9349	0.9349
HKYIG	41	-2507.0757	5274.4660	0.0391	0.9740
HKYI	40	-2510.6578	5275.2812	0.0260	1.0000
GTRG	44	-2504.8829	5289.1279	0.0000	1.0000
GTRIG	45	-2504.8829	5295.4770	0.0000	1.0000
K2PG	37	-2530.3509	5295.6200	0.0000	1.0000
GTRI	44	-2508.5705	5296.5030	0.0000	1.0000
K2PIG	38	-2530.3509	5301.9692	0.0000	1.0000
K2PI	37	-2534.1170	5303.1522	0.0000	1.0000
HKY	39	-2535.1829	5317.9822	0.0000	1.0000
SYMG	41	-2529.2881	5318.8910	0.0000	1.0000
SYMIG	42	-2529.2881	5325.2401	0.0000	1.0000
SYMI	41	-2533.1383	5326.5913	0.0000	1.0000
GTR	43	-2532.6981	5338.4091	0.0000	1.0000
K2P	36	-2559.8293	5348.2276	0.0000	1.0000
SYM	40	-2558.4785	5370.9226	0.0000	1.0000
F81G	39	-2586.5287	5420.6738	0.0000	1.0000
F81IG	40	-2586.5286	5427.0228	0.0000	1.0000
F81I	39	-2590.0735	5427.7634	0.0000	1.0000
JC69G	36	-2604.0044	5436.5778	0.0000	1.0000
JC69IG	37	-2604.0044	5442.9269	0.0000	1.0000
JC69I	36	-2607.6531	5443.8752	0.0000	1.0000
F81	38	-2613.5469	5468.3610	0.0000	1.0000
JC69	35	-2631.8459	5485.9116	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/chr8_3308.phylip = 1 - 572;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_3308.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/chr8_3308.phylip = 1 - 572;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_3308.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/chr8_3308.phylip = 1 - 572;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_3308.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

