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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	40	-560.8036	1220.9013	0.6562	0.6562
HKYIG	41	-560.4468	1223.2722	0.2005	0.8568
HKYI	40	-562.5555	1224.4052	0.1138	0.9706
GTRG	44	-558.2220	1228.2995	0.0162	0.9868
GTRI	44	-558.9191	1229.6936	0.0081	0.9949
GTRIG	45	-557.7622	1230.6152	0.0051	1.0000
K2PG	37	-578.0134	1246.2811	0.0000	1.0000
K2PI	37	-578.4759	1247.2062	0.0000	1.0000
F81G	39	-575.5538	1247.3532	0.0000	1.0000
SYMG	41	-572.6793	1247.7372	0.0000	1.0000
SYMI	41	-572.9084	1248.1954	0.0000	1.0000
K2PIG	38	-578.0705	1249.3735	0.0000	1.0000
SYMIG	42	-572.3994	1250.2988	0.0000	1.0000
F81IG	40	-576.3228	1251.9398	0.0000	1.0000
F81I	39	-578.9836	1254.2127	0.0000	1.0000
HKY	39	-585.4170	1267.0796	0.0000	1.0000
GTR	43	-580.2702	1269.1991	0.0000	1.0000
JC69G	36	-594.6723	1276.6550	0.0000	1.0000
JC69IG	37	-595.5315	1281.3173	0.0000	1.0000
JC69I	36	-597.2384	1281.7872	0.0000	1.0000
SYM	40	-592.2621	1283.8184	0.0000	1.0000
F81	38	-598.6941	1290.6207	0.0000	1.0000
K2P	36	-602.1245	1291.5593	0.0000	1.0000
JC69	35	-617.9573	1320.3146	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-560.8036	1201.6072	0.4540	0.4540
HKYIG	41	-560.4468	1202.8935	0.2386	0.6926
GTRG	44	-558.2220	1204.4440	0.1099	0.8025
HKYI	40	-562.5555	1205.1110	0.0787	0.8812
GTRIG	45	-557.7622	1205.5243	0.0640	0.9453
GTRI	44	-558.9191	1205.8381	0.0547	1.0000
SYMG	41	-572.6793	1227.3585	0.0000	1.0000
SYMI	41	-572.9084	1227.8167	0.0000	1.0000
SYMIG	42	-572.3994	1228.7988	0.0000	1.0000
F81G	39	-575.5538	1229.1075	0.0000	1.0000
K2PG	37	-578.0134	1230.0267	0.0000	1.0000
K2PI	37	-578.4759	1230.9518	0.0000	1.0000
K2PIG	38	-578.0705	1232.1409	0.0000	1.0000
F81IG	40	-576.3228	1232.6457	0.0000	1.0000
F81I	39	-578.9836	1235.9671	0.0000	1.0000
GTR	43	-580.2702	1246.5404	0.0000	1.0000
HKY	39	-585.4170	1248.8340	0.0000	1.0000
JC69G	36	-594.6723	1261.3447	0.0000	1.0000
SYM	40	-592.2621	1264.5243	0.0000	1.0000
JC69IG	37	-595.5315	1265.0630	0.0000	1.0000
JC69I	36	-597.2384	1266.4769	0.0000	1.0000
F81	38	-598.6941	1273.3881	0.0000	1.0000
K2P	36	-602.1245	1276.2490	0.0000	1.0000
JC69	35	-617.9573	1305.9146	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-560.8036	1335.6815	0.7859	0.7859
HKYI	40	-562.5555	1339.1854	0.1363	0.9222
HKYIG	41	-560.4468	1340.3197	0.0773	0.9995
GTRG	44	-558.2220	1351.9258	0.0002	0.9997
GTRI	44	-558.9191	1353.3199	0.0001	0.9998
K2PG	37	-578.0134	1354.0455	0.0001	0.9999
K2PI	37	-578.4759	1354.9706	0.0001	1.0000
GTRIG	45	-557.7622	1356.3580	0.0000	1.0000
K2PIG	38	-578.0705	1359.5115	0.0000	1.0000
F81G	39	-575.5538	1359.8300	0.0000	1.0000
SYMG	41	-572.6793	1364.7847	0.0000	1.0000
SYMI	41	-572.9084	1365.2429	0.0000	1.0000
F81I	39	-578.9836	1366.6896	0.0000	1.0000
F81IG	40	-576.3228	1366.7200	0.0000	1.0000
SYMIG	42	-572.3994	1369.5768	0.0000	1.0000
HKY	39	-585.4170	1379.5565	0.0000	1.0000
JC69G	36	-594.6723	1382.0116	0.0000	1.0000
JC69I	36	-597.2384	1387.1438	0.0000	1.0000
JC69IG	37	-595.5315	1389.0817	0.0000	1.0000
GTR	43	-580.2702	1390.6703	0.0000	1.0000
K2P	36	-602.1245	1396.9159	0.0000	1.0000
SYM	40	-592.2621	1398.5986	0.0000	1.0000
F81	38	-598.6941	1400.7587	0.0000	1.0000
JC69	35	-617.9573	1423.2296	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/chr1_32308.phylip = 1 - 211;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32308.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 (HKYG)]

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

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


-------------------------------------------------------------
End of Output

