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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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	-1182.8493	2453.6404	0.4019	0.4019
HKYI	40	-1183.2590	2454.4599	0.2668	0.6687
GTRG	44	-1178.8963	2455.4747	0.1606	0.8293
GTRIG	45	-1178.3907	2456.9284	0.0777	0.9070
HKYIG	41	-1183.7462	2457.8516	0.0489	0.9559
GTRI	44	-1180.1894	2458.0610	0.0441	1.0000
HKY	39	-1204.0212	2493.5787	0.0000	1.0000
GTR	43	-1202.7622	2500.7537	0.0000	1.0000
K2PIG	38	-1209.6934	2502.5289	0.0000	1.0000
K2PG	37	-1210.9805	2502.7206	0.0000	1.0000
K2PI	37	-1211.7247	2504.2090	0.0000	1.0000
SYMG	41	-1209.8246	2510.0083	0.0000	1.0000
SYMI	41	-1210.2707	2510.9006	0.0000	1.0000
SYMIG	42	-1210.2940	2513.3764	0.0000	1.0000
F81G	39	-1227.8837	2541.3036	0.0000	1.0000
F81IG	40	-1227.8869	2543.7156	0.0000	1.0000
F81I	39	-1229.7282	2544.9926	0.0000	1.0000
K2P	36	-1234.5274	2547.4433	0.0000	1.0000
SYM	40	-1231.9561	2551.8540	0.0000	1.0000
JC69G	36	-1251.8858	2582.1601	0.0000	1.0000
F81	38	-1251.7481	2586.6383	0.0000	1.0000
JC69IG	37	-1254.3517	2589.4630	0.0000	1.0000
JC69I	36	-1256.0159	2590.4202	0.0000	1.0000
JC69	35	-1272.4380	2620.9048	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1182.8493	2445.6985	0.2769	0.2769
GTRG	44	-1178.8963	2445.7926	0.2642	0.5410
HKYI	40	-1183.2590	2446.5180	0.1838	0.7248
GTRIG	45	-1178.3907	2446.7814	0.1611	0.8860
GTRI	44	-1180.1894	2448.3788	0.0725	0.9585
HKYIG	41	-1183.7462	2449.4924	0.0415	1.0000
HKY	39	-1204.0212	2486.0424	0.0000	1.0000
GTR	43	-1202.7622	2491.5244	0.0000	1.0000
K2PIG	38	-1209.6934	2495.3868	0.0000	1.0000
K2PG	37	-1210.9805	2495.9610	0.0000	1.0000
K2PI	37	-1211.7247	2497.4494	0.0000	1.0000
SYMG	41	-1209.8246	2501.6491	0.0000	1.0000
SYMI	41	-1210.2707	2502.5413	0.0000	1.0000
SYMIG	42	-1210.2940	2504.5881	0.0000	1.0000
F81G	39	-1227.8837	2533.7674	0.0000	1.0000
F81IG	40	-1227.8869	2535.7738	0.0000	1.0000
F81I	39	-1229.7282	2537.4564	0.0000	1.0000
K2P	36	-1234.5274	2541.0549	0.0000	1.0000
SYM	40	-1231.9561	2543.9121	0.0000	1.0000
JC69G	36	-1251.8858	2575.7717	0.0000	1.0000
F81	38	-1251.7481	2579.4961	0.0000	1.0000
JC69IG	37	-1254.3517	2582.7034	0.0000	1.0000
JC69I	36	-1256.0159	2584.0317	0.0000	1.0000
JC69	35	-1272.4380	2614.8761	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1182.8493	2610.4224	0.5941	0.5941
HKYI	40	-1183.2590	2611.2419	0.3944	0.9884
HKYIG	41	-1183.7462	2618.3343	0.0114	0.9998
GTRG	44	-1178.8963	2626.9889	0.0002	0.9999
GTRI	44	-1180.1894	2629.5751	0.0000	1.0000
GTRIG	45	-1178.3907	2632.0958	0.0000	1.0000
HKY	39	-1204.0212	2646.6482	0.0000	1.0000
K2PG	37	-1210.9805	2648.3306	0.0000	1.0000
K2PI	37	-1211.7247	2649.8190	0.0000	1.0000
K2PIG	38	-1209.6934	2651.8745	0.0000	1.0000
GTR	43	-1202.7622	2668.6026	0.0000	1.0000
SYMG	41	-1209.8246	2670.4911	0.0000	1.0000
SYMI	41	-1210.2707	2671.3833	0.0000	1.0000
SYMIG	42	-1210.2940	2677.5482	0.0000	1.0000
K2P	36	-1234.5274	2689.3064	0.0000	1.0000
F81G	39	-1227.8837	2694.3732	0.0000	1.0000
F81I	39	-1229.7282	2698.0622	0.0000	1.0000
F81IG	40	-1227.8869	2700.4976	0.0000	1.0000
SYM	40	-1231.9561	2708.6360	0.0000	1.0000
JC69G	36	-1251.8858	2724.0232	0.0000	1.0000
JC69I	36	-1256.0159	2732.2832	0.0000	1.0000
JC69IG	37	-1254.3517	2735.0730	0.0000	1.0000
F81	38	-1251.7481	2735.9838	0.0000	1.0000
JC69	35	-1272.4380	2759.0095	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/chr3_6248.phylip = 1 - 454;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_6248.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/chr3_6248.phylip = 1 - 454;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_6248.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 (HKYG)]

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


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

