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

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


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	-1291.7377	2679.7082	0.5302	0.5302
GTRIG	45	-1291.6241	2681.8733	0.1796	0.7098
GTRI	44	-1293.1980	2682.6288	0.1231	0.8329
HKYG	40	-1298.4427	2683.6483	0.0739	0.9069
HKYIG	41	-1298.2562	2685.6282	0.0275	0.9344
SYMG	41	-1298.3148	2685.7452	0.0259	0.9603
HKYI	40	-1299.5936	2685.9500	0.0234	0.9837
SYMIG	42	-1298.2113	2687.9008	0.0088	0.9925
SYMI	41	-1299.7319	2688.5795	0.0063	0.9988
K2PG	37	-1306.5835	2692.9293	0.0007	0.9995
K2PIG	38	-1306.4001	2694.8865	0.0003	0.9998
K2PI	37	-1307.6619	2695.0861	0.0002	1.0000
GTR	43	-1305.4502	2704.7511	0.0000	1.0000
SYM	40	-1311.4182	2709.5994	0.0000	1.0000
F81G	39	-1315.1006	2714.6209	0.0000	1.0000
HKY	39	-1315.2398	2714.8993	0.0000	1.0000
F81IG	40	-1314.9242	2716.6112	0.0000	1.0000
F81I	39	-1316.1886	2716.7969	0.0000	1.0000
K2P	36	-1323.5609	2724.5697	0.0000	1.0000
JC69G	36	-1323.6453	2724.7384	0.0000	1.0000
JC69IG	37	-1323.4996	2726.7615	0.0000	1.0000
JC69I	36	-1324.7149	2726.8777	0.0000	1.0000
F81	38	-1332.2021	2746.4904	0.0000	1.0000
JC69	35	-1339.9334	2755.0097	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1291.7377	2671.4753	0.5529	0.5529
GTRIG	45	-1291.6241	2673.2483	0.2278	0.7807
GTRI	44	-1293.1980	2674.3960	0.1284	0.9091
HKYG	40	-1298.4427	2676.8855	0.0370	0.9460
HKYIG	41	-1298.2562	2678.5125	0.0164	0.9624
SYMG	41	-1298.3148	2678.6295	0.0155	0.9779
HKYI	40	-1299.5936	2679.1872	0.0117	0.9896
SYMIG	42	-1298.2113	2680.4225	0.0063	0.9959
SYMI	41	-1299.7319	2681.4638	0.0037	0.9996
K2PG	37	-1306.5835	2687.1670	0.0002	0.9998
K2PIG	38	-1306.4001	2688.8002	0.0001	0.9999
K2PI	37	-1307.6619	2689.3238	0.0001	1.0000
GTR	43	-1305.4502	2696.9004	0.0000	1.0000
SYM	40	-1311.4182	2702.8365	0.0000	1.0000
F81G	39	-1315.1006	2708.2011	0.0000	1.0000
HKY	39	-1315.2398	2708.4795	0.0000	1.0000
F81IG	40	-1314.9242	2709.8483	0.0000	1.0000
F81I	39	-1316.1886	2710.3771	0.0000	1.0000
K2P	36	-1323.5609	2719.1218	0.0000	1.0000
JC69G	36	-1323.6453	2719.2905	0.0000	1.0000
JC69IG	37	-1323.4996	2720.9992	0.0000	1.0000
JC69I	36	-1324.7149	2721.4298	0.0000	1.0000
F81	38	-1332.2021	2740.4042	0.0000	1.0000
JC69	35	-1339.9334	2749.8668	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1306.5835	2844.9831	0.5552	0.5552
K2PI	37	-1307.6619	2847.1399	0.1888	0.7440
HKYG	40	-1298.4427	2847.4975	0.1579	0.9020
HKYI	40	-1299.5936	2849.7992	0.0500	0.9519
K2PIG	38	-1306.4001	2850.8817	0.0291	0.9810
HKYIG	41	-1298.2562	2853.3898	0.0083	0.9893
SYMG	41	-1298.3148	2853.5069	0.0078	0.9971
SYMI	41	-1299.7319	2856.3412	0.0019	0.9990
GTRG	44	-1291.7377	2859.1486	0.0005	0.9995
SYMIG	42	-1298.2113	2859.5652	0.0004	0.9999
GTRI	44	-1293.1980	2862.0692	0.0001	1.0000
GTRIG	45	-1291.6241	2865.1868	0.0000	1.0000
K2P	36	-1323.5609	2872.6726	0.0000	1.0000
JC69G	36	-1323.6453	2872.8414	0.0000	1.0000
SYM	40	-1311.4182	2873.4485	0.0000	1.0000
F81G	39	-1315.1006	2874.5479	0.0000	1.0000
HKY	39	-1315.2398	2874.8263	0.0000	1.0000
JC69I	36	-1324.7149	2874.9807	0.0000	1.0000
F81I	39	-1316.1886	2876.7239	0.0000	1.0000
JC69IG	37	-1323.4996	2878.8153	0.0000	1.0000
GTR	43	-1305.4502	2880.3084	0.0000	1.0000
F81IG	40	-1314.9242	2880.4603	0.0000	1.0000
JC69	35	-1339.9334	2899.1524	0.0000	1.0000
F81	38	-1332.2021	2902.4856	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/chr2_18144.phylip = 1 - 526;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18144.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/chr2_18144.phylip = 1 - 526;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18144.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 (K2PG)]

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


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

