
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:54 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_17721.phylip" (19 taxa, 426 characters)

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	37	-1081.1715	2243.5904	0.3157	0.3157
K2PI	37	-1081.2704	2243.7882	0.2860	0.6016
HKYI	40	-1078.5843	2245.6880	0.1106	0.7123
HKYG	40	-1078.5997	2245.7188	0.1089	0.8212
K2PIG	38	-1081.0792	2245.8173	0.1037	0.9248
HKYIG	41	-1078.5154	2247.9995	0.0348	0.9597
SYMI	41	-1079.3941	2249.7569	0.0145	0.9741
SYMG	41	-1079.5463	2250.0614	0.0124	0.9865
GTRI	44	-1076.8186	2252.0309	0.0046	0.9912
GTRG	44	-1076.9806	2252.3548	0.0039	0.9951
SYMIG	42	-1079.6082	2252.6472	0.0034	0.9985
GTRIG	45	-1076.9922	2254.8790	0.0011	0.9997
K2P	36	-1089.6009	2258.0502	0.0002	0.9999
HKY	39	-1086.8874	2259.8577	0.0001	1.0000
SYM	40	-1087.8198	2264.1590	0.0000	1.0000
GTR	43	-1085.2808	2266.4673	0.0000	1.0000
JC69G	36	-1093.8138	2266.4759	0.0000	1.0000
JC69I	36	-1093.9109	2266.6701	0.0000	1.0000
F81I	39	-1091.0848	2268.2525	0.0000	1.0000
F81G	39	-1091.1121	2268.3071	0.0000	1.0000
JC69IG	37	-1093.7394	2268.7263	0.0000	1.0000
F81IG	40	-1091.0437	2270.6070	0.0000	1.0000
JC69	35	-1101.5526	2279.5668	0.0000	1.0000
F81	38	-1098.7441	2281.1470	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-1081.1715	2236.3430	0.2341	0.2341
K2PI	37	-1081.2704	2236.5408	0.2121	0.4462
HKYI	40	-1078.5843	2237.1685	0.1549	0.6011
HKYG	40	-1078.5997	2237.1993	0.1526	0.7537
K2PIG	38	-1081.0792	2238.1584	0.0945	0.8482
HKYIG	41	-1078.5154	2239.0307	0.0611	0.9092
SYMI	41	-1079.3941	2240.7881	0.0254	0.9346
SYMG	41	-1079.5463	2241.0927	0.0218	0.9564
GTRI	44	-1076.8186	2241.6372	0.0166	0.9730
GTRG	44	-1076.9806	2241.9611	0.0141	0.9871
SYMIG	42	-1079.6082	2243.2164	0.0075	0.9946
GTRIG	45	-1076.9922	2243.9843	0.0051	0.9997
K2P	36	-1089.6009	2251.2018	0.0001	0.9999
HKY	39	-1086.8874	2251.7748	0.0001	1.0000
SYM	40	-1087.8198	2255.6395	0.0000	1.0000
GTR	43	-1085.2808	2256.5616	0.0000	1.0000
JC69G	36	-1093.8138	2259.6276	0.0000	1.0000
JC69I	36	-1093.9109	2259.8218	0.0000	1.0000
F81I	39	-1091.0848	2260.1696	0.0000	1.0000
F81G	39	-1091.1121	2260.2242	0.0000	1.0000
JC69IG	37	-1093.7394	2261.4789	0.0000	1.0000
F81IG	40	-1091.0437	2262.0875	0.0000	1.0000
JC69	35	-1101.5526	2273.1052	0.0000	1.0000
F81	38	-1098.7441	2273.4881	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1081.1715	2386.3572	0.5085	0.5085
K2PI	37	-1081.2704	2386.5551	0.4606	0.9690
K2PIG	38	-1081.0792	2392.2271	0.0270	0.9960
K2P	36	-1089.6009	2397.1617	0.0023	0.9983
HKYI	40	-1078.5843	2399.3461	0.0008	0.9991
HKYG	40	-1078.5997	2399.3769	0.0008	0.9999
HKYIG	41	-1078.5154	2405.2628	0.0000	0.9999
JC69G	36	-1093.8138	2405.5874	0.0000	0.9999
JC69I	36	-1093.9109	2405.7816	0.0000	1.0000
SYMI	41	-1079.3941	2407.0201	0.0000	1.0000
SYMG	41	-1079.5463	2407.3247	0.0000	1.0000
HKY	39	-1086.8874	2409.8979	0.0000	1.0000
JC69IG	37	-1093.7394	2411.4931	0.0000	1.0000
SYMIG	42	-1079.6082	2413.5029	0.0000	1.0000
JC69	35	-1101.5526	2415.0106	0.0000	1.0000
SYM	40	-1087.8198	2417.8171	0.0000	1.0000
F81I	39	-1091.0848	2418.2927	0.0000	1.0000
F81G	39	-1091.1121	2418.3473	0.0000	1.0000
GTRI	44	-1076.8186	2420.0326	0.0000	1.0000
GTRG	44	-1076.9806	2420.3565	0.0000	1.0000
F81IG	40	-1091.0437	2424.2651	0.0000	1.0000
GTRIG	45	-1076.9922	2426.4341	0.0000	1.0000
F81	38	-1098.7441	2427.5568	0.0000	1.0000
GTR	43	-1085.2808	2430.9025	0.0000	1.0000



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

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

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


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

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


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

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


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

