
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:50 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/chr2_23124.phylip" (19 taxa, 459 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	-1459.4311	3006.7090	0.6457	0.6457
GTRG	44	-1456.1132	3009.7916	0.1383	0.7840
HKYIG	41	-1460.2004	3010.6598	0.0896	0.8736
HKYI	40	-1461.4253	3010.6975	0.0879	0.9615
GTRIG	45	-1456.8028	3013.6298	0.0203	0.9817
GTRI	44	-1458.1462	3013.8576	0.0181	0.9998
K2PG	37	-1471.9386	3024.5565	0.0001	0.9999
SYMG	41	-1468.4746	3027.2082	0.0000	1.0000
SYMI	41	-1468.7477	3027.7544	0.0000	1.0000
K2PI	37	-1473.8225	3028.3243	0.0000	1.0000
K2PIG	38	-1472.9691	3028.9953	0.0000	1.0000
SYMIG	42	-1468.5671	3029.8170	0.0000	1.0000
HKY	39	-1482.8070	3051.0604	0.0000	1.0000
F81G	39	-1483.3859	3052.2181	0.0000	1.0000
GTR	43	-1480.2859	3055.6899	0.0000	1.0000
F81I	39	-1485.6203	3056.6870	0.0000	1.0000
F81IG	40	-1484.8232	3057.4932	0.0000	1.0000
JC69G	36	-1493.6891	3065.6911	0.0000	1.0000
K2P	36	-1494.8044	3067.9216	0.0000	1.0000
JC69I	36	-1495.9521	3070.2170	0.0000	1.0000
JC69IG	37	-1495.2988	3071.2769	0.0000	1.0000
SYM	40	-1492.2861	3072.4190	0.0000	1.0000
F81	38	-1504.9486	3092.9544	0.0000	1.0000
JC69	35	-1515.1084	3106.1742	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1459.4311	2998.8621	0.5071	0.5071
GTRG	44	-1456.1132	3000.2264	0.2564	0.7635
HKYIG	41	-1460.2004	3002.4008	0.0864	0.8500
HKYI	40	-1461.4253	3002.8506	0.0690	0.9190
GTRIG	45	-1456.8028	3003.6056	0.0473	0.9663
GTRI	44	-1458.1462	3004.2924	0.0336	0.9999
K2PG	37	-1471.9386	3017.8771	0.0000	0.9999
SYMG	41	-1468.4746	3018.9492	0.0000	1.0000
SYMI	41	-1468.7477	3019.4954	0.0000	1.0000
SYMIG	42	-1468.5671	3021.1343	0.0000	1.0000
K2PI	37	-1473.8225	3021.6450	0.0000	1.0000
K2PIG	38	-1472.9691	3021.9381	0.0000	1.0000
HKY	39	-1482.8070	3043.6141	0.0000	1.0000
F81G	39	-1483.3859	3044.7718	0.0000	1.0000
GTR	43	-1480.2859	3046.5718	0.0000	1.0000
F81I	39	-1485.6203	3049.2407	0.0000	1.0000
F81IG	40	-1484.8232	3049.6463	0.0000	1.0000
JC69G	36	-1493.6891	3059.3783	0.0000	1.0000
K2P	36	-1494.8044	3061.6088	0.0000	1.0000
JC69I	36	-1495.9521	3063.9042	0.0000	1.0000
SYM	40	-1492.2861	3064.5721	0.0000	1.0000
JC69IG	37	-1495.2988	3064.5976	0.0000	1.0000
F81	38	-1504.9486	3085.8973	0.0000	1.0000
JC69	35	-1515.1084	3100.2168	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1459.4311	3164.0241	0.8330	0.8330
HKYI	40	-1461.4253	3168.0126	0.1134	0.9464
K2PG	37	-1471.9386	3170.6520	0.0303	0.9767
HKYIG	41	-1460.2004	3171.6918	0.0180	0.9948
K2PI	37	-1473.8225	3174.4198	0.0046	0.9994
K2PIG	38	-1472.9691	3178.8420	0.0005	0.9999
GTRG	44	-1456.1132	3181.9046	0.0001	1.0000
GTRI	44	-1458.1462	3185.9706	0.0000	1.0000
SYMG	41	-1468.4746	3188.2403	0.0000	1.0000
SYMI	41	-1468.7477	3188.7865	0.0000	1.0000
GTRIG	45	-1456.8028	3189.4128	0.0000	1.0000
SYMIG	42	-1468.5671	3194.5544	0.0000	1.0000
HKY	39	-1482.8070	3204.6470	0.0000	1.0000
F81G	39	-1483.3859	3205.8047	0.0000	1.0000
JC69G	36	-1493.6891	3208.0241	0.0000	1.0000
K2P	36	-1494.8044	3210.2546	0.0000	1.0000
F81I	39	-1485.6203	3210.2736	0.0000	1.0000
JC69I	36	-1495.9521	3212.5500	0.0000	1.0000
F81IG	40	-1484.8232	3214.8083	0.0000	1.0000
JC69IG	37	-1495.2988	3217.3725	0.0000	1.0000
GTR	43	-1480.2859	3224.1209	0.0000	1.0000
SYM	40	-1492.2861	3229.7341	0.0000	1.0000
F81	38	-1504.9486	3242.8012	0.0000	1.0000
JC69	35	-1515.1084	3244.7335	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/chr2_23124.phylip = 1 - 459;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23124.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/chr2_23124.phylip = 1 - 459;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23124.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/chr2_23124.phylip = 1 - 459;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23124.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

