
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:18 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_24903.phylip" (19 taxa, 523 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	-1555.9534	3198.7118	0.9693	0.9693
HKYIG	41	-1558.7388	3206.6377	0.0184	0.9877
HKYI	40	-1560.9109	3208.6267	0.0068	0.9945
GTRG	44	-1557.0672	3210.4189	0.0028	0.9973
GTRIG	45	-1556.5775	3211.8342	0.0014	0.9987
K2PG	37	-1566.5707	3212.9393	0.0008	0.9995
GTRI	44	-1558.7609	3213.8063	0.0005	1.0000
K2PIG	38	-1569.3230	3220.7700	0.0000	1.0000
K2PI	37	-1571.3827	3222.5633	0.0000	1.0000
SYMG	41	-1567.1330	3223.4260	0.0000	1.0000
SYMIG	42	-1566.4946	3224.5142	0.0000	1.0000
SYMI	41	-1568.8418	3226.8437	0.0000	1.0000
F81G	39	-1578.5389	3241.5374	0.0000	1.0000
F81IG	40	-1578.0830	3242.9709	0.0000	1.0000
HKY	39	-1579.8753	3244.2102	0.0000	1.0000
F81I	39	-1580.0258	3244.5113	0.0000	1.0000
JC69G	36	-1588.7369	3254.9553	0.0000	1.0000
JC69IG	37	-1588.2377	3256.2734	0.0000	1.0000
GTR	43	-1582.0074	3257.9146	0.0000	1.0000
JC69I	36	-1590.2353	3257.9521	0.0000	1.0000
K2P	36	-1591.6887	3260.8589	0.0000	1.0000
SYM	40	-1594.1089	3275.0228	0.0000	1.0000
F81	38	-1604.7653	3291.6545	0.0000	1.0000
JC69	35	-1615.4959	3306.1664	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1555.9534	3191.9068	0.9606	0.9606
HKYIG	41	-1558.7388	3199.4776	0.0218	0.9824
HKYI	40	-1560.9109	3201.8217	0.0068	0.9892
GTRG	44	-1557.0672	3202.1344	0.0058	0.9950
GTRIG	45	-1556.5775	3203.1550	0.0035	0.9984
GTRI	44	-1558.7609	3205.5217	0.0011	0.9995
K2PG	37	-1566.5707	3207.1413	0.0005	1.0000
K2PIG	38	-1569.3230	3214.6461	0.0000	1.0000
SYMG	41	-1567.1330	3216.2660	0.0000	1.0000
K2PI	37	-1571.3827	3216.7653	0.0000	1.0000
SYMIG	42	-1566.4946	3216.9892	0.0000	1.0000
SYMI	41	-1568.8418	3219.6836	0.0000	1.0000
F81G	39	-1578.5389	3235.0778	0.0000	1.0000
F81IG	40	-1578.0830	3236.1659	0.0000	1.0000
HKY	39	-1579.8753	3237.7506	0.0000	1.0000
F81I	39	-1580.0258	3238.0516	0.0000	1.0000
JC69G	36	-1588.7369	3249.4739	0.0000	1.0000
GTR	43	-1582.0074	3250.0148	0.0000	1.0000
JC69IG	37	-1588.2377	3250.4755	0.0000	1.0000
JC69I	36	-1590.2353	3252.4706	0.0000	1.0000
K2P	36	-1591.6887	3255.3774	0.0000	1.0000
SYM	40	-1594.1089	3268.2178	0.0000	1.0000
F81	38	-1604.7653	3285.5305	0.0000	1.0000
JC69	35	-1615.4959	3300.9918	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1555.9534	3362.2901	0.7658	0.7658
K2PG	37	-1566.5707	3364.7459	0.2243	0.9901
HKYI	40	-1560.9109	3372.2050	0.0054	0.9955
HKYIG	41	-1558.7388	3374.1204	0.0021	0.9975
K2PI	37	-1571.3827	3374.3699	0.0018	0.9994
K2PIG	38	-1569.3230	3376.5102	0.0006	1.0000
GTRG	44	-1557.0672	3389.5560	0.0000	1.0000
SYMG	41	-1567.1330	3390.9088	0.0000	1.0000
GTRI	44	-1558.7609	3392.9433	0.0000	1.0000
SYMI	41	-1568.8418	3394.3265	0.0000	1.0000
GTRIG	45	-1556.5775	3394.8362	0.0000	1.0000
SYMIG	42	-1566.4946	3395.8916	0.0000	1.0000
F81G	39	-1578.5389	3401.2014	0.0000	1.0000
JC69G	36	-1588.7369	3402.8188	0.0000	1.0000
HKY	39	-1579.8753	3403.8743	0.0000	1.0000
F81I	39	-1580.0258	3404.1753	0.0000	1.0000
JC69I	36	-1590.2353	3405.8156	0.0000	1.0000
F81IG	40	-1578.0830	3406.5492	0.0000	1.0000
JC69IG	37	-1588.2377	3408.0800	0.0000	1.0000
K2P	36	-1591.6887	3408.7224	0.0000	1.0000
GTR	43	-1582.0074	3433.1768	0.0000	1.0000
SYM	40	-1594.1089	3438.6011	0.0000	1.0000
F81	38	-1604.7653	3447.3946	0.0000	1.0000
JC69	35	-1615.4959	3450.0772	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_24903.phylip = 1 - 523;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24903.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_24903.phylip = 1 - 523;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24903.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_24903.phylip = 1 - 523;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24903.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

