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

Minimum AIC  model: K2PI
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	37	-630.6786	1351.1550	0.6779	0.6779
K2PG	37	-631.9796	1353.7570	0.1846	0.8625
K2PIG	38	-631.4924	1355.7306	0.0688	0.9313
HKYG	40	-629.3873	1357.5175	0.0282	0.9594
HKYI	40	-629.8077	1358.3582	0.0185	0.9779
SYMG	41	-628.9605	1359.7140	0.0094	0.9873
SYMI	41	-629.1923	1360.1776	0.0074	0.9948
HKYIG	41	-630.5220	1362.8371	0.0020	0.9967
K2P	36	-638.2818	1363.4464	0.0015	0.9982
SYMIG	42	-629.4367	1363.7521	0.0012	0.9994
GTRG	44	-627.8634	1366.8846	0.0003	0.9997
GTRI	44	-628.1915	1367.5410	0.0002	0.9999
SYM	40	-635.7417	1370.2262	0.0000	0.9999
HKY	39	-637.4009	1370.5292	0.0000	1.0000
GTRIG	45	-628.3521	1371.0571	0.0000	1.0000
GTR	43	-634.8067	1377.6133	0.0000	1.0000
JC69G	36	-648.3458	1383.5742	0.0000	1.0000
JC69I	36	-648.6867	1384.2560	0.0000	1.0000
JC69IG	37	-648.1857	1386.1692	0.0000	1.0000
F81G	39	-647.0400	1389.8073	0.0000	1.0000
F81I	39	-647.3713	1390.4698	0.0000	1.0000
F81IG	40	-646.8885	1392.5199	0.0000	1.0000
JC69	35	-655.9856	1395.9712	0.0000	1.0000
F81	38	-654.7240	1402.1938	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	37	-630.6786	1335.3572	0.5013	0.5013
K2PG	37	-631.9796	1337.9593	0.1365	0.6378
HKYG	40	-629.3873	1338.7747	0.0908	0.7286
K2PIG	38	-631.4924	1338.9848	0.0817	0.8104
HKYI	40	-629.8077	1339.6153	0.0596	0.8700
SYMG	41	-628.9605	1339.9209	0.0512	0.9212
SYMI	41	-629.1923	1340.3845	0.0406	0.9618
SYMIG	42	-629.4367	1342.8735	0.0117	0.9735
HKYIG	41	-630.5220	1343.0440	0.0107	0.9842
GTRG	44	-627.8634	1343.7267	0.0076	0.9918
GTRI	44	-628.1915	1344.3831	0.0055	0.9973
GTRIG	45	-628.3521	1346.7042	0.0017	0.9991
K2P	36	-638.2818	1348.5637	0.0007	0.9997
SYM	40	-635.7417	1351.4834	0.0002	0.9999
HKY	39	-637.4009	1352.8019	0.0001	1.0000
GTR	43	-634.8067	1355.6133	0.0000	1.0000
JC69G	36	-648.3458	1368.6916	0.0000	1.0000
JC69I	36	-648.6867	1369.3734	0.0000	1.0000
JC69IG	37	-648.1857	1370.3715	0.0000	1.0000
F81G	39	-647.0400	1372.0800	0.0000	1.0000
F81I	39	-647.3713	1372.7425	0.0000	1.0000
F81IG	40	-646.8885	1373.7771	0.0000	1.0000
JC69	35	-655.9856	1381.9712	0.0000	1.0000
F81	38	-654.7240	1385.4480	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-630.6786	1460.2425	0.7623	0.7623
K2PG	37	-631.9796	1462.8446	0.2075	0.9698
K2PIG	38	-631.4924	1467.2454	0.0230	0.9928
K2P	36	-638.2818	1470.0737	0.0056	0.9984
HKYG	40	-629.3873	1473.7858	0.0009	0.9992
HKYI	40	-629.8077	1474.6265	0.0006	0.9998
SYMG	41	-628.9605	1478.3074	0.0001	0.9999
SYMI	41	-629.1923	1478.7710	0.0001	1.0000
HKYIG	41	-630.5220	1481.4304	0.0000	1.0000
HKY	39	-637.4009	1484.4377	0.0000	1.0000
SYMIG	42	-629.4367	1484.6352	0.0000	1.0000
SYM	40	-635.7417	1486.4945	0.0000	1.0000
JC69G	36	-648.3458	1490.2016	0.0000	1.0000
JC69I	36	-648.6867	1490.8834	0.0000	1.0000
GTRG	44	-627.8634	1492.2390	0.0000	1.0000
GTRI	44	-628.1915	1492.8953	0.0000	1.0000
JC69IG	37	-648.1857	1495.2568	0.0000	1.0000
GTRIG	45	-628.3521	1498.5917	0.0000	1.0000
JC69	35	-655.9856	1500.1060	0.0000	1.0000
GTR	43	-634.8067	1500.7503	0.0000	1.0000
F81G	39	-647.0400	1503.7159	0.0000	1.0000
F81I	39	-647.3713	1504.3784	0.0000	1.0000
F81IG	40	-646.8885	1508.7882	0.0000	1.0000
F81	38	-654.7240	1513.7086	0.0000	1.0000



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

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

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


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

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


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

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


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

