
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:13 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/chr7_4777.phylip" (19 taxa, 562 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	-2240.3487	4566.9930	0.4895	0.4895
HKYI	40	-2240.9622	4568.2200	0.2650	0.7545
HKYIG	41	-2240.0956	4568.8143	0.1969	0.9514
GTRG	44	-2238.6807	4573.0210	0.0240	0.9754
GTRI	44	-2239.1796	4574.0188	0.0146	0.9900
GTRIG	45	-2238.3789	4574.7811	0.0100	1.0000
K2PG	37	-2255.6703	4590.7070	0.0000	1.0000
K2PI	37	-2256.5124	4592.3913	0.0000	1.0000
K2PIG	38	-2255.4763	4592.6200	0.0000	1.0000
SYMG	41	-2254.5132	4597.6496	0.0000	1.0000
SYMI	41	-2255.1481	4598.9193	0.0000	1.0000
SYMIG	42	-2254.2708	4599.5011	0.0000	1.0000
HKY	39	-2272.3551	4628.6872	0.0000	1.0000
GTR	43	-2270.2540	4633.8130	0.0000	1.0000
K2P	36	-2289.3734	4655.8211	0.0000	1.0000
SYM	40	-2287.8447	4661.9849	0.0000	1.0000
F81G	39	-2306.7990	4697.5751	0.0000	1.0000
F81I	39	-2307.1375	4698.2520	0.0000	1.0000
F81IG	40	-2306.4858	4699.2671	0.0000	1.0000
JC69G	36	-2319.8898	4716.8538	0.0000	1.0000
JC69I	36	-2320.3204	4717.7151	0.0000	1.0000
JC69IG	37	-2319.6034	4718.5732	0.0000	1.0000
F81	38	-2337.3773	4756.4219	0.0000	1.0000
JC69	35	-2351.3015	4777.3938	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2240.3487	4560.6975	0.4505	0.4505
HKYI	40	-2240.9622	4561.9245	0.2439	0.6944
HKYIG	41	-2240.0956	4562.1912	0.2135	0.9079
GTRG	44	-2238.6807	4565.3614	0.0437	0.9517
GTRI	44	-2239.1796	4566.3592	0.0266	0.9782
GTRIG	45	-2238.3789	4566.7579	0.0218	1.0000
K2PG	37	-2255.6703	4585.3406	0.0000	1.0000
K2PIG	38	-2255.4763	4586.9527	0.0000	1.0000
K2PI	37	-2256.5124	4587.0248	0.0000	1.0000
SYMG	41	-2254.5132	4591.0265	0.0000	1.0000
SYMI	41	-2255.1481	4592.2962	0.0000	1.0000
SYMIG	42	-2254.2708	4592.5415	0.0000	1.0000
HKY	39	-2272.3551	4622.7102	0.0000	1.0000
GTR	43	-2270.2540	4626.5079	0.0000	1.0000
K2P	36	-2289.3734	4650.7469	0.0000	1.0000
SYM	40	-2287.8447	4655.6894	0.0000	1.0000
F81G	39	-2306.7990	4691.5980	0.0000	1.0000
F81I	39	-2307.1375	4692.2749	0.0000	1.0000
F81IG	40	-2306.4858	4692.9716	0.0000	1.0000
JC69G	36	-2319.8898	4711.7795	0.0000	1.0000
JC69I	36	-2320.3204	4712.6408	0.0000	1.0000
JC69IG	37	-2319.6034	4713.2067	0.0000	1.0000
F81	38	-2337.3773	4750.7546	0.0000	1.0000
JC69	35	-2351.3015	4772.6029	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2240.3487	4733.9575	0.6249	0.6249
HKYI	40	-2240.9622	4735.1845	0.3384	0.9633
HKYIG	41	-2240.0956	4739.7828	0.0340	0.9972
K2PG	37	-2255.6703	4745.6061	0.0018	0.9991
K2PI	37	-2256.5124	4747.2904	0.0008	0.9999
K2PIG	38	-2255.4763	4751.5497	0.0001	1.0000
GTRG	44	-2238.6807	4755.9475	0.0000	1.0000
GTRI	44	-2239.1796	4756.9453	0.0000	1.0000
GTRIG	45	-2238.3789	4761.6755	0.0000	1.0000
SYMG	41	-2254.5132	4768.6181	0.0000	1.0000
SYMI	41	-2255.1481	4769.8878	0.0000	1.0000
SYMIG	42	-2254.2708	4774.4646	0.0000	1.0000
HKY	39	-2272.3551	4791.6388	0.0000	1.0000
K2P	36	-2289.3734	4806.6809	0.0000	1.0000
GTR	43	-2270.2540	4812.7625	0.0000	1.0000
SYM	40	-2287.8447	4828.9494	0.0000	1.0000
F81G	39	-2306.7990	4860.5266	0.0000	1.0000
F81I	39	-2307.1375	4861.2035	0.0000	1.0000
F81IG	40	-2306.4858	4866.2316	0.0000	1.0000
JC69G	36	-2319.8898	4867.7136	0.0000	1.0000
JC69I	36	-2320.3204	4868.5748	0.0000	1.0000
JC69IG	37	-2319.6034	4873.4723	0.0000	1.0000
F81	38	-2337.3773	4915.3517	0.0000	1.0000
JC69	35	-2351.3015	4924.2055	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/chr7_4777.phylip = 1 - 562;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4777.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/chr7_4777.phylip = 1 - 562;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4777.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/chr7_4777.phylip = 1 - 562;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4777.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

