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

Minimum AIC  model: K2PG
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	-564.8278	1218.3781	0.4087	0.4087
K2PI	37	-564.9968	1218.7161	0.3452	0.7539
K2PIG	38	-564.6727	1220.9453	0.1132	0.8671
SYMG	41	-561.1628	1222.7428	0.0461	0.9132
SYMI	41	-561.3280	1223.0731	0.0391	0.9523
SYMIG	42	-560.9525	1225.3243	0.0127	0.9650
HKYG	40	-564.0053	1225.4574	0.0119	0.9768
HKYI	40	-564.1809	1225.8087	0.0100	0.9868
K2P	36	-570.3144	1226.5039	0.0070	0.9938
HKYIG	41	-563.8539	1228.1250	0.0031	0.9970
GTRG	44	-560.4012	1230.3241	0.0010	0.9980
GTRI	44	-560.5680	1230.6578	0.0009	0.9989
SYM	40	-567.1743	1231.7954	0.0005	0.9994
GTRIG	45	-560.1897	1233.0024	0.0003	0.9996
HKY	39	-569.3879	1233.2837	0.0002	0.9999
JC69G	36	-575.2994	1236.4738	0.0000	0.9999
JC69I	36	-575.4625	1236.8000	0.0000	1.0000
JC69IG	37	-575.1643	1239.0511	0.0000	1.0000
GTR	43	-566.3415	1239.1370	0.0000	1.0000
F81G	39	-574.2733	1243.0545	0.0000	1.0000
F81I	39	-574.4414	1243.3906	0.0000	1.0000
JC69	35	-580.6844	1244.4259	0.0000	1.0000
F81IG	40	-574.1349	1245.7166	0.0000	1.0000
F81	38	-579.6522	1250.9044	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-564.8278	1203.6556	0.2288	0.2288
K2PI	37	-564.9968	1203.9935	0.1932	0.4220
SYMG	41	-561.1628	1204.3256	0.1636	0.5856
SYMI	41	-561.3280	1204.6559	0.1387	0.7243
K2PIG	38	-564.6727	1205.3453	0.0983	0.8226
SYMIG	42	-560.9525	1205.9050	0.0743	0.8969
HKYG	40	-564.0053	1208.0106	0.0259	0.9228
HKYI	40	-564.1809	1208.3619	0.0217	0.9446
GTRG	44	-560.4012	1208.8023	0.0175	0.9620
GTRI	44	-560.5680	1209.1361	0.0148	0.9768
HKYIG	41	-563.8539	1209.7079	0.0111	0.9879
GTRIG	45	-560.1897	1210.3794	0.0079	0.9958
K2P	36	-570.3144	1212.6289	0.0026	0.9984
SYM	40	-567.1743	1214.3486	0.0011	0.9995
HKY	39	-569.3879	1216.7758	0.0003	0.9998
GTR	43	-566.3415	1218.6829	0.0001	1.0000
JC69G	36	-575.2994	1222.5988	0.0000	1.0000
JC69I	36	-575.4625	1222.9250	0.0000	1.0000
JC69IG	37	-575.1643	1224.3286	0.0000	1.0000
F81G	39	-574.2733	1226.5465	0.0000	1.0000
F81I	39	-574.4414	1226.8827	0.0000	1.0000
F81IG	40	-574.1349	1228.2698	0.0000	1.0000
JC69	35	-580.6844	1231.3689	0.0000	1.0000
F81	38	-579.6522	1235.3044	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-564.8278	1330.7033	0.5030	0.5030
K2PI	37	-564.9968	1331.0413	0.4248	0.9278
K2PIG	38	-564.6727	1335.8268	0.0388	0.9666
K2P	36	-570.3144	1336.2429	0.0315	0.9982
SYMG	41	-561.1628	1345.1082	0.0004	0.9986
HKYG	40	-564.0053	1345.3595	0.0003	0.9989
SYMI	41	-561.3280	1345.4385	0.0003	0.9992
HKYI	40	-564.1809	1345.7107	0.0003	0.9995
JC69G	36	-575.2994	1346.2128	0.0002	0.9997
JC69I	36	-575.4625	1346.5390	0.0002	0.9999
SYMIG	42	-560.9525	1350.1213	0.0000	0.9999
HKYIG	41	-563.8539	1350.4905	0.0000	0.9999
HKY	39	-569.3879	1350.6910	0.0000	1.0000
JC69IG	37	-575.1643	1351.3763	0.0000	1.0000
JC69	35	-580.6844	1351.5491	0.0000	1.0000
SYM	40	-567.1743	1351.6975	0.0000	1.0000
GTRG	44	-560.4012	1359.8861	0.0000	1.0000
GTRI	44	-560.5680	1360.2198	0.0000	1.0000
F81G	39	-574.2733	1360.4617	0.0000	1.0000
F81I	39	-574.4414	1360.7979	0.0000	1.0000
GTRIG	45	-560.1897	1364.8969	0.0000	1.0000
F81IG	40	-574.1349	1365.6187	0.0000	1.0000
F81	38	-579.6522	1365.7858	0.0000	1.0000
GTR	43	-566.3415	1366.3329	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14664.phylip = 1 - 229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14664.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 AICc model (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14664.phylip = 1 - 229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14664.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/chr5_14664.phylip = 1 - 229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14664.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

