
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:48 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_25863.phylip" (19 taxa, 657 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	-2091.8884	4269.1016	0.5587	0.5587
HKYIG	41	-2091.7662	4271.1324	0.2024	0.7610
HKYI	40	-2092.9270	4271.1786	0.1978	0.9588
GTRG	44	-2090.4834	4275.4374	0.0235	0.9823
GTRI	44	-2091.3479	4277.1664	0.0099	0.9922
GTRIG	45	-2090.4345	4277.6448	0.0078	1.0000
K2PG	37	-2114.6100	4307.7628	0.0000	1.0000
K2PIG	38	-2114.4153	4309.6268	0.0000	1.0000
SYMIG	42	-2112.4347	4314.7520	0.0000	1.0000
SYMI	41	-2113.9443	4315.4887	0.0000	1.0000
K2PI	37	-2119.4351	4317.4131	0.0000	1.0000
SYMG	41	-2116.1640	4319.9280	0.0000	1.0000
F81G	39	-2126.0468	4335.1503	0.0000	1.0000
F81I	39	-2126.7230	4336.5027	0.0000	1.0000
F81IG	40	-2126.0267	4337.3782	0.0000	1.0000
JC69IG	37	-2145.8324	4370.2076	0.0000	1.0000
GTR	43	-2139.3654	4370.9037	0.0000	1.0000
HKY	39	-2144.9547	4372.9662	0.0000	1.0000
JC69G	36	-2149.0571	4374.4110	0.0000	1.0000
JC69I	36	-2149.9541	4376.2050	0.0000	1.0000
K2P	36	-2164.4325	4405.1619	0.0000	1.0000
SYM	40	-2162.4948	4410.3142	0.0000	1.0000
F81	38	-2175.0012	4430.7984	0.0000	1.0000
JC69	35	-2192.6097	4459.2773	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2091.8884	4263.7769	0.5250	0.5250
HKYIG	41	-2091.7662	4265.5324	0.2183	0.7433
HKYI	40	-2092.9270	4265.8539	0.1859	0.9292
GTRG	44	-2090.4834	4268.9669	0.0392	0.9683
GTRI	44	-2091.3479	4270.6958	0.0165	0.9849
GTRIG	45	-2090.4345	4270.8690	0.0151	1.0000
K2PG	37	-2114.6100	4303.2200	0.0000	1.0000
K2PIG	38	-2114.4153	4304.8307	0.0000	1.0000
SYMIG	42	-2112.4347	4308.8693	0.0000	1.0000
SYMI	41	-2113.9443	4309.8887	0.0000	1.0000
K2PI	37	-2119.4351	4312.8702	0.0000	1.0000
SYMG	41	-2116.1640	4314.3280	0.0000	1.0000
F81G	39	-2126.0468	4330.0936	0.0000	1.0000
F81I	39	-2126.7230	4331.4460	0.0000	1.0000
F81IG	40	-2126.0267	4332.0535	0.0000	1.0000
GTR	43	-2139.3654	4364.7308	0.0000	1.0000
JC69IG	37	-2145.8324	4365.6648	0.0000	1.0000
HKY	39	-2144.9547	4367.9095	0.0000	1.0000
JC69G	36	-2149.0571	4370.1143	0.0000	1.0000
JC69I	36	-2149.9541	4371.9082	0.0000	1.0000
K2P	36	-2164.4325	4400.8651	0.0000	1.0000
SYM	40	-2162.4948	4404.9895	0.0000	1.0000
F81	38	-2175.0012	4426.0023	0.0000	1.0000
JC69	35	-2192.6097	4455.2193	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2091.8884	4443.2842	0.7153	0.7153
HKYI	40	-2092.9270	4445.3613	0.2532	0.9685
HKYIG	41	-2091.7662	4449.5274	0.0315	1.0000
GTRG	44	-2090.4834	4466.4250	0.0000	1.0000
GTRI	44	-2091.3479	4468.1539	0.0000	1.0000
K2PG	37	-2114.6100	4469.2643	0.0000	1.0000
GTRIG	45	-2090.4345	4472.8148	0.0000	1.0000
K2PIG	38	-2114.4153	4475.3627	0.0000	1.0000
K2PI	37	-2119.4351	4478.9145	0.0000	1.0000
SYMI	41	-2113.9443	4493.8837	0.0000	1.0000
SYMIG	42	-2112.4347	4497.3520	0.0000	1.0000
SYMG	41	-2116.1640	4498.3231	0.0000	1.0000
F81G	39	-2126.0468	4505.1133	0.0000	1.0000
F81I	39	-2126.7230	4506.4657	0.0000	1.0000
F81IG	40	-2126.0267	4511.5608	0.0000	1.0000
JC69G	36	-2149.0571	4531.6709	0.0000	1.0000
JC69IG	37	-2145.8324	4531.7091	0.0000	1.0000
JC69I	36	-2149.9541	4533.4648	0.0000	1.0000
HKY	39	-2144.9547	4542.9292	0.0000	1.0000
GTR	43	-2139.3654	4557.7012	0.0000	1.0000
K2P	36	-2164.4325	4562.4217	0.0000	1.0000
SYM	40	-2162.4948	4584.4969	0.0000	1.0000
F81	38	-2175.0012	4596.5343	0.0000	1.0000
JC69	35	-2192.6097	4612.2883	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_25863.phylip = 1 - 657;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_25863.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_25863.phylip = 1 - 657;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_25863.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/chr2_25863.phylip = 1 - 657;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_25863.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

