
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:38 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_12015.phylip" (19 taxa, 429 characters)

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	41	-1206.4061	2503.7114	0.3947	0.3947
K2PIG	38	-1210.8179	2505.2358	0.1842	0.5789
HKYG	40	-1208.8984	2506.2504	0.1109	0.6899
HKYI	40	-1209.2590	2506.9717	0.0773	0.7672
GTRIG	45	-1203.1013	2507.0121	0.0758	0.8430
K2PG	37	-1213.2977	2507.7872	0.0514	0.8944
K2PI	37	-1213.6187	2508.4292	0.0373	0.9317
GTRG	44	-1205.1607	2508.6339	0.0337	0.9654
GTRI	44	-1205.8566	2510.0257	0.0168	0.9822
SYMIG	42	-1208.7711	2510.8996	0.0108	0.9931
SYMG	41	-1210.9026	2512.7044	0.0044	0.9975
SYMI	41	-1211.4492	2513.7977	0.0025	1.0000
F81IG	40	-1228.0765	2544.6065	0.0000	1.0000
JC69IG	37	-1231.8052	2544.8022	0.0000	1.0000
F81G	39	-1230.1325	2546.2856	0.0000	1.0000
JC69G	36	-1233.8432	2546.4823	0.0000	1.0000
F81I	39	-1230.5319	2547.0844	0.0000	1.0000
JC69I	36	-1234.2028	2547.2015	0.0000	1.0000
GTR	43	-1236.4308	2568.6903	0.0000	1.0000
HKY	39	-1242.7913	2571.6032	0.0000	1.0000
SYM	40	-1242.2679	2572.9895	0.0000	1.0000
K2P	36	-1247.1163	2573.0285	0.0000	1.0000
F81	38	-1262.9710	2609.5419	0.0000	1.0000
JC69	35	-1266.9003	2610.2128	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-1206.4061	2494.8122	0.3950	0.3950
GTRIG	45	-1203.1013	2496.2027	0.1971	0.5921
K2PIG	38	-1210.8179	2497.6358	0.0963	0.6884
HKYG	40	-1208.8984	2497.7968	0.0888	0.7772
GTRG	44	-1205.1607	2498.3214	0.0683	0.8456
HKYI	40	-1209.2590	2498.5181	0.0619	0.9075
GTRI	44	-1205.8566	2499.7132	0.0341	0.9416
K2PG	37	-1213.2977	2500.5954	0.0219	0.9635
K2PI	37	-1213.6187	2501.2374	0.0159	0.9794
SYMIG	42	-1208.7711	2501.5421	0.0137	0.9930
SYMG	41	-1210.9026	2503.8052	0.0044	0.9975
SYMI	41	-1211.4492	2504.8984	0.0025	1.0000
F81IG	40	-1228.0765	2536.1529	0.0000	1.0000
JC69IG	37	-1231.8052	2537.6103	0.0000	1.0000
F81G	39	-1230.1325	2538.2651	0.0000	1.0000
F81I	39	-1230.5319	2539.0638	0.0000	1.0000
JC69G	36	-1233.8432	2539.6864	0.0000	1.0000
JC69I	36	-1234.2028	2540.4056	0.0000	1.0000
GTR	43	-1236.4308	2558.8617	0.0000	1.0000
HKY	39	-1242.7913	2563.5826	0.0000	1.0000
SYM	40	-1242.2679	2564.5359	0.0000	1.0000
K2P	36	-1247.1163	2566.2326	0.0000	1.0000
F81	38	-1262.9710	2601.9419	0.0000	1.0000
JC69	35	-1266.9003	2603.8006	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1213.2977	2650.8693	0.4305	0.4305
K2PI	37	-1213.6187	2651.5113	0.3123	0.7428
K2PIG	38	-1210.8179	2651.9712	0.2481	0.9910
HKYG	40	-1208.8984	2660.2551	0.0039	0.9949
HKYI	40	-1209.2590	2660.9764	0.0027	0.9976
HKYIG	41	-1206.4061	2661.3319	0.0023	0.9999
SYMG	41	-1210.9026	2670.3249	0.0000	1.0000
SYMI	41	-1211.4492	2671.4182	0.0000	1.0000
SYMIG	42	-1208.7711	2672.1233	0.0000	1.0000
GTRG	44	-1205.1607	2677.0255	0.0000	1.0000
GTRI	44	-1205.8566	2678.4173	0.0000	1.0000
GTRIG	45	-1203.1013	2678.9683	0.0000	1.0000
JC69G	36	-1233.8432	2685.8988	0.0000	1.0000
JC69I	36	-1234.2028	2686.6181	0.0000	1.0000
JC69IG	37	-1231.8052	2687.8842	0.0000	1.0000
F81G	39	-1230.1325	2696.6619	0.0000	1.0000
F81I	39	-1230.5319	2697.4607	0.0000	1.0000
F81IG	40	-1228.0765	2698.6112	0.0000	1.0000
K2P	36	-1247.1163	2712.4450	0.0000	1.0000
HKY	39	-1242.7913	2721.9794	0.0000	1.0000
SYM	40	-1242.2679	2726.9942	0.0000	1.0000
GTR	43	-1236.4308	2733.5043	0.0000	1.0000
JC69	35	-1266.9003	2745.9516	0.0000	1.0000
F81	38	-1262.9710	2756.2773	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12015.phylip = 1 - 429;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12015.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12015.phylip = 1 - 429;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12015.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.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/chr5_12015.phylip = 1 - 429;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12015.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

