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

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: K2P


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
HKY	39	-617.6526	1329.3876	0.5213	0.5213
HKYI	40	-617.3187	1331.6321	0.1697	0.6910
HKYG	40	-617.4230	1331.8409	0.1529	0.8439
HKYIG	41	-617.3230	1334.5834	0.0388	0.8827
K2P	36	-624.6256	1334.7741	0.0353	0.9180
GTR	43	-614.6278	1335.1713	0.0289	0.9469
K2PI	37	-623.9357	1336.2184	0.0171	0.9640
K2PG	37	-624.0803	1336.5075	0.0148	0.9788
GTRI	44	-614.4738	1337.9000	0.0074	0.9862
GTRG	44	-614.5416	1338.0355	0.0069	0.9931
K2PIG	38	-623.9400	1339.0800	0.0041	0.9972
GTRIG	45	-614.4781	1340.9775	0.0016	0.9988
SYM	40	-622.8639	1342.7225	0.0007	0.9995
SYMI	41	-622.3788	1344.6951	0.0002	0.9997
SYMG	41	-622.4993	1344.9361	0.0002	0.9999
SYMIG	42	-622.3835	1347.6780	0.0001	1.0000
F81	38	-632.3217	1355.8434	0.0000	1.0000
F81I	39	-631.8418	1357.7662	0.0000	1.0000
F81G	39	-631.9698	1358.0220	0.0000	1.0000
JC69	35	-638.5804	1359.8881	0.0000	1.0000
F81IG	40	-631.8459	1360.6866	0.0000	1.0000
JC69I	36	-637.8976	1361.3180	0.0000	1.0000
JC69G	36	-638.0418	1361.6064	0.0000	1.0000
JC69IG	37	-637.9025	1364.1520	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	39	-617.6526	1313.3052	0.3369	0.3369
HKYI	40	-617.3187	1314.6373	0.1731	0.5099
HKYG	40	-617.4230	1314.8461	0.1559	0.6658
GTR	43	-614.6278	1315.2555	0.1270	0.7929
HKYIG	41	-617.3230	1316.6459	0.0634	0.8563
GTRI	44	-614.4738	1316.9476	0.0545	0.9108
GTRG	44	-614.5416	1317.0831	0.0509	0.9617
GTRIG	45	-614.4781	1318.9563	0.0200	0.9817
K2P	36	-624.6256	1321.2513	0.0063	0.9880
K2PI	37	-623.9357	1321.8715	0.0046	0.9927
K2PG	37	-624.0803	1322.1605	0.0040	0.9967
K2PIG	38	-623.9400	1323.8800	0.0017	0.9984
SYM	40	-622.8639	1325.7277	0.0007	0.9991
SYMI	41	-622.3788	1326.7576	0.0004	0.9995
SYMG	41	-622.4993	1326.9986	0.0004	0.9999
SYMIG	42	-622.3835	1328.7670	0.0001	1.0000
F81	38	-632.3217	1340.6434	0.0000	1.0000
F81I	39	-631.8418	1341.6837	0.0000	1.0000
F81G	39	-631.9698	1341.9395	0.0000	1.0000
F81IG	40	-631.8459	1343.6918	0.0000	1.0000
JC69	35	-638.5804	1347.1608	0.0000	1.0000
JC69I	36	-637.8976	1347.7951	0.0000	1.0000
JC69G	36	-638.0418	1348.0836	0.0000	1.0000
JC69IG	37	-637.9025	1349.8051	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	36	-624.6256	1445.6429	0.6236	0.6236
HKY	39	-617.6526	1448.0627	0.1860	0.8096
K2PI	37	-623.9357	1449.7184	0.0813	0.8908
K2PG	37	-624.0803	1450.0074	0.0703	0.9612
HKYI	40	-617.3187	1452.8501	0.0170	0.9781
HKYG	40	-617.4230	1453.0589	0.0153	0.9934
K2PIG	38	-623.9400	1455.1822	0.0053	0.9987
HKYIG	41	-617.3230	1458.3141	0.0011	0.9998
GTR	43	-614.6278	1463.8343	0.0001	0.9999
SYM	40	-622.8639	1463.9406	0.0001	1.0000
JC69	35	-638.5804	1468.0970	0.0000	1.0000
SYMI	41	-622.3788	1468.4258	0.0000	1.0000
SYMG	41	-622.4993	1468.6667	0.0000	1.0000
GTRI	44	-614.4738	1468.9817	0.0000	1.0000
GTRG	44	-614.5416	1469.1173	0.0000	1.0000
F81	38	-632.3217	1471.9456	0.0000	1.0000
JC69I	36	-637.8976	1472.1867	0.0000	1.0000
JC69G	36	-638.0418	1472.4752	0.0000	1.0000
SYMIG	42	-622.3835	1473.8905	0.0000	1.0000
GTRIG	45	-614.4781	1474.4457	0.0000	1.0000
F81I	39	-631.8418	1476.4412	0.0000	1.0000
F81G	39	-631.9698	1476.6971	0.0000	1.0000
JC69IG	37	-637.9025	1477.6519	0.0000	1.0000
F81IG	40	-631.8459	1481.9047	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_1574.phylip = 1 - 234;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_1574.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


[Mrbayes block for the best AICc model (HKY)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_1574.phylip = 1 - 234;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_1574.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


[Mrbayes block for the best BIC model (K2P)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_1574.phylip = 1 - 234;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_1574.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

