
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:37 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/chr5_11797.phylip" (19 taxa, 644 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	-1401.5977	2888.6349	0.4614	0.4614
HKYIG	41	-1400.9716	2889.6641	0.2758	0.7372
HKYI	40	-1402.1623	2889.7641	0.2623	0.9995
HKY	39	-1410.5776	2904.3208	0.0002	0.9997
GTRG	44	-1405.3249	2905.2609	0.0001	0.9998
GTRIG	45	-1404.2822	2905.4874	0.0001	0.9999
GTRI	44	-1406.0193	2906.6495	0.0001	1.0000
GTR	43	-1415.7980	2923.9026	0.0000	1.0000
F81G	39	-1423.6617	2930.4890	0.0000	1.0000
F81IG	40	-1422.6020	2930.6436	0.0000	1.0000
F81I	39	-1424.2497	2931.6649	0.0000	1.0000
F81	38	-1434.3453	2949.5897	0.0000	1.0000
K2PG	37	-1438.3714	2955.3831	0.0000	1.0000
K2PI	37	-1438.7940	2956.2282	0.0000	1.0000
K2PIG	38	-1437.6801	2956.2594	0.0000	1.0000
SYMG	41	-1437.1177	2961.9563	0.0000	1.0000
SYMI	41	-1437.5063	2962.7336	0.0000	1.0000
SYMIG	42	-1436.4285	2962.8669	0.0000	1.0000
K2P	36	-1448.0321	2972.4530	0.0000	1.0000
SYM	40	-1446.9037	2979.2469	0.0000	1.0000
JC69G	36	-1455.4450	2987.2788	0.0000	1.0000
JC69I	36	-1455.8923	2988.1734	0.0000	1.0000
JC69IG	37	-1454.7691	2988.1785	0.0000	1.0000
JC69	35	-1464.9637	3004.0721	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1401.5977	2883.1954	0.4428	0.4428
HKYIG	41	-1400.9716	2883.9432	0.3047	0.7476
HKYI	40	-1402.1623	2884.3247	0.2518	0.9994
GTRIG	45	-1404.2822	2898.5643	0.0002	0.9996
GTRG	44	-1405.3249	2898.6499	0.0002	0.9998
HKY	39	-1410.5776	2899.1552	0.0002	0.9999
GTRI	44	-1406.0193	2900.0385	0.0001	1.0000
GTR	43	-1415.7980	2917.5959	0.0000	1.0000
F81IG	40	-1422.6020	2925.2041	0.0000	1.0000
F81G	39	-1423.6617	2925.3234	0.0000	1.0000
F81I	39	-1424.2497	2926.4993	0.0000	1.0000
F81	38	-1434.3453	2944.6906	0.0000	1.0000
K2PG	37	-1438.3714	2950.7428	0.0000	1.0000
K2PIG	38	-1437.6801	2951.3602	0.0000	1.0000
K2PI	37	-1438.7940	2951.5879	0.0000	1.0000
SYMG	41	-1437.1177	2956.2354	0.0000	1.0000
SYMIG	42	-1436.4285	2956.8569	0.0000	1.0000
SYMI	41	-1437.5063	2957.0126	0.0000	1.0000
K2P	36	-1448.0321	2968.0642	0.0000	1.0000
SYM	40	-1446.9037	2973.8075	0.0000	1.0000
JC69G	36	-1455.4450	2982.8900	0.0000	1.0000
JC69IG	37	-1454.7691	2983.5383	0.0000	1.0000
JC69I	36	-1455.8923	2983.7846	0.0000	1.0000
JC69	35	-1464.9637	2999.9273	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1401.5977	3061.9034	0.6077	0.6077
HKYI	40	-1402.1623	3063.0326	0.3455	0.9533
HKYIG	41	-1400.9716	3067.1188	0.0448	0.9981
HKY	39	-1410.5776	3073.3955	0.0019	1.0000
GTRG	44	-1405.3249	3095.2286	0.0000	1.0000
GTRI	44	-1406.0193	3096.6172	0.0000	1.0000
F81G	39	-1423.6617	3099.5637	0.0000	1.0000
GTRIG	45	-1404.2822	3099.6107	0.0000	1.0000
F81I	39	-1424.2497	3100.7396	0.0000	1.0000
F81IG	40	-1422.6020	3103.9120	0.0000	1.0000
GTR	43	-1415.7980	3109.7070	0.0000	1.0000
F81	38	-1434.3453	3114.4631	0.0000	1.0000
K2PG	37	-1438.3714	3116.0477	0.0000	1.0000
K2PI	37	-1438.7940	3116.8928	0.0000	1.0000
K2PIG	38	-1437.6801	3121.1328	0.0000	1.0000
K2P	36	-1448.0321	3128.9013	0.0000	1.0000
SYMG	41	-1437.1177	3139.4110	0.0000	1.0000
SYMI	41	-1437.5063	3140.1883	0.0000	1.0000
JC69G	36	-1455.4450	3143.7272	0.0000	1.0000
SYMIG	42	-1436.4285	3144.5002	0.0000	1.0000
JC69I	36	-1455.8923	3144.6218	0.0000	1.0000
JC69IG	37	-1454.7691	3148.8431	0.0000	1.0000
SYM	40	-1446.9037	3152.5154	0.0000	1.0000
JC69	35	-1464.9637	3156.2968	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/chr5_11797.phylip = 1 - 644;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11797.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/chr5_11797.phylip = 1 - 644;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11797.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/chr5_11797.phylip = 1 - 644;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11797.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

