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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: SYMIG


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
GTRIG	45	-2511.7817	5118.8440	0.9900	0.9900
SYMIG	42	-2520.2038	5128.9971	0.0062	0.9962
GTRG	44	-2518.4924	5130.0294	0.0037	0.9999
GTRI	44	-2522.3744	5137.7933	0.0001	1.0000
SYMG	41	-2526.2589	5138.8883	0.0000	1.0000
SYMI	41	-2529.5205	5145.4116	0.0000	1.0000
HKYIG	41	-2532.9375	5152.2456	0.0000	1.0000
K2PIG	38	-2536.7721	5153.2913	0.0000	1.0000
K2PG	37	-2540.7486	5159.0477	0.0000	1.0000
HKYG	40	-2538.6617	5161.4805	0.0000	1.0000
HKYI	40	-2541.7777	5167.7125	0.0000	1.0000
K2PI	37	-2545.4429	5168.4364	0.0000	1.0000
JC69IG	37	-2578.7126	5234.9758	0.0000	1.0000
F81IG	40	-2575.5892	5235.3355	0.0000	1.0000
JC69G	36	-2583.9315	5243.2224	0.0000	1.0000
F81G	39	-2580.8464	5243.6422	0.0000	1.0000
JC69I	36	-2586.9897	5249.3389	0.0000	1.0000
F81I	39	-2583.9641	5249.8775	0.0000	1.0000
GTR	43	-2596.5417	5283.8976	0.0000	1.0000
SYM	40	-2603.3017	5290.7605	0.0000	1.0000
HKY	39	-2612.7497	5307.4488	0.0000	1.0000
K2P	36	-2616.5416	5308.4425	0.0000	1.0000
F81	38	-2653.7368	5387.2208	0.0000	1.0000
JC69	35	-2657.0660	5387.3058	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-2511.7817	5113.5634	0.9922	0.9922
SYMIG	42	-2520.2038	5124.4075	0.0044	0.9966
GTRG	44	-2518.4924	5124.9849	0.0033	0.9999
GTRI	44	-2522.3744	5132.7487	0.0001	1.0000
SYMG	41	-2526.2589	5134.5177	0.0000	1.0000
SYMI	41	-2529.5205	5141.0410	0.0000	1.0000
HKYIG	41	-2532.9375	5147.8750	0.0000	1.0000
K2PIG	38	-2536.7721	5149.5442	0.0000	1.0000
K2PG	37	-2540.7486	5155.4972	0.0000	1.0000
HKYG	40	-2538.6617	5157.3234	0.0000	1.0000
HKYI	40	-2541.7777	5163.5553	0.0000	1.0000
K2PI	37	-2545.4429	5164.8859	0.0000	1.0000
F81IG	40	-2575.5892	5231.1783	0.0000	1.0000
JC69IG	37	-2578.7126	5231.4253	0.0000	1.0000
F81G	39	-2580.8464	5239.6929	0.0000	1.0000
JC69G	36	-2583.9315	5239.8630	0.0000	1.0000
F81I	39	-2583.9641	5245.9282	0.0000	1.0000
JC69I	36	-2586.9897	5245.9795	0.0000	1.0000
GTR	43	-2596.5417	5279.0834	0.0000	1.0000
SYM	40	-2603.3017	5286.6033	0.0000	1.0000
HKY	39	-2612.7497	5303.4995	0.0000	1.0000
K2P	36	-2616.5416	5305.0831	0.0000	1.0000
F81	38	-2653.7368	5383.4736	0.0000	1.0000
JC69	35	-2657.0660	5384.1320	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMIG	42	-2520.2038	5322.7074	0.7490	0.7490
GTRIG	45	-2511.7817	5326.0275	0.1424	0.8914
SYMG	41	-2526.2589	5328.0962	0.0506	0.9421
K2PIG	38	-2536.7721	5328.9584	0.0329	0.9750
K2PG	37	-2540.7486	5330.1900	0.0178	0.9927
GTRG	44	-2518.4924	5332.7276	0.0050	0.9977
SYMI	41	-2529.5205	5334.6195	0.0019	0.9997
K2PI	37	-2545.4429	5339.5786	0.0002	0.9998
GTRI	44	-2522.3744	5340.4915	0.0001	0.9999
HKYIG	41	-2532.9375	5341.4535	0.0001	1.0000
HKYG	40	-2538.6617	5346.1804	0.0000	1.0000
HKYI	40	-2541.7777	5352.4124	0.0000	1.0000
JC69IG	37	-2578.7126	5406.1180	0.0000	1.0000
JC69G	36	-2583.9315	5409.8343	0.0000	1.0000
JC69I	36	-2586.9897	5415.9508	0.0000	1.0000
F81IG	40	-2575.5892	5420.0353	0.0000	1.0000
F81G	39	-2580.8464	5423.8285	0.0000	1.0000
F81I	39	-2583.9641	5430.0638	0.0000	1.0000
K2P	36	-2616.5416	5475.0544	0.0000	1.0000
SYM	40	-2603.3017	5475.4603	0.0000	1.0000
GTR	43	-2596.5417	5482.1047	0.0000	1.0000
HKY	39	-2612.7497	5487.6351	0.0000	1.0000
JC69	35	-2657.0660	5549.3819	0.0000	1.0000
F81	38	-2653.7368	5562.8878	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_13070.phylip = 1 - 830;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_13070.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_13070.phylip = 1 - 830;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_13070.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_13070.phylip = 1 - 830;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_13070.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

