
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:22 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/chr1_2930.phylip" (19 taxa, 972 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	-5215.9514	10524.1746	0.8818	0.8818
GTRIG	45	-5216.8734	10528.2177	0.1168	0.9985
HKYIG	41	-5225.8778	10537.4589	0.0012	0.9997
GTRI	44	-5224.3374	10540.9467	0.0002	0.9999
HKYG	40	-5229.5438	10542.6106	0.0001	1.0000
SYMIG	42	-5229.9293	10547.7467	0.0000	1.0000
K2PG	37	-5236.2393	10549.4892	0.0000	1.0000
HKYI	40	-5233.0767	10549.6765	0.0000	1.0000
SYMG	41	-5233.1783	10552.0598	0.0000	1.0000
K2PIG	38	-5237.6150	10554.4069	0.0000	1.0000
K2PI	37	-5244.3728	10565.7563	0.0000	1.0000
SYMI	41	-5240.9459	10567.5949	0.0000	1.0000
GTR	43	-5357.0551	10804.1877	0.0000	1.0000
HKY	39	-5364.7417	10810.8310	0.0000	1.0000
F81IG	40	-5366.1686	10815.8603	0.0000	1.0000
F81G	39	-5368.8452	10819.0381	0.0000	1.0000
K2P	36	-5375.3077	10825.4647	0.0000	1.0000
SYM	40	-5372.3861	10828.2953	0.0000	1.0000
F81I	39	-5376.0018	10833.3512	0.0000	1.0000
JC69IG	37	-5379.6421	10836.2950	0.0000	1.0000
JC69G	36	-5382.0643	10838.9777	0.0000	1.0000
JC69I	36	-5390.0596	10854.9685	0.0000	1.0000
F81	38	-5494.5029	11068.1826	0.0000	1.0000
JC69	35	-5508.6416	11089.9755	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-5215.9514	10519.9028	0.8714	0.8714
GTRIG	45	-5216.8734	10523.7468	0.1275	0.9989
HKYIG	41	-5225.8778	10533.7557	0.0009	0.9997
GTRI	44	-5224.3374	10536.6749	0.0002	0.9999
HKYG	40	-5229.5438	10539.0875	0.0001	1.0000
SYMIG	42	-5229.9293	10543.8586	0.0000	1.0000
HKYI	40	-5233.0767	10546.1534	0.0000	1.0000
K2PG	37	-5236.2393	10546.4785	0.0000	1.0000
SYMG	41	-5233.1783	10548.3566	0.0000	1.0000
K2PIG	38	-5237.6150	10551.2301	0.0000	1.0000
K2PI	37	-5244.3728	10562.7456	0.0000	1.0000
SYMI	41	-5240.9459	10563.8917	0.0000	1.0000
GTR	43	-5357.0551	10800.1101	0.0000	1.0000
HKY	39	-5364.7417	10807.4833	0.0000	1.0000
F81IG	40	-5366.1686	10812.3372	0.0000	1.0000
F81G	39	-5368.8452	10815.6904	0.0000	1.0000
K2P	36	-5375.3077	10822.6155	0.0000	1.0000
SYM	40	-5372.3861	10824.7722	0.0000	1.0000
F81I	39	-5376.0018	10830.0035	0.0000	1.0000
JC69IG	37	-5379.6421	10833.2843	0.0000	1.0000
JC69G	36	-5382.0643	10836.1285	0.0000	1.0000
JC69I	36	-5390.0596	10852.1193	0.0000	1.0000
F81	38	-5494.5029	11065.0058	0.0000	1.0000
JC69	35	-5508.6416	11087.2832	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-5236.2393	10727.0147	0.9155	0.9155
HKYIG	41	-5225.8778	10733.8093	0.0306	0.9462
HKYG	40	-5229.5438	10734.2618	0.0244	0.9706
GTRG	44	-5215.9514	10734.5945	0.0207	0.9913
K2PIG	38	-5237.6150	10736.6456	0.0074	0.9987
HKYI	40	-5233.0767	10741.3277	0.0007	0.9994
K2PI	37	-5244.3728	10743.2818	0.0003	0.9997
GTRIG	45	-5216.8734	10743.3179	0.0003	1.0000
SYMG	41	-5233.1783	10748.4102	0.0000	1.0000
SYMIG	42	-5229.9293	10748.7916	0.0000	1.0000
GTRI	44	-5224.3374	10751.3666	0.0000	1.0000
SYMI	41	-5240.9459	10763.9453	0.0000	1.0000
HKY	39	-5364.7417	10997.7782	0.0000	1.0000
K2P	36	-5375.3077	10998.2723	0.0000	1.0000
F81G	39	-5368.8452	11005.9853	0.0000	1.0000
F81IG	40	-5366.1686	11007.5114	0.0000	1.0000
GTR	43	-5357.0551	11009.9224	0.0000	1.0000
JC69G	36	-5382.0643	11011.7853	0.0000	1.0000
JC69IG	37	-5379.6421	11013.8205	0.0000	1.0000
SYM	40	-5372.3861	11019.9464	0.0000	1.0000
F81I	39	-5376.0018	11020.2984	0.0000	1.0000
JC69I	36	-5390.0596	11027.7761	0.0000	1.0000
F81	38	-5494.5029	11250.4213	0.0000	1.0000
JC69	35	-5508.6416	11258.0607	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/chr1_2930.phylip = 1 - 972;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_2930.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/chr1_2930.phylip = 1 - 972;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_2930.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/chr1_2930.phylip = 1 - 972;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_2930.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

