
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:51 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/chr2_26586.phylip" (19 taxa, 330 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	-1196.9499	2485.2493	0.4467	0.4467
HKYI	40	-1197.2879	2485.9253	0.3186	0.7653
HKYIG	41	-1196.6140	2487.1863	0.1696	0.9349
K2PG	37	-1203.6683	2490.9668	0.0256	0.9605
K2PI	37	-1204.5238	2492.6778	0.0109	0.9714
K2PIG	38	-1203.4084	2493.0024	0.0093	0.9806
GTRG	44	-1195.6543	2493.2034	0.0084	0.9890
GTRI	44	-1196.1447	2494.1841	0.0051	0.9941
GTRIG	45	-1195.1587	2494.8949	0.0036	0.9977
SYMG	41	-1201.4719	2496.9022	0.0013	0.9991
SYMIG	42	-1200.9797	2498.5447	0.0006	0.9996
SYMI	41	-1202.7633	2499.4850	0.0004	1.0000
F81G	39	-1212.5785	2513.9157	0.0000	1.0000
F81I	39	-1213.0739	2514.9064	0.0000	1.0000
F81IG	40	-1212.7523	2516.8540	0.0000	1.0000
JC69G	36	-1218.4943	2518.0807	0.0000	1.0000
JC69I	36	-1219.2305	2519.5531	0.0000	1.0000
JC69IG	37	-1218.2947	2520.2195	0.0000	1.0000
HKY	39	-1246.8597	2582.4780	0.0000	1.0000
SYM	40	-1245.6984	2582.7463	0.0000	1.0000
GTR	43	-1241.8580	2582.9468	0.0000	1.0000
K2P	36	-1251.5922	2584.2765	0.0000	1.0000
F81	38	-1260.1441	2606.4738	0.0000	1.0000
JC69	35	-1264.3035	2607.1785	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1196.9499	2473.8998	0.4115	0.4115
HKYI	40	-1197.2879	2474.5759	0.2935	0.7050
HKYIG	41	-1196.6140	2475.2279	0.2118	0.9169
GTRG	44	-1195.6543	2479.3086	0.0275	0.9444
GTRI	44	-1196.1447	2480.2894	0.0169	0.9613
GTRIG	45	-1195.1587	2480.3174	0.0166	0.9779
K2PG	37	-1203.6683	2481.3366	0.0100	0.9879
K2PIG	38	-1203.4084	2482.8169	0.0048	0.9927
K2PI	37	-1204.5238	2483.0477	0.0042	0.9969
SYMG	41	-1201.4719	2484.9438	0.0016	0.9986
SYMIG	42	-1200.9797	2485.9594	0.0010	0.9995
SYMI	41	-1202.7633	2487.5267	0.0005	1.0000
F81G	39	-1212.5785	2503.1571	0.0000	1.0000
F81I	39	-1213.0739	2504.1478	0.0000	1.0000
F81IG	40	-1212.7523	2505.5046	0.0000	1.0000
JC69G	36	-1218.4943	2508.9885	0.0000	1.0000
JC69I	36	-1219.2305	2510.4609	0.0000	1.0000
JC69IG	37	-1218.2947	2510.5894	0.0000	1.0000
GTR	43	-1241.8580	2569.7161	0.0000	1.0000
SYM	40	-1245.6984	2571.3969	0.0000	1.0000
HKY	39	-1246.8597	2571.7194	0.0000	1.0000
K2P	36	-1251.5922	2575.1843	0.0000	1.0000
F81	38	-1260.1441	2596.2883	0.0000	1.0000
JC69	35	-1264.3035	2598.6071	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1203.6683	2621.9031	0.5735	0.5735
K2PI	37	-1204.5238	2623.6141	0.2438	0.8173
HKYG	40	-1196.9499	2625.8635	0.0792	0.8964
HKYI	40	-1197.2879	2626.5396	0.0565	0.9529
K2PIG	38	-1203.4084	2627.1824	0.0409	0.9938
HKYIG	41	-1196.6140	2630.9907	0.0061	0.9999
SYMG	41	-1201.4719	2640.7066	0.0000	1.0000
SYMI	41	-1202.7633	2643.2895	0.0000	1.0000
SYMIG	42	-1200.9797	2645.5213	0.0000	1.0000
JC69G	36	-1218.4943	2645.7559	0.0000	1.0000
GTRG	44	-1195.6543	2646.4687	0.0000	1.0000
JC69I	36	-1219.2305	2647.2283	0.0000	1.0000
GTRI	44	-1196.1447	2647.4494	0.0000	1.0000
JC69IG	37	-1218.2947	2651.1558	0.0000	1.0000
GTRIG	45	-1195.1587	2651.2766	0.0000	1.0000
F81G	39	-1212.5785	2651.3217	0.0000	1.0000
F81I	39	-1213.0739	2652.3124	0.0000	1.0000
F81IG	40	-1212.7523	2657.4683	0.0000	1.0000
K2P	36	-1251.5922	2711.9517	0.0000	1.0000
HKY	39	-1246.8597	2719.8840	0.0000	1.0000
SYM	40	-1245.6984	2723.3606	0.0000	1.0000
JC69	35	-1264.3035	2731.5753	0.0000	1.0000
GTR	43	-1241.8580	2733.0771	0.0000	1.0000
F81	38	-1260.1441	2740.6538	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/chr2_26586.phylip = 1 - 330;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26586.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/chr2_26586.phylip = 1 - 330;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26586.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/chr2_26586.phylip = 1 - 330;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26586.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

