
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:28 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/chr6_3228.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	-1513.7523	3115.3515	0.4420	0.4420
HKYIG	41	-1512.8110	3115.8811	0.3392	0.7813
HKYI	40	-1515.0470	3117.9408	0.1211	0.9024
GTRG	44	-1511.0764	3119.7180	0.0498	0.9522
GTRIG	45	-1510.1364	3120.2970	0.0373	0.9895
GTRI	44	-1512.6315	3122.8282	0.0105	1.0000
K2PG	37	-1528.7325	3138.1444	0.0000	1.0000
K2PIG	38	-1527.7336	3138.5243	0.0000	1.0000
K2PI	37	-1529.9107	3140.5008	0.0000	1.0000
SYMG	41	-1527.5461	3145.3512	0.0000	1.0000
SYMIG	42	-1526.5224	3145.7275	0.0000	1.0000
SYMI	41	-1528.8731	3148.0051	0.0000	1.0000
F81G	39	-1550.9802	3187.4068	0.0000	1.0000
F81IG	40	-1550.2713	3188.3894	0.0000	1.0000
HKY	39	-1551.6224	3188.6910	0.0000	1.0000
F81I	39	-1551.9251	3189.2965	0.0000	1.0000
GTR	43	-1547.5134	3190.1450	0.0000	1.0000
JC69G	36	-1564.8583	3208.0294	0.0000	1.0000
JC69IG	37	-1564.1282	3208.9358	0.0000	1.0000
JC69I	36	-1565.8077	3209.9282	0.0000	1.0000
K2P	36	-1569.1076	3216.5280	0.0000	1.0000
SYM	40	-1566.8090	3221.4649	0.0000	1.0000
F81	38	-1587.4728	3258.0027	0.0000	1.0000
JC69	35	-1602.3812	3280.7198	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1513.7523	3107.5047	0.3585	0.3585
HKYIG	41	-1512.8110	3107.6221	0.3380	0.6965
HKYI	40	-1515.0470	3110.0940	0.0982	0.7947
GTRG	44	-1511.0764	3110.1528	0.0954	0.8900
GTRIG	45	-1510.1364	3110.2728	0.0898	0.9799
GTRI	44	-1512.6315	3113.2630	0.0201	1.0000
K2PG	37	-1528.7325	3131.4651	0.0000	1.0000
K2PIG	38	-1527.7336	3131.4672	0.0000	1.0000
K2PI	37	-1529.9107	3133.8215	0.0000	1.0000
SYMIG	42	-1526.5224	3137.0448	0.0000	1.0000
SYMG	41	-1527.5461	3137.0922	0.0000	1.0000
SYMI	41	-1528.8731	3139.7461	0.0000	1.0000
F81G	39	-1550.9802	3179.9605	0.0000	1.0000
F81IG	40	-1550.2713	3180.5425	0.0000	1.0000
GTR	43	-1547.5134	3181.0269	0.0000	1.0000
HKY	39	-1551.6224	3181.2447	0.0000	1.0000
F81I	39	-1551.9251	3181.8502	0.0000	1.0000
JC69G	36	-1564.8583	3201.7166	0.0000	1.0000
JC69IG	37	-1564.1282	3202.2565	0.0000	1.0000
JC69I	36	-1565.8077	3203.6154	0.0000	1.0000
K2P	36	-1569.1076	3210.2152	0.0000	1.0000
SYM	40	-1566.8090	3213.6180	0.0000	1.0000
F81	38	-1587.4728	3250.9455	0.0000	1.0000
JC69	35	-1602.3812	3274.7623	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1513.7523	3272.6667	0.7152	0.7152
HKYI	40	-1515.0470	3275.2560	0.1960	0.9112
HKYIG	41	-1512.8110	3276.9131	0.0856	0.9968
K2PG	37	-1528.7325	3284.2399	0.0022	0.9990
K2PI	37	-1529.9107	3286.5963	0.0007	0.9997
K2PIG	38	-1527.7336	3288.3711	0.0003	0.9999
GTRG	44	-1511.0764	3291.8310	0.0000	1.0000
GTRI	44	-1512.6315	3294.9412	0.0000	1.0000
GTRIG	45	-1510.1364	3296.0801	0.0000	1.0000
SYMG	41	-1527.5461	3306.3833	0.0000	1.0000
SYMI	41	-1528.8731	3309.0372	0.0000	1.0000
SYMIG	42	-1526.5224	3310.4649	0.0000	1.0000
F81G	39	-1550.9802	3340.9934	0.0000	1.0000
HKY	39	-1551.6224	3342.2777	0.0000	1.0000
F81I	39	-1551.9251	3342.8831	0.0000	1.0000
F81IG	40	-1550.2713	3345.7045	0.0000	1.0000
JC69G	36	-1564.8583	3350.3624	0.0000	1.0000
JC69I	36	-1565.8077	3352.2612	0.0000	1.0000
JC69IG	37	-1564.1282	3355.0313	0.0000	1.0000
GTR	43	-1547.5134	3358.5761	0.0000	1.0000
K2P	36	-1569.1076	3358.8610	0.0000	1.0000
SYM	40	-1566.8090	3378.7800	0.0000	1.0000
F81	38	-1587.4728	3407.8494	0.0000	1.0000
JC69	35	-1602.3812	3419.2791	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/chr6_3228.phylip = 1 - 459;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3228.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/chr6_3228.phylip = 1 - 459;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3228.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/chr6_3228.phylip = 1 - 459;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3228.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

