
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:20 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/chr17_2957.phylip" (19 taxa, 476 characters)

Minimum AIC  model: K2P
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	36	-1276.5427	2631.1537	0.4275	0.4275
K2PG	37	-1275.9824	2632.3849	0.2310	0.6585
K2PI	37	-1276.2151	2632.8502	0.1830	0.8415
K2PIG	38	-1275.9822	2634.7470	0.0709	0.9124
HKY	39	-1275.5652	2636.2863	0.0328	0.9453
HKYG	40	-1275.0001	2637.5405	0.0175	0.9628
HKYI	40	-1275.2348	2638.0098	0.0139	0.9767
SYM	40	-1275.8020	2639.1443	0.0079	0.9846
HKYIG	41	-1274.9999	2639.9353	0.0053	0.9899
SYMG	41	-1275.2106	2640.3568	0.0043	0.9941
SYMI	41	-1275.4449	2640.8252	0.0034	0.9975
SYMIG	42	-1275.2103	2642.7625	0.0013	0.9988
GTR	43	-1274.8527	2644.4647	0.0006	0.9994
GTRG	44	-1274.2458	2645.6795	0.0003	0.9997
GTRI	44	-1274.4837	2646.1554	0.0002	0.9999
GTRIG	45	-1274.2455	2648.1189	0.0001	1.0000
JC69	35	-1307.4460	2690.6192	0.0000	1.0000
JC69G	36	-1307.1221	2692.3126	0.0000	1.0000
JC69I	36	-1307.3090	2692.6862	0.0000	1.0000
JC69IG	37	-1307.1218	2694.6638	0.0000	1.0000
F81	38	-1306.6258	2696.0343	0.0000	1.0000
F81G	39	-1306.3040	2697.7640	0.0000	1.0000
F81I	39	-1306.4919	2698.1397	0.0000	1.0000
F81IG	40	-1306.3037	2700.1476	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2P	36	-1276.5427	2625.0854	0.3553	0.3553
K2PG	37	-1275.9824	2625.9648	0.2289	0.5842
K2PI	37	-1276.2151	2626.4301	0.1814	0.7655
K2PIG	38	-1275.9822	2627.9644	0.0842	0.8498
HKY	39	-1275.5652	2629.1304	0.0470	0.8968
HKYG	40	-1275.0001	2630.0003	0.0304	0.9272
HKYI	40	-1275.2348	2630.4696	0.0241	0.9513
SYM	40	-1275.8020	2631.6041	0.0136	0.9649
HKYIG	41	-1274.9999	2631.9998	0.0112	0.9761
SYMG	41	-1275.2106	2632.4213	0.0091	0.9852
SYMI	41	-1275.4449	2632.8897	0.0072	0.9924
SYMIG	42	-1275.2103	2634.4207	0.0033	0.9957
GTR	43	-1274.8527	2635.7054	0.0018	0.9974
GTRG	44	-1274.2458	2636.4915	0.0012	0.9986
GTRI	44	-1274.4837	2636.9674	0.0009	0.9996
GTRIG	45	-1274.2455	2638.4910	0.0004	1.0000
JC69	35	-1307.4460	2684.8919	0.0000	1.0000
JC69G	36	-1307.1221	2686.2442	0.0000	1.0000
JC69I	36	-1307.3090	2686.6179	0.0000	1.0000
JC69IG	37	-1307.1218	2688.2437	0.0000	1.0000
F81	38	-1306.6258	2689.2517	0.0000	1.0000
F81G	39	-1306.3040	2690.6080	0.0000	1.0000
F81I	39	-1306.4919	2690.9837	0.0000	1.0000
F81IG	40	-1306.3037	2692.6074	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	36	-1276.5427	2775.0404	0.8712	0.8712
K2PG	37	-1275.9824	2780.0853	0.0699	0.9411
K2PI	37	-1276.2151	2780.5506	0.0554	0.9965
K2PIG	38	-1275.9822	2786.2502	0.0032	0.9997
HKY	39	-1275.5652	2791.5817	0.0002	1.0000
HKYG	40	-1275.0001	2796.6170	0.0000	1.0000
HKYI	40	-1275.2348	2797.0863	0.0000	1.0000
SYM	40	-1275.8020	2798.2208	0.0000	1.0000
HKYIG	41	-1274.9999	2802.7820	0.0000	1.0000
SYMG	41	-1275.2106	2803.2034	0.0000	1.0000
SYMI	41	-1275.4449	2803.6719	0.0000	1.0000
SYMIG	42	-1275.2103	2809.3682	0.0000	1.0000
GTR	43	-1274.8527	2814.8184	0.0000	1.0000
GTRG	44	-1274.2458	2819.7699	0.0000	1.0000
GTRI	44	-1274.4837	2820.2458	0.0000	1.0000
GTRIG	45	-1274.2455	2825.9348	0.0000	1.0000
JC69	35	-1307.4460	2830.6815	0.0000	1.0000
JC69G	36	-1307.1221	2836.1993	0.0000	1.0000
JC69I	36	-1307.3090	2836.5729	0.0000	1.0000
JC69IG	37	-1307.1218	2842.3641	0.0000	1.0000
F81	38	-1306.6258	2847.5376	0.0000	1.0000
F81G	39	-1306.3040	2853.0593	0.0000	1.0000
F81I	39	-1306.4919	2853.4350	0.0000	1.0000
F81IG	40	-1306.3037	2859.2241	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2957.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2957.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2957.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2957.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2957.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2957.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

