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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	44	-1562.6007	3220.9967	0.5551	0.5551
GTRI	44	-1563.4419	3222.6791	0.2394	0.7945
GTRIG	45	-1562.4093	3222.9844	0.2055	1.0000
HKYG	40	-1577.2953	3240.9969	0.0000	1.0000
HKYI	40	-1577.8466	3242.0995	0.0000	1.0000
HKYIG	41	-1577.1218	3242.9832	0.0000	1.0000
SYMG	41	-1583.7267	3256.1931	0.0000	1.0000
SYMI	41	-1584.1457	3257.0311	0.0000	1.0000
SYMIG	42	-1583.6560	3258.3943	0.0000	1.0000
K2PG	37	-1591.6074	3262.6750	0.0000	1.0000
K2PI	37	-1592.1282	3263.7167	0.0000	1.0000
K2PIG	38	-1591.4318	3264.6302	0.0000	1.0000
GTR	43	-1591.4060	3276.2461	0.0000	1.0000
HKY	39	-1605.9879	3296.0576	0.0000	1.0000
F81G	39	-1607.3210	3298.7239	0.0000	1.0000
F81I	39	-1607.7302	3299.5423	0.0000	1.0000
F81IG	40	-1607.2677	3300.9416	0.0000	1.0000
SYM	40	-1612.4893	3311.3849	0.0000	1.0000
JC69G	36	-1618.9449	3315.0525	0.0000	1.0000
JC69I	36	-1619.3700	3315.9028	0.0000	1.0000
JC69IG	37	-1618.8702	3317.2006	0.0000	1.0000
K2P	36	-1621.4061	3319.9750	0.0000	1.0000
F81	38	-1635.3561	3352.4787	0.0000	1.0000
JC69	35	-1647.8999	3370.6741	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1562.6007	3213.2014	0.5329	0.5329
GTRIG	45	-1562.4093	3214.8187	0.2374	0.7702
GTRI	44	-1563.4419	3214.8839	0.2298	1.0000
HKYG	40	-1577.2953	3234.5907	0.0000	1.0000
HKYI	40	-1577.8466	3235.6933	0.0000	1.0000
HKYIG	41	-1577.1218	3236.2435	0.0000	1.0000
SYMG	41	-1583.7267	3249.4534	0.0000	1.0000
SYMI	41	-1584.1457	3250.2914	0.0000	1.0000
SYMIG	42	-1583.6560	3251.3119	0.0000	1.0000
K2PG	37	-1591.6074	3257.2148	0.0000	1.0000
K2PI	37	-1592.1282	3258.2565	0.0000	1.0000
K2PIG	38	-1591.4318	3258.8637	0.0000	1.0000
GTR	43	-1591.4060	3268.8119	0.0000	1.0000
HKY	39	-1605.9879	3289.9758	0.0000	1.0000
F81G	39	-1607.3210	3292.6420	0.0000	1.0000
F81I	39	-1607.7302	3293.4604	0.0000	1.0000
F81IG	40	-1607.2677	3294.5353	0.0000	1.0000
SYM	40	-1612.4893	3304.9787	0.0000	1.0000
JC69G	36	-1618.9449	3309.8897	0.0000	1.0000
JC69I	36	-1619.3700	3310.7400	0.0000	1.0000
JC69IG	37	-1618.8702	3311.7404	0.0000	1.0000
K2P	36	-1621.4061	3314.8122	0.0000	1.0000
F81	38	-1635.3561	3346.7122	0.0000	1.0000
JC69	35	-1647.8999	3365.7998	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1562.6007	3403.0771	0.5914	0.5914
GTRI	44	-1563.4419	3404.7596	0.2550	0.8464
HKYG	40	-1577.2953	3407.2050	0.0751	0.9215
HKYI	40	-1577.8466	3408.3076	0.0433	0.9648
GTRIG	45	-1562.4093	3409.0098	0.0305	0.9952
HKYIG	41	-1577.1218	3413.1732	0.0038	0.9990
K2PG	37	-1591.6074	3416.8830	0.0006	0.9996
K2PI	37	-1592.1282	3417.9247	0.0004	1.0000
K2PIG	38	-1591.4318	3422.8473	0.0000	1.0000
SYMG	41	-1583.7267	3426.3831	0.0000	1.0000
SYMI	41	-1584.1457	3427.2211	0.0000	1.0000
SYMIG	42	-1583.6560	3432.5570	0.0000	1.0000
GTR	43	-1591.4060	3454.3723	0.0000	1.0000
HKY	39	-1605.9879	3458.2747	0.0000	1.0000
F81G	39	-1607.3210	3460.9409	0.0000	1.0000
F81I	39	-1607.7302	3461.7594	0.0000	1.0000
JC69G	36	-1618.9449	3465.2426	0.0000	1.0000
JC69I	36	-1619.3700	3466.0929	0.0000	1.0000
F81IG	40	-1607.2677	3467.1497	0.0000	1.0000
K2P	36	-1621.4061	3470.1651	0.0000	1.0000
JC69IG	37	-1618.8702	3471.4086	0.0000	1.0000
SYM	40	-1612.4893	3477.5930	0.0000	1.0000
F81	38	-1635.3561	3510.6958	0.0000	1.0000
JC69	35	-1647.8999	3516.8373	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2142.phylip = 1 - 553;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2142.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2142.phylip = 1 - 553;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2142.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2142.phylip = 1 - 553;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr15_2142.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

