
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:20 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/chr1_32397.phylip" (19 taxa, 411 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	-1240.7208	2570.3065	0.4508	0.4508
HKYI	40	-1240.8531	2570.5711	0.3950	0.8458
HKYIG	41	-1240.7215	2572.7763	0.1311	0.9769
GTRG	44	-1239.4189	2577.6576	0.0114	0.9883
GTRI	44	-1239.7277	2578.2751	0.0084	0.9967
GTRIG	45	-1239.4102	2580.1629	0.0033	1.0000
K2PI	37	-1255.7975	2593.1338	0.0000	1.0000
K2PG	37	-1255.8471	2593.2331	0.0000	1.0000
K2PIG	38	-1255.8469	2595.6615	0.0000	1.0000
SYMG	41	-1253.9349	2599.2032	0.0000	1.0000
SYMI	41	-1254.0497	2599.4328	0.0000	1.0000
SYMIG	42	-1253.9352	2601.6856	0.0000	1.0000
F81I	39	-1274.9938	2636.3974	0.0000	1.0000
F81G	39	-1275.1230	2636.6556	0.0000	1.0000
HKY	39	-1275.4106	2637.2310	0.0000	1.0000
F81IG	40	-1275.1228	2639.1104	0.0000	1.0000
GTR	43	-1273.9678	2644.2462	0.0000	1.0000
JC69I	36	-1287.6923	2654.5077	0.0000	1.0000
JC69G	36	-1287.8722	2654.8675	0.0000	1.0000
JC69IG	37	-1287.8720	2657.2829	0.0000	1.0000
K2P	36	-1290.7388	2660.6006	0.0000	1.0000
SYM	40	-1289.0182	2666.9013	0.0000	1.0000
F81	38	-1308.0452	2700.0581	0.0000	1.0000
JC69	35	-1320.6837	2718.0873	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1240.7208	2561.4417	0.4192	0.4192
HKYI	40	-1240.8531	2561.7062	0.3673	0.7865
HKYIG	41	-1240.7215	2563.4430	0.1541	0.9406
GTRG	44	-1239.4189	2566.8379	0.0282	0.9688
GTRI	44	-1239.7277	2567.4554	0.0207	0.9895
GTRIG	45	-1239.4102	2568.8204	0.0105	1.0000
K2PI	37	-1255.7975	2585.5950	0.0000	1.0000
K2PG	37	-1255.8471	2585.6942	0.0000	1.0000
K2PIG	38	-1255.8469	2587.6938	0.0000	1.0000
SYMG	41	-1253.9349	2589.8699	0.0000	1.0000
SYMI	41	-1254.0497	2590.0995	0.0000	1.0000
SYMIG	42	-1253.9352	2591.8704	0.0000	1.0000
F81I	39	-1274.9938	2627.9877	0.0000	1.0000
F81G	39	-1275.1230	2628.2459	0.0000	1.0000
HKY	39	-1275.4106	2628.8213	0.0000	1.0000
F81IG	40	-1275.1228	2630.2456	0.0000	1.0000
GTR	43	-1273.9678	2633.9356	0.0000	1.0000
JC69I	36	-1287.6923	2647.3847	0.0000	1.0000
JC69G	36	-1287.8722	2647.7445	0.0000	1.0000
JC69IG	37	-1287.8720	2649.7440	0.0000	1.0000
K2P	36	-1290.7388	2653.4776	0.0000	1.0000
SYM	40	-1289.0182	2658.0364	0.0000	1.0000
F81	38	-1308.0452	2692.0904	0.0000	1.0000
JC69	35	-1320.6837	2711.3673	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1240.7208	2722.1854	0.5181	0.5181
HKYI	40	-1240.8531	2722.4499	0.4539	0.9720
HKYIG	41	-1240.7215	2728.2053	0.0255	0.9975
K2PI	37	-1255.7975	2734.2829	0.0012	0.9988
K2PG	37	-1255.8471	2734.3822	0.0012	0.9999
K2PIG	38	-1255.8469	2740.4003	0.0001	1.0000
GTRG	44	-1239.4189	2743.6560	0.0000	1.0000
GTRI	44	-1239.7277	2744.2735	0.0000	1.0000
GTRIG	45	-1239.4102	2749.6571	0.0000	1.0000
SYMG	41	-1253.9349	2754.6322	0.0000	1.0000
SYMI	41	-1254.0497	2754.8618	0.0000	1.0000
SYMIG	42	-1253.9352	2760.6513	0.0000	1.0000
F81I	39	-1274.9938	2784.7128	0.0000	1.0000
F81G	39	-1275.1230	2784.9711	0.0000	1.0000
HKY	39	-1275.4106	2785.5464	0.0000	1.0000
F81IG	40	-1275.1228	2790.9893	0.0000	1.0000
JC69I	36	-1287.6923	2792.0540	0.0000	1.0000
JC69G	36	-1287.8722	2792.4138	0.0000	1.0000
K2P	36	-1290.7388	2798.1470	0.0000	1.0000
JC69IG	37	-1287.8720	2798.4319	0.0000	1.0000
GTR	43	-1273.9678	2806.7351	0.0000	1.0000
SYM	40	-1289.0182	2818.7801	0.0000	1.0000
F81	38	-1308.0452	2844.7969	0.0000	1.0000
JC69	35	-1320.6837	2852.0181	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/chr1_32397.phylip = 1 - 411;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32397.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/chr1_32397.phylip = 1 - 411;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32397.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/chr1_32397.phylip = 1 - 411;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32397.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

