
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:41 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/chr7_10512.phylip" (19 taxa, 699 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-2482.3516	5058.7582	0.6682	0.6682
GTRIG	45	-2481.9210	5060.1819	0.3279	0.9961
HKYG	40	-2492.6460	5070.2769	0.0021	0.9982
GTRI	44	-2488.5837	5071.2225	0.0013	0.9995
HKYIG	41	-2493.2034	5073.6487	0.0004	0.9999
SYMG	41	-2495.2060	5077.6541	0.0001	1.0000
SYMIG	42	-2494.7390	5078.9842	0.0000	1.0000
K2PG	37	-2501.2270	5080.7082	0.0000	1.0000
K2PIG	38	-2501.6429	5083.7767	0.0000	1.0000
HKYI	40	-2500.6502	5086.2852	0.0000	1.0000
SYMI	41	-2501.5662	5090.3745	0.0000	1.0000
K2PI	37	-2508.4652	5095.1846	0.0000	1.0000
GTR	43	-2527.7400	5147.2571	0.0000	1.0000
HKY	39	-2538.6283	5159.9910	0.0000	1.0000
SYM	40	-2540.3329	5165.6507	0.0000	1.0000
K2P	36	-2547.6676	5171.3593	0.0000	1.0000
F81G	39	-2562.3381	5207.4107	0.0000	1.0000
F81IG	40	-2561.9567	5208.8982	0.0000	1.0000
JC69G	36	-2568.1923	5212.4087	0.0000	1.0000
JC69IG	37	-2567.7946	5213.8433	0.0000	1.0000
F81I	39	-2568.4869	5219.7083	0.0000	1.0000
JC69I	36	-2574.5418	5225.1077	0.0000	1.0000
F81	38	-2604.4949	5289.4806	0.0000	1.0000
JC69	35	-2610.9399	5295.6807	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2482.3516	5052.7032	0.6369	0.6369
GTRIG	45	-2481.9210	5053.8419	0.3604	0.9973
GTRI	44	-2488.5837	5065.1674	0.0013	0.9985
HKYG	40	-2492.6460	5065.2921	0.0012	0.9997
HKYIG	41	-2493.2034	5068.4067	0.0002	0.9999
SYMG	41	-2495.2060	5072.4121	0.0000	1.0000
SYMIG	42	-2494.7390	5073.4781	0.0000	1.0000
K2PG	37	-2501.2270	5076.4541	0.0000	1.0000
K2PIG	38	-2501.6429	5079.2858	0.0000	1.0000
HKYI	40	-2500.6502	5081.3004	0.0000	1.0000
SYMI	41	-2501.5662	5085.1325	0.0000	1.0000
K2PI	37	-2508.4652	5090.9304	0.0000	1.0000
GTR	43	-2527.7400	5141.4800	0.0000	1.0000
HKY	39	-2538.6283	5155.2565	0.0000	1.0000
SYM	40	-2540.3329	5160.6659	0.0000	1.0000
K2P	36	-2547.6676	5167.3351	0.0000	1.0000
F81G	39	-2562.3381	5202.6763	0.0000	1.0000
F81IG	40	-2561.9567	5203.9134	0.0000	1.0000
JC69G	36	-2568.1923	5208.3846	0.0000	1.0000
JC69IG	37	-2567.7946	5209.5891	0.0000	1.0000
F81I	39	-2568.4869	5214.9739	0.0000	1.0000
JC69I	36	-2574.5418	5221.0835	0.0000	1.0000
F81	38	-2604.4949	5284.9897	0.0000	1.0000
JC69	35	-2610.9399	5291.8798	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-2501.2270	5244.7911	0.7464	0.7464
HKYG	40	-2492.6460	5247.2781	0.2152	0.9616
K2PIG	38	-2501.6429	5252.1726	0.0186	0.9803
GTRG	44	-2482.3516	5252.8878	0.0130	0.9933
HKYIG	41	-2493.2034	5254.9424	0.0047	0.9979
GTRIG	45	-2481.9210	5258.5762	0.0008	0.9987
SYMG	41	-2495.2060	5258.9478	0.0006	0.9993
K2PI	37	-2508.4652	5259.2675	0.0005	0.9999
HKYI	40	-2500.6502	5263.2864	0.0001	0.9999
SYMIG	42	-2494.7390	5264.5634	0.0000	1.0000
GTRI	44	-2488.5837	5265.3521	0.0000	1.0000
SYMI	41	-2501.5662	5271.6682	0.0000	1.0000
K2P	36	-2547.6676	5331.1226	0.0000	1.0000
HKY	39	-2538.6283	5332.6929	0.0000	1.0000
GTR	43	-2527.7400	5337.1150	0.0000	1.0000
SYM	40	-2540.3329	5342.6519	0.0000	1.0000
JC69G	36	-2568.1923	5372.1720	0.0000	1.0000
JC69IG	37	-2567.7946	5377.9262	0.0000	1.0000
F81G	39	-2562.3381	5380.1126	0.0000	1.0000
JC69I	36	-2574.5418	5384.8709	0.0000	1.0000
F81IG	40	-2561.9567	5385.8994	0.0000	1.0000
F81I	39	-2568.4869	5392.4103	0.0000	1.0000
JC69	35	-2610.9399	5451.1176	0.0000	1.0000
F81	38	-2604.4949	5457.8764	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10512.phylip = 1 - 699;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10512.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10512.phylip = 1 - 699;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10512.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
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/chr7_10512.phylip = 1 - 699;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10512.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

