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

Minimum AIC  model: HKYG
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	-447.7189	993.4719	0.3597	0.3597
K2PI	37	-448.0745	994.1831	0.2521	0.6117
K2P	36	-450.6103	995.7970	0.1125	0.7242
HKYG	40	-444.0568	996.8855	0.0653	0.7895
K2PIG	38	-447.6985	996.9487	0.0632	0.8527
HKYI	40	-444.4097	997.5914	0.0459	0.8986
HKY	39	-446.7951	998.7207	0.0261	0.9246
SYMG	41	-443.1359	998.7496	0.0257	0.9503
SYMI	41	-443.5700	999.6180	0.0166	0.9670
SYM	40	-445.6066	999.9852	0.0139	0.9808
HKYIG	41	-444.0431	1000.5641	0.0104	0.9912
SYMIG	42	-443.1311	1002.5123	0.0039	0.9951
GTRG	44	-440.5339	1005.0678	0.0011	0.9962
GTRI	44	-440.9147	1005.8295	0.0007	0.9970
GTR	43	-442.9321	1005.9543	0.0007	0.9977
JC69G	36	-455.7490	1006.0744	0.0007	0.9983
JC69I	36	-456.0633	1006.7029	0.0005	0.9988
F81G	39	-451.2468	1007.6241	0.0003	0.9991
JC69	35	-458.4482	1008.0728	0.0002	0.9993
F81I	39	-451.5358	1008.2021	0.0002	0.9996
GTRIG	45	-440.5349	1009.0514	0.0001	0.9997
JC69IG	37	-455.7393	1009.5128	0.0001	0.9998
F81	38	-454.0800	1009.7118	0.0001	1.0000
F81IG	40	-451.2361	1011.2442	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-444.0568	968.1136	0.1448	0.1448
SYMG	41	-443.1359	968.2717	0.1338	0.2785
HKYI	40	-444.4097	968.8195	0.1017	0.3802
GTRG	44	-440.5339	969.0678	0.0898	0.4701
SYMI	41	-443.5700	969.1401	0.0866	0.5567
K2PG	37	-447.7189	969.4378	0.0747	0.6314
GTRI	44	-440.9147	969.8295	0.0614	0.6927
HKYIG	41	-444.0431	970.0863	0.0540	0.7467
K2PI	37	-448.0745	970.1489	0.0523	0.7990
SYMIG	42	-443.1311	970.2623	0.0494	0.8485
GTRIG	45	-440.5349	971.0697	0.0330	0.8815
SYM	40	-445.6066	971.2132	0.0307	0.9122
K2PIG	38	-447.6985	971.3970	0.0280	0.9403
HKY	39	-446.7951	971.5903	0.0255	0.9657
GTR	43	-442.9321	971.8642	0.0222	0.9879
K2P	36	-450.6103	973.2207	0.0113	0.9992
F81G	39	-451.2468	980.4937	0.0003	0.9995
F81I	39	-451.5358	981.0717	0.0002	0.9997
F81IG	40	-451.2361	982.4723	0.0001	0.9998
JC69G	36	-455.7490	983.4981	0.0001	0.9999
JC69I	36	-456.0633	984.1266	0.0000	0.9999
F81	38	-454.0800	984.1601	0.0000	1.0000
JC69IG	37	-455.7393	985.4786	0.0000	1.0000
JC69	35	-458.4482	986.8963	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-447.7189	1082.0445	0.3924	0.3924
K2PI	37	-448.0745	1082.7557	0.2750	0.6674
K2P	36	-450.6103	1082.7840	0.2711	0.9385
K2PIG	38	-447.6985	1087.0471	0.0322	0.9706
HKYG	40	-444.0568	1089.8506	0.0079	0.9785
HKY	39	-446.7951	1090.2838	0.0064	0.9849
HKYI	40	-444.4097	1090.5565	0.0056	0.9905
SYM	40	-445.6066	1092.9502	0.0017	0.9922
SYMG	41	-443.1359	1093.0521	0.0016	0.9938
JC69G	36	-455.7490	1093.0614	0.0016	0.9953
JC69	35	-458.4482	1093.4162	0.0013	0.9967
JC69I	36	-456.0633	1093.6899	0.0012	0.9978
SYMI	41	-443.5700	1093.9205	0.0010	0.9989
HKYIG	41	-444.0431	1094.8667	0.0006	0.9995
JC69IG	37	-455.7393	1098.0853	0.0001	0.9996
SYMIG	42	-443.1311	1098.0861	0.0001	0.9998
F81G	39	-451.2468	1099.1873	0.0001	0.9999
F81I	39	-451.5358	1099.7653	0.0001	0.9999
F81	38	-454.0800	1099.8102	0.0001	1.0000
GTR	43	-442.9321	1102.7315	0.0000	1.0000
GTRG	44	-440.5339	1102.9785	0.0000	1.0000
GTRI	44	-440.9147	1103.7402	0.0000	1.0000
F81IG	40	-451.2361	1104.2093	0.0000	1.0000
GTRIG	45	-440.5349	1108.0239	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5857.phylip = 1 - 155;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5857.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) 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_5857.phylip = 1 - 155;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5857.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_5857.phylip = 1 - 155;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5857.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

