
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:23 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/chrZ_7833.phylip" (19 taxa, 464 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-1382.3398	2852.4336	0.4023	0.4023
HKYG	40	-1382.5402	2852.8346	0.3292	0.7316
HKYIG	41	-1382.5402	2855.2415	0.0988	0.8304
GTRI	44	-1378.9214	2855.2938	0.0963	0.9267
GTRG	44	-1379.4510	2856.3531	0.0567	0.9834
GTRIG	45	-1379.4511	2858.8064	0.0166	1.0000
K2PI	37	-1400.1578	2880.9166	0.0000	1.0000
K2PG	37	-1400.2390	2881.0790	0.0000	1.0000
K2PIG	38	-1400.2390	2883.4521	0.0000	1.0000
SYMI	41	-1399.5176	2889.1964	0.0000	1.0000
SYMG	41	-1399.7145	2889.5901	0.0000	1.0000
SYMIG	42	-1399.7144	2892.0084	0.0000	1.0000
HKY	39	-1408.8974	2903.1534	0.0000	1.0000
GTR	43	-1404.1999	2903.4093	0.0000	1.0000
F81I	39	-1413.2640	2911.8866	0.0000	1.0000
F81G	39	-1413.4403	2912.2391	0.0000	1.0000
F81IG	40	-1413.4402	2914.6345	0.0000	1.0000
K2P	36	-1429.6476	2937.5340	0.0000	1.0000
JC69I	36	-1430.4425	2939.1238	0.0000	1.0000
JC69G	36	-1430.5192	2939.2772	0.0000	1.0000
JC69IG	37	-1430.5190	2941.6389	0.0000	1.0000
SYM	40	-1428.5241	2944.8023	0.0000	1.0000
F81	38	-1440.4617	2963.8976	0.0000	1.0000
JC69	35	-1459.1789	2994.2458	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1382.3398	2844.6795	0.3196	0.3196
HKYG	40	-1382.5402	2845.0805	0.2616	0.5812
GTRI	44	-1378.9214	2845.8427	0.1787	0.7599
GTRG	44	-1379.4510	2846.9020	0.1052	0.8651
HKYIG	41	-1382.5402	2847.0804	0.0962	0.9613
GTRIG	45	-1379.4511	2848.9021	0.0387	1.0000
K2PI	37	-1400.1578	2874.3157	0.0000	1.0000
K2PG	37	-1400.2390	2874.4781	0.0000	1.0000
K2PIG	38	-1400.2390	2876.4780	0.0000	1.0000
SYMI	41	-1399.5176	2881.0352	0.0000	1.0000
SYMG	41	-1399.7145	2881.4289	0.0000	1.0000
SYMIG	42	-1399.7144	2883.4288	0.0000	1.0000
GTR	43	-1404.1999	2894.3997	0.0000	1.0000
HKY	39	-1408.8974	2895.7949	0.0000	1.0000
F81I	39	-1413.2640	2904.5281	0.0000	1.0000
F81G	39	-1413.4403	2904.8806	0.0000	1.0000
F81IG	40	-1413.4402	2906.8804	0.0000	1.0000
K2P	36	-1429.6476	2931.2951	0.0000	1.0000
JC69I	36	-1430.4425	2932.8850	0.0000	1.0000
JC69G	36	-1430.5192	2933.0383	0.0000	1.0000
JC69IG	37	-1430.5190	2935.0379	0.0000	1.0000
SYM	40	-1428.5241	2937.0482	0.0000	1.0000
F81	38	-1440.4617	2956.9234	0.0000	1.0000
JC69	35	-1459.1789	2988.3579	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1382.3398	3010.2749	0.5385	0.5385
HKYG	40	-1382.5402	3010.6758	0.4407	0.9792
HKYIG	41	-1382.5402	3016.8156	0.0205	0.9997
K2PI	37	-1400.1578	3027.4914	0.0001	0.9998
K2PG	37	-1400.2390	3027.6538	0.0001	0.9999
GTRI	44	-1378.9214	3027.9977	0.0001	0.9999
GTRG	44	-1379.4510	3029.0569	0.0000	1.0000
K2PIG	38	-1400.2390	3033.7936	0.0000	1.0000
GTRIG	45	-1379.4511	3035.1969	0.0000	1.0000
SYMI	41	-1399.5176	3050.7705	0.0000	1.0000
SYMG	41	-1399.7145	3051.1642	0.0000	1.0000
HKY	39	-1408.8974	3057.2504	0.0000	1.0000
SYMIG	42	-1399.7144	3057.3040	0.0000	1.0000
F81I	39	-1413.2640	3065.9836	0.0000	1.0000
F81G	39	-1413.4403	3066.3361	0.0000	1.0000
GTR	43	-1404.1999	3072.4148	0.0000	1.0000
F81IG	40	-1413.4402	3072.4757	0.0000	1.0000
K2P	36	-1429.6476	3080.3310	0.0000	1.0000
JC69I	36	-1430.4425	3081.9208	0.0000	1.0000
JC69G	36	-1430.5192	3082.0742	0.0000	1.0000
JC69IG	37	-1430.5190	3088.2136	0.0000	1.0000
SYM	40	-1428.5241	3102.6436	0.0000	1.0000
F81	38	-1440.4617	3114.2391	0.0000	1.0000
JC69	35	-1459.1789	3133.2539	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7833.phylip = 1 - 464;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7833.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7833.phylip = 1 - 464;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7833.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7833.phylip = 1 - 464;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7833.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

