
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:11 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_8698.phylip" (19 taxa, 626 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	-1615.7967	3317.2002	0.5448	0.5448
HKYI	40	-1616.6151	3318.8371	0.2403	0.7851
HKYIG	41	-1615.7670	3319.4313	0.1786	0.9637
GTRG	44	-1614.4840	3323.7839	0.0203	0.9839
GTRI	44	-1615.3561	3325.5279	0.0085	0.9924
GTRIG	45	-1614.4569	3326.0518	0.0065	0.9989
HKY	39	-1623.2214	3329.7670	0.0010	0.9999
GTR	43	-1621.6133	3335.7284	0.0001	1.0000
K2PG	37	-1654.6195	3388.0213	0.0000	1.0000
K2PI	37	-1655.3576	3389.4975	0.0000	1.0000
K2PIG	38	-1654.5915	3390.2324	0.0000	1.0000
SYMG	41	-1651.8051	3391.5074	0.0000	1.0000
SYMI	41	-1652.5994	3393.0962	0.0000	1.0000
SYMIG	42	-1651.7740	3393.7436	0.0000	1.0000
F81G	39	-1657.9301	3399.1844	0.0000	1.0000
F81I	39	-1658.3675	3400.0593	0.0000	1.0000
K2P	36	-1662.2506	3401.0240	0.0000	1.0000
F81IG	40	-1657.9250	3401.4568	0.0000	1.0000
SYM	40	-1659.2162	3404.0393	0.0000	1.0000
F81	38	-1664.7077	3410.4647	0.0000	1.0000
JC69G	36	-1691.1619	3458.8467	0.0000	1.0000
JC69I	36	-1691.6144	3459.7516	0.0000	1.0000
JC69IG	37	-1691.1599	3461.1020	0.0000	1.0000
JC69	35	-1698.1806	3470.6324	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1615.7967	3311.5933	0.5144	0.5144
HKYI	40	-1616.6151	3313.2302	0.2269	0.7413
HKYIG	41	-1615.7670	3313.5340	0.1949	0.9362
GTRG	44	-1614.4840	3316.9680	0.0350	0.9712
GTRI	44	-1615.3561	3318.7121	0.0146	0.9859
GTRIG	45	-1614.4569	3318.9138	0.0132	0.9991
HKY	39	-1623.2214	3324.4427	0.0008	0.9999
GTR	43	-1621.6133	3329.2266	0.0001	1.0000
K2PG	37	-1654.6195	3383.2390	0.0000	1.0000
K2PI	37	-1655.3576	3384.7152	0.0000	1.0000
K2PIG	38	-1654.5915	3385.1830	0.0000	1.0000
SYMG	41	-1651.8051	3385.6102	0.0000	1.0000
SYMI	41	-1652.5994	3387.1989	0.0000	1.0000
SYMIG	42	-1651.7740	3387.5481	0.0000	1.0000
F81G	39	-1657.9301	3393.8601	0.0000	1.0000
F81I	39	-1658.3675	3394.7351	0.0000	1.0000
F81IG	40	-1657.9250	3395.8499	0.0000	1.0000
K2P	36	-1662.2506	3396.5011	0.0000	1.0000
SYM	40	-1659.2162	3398.4325	0.0000	1.0000
F81	38	-1664.7077	3405.4153	0.0000	1.0000
JC69G	36	-1691.1619	3454.3238	0.0000	1.0000
JC69I	36	-1691.6144	3455.2287	0.0000	1.0000
JC69IG	37	-1691.1599	3456.3197	0.0000	1.0000
JC69	35	-1698.1806	3466.3612	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1615.7967	3489.1674	0.6679	0.6679
HKYI	40	-1616.6151	3490.8042	0.2946	0.9625
HKYIG	41	-1615.7670	3495.5474	0.0275	0.9900
HKY	39	-1623.2214	3497.5774	0.0100	1.0000
GTRG	44	-1614.4840	3512.2994	0.0000	1.0000
GTRI	44	-1615.3561	3514.0435	0.0000	1.0000
GTRIG	45	-1614.4569	3518.6846	0.0000	1.0000
GTR	43	-1621.6133	3520.1187	0.0000	1.0000
K2PG	37	-1654.6195	3547.4950	0.0000	1.0000
K2PI	37	-1655.3576	3548.9711	0.0000	1.0000
K2PIG	38	-1654.5915	3553.8783	0.0000	1.0000
K2P	36	-1662.2506	3556.3177	0.0000	1.0000
F81G	39	-1657.9301	3566.9948	0.0000	1.0000
SYMG	41	-1651.8051	3567.6235	0.0000	1.0000
F81I	39	-1658.3675	3567.8697	0.0000	1.0000
SYMI	41	-1652.5994	3569.2123	0.0000	1.0000
F81IG	40	-1657.9250	3573.4239	0.0000	1.0000
SYMIG	42	-1651.7740	3574.0008	0.0000	1.0000
F81	38	-1664.7077	3574.1107	0.0000	1.0000
SYM	40	-1659.2162	3576.0065	0.0000	1.0000
JC69G	36	-1691.1619	3614.1404	0.0000	1.0000
JC69I	36	-1691.6144	3615.0453	0.0000	1.0000
JC69IG	37	-1691.1599	3620.5757	0.0000	1.0000
JC69	35	-1698.1806	3621.7385	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_8698.phylip = 1 - 626;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8698.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_8698.phylip = 1 - 626;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8698.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_8698.phylip = 1 - 626;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8698.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

