
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:33 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/chr18_3214.phylip" (19 taxa, 435 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	-1696.7535	3491.6608	0.7071	0.7071
GTRIG	45	-1696.6411	3493.9249	0.2280	0.9351
HKYG	40	-1705.4383	3499.2015	0.0163	0.9514
SYMG	41	-1704.4472	3499.6578	0.0130	0.9643
K2PG	37	-1709.3276	3499.7383	0.0125	0.9768
HKYIG	41	-1704.9446	3500.6526	0.0079	0.9847
GTRI	44	-1701.6152	3501.3841	0.0055	0.9902
K2PIG	38	-1709.1187	3501.7223	0.0046	0.9948
SYMIG	42	-1704.2992	3501.8127	0.0044	0.9992
HKYI	40	-1709.2681	3506.8610	0.0004	0.9995
K2PI	37	-1713.2266	3507.5363	0.0003	0.9998
SYMI	41	-1708.6021	3507.9677	0.0002	1.0000
GTR	43	-1725.3909	3546.4596	0.0000	1.0000
HKY	39	-1734.4663	3554.8314	0.0000	1.0000
SYM	40	-1734.2157	3556.7563	0.0000	1.0000
K2P	36	-1740.0393	3558.7720	0.0000	1.0000
JC69G	36	-1746.8618	3572.4171	0.0000	1.0000
F81G	39	-1743.6376	3573.1739	0.0000	1.0000
JC69IG	37	-1746.7415	3574.5660	0.0000	1.0000
F81IG	40	-1743.4621	3575.2490	0.0000	1.0000
JC69I	36	-1750.5188	3579.7310	0.0000	1.0000
F81I	39	-1747.1326	3580.1640	0.0000	1.0000
F81	38	-1771.5778	3626.6405	0.0000	1.0000
JC69	35	-1775.5557	3627.4273	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1696.7535	3481.5069	0.6882	0.6882
GTRIG	45	-1696.6411	3483.2822	0.2833	0.9714
HKYG	40	-1705.4383	3490.8766	0.0064	0.9778
SYMG	41	-1704.4472	3490.8944	0.0063	0.9841
GTRI	44	-1701.6152	3491.2303	0.0053	0.9894
HKYIG	41	-1704.9446	3491.8892	0.0038	0.9932
SYMIG	42	-1704.2992	3492.5985	0.0027	0.9959
K2PG	37	-1709.3276	3492.6552	0.0026	0.9985
K2PIG	38	-1709.1187	3494.2375	0.0012	0.9997
HKYI	40	-1709.2681	3498.5361	0.0001	0.9998
SYMI	41	-1708.6021	3499.2043	0.0001	0.9999
K2PI	37	-1713.2266	3500.4532	0.0001	1.0000
GTR	43	-1725.3909	3536.7819	0.0000	1.0000
HKY	39	-1734.4663	3546.9327	0.0000	1.0000
SYM	40	-1734.2157	3548.4314	0.0000	1.0000
K2P	36	-1740.0393	3552.0785	0.0000	1.0000
F81G	39	-1743.6376	3565.2752	0.0000	1.0000
JC69G	36	-1746.8618	3565.7237	0.0000	1.0000
F81IG	40	-1743.4621	3566.9242	0.0000	1.0000
JC69IG	37	-1746.7415	3567.4829	0.0000	1.0000
F81I	39	-1747.1326	3572.2652	0.0000	1.0000
JC69I	36	-1750.5188	3573.0375	0.0000	1.0000
F81	38	-1771.5778	3619.1557	0.0000	1.0000
JC69	35	-1775.5557	3621.1115	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1709.3276	3643.4430	0.9206	0.9206
K2PIG	38	-1709.1187	3649.1006	0.0544	0.9750
K2PI	37	-1713.2266	3651.2410	0.0187	0.9937
HKYG	40	-1705.4383	3653.8905	0.0050	0.9987
SYMG	41	-1704.4472	3657.9836	0.0006	0.9993
HKYIG	41	-1704.9446	3658.9784	0.0004	0.9997
GTRG	44	-1696.7535	3660.8222	0.0002	0.9998
HKYI	40	-1709.2681	3661.5500	0.0001	0.9999
SYMIG	42	-1704.2992	3663.7630	0.0000	1.0000
SYMI	41	-1708.6021	3666.2935	0.0000	1.0000
GTRIG	45	-1696.6411	3666.6728	0.0000	1.0000
GTRI	44	-1701.6152	3670.5455	0.0000	1.0000
K2P	36	-1740.0393	3698.7910	0.0000	1.0000
HKY	39	-1734.4663	3705.8712	0.0000	1.0000
SYM	40	-1734.2157	3711.4453	0.0000	1.0000
GTR	43	-1725.3909	3712.0218	0.0000	1.0000
JC69G	36	-1746.8618	3712.4361	0.0000	1.0000
JC69IG	37	-1746.7415	3718.2707	0.0000	1.0000
JC69I	36	-1750.5188	3719.7500	0.0000	1.0000
F81G	39	-1743.6376	3724.2137	0.0000	1.0000
F81IG	40	-1743.4621	3729.9380	0.0000	1.0000
F81I	39	-1747.1326	3731.2037	0.0000	1.0000
JC69	35	-1775.5557	3763.7486	0.0000	1.0000
F81	38	-1771.5778	3774.0188	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/chr18_3214.phylip = 1 - 435;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr18_3214.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/chr18_3214.phylip = 1 - 435;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr18_3214.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/chr18_3214.phylip = 1 - 435;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr18_3214.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

