
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:32 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/chr23_16.phylip" (19 taxa, 736 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	-2509.3635	5103.4464	0.3560	0.3560
HKYIG	41	-2508.8233	5104.6092	0.1990	0.5550
SYMG	41	-2509.5489	5106.0603	0.0964	0.6514
GTRG	44	-2506.2128	5106.1565	0.0918	0.7432
K2PG	37	-2514.0903	5106.2092	0.0894	0.8326
SYMIG	42	-2508.9969	5107.2060	0.0543	0.8870
K2PIG	38	-2513.5434	5107.3394	0.0508	0.9378
GTRIG	45	-2505.6877	5107.3753	0.0499	0.9877
HKYI	40	-2513.2838	5111.2871	0.0071	0.9948
K2PI	37	-2517.8963	5113.8212	0.0020	0.9968
SYMI	41	-2513.5872	5114.1369	0.0017	0.9985
GTRI	44	-2510.3147	5114.3602	0.0015	1.0000
HKY	39	-2548.6922	5179.8672	0.0000	1.0000
GTR	43	-2545.7895	5183.0471	0.0000	1.0000
K2P	36	-2553.7544	5183.3199	0.0000	1.0000
SYM	40	-2549.7014	5184.1222	0.0000	1.0000
F81G	39	-2558.4313	5199.3454	0.0000	1.0000
F81IG	40	-2558.4431	5201.6056	0.0000	1.0000
JC69G	36	-2564.1668	5204.1448	0.0000	1.0000
F81I	39	-2561.7620	5206.0068	0.0000	1.0000
JC69IG	37	-2564.1845	5206.3976	0.0000	1.0000
JC69I	36	-2567.7580	5211.3272	0.0000	1.0000
F81	38	-2594.3322	5268.9169	0.0000	1.0000
JC69	35	-2600.0609	5273.7219	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2509.3635	5098.7270	0.3171	0.3171
HKYIG	41	-2508.8233	5099.6467	0.2002	0.5174
GTRG	44	-2506.2128	5100.4256	0.1356	0.6530
SYMG	41	-2509.5489	5101.0977	0.0969	0.7500
GTRIG	45	-2505.6877	5101.3753	0.0844	0.8343
SYMIG	42	-2508.9969	5101.9938	0.0619	0.8962
K2PG	37	-2514.0903	5102.1805	0.0564	0.9527
K2PIG	38	-2513.5434	5103.0869	0.0359	0.9885
HKYI	40	-2513.2838	5106.5677	0.0063	0.9948
GTRI	44	-2510.3147	5108.6294	0.0022	0.9970
SYMI	41	-2513.5872	5109.1744	0.0017	0.9987
K2PI	37	-2517.8963	5109.7925	0.0013	1.0000
HKY	39	-2548.6922	5175.3845	0.0000	1.0000
GTR	43	-2545.7895	5177.5789	0.0000	1.0000
SYM	40	-2549.7014	5179.4028	0.0000	1.0000
K2P	36	-2553.7544	5179.5087	0.0000	1.0000
F81G	39	-2558.4313	5194.8626	0.0000	1.0000
F81IG	40	-2558.4431	5196.8862	0.0000	1.0000
JC69G	36	-2564.1668	5200.3337	0.0000	1.0000
F81I	39	-2561.7620	5201.5241	0.0000	1.0000
JC69IG	37	-2564.1845	5202.3689	0.0000	1.0000
JC69I	36	-2567.7580	5207.5160	0.0000	1.0000
F81	38	-2594.3322	5264.6644	0.0000	1.0000
JC69	35	-2600.0609	5270.1219	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-2514.0903	5272.4261	0.9156	0.9156
K2PIG	38	-2513.5434	5277.9336	0.0583	0.9739
K2PI	37	-2517.8963	5280.0381	0.0204	0.9942
HKYG	40	-2509.3635	5282.7762	0.0052	0.9994
HKYIG	41	-2508.8233	5288.2971	0.0003	0.9997
SYMG	41	-2509.5489	5289.7482	0.0002	0.9999
HKYI	40	-2513.2838	5290.6169	0.0001	1.0000
SYMIG	42	-2508.9969	5295.2455	0.0000	1.0000
SYMI	41	-2513.5872	5297.8248	0.0000	1.0000
GTRG	44	-2506.2128	5302.8798	0.0000	1.0000
GTRIG	45	-2505.6877	5308.4307	0.0000	1.0000
GTRI	44	-2510.3147	5311.0835	0.0000	1.0000
K2P	36	-2553.7544	5345.1530	0.0000	1.0000
HKY	39	-2548.6922	5354.8324	0.0000	1.0000
SYM	40	-2549.7014	5363.4520	0.0000	1.0000
JC69G	36	-2564.1668	5365.9780	0.0000	1.0000
JC69IG	37	-2564.1845	5372.6144	0.0000	1.0000
JC69I	36	-2567.7580	5373.1603	0.0000	1.0000
F81G	39	-2558.4313	5374.3106	0.0000	1.0000
GTR	43	-2545.7895	5375.4318	0.0000	1.0000
F81IG	40	-2558.4431	5380.9354	0.0000	1.0000
F81I	39	-2561.7620	5380.9721	0.0000	1.0000
JC69	35	-2600.0609	5431.1649	0.0000	1.0000
F81	38	-2594.3322	5439.5112	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/chr23_16.phylip = 1 - 736;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr23_16.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/chr23_16.phylip = 1 - 736;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr23_16.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/chr23_16.phylip = 1 - 736;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr23_16.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

