
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:56 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_7954.phylip" (19 taxa, 385 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	-1319.6628	2728.8604	0.7705	0.7705
HKYI	40	-1321.9859	2733.5066	0.0755	0.8460
HKYIG	41	-1320.8806	2733.8021	0.0651	0.9111
GTRG	44	-1317.3387	2734.3246	0.0501	0.9613
GTRI	44	-1318.2675	2736.1820	0.0198	0.9811
GTRIG	45	-1317.0307	2736.2738	0.0189	1.0000
K2PG	37	-1334.5569	2751.2175	0.0000	1.0000
K2PIG	38	-1335.8409	2756.2484	0.0000	1.0000
K2PI	37	-1337.3512	2756.8061	0.0000	1.0000
SYMG	41	-1334.1051	2760.2510	0.0000	1.0000
SYMIG	42	-1334.9947	2764.5507	0.0000	1.0000
SYMI	41	-1336.4178	2764.8764	0.0000	1.0000
HKY	39	-1357.1037	2801.2510	0.0000	1.0000
F81G	39	-1358.5659	2804.1754	0.0000	1.0000
GTR	43	-1353.5915	2804.2798	0.0000	1.0000
F81IG	40	-1357.6316	2804.7981	0.0000	1.0000
F81I	39	-1358.9019	2804.8473	0.0000	1.0000
JC69G	36	-1370.4591	2820.5734	0.0000	1.0000
JC69IG	37	-1369.9660	2822.0357	0.0000	1.0000
JC69I	36	-1371.4443	2822.5437	0.0000	1.0000
K2P	36	-1371.6720	2822.9991	0.0000	1.0000
SYM	40	-1370.9331	2831.4011	0.0000	1.0000
F81	38	-1392.0191	2868.6046	0.0000	1.0000
JC69	35	-1404.7948	2886.8103	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1319.6628	2719.3256	0.6404	0.6404
GTRG	44	-1317.3387	2722.6775	0.1198	0.7602
HKYIG	41	-1320.8806	2723.7613	0.0697	0.8299
HKYI	40	-1321.9859	2723.9717	0.0627	0.8927
GTRIG	45	-1317.0307	2724.0614	0.0600	0.9527
GTRI	44	-1318.2675	2724.5350	0.0473	1.0000
K2PG	37	-1334.5569	2743.1138	0.0000	1.0000
K2PIG	38	-1335.8409	2747.6819	0.0000	1.0000
K2PI	37	-1337.3512	2748.7024	0.0000	1.0000
SYMG	41	-1334.1051	2750.2102	0.0000	1.0000
SYMIG	42	-1334.9947	2753.9893	0.0000	1.0000
SYMI	41	-1336.4178	2754.8356	0.0000	1.0000
HKY	39	-1357.1037	2792.2075	0.0000	1.0000
GTR	43	-1353.5915	2793.1830	0.0000	1.0000
F81G	39	-1358.5659	2795.1319	0.0000	1.0000
F81IG	40	-1357.6316	2795.2632	0.0000	1.0000
F81I	39	-1358.9019	2795.8038	0.0000	1.0000
JC69G	36	-1370.4591	2812.9183	0.0000	1.0000
JC69IG	37	-1369.9660	2813.9319	0.0000	1.0000
JC69I	36	-1371.4443	2814.8886	0.0000	1.0000
K2P	36	-1371.6720	2815.3440	0.0000	1.0000
SYM	40	-1370.9331	2821.8662	0.0000	1.0000
F81	38	-1392.0191	2860.0382	0.0000	1.0000
JC69	35	-1404.7948	2879.5897	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1319.6628	2877.4553	0.8961	0.8961
HKYI	40	-1321.9859	2882.1015	0.0878	0.9839
HKYIG	41	-1320.8806	2885.8442	0.0135	0.9974
K2PG	37	-1334.5569	2889.3838	0.0023	0.9997
K2PI	37	-1337.3512	2894.9724	0.0001	0.9999
GTRG	44	-1317.3387	2896.6202	0.0001	0.9999
K2PIG	38	-1335.8409	2897.9051	0.0000	1.0000
GTRI	44	-1318.2675	2898.4777	0.0000	1.0000
GTRIG	45	-1317.0307	2901.9574	0.0000	1.0000
SYMG	41	-1334.1051	2912.2932	0.0000	1.0000
SYMI	41	-1336.4178	2916.9186	0.0000	1.0000
SYMIG	42	-1334.9947	2920.0255	0.0000	1.0000
HKY	39	-1357.1037	2946.3840	0.0000	1.0000
F81G	39	-1358.5659	2949.3084	0.0000	1.0000
F81I	39	-1358.9019	2949.9803	0.0000	1.0000
F81IG	40	-1357.6316	2953.3929	0.0000	1.0000
JC69G	36	-1370.4591	2955.2350	0.0000	1.0000
JC69I	36	-1371.4443	2957.2053	0.0000	1.0000
K2P	36	-1371.6720	2957.6607	0.0000	1.0000
JC69IG	37	-1369.9660	2960.2019	0.0000	1.0000
GTR	43	-1353.5915	2963.1725	0.0000	1.0000
SYM	40	-1370.9331	2979.9960	0.0000	1.0000
F81	38	-1392.0191	3010.2614	0.0000	1.0000
JC69	35	-1404.7948	3017.9532	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/chr2_7954.phylip = 1 - 385;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7954.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/chr2_7954.phylip = 1 - 385;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7954.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/chr2_7954.phylip = 1 - 385;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7954.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

