
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:03 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_31769.phylip" (19 taxa, 342 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	-1280.8190	2652.5350	0.5886	0.5886
HKYIG	41	-1280.5593	2654.5986	0.2097	0.7983
HKYI	40	-1282.3155	2655.5280	0.1318	0.9301
GTRG	44	-1278.2189	2657.7712	0.0429	0.9730
GTRIG	45	-1277.9367	2659.8600	0.0151	0.9881
GTRI	44	-1279.5035	2660.3404	0.0119	1.0000
HKY	39	-1299.7855	2687.9020	0.0000	1.0000
GTR	43	-1296.7132	2692.1244	0.0000	1.0000
F81G	39	-1303.3783	2695.0877	0.0000	1.0000
K2PG	37	-1306.2402	2695.7304	0.0000	1.0000
F81I	39	-1304.2645	2696.8600	0.0000	1.0000
F81IG	40	-1303.1075	2697.1119	0.0000	1.0000
K2PIG	38	-1305.9373	2697.6568	0.0000	1.0000
K2PI	37	-1307.3417	2697.9335	0.0000	1.0000
SYMG	41	-1304.9722	2703.4244	0.0000	1.0000
SYMIG	42	-1304.6808	2705.4418	0.0000	1.0000
SYMI	41	-1306.1013	2705.6825	0.0000	1.0000
F81	38	-1320.8794	2727.5410	0.0000	1.0000
K2P	36	-1325.3830	2731.5004	0.0000	1.0000
JC69G	36	-1326.4736	2733.6816	0.0000	1.0000
JC69IG	37	-1326.2649	2735.7797	0.0000	1.0000
JC69I	36	-1327.5376	2735.8095	0.0000	1.0000
SYM	40	-1323.8362	2738.5693	0.0000	1.0000
JC69	35	-1343.7109	2765.6570	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1280.8190	2641.6380	0.4682	0.4682
HKYIG	41	-1280.5593	2643.1186	0.2233	0.6915
GTRG	44	-1278.2189	2644.4379	0.1155	0.8069
HKYI	40	-1282.3155	2644.6310	0.1048	0.9117
GTRIG	45	-1277.9367	2645.8735	0.0563	0.9680
GTRI	44	-1279.5035	2647.0071	0.0320	1.0000
HKY	39	-1299.7855	2677.5709	0.0000	1.0000
GTR	43	-1296.7132	2679.4265	0.0000	1.0000
F81G	39	-1303.3783	2684.7565	0.0000	1.0000
F81IG	40	-1303.1075	2686.2149	0.0000	1.0000
K2PG	37	-1306.2402	2686.4804	0.0000	1.0000
F81I	39	-1304.2645	2686.5289	0.0000	1.0000
K2PIG	38	-1305.9373	2687.8746	0.0000	1.0000
K2PI	37	-1307.3417	2688.6835	0.0000	1.0000
SYMG	41	-1304.9722	2691.9444	0.0000	1.0000
SYMIG	42	-1304.6808	2693.3615	0.0000	1.0000
SYMI	41	-1306.1013	2694.2025	0.0000	1.0000
F81	38	-1320.8794	2717.7588	0.0000	1.0000
K2P	36	-1325.3830	2722.7660	0.0000	1.0000
JC69G	36	-1326.4736	2724.9472	0.0000	1.0000
JC69IG	37	-1326.2649	2726.5297	0.0000	1.0000
JC69I	36	-1327.5376	2727.0751	0.0000	1.0000
SYM	40	-1323.8362	2727.6723	0.0000	1.0000
JC69	35	-1343.7109	2757.4217	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1280.8190	2795.0304	0.7727	0.7727
HKYI	40	-1282.3155	2798.0235	0.1730	0.9457
HKYIG	41	-1280.5593	2800.3458	0.0542	0.9999
GTRG	44	-1278.2189	2813.1695	0.0001	1.0000
GTRI	44	-1279.5035	2815.7388	0.0000	1.0000
GTRIG	45	-1277.9367	2818.4400	0.0000	1.0000
HKY	39	-1299.7855	2827.1285	0.0000	1.0000
K2PG	37	-1306.2402	2828.3684	0.0000	1.0000
K2PI	37	-1307.3417	2830.5715	0.0000	1.0000
K2PIG	38	-1305.9373	2833.5974	0.0000	1.0000
F81G	39	-1303.3783	2834.3142	0.0000	1.0000
F81I	39	-1304.2645	2836.0865	0.0000	1.0000
F81IG	40	-1303.1075	2839.6073	0.0000	1.0000
GTR	43	-1296.7132	2844.3233	0.0000	1.0000
SYMG	41	-1304.9722	2849.1717	0.0000	1.0000
SYMI	41	-1306.1013	2851.4298	0.0000	1.0000
SYMIG	42	-1304.6808	2854.4236	0.0000	1.0000
K2P	36	-1325.3830	2860.8192	0.0000	1.0000
JC69G	36	-1326.4736	2863.0004	0.0000	1.0000
F81	38	-1320.8794	2863.4816	0.0000	1.0000
JC69I	36	-1327.5376	2865.1283	0.0000	1.0000
JC69IG	37	-1326.2649	2868.4177	0.0000	1.0000
SYM	40	-1323.8362	2881.0647	0.0000	1.0000
JC69	35	-1343.7109	2891.6401	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_31769.phylip = 1 - 342;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_31769.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_31769.phylip = 1 - 342;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_31769.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_31769.phylip = 1 - 342;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_31769.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

