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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	40	-1180.2079	2448.7831	0.4083	0.4083
HKYI	40	-1180.6510	2449.6694	0.2621	0.6705
HKYIG	41	-1179.8059	2450.4201	0.1801	0.8506
K2PG	37	-1186.1388	2453.3966	0.0407	0.8912
GTRG	44	-1177.8885	2453.9832	0.0303	0.9216
K2PI	37	-1186.5653	2454.2495	0.0265	0.9481
K2PIG	38	-1185.7694	2455.0615	0.0177	0.9658
GTRI	44	-1178.5194	2455.2450	0.0161	0.9819
GTRIG	45	-1177.4153	2455.5283	0.0140	0.9959
SYMG	41	-1184.2895	2459.3871	0.0020	0.9980
SYMI	41	-1184.9219	2460.6520	0.0011	0.9991
SYMIG	42	-1183.8253	2460.9122	0.0009	1.0000
HKY	39	-1196.7645	2479.4680	0.0000	1.0000
GTR	43	-1194.2091	2484.1457	0.0000	1.0000
K2P	36	-1202.7623	2484.2518	0.0000	1.0000
SYM	40	-1201.0161	2490.3996	0.0000	1.0000
F81G	39	-1210.4327	2506.8044	0.0000	1.0000
F81I	39	-1210.4981	2506.9352	0.0000	1.0000
F81IG	40	-1210.4259	2509.2192	0.0000	1.0000
JC69G	36	-1219.0578	2516.8429	0.0000	1.0000
JC69I	36	-1219.1148	2516.9568	0.0000	1.0000
JC69IG	37	-1219.0494	2519.2179	0.0000	1.0000
F81	38	-1225.9164	2535.3557	0.0000	1.0000
JC69	35	-1234.2552	2544.8580	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1180.2079	2440.4157	0.3681	0.3681
HKYI	40	-1180.6510	2441.3020	0.2363	0.6044
HKYIG	41	-1179.8059	2441.6119	0.2024	0.8067
GTRG	44	-1177.8885	2443.7771	0.0686	0.8753
GTRIG	45	-1177.4153	2444.8306	0.0405	0.9158
GTRI	44	-1178.5194	2445.0388	0.0365	0.9523
K2PG	37	-1186.1388	2446.2776	0.0196	0.9719
K2PI	37	-1186.5653	2447.1305	0.0128	0.9847
K2PIG	38	-1185.7694	2447.5387	0.0105	0.9952
SYMG	41	-1184.2895	2450.5789	0.0023	0.9974
SYMIG	42	-1183.8253	2451.6507	0.0013	0.9988
SYMI	41	-1184.9219	2451.8438	0.0012	1.0000
HKY	39	-1196.7645	2471.5290	0.0000	1.0000
GTR	43	-1194.2091	2474.4182	0.0000	1.0000
K2P	36	-1202.7623	2477.5246	0.0000	1.0000
SYM	40	-1201.0161	2482.0323	0.0000	1.0000
F81G	39	-1210.4327	2498.8655	0.0000	1.0000
F81I	39	-1210.4981	2498.9963	0.0000	1.0000
F81IG	40	-1210.4259	2500.8518	0.0000	1.0000
JC69G	36	-1219.0578	2510.1156	0.0000	1.0000
JC69I	36	-1219.1148	2510.2296	0.0000	1.0000
JC69IG	37	-1219.0494	2512.0989	0.0000	1.0000
F81	38	-1225.9164	2527.8328	0.0000	1.0000
JC69	35	-1234.2552	2538.5104	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1186.1388	2596.8949	0.5574	0.5574
K2PI	37	-1186.5653	2597.7478	0.3639	0.9213
K2PIG	38	-1185.7694	2602.2267	0.0388	0.9600
HKYG	40	-1180.2079	2603.2452	0.0233	0.9833
HKYI	40	-1180.6510	2604.1315	0.0150	0.9983
HKYIG	41	-1179.8059	2608.5121	0.0017	1.0000
SYMG	41	-1184.2895	2617.4792	0.0000	1.0000
SYMI	41	-1184.9219	2618.7441	0.0000	1.0000
SYMIG	42	-1183.8253	2622.6217	0.0000	1.0000
GTRG	44	-1177.8885	2622.8895	0.0000	1.0000
K2P	36	-1202.7623	2624.0711	0.0000	1.0000
GTRI	44	-1178.5194	2624.1513	0.0000	1.0000
GTRIG	45	-1177.4153	2628.0138	0.0000	1.0000
HKY	39	-1196.7645	2630.2878	0.0000	1.0000
SYM	40	-1201.0161	2644.8618	0.0000	1.0000
GTR	43	-1194.2091	2649.4599	0.0000	1.0000
JC69G	36	-1219.0578	2656.6622	0.0000	1.0000
JC69I	36	-1219.1148	2656.7761	0.0000	1.0000
F81G	39	-1210.4327	2657.6243	0.0000	1.0000
F81I	39	-1210.4981	2657.7551	0.0000	1.0000
JC69IG	37	-1219.0494	2662.7162	0.0000	1.0000
F81IG	40	-1210.4259	2663.6813	0.0000	1.0000
JC69	35	-1234.2552	2680.9863	0.0000	1.0000
F81	38	-1225.9164	2682.5209	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_12347.phylip = 1 - 433;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12347.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_12347.phylip = 1 - 433;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12347.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 (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12347.phylip = 1 - 433;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_12347.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

