
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:34 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_14426.phylip" (19 taxa, 476 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-1481.3149	3050.1700	0.4634	0.4634
HKYG	40	-1481.4525	3050.4452	0.4039	0.8673
HKYIG	41	-1481.4523	3052.8401	0.1219	0.9892
GTRG	44	-1480.9514	3059.0907	0.0054	0.9946
GTRI	44	-1481.5012	3060.1902	0.0031	0.9977
GTRIG	45	-1480.5631	3060.7541	0.0023	1.0000
HKY	39	-1507.4479	3100.0518	0.0000	1.0000
K2PI	37	-1511.1106	3102.6413	0.0000	1.0000
K2PG	37	-1511.5428	3103.5058	0.0000	1.0000
K2PIG	38	-1511.5427	3105.8681	0.0000	1.0000
GTR	43	-1508.4196	3111.5985	0.0000	1.0000
SYMG	41	-1512.8917	3115.7188	0.0000	1.0000
SYMI	41	-1513.1436	3116.2228	0.0000	1.0000
SYMIG	42	-1512.6008	3117.5434	0.0000	1.0000
F81G	39	-1517.9832	3121.1223	0.0000	1.0000
F81I	39	-1518.6482	3122.4523	0.0000	1.0000
F81IG	40	-1517.6946	3122.9295	0.0000	1.0000
K2P	36	-1539.0198	3156.1079	0.0000	1.0000
JC69G	36	-1543.9268	3165.9219	0.0000	1.0000
JC69I	36	-1544.4591	3166.9865	0.0000	1.0000
JC69IG	37	-1543.6780	3167.7762	0.0000	1.0000
SYM	40	-1542.1757	3171.8916	0.0000	1.0000
F81	38	-1544.6357	3172.0540	0.0000	1.0000
JC69	35	-1571.4964	3218.7202	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1481.3149	3042.6298	0.4449	0.4449
HKYG	40	-1481.4525	3042.9050	0.3877	0.8325
HKYIG	41	-1481.4523	3044.9047	0.1426	0.9752
GTRG	44	-1480.9514	3049.9027	0.0117	0.9869
GTRI	44	-1481.5012	3051.0023	0.0068	0.9936
GTRIG	45	-1480.5631	3051.1262	0.0064	1.0000
HKY	39	-1507.4479	3092.8959	0.0000	1.0000
K2PI	37	-1511.1106	3096.2212	0.0000	1.0000
K2PG	37	-1511.5428	3097.0857	0.0000	1.0000
K2PIG	38	-1511.5427	3099.0855	0.0000	1.0000
GTR	43	-1508.4196	3102.8392	0.0000	1.0000
SYMG	41	-1512.8917	3107.7833	0.0000	1.0000
SYMI	41	-1513.1436	3108.2873	0.0000	1.0000
SYMIG	42	-1512.6008	3109.2016	0.0000	1.0000
F81G	39	-1517.9832	3113.9663	0.0000	1.0000
F81I	39	-1518.6482	3115.2964	0.0000	1.0000
F81IG	40	-1517.6946	3115.3892	0.0000	1.0000
K2P	36	-1539.0198	3150.0396	0.0000	1.0000
JC69G	36	-1543.9268	3159.8536	0.0000	1.0000
JC69I	36	-1544.4591	3160.9181	0.0000	1.0000
JC69IG	37	-1543.6780	3161.3561	0.0000	1.0000
SYM	40	-1542.1757	3164.3514	0.0000	1.0000
F81	38	-1544.6357	3165.2714	0.0000	1.0000
JC69	35	-1571.4964	3212.9929	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1481.3149	3209.2465	0.5232	0.5232
HKYG	40	-1481.4525	3209.5217	0.4559	0.9791
HKYIG	41	-1481.4523	3215.6868	0.0209	1.0000
GTRG	44	-1480.9514	3233.1811	0.0000	1.0000
GTRI	44	-1481.5012	3234.2807	0.0000	1.0000
GTRIG	45	-1480.5631	3238.5700	0.0000	1.0000
K2PI	37	-1511.1106	3250.3416	0.0000	1.0000
K2PG	37	-1511.5428	3251.2061	0.0000	1.0000
HKY	39	-1507.4479	3255.3472	0.0000	1.0000
K2PIG	38	-1511.5427	3257.3714	0.0000	1.0000
F81G	39	-1517.9832	3276.4176	0.0000	1.0000
F81I	39	-1518.6482	3277.7477	0.0000	1.0000
SYMG	41	-1512.8917	3278.5655	0.0000	1.0000
SYMI	41	-1513.1436	3279.0694	0.0000	1.0000
GTR	43	-1508.4196	3281.9522	0.0000	1.0000
F81IG	40	-1517.6946	3282.0060	0.0000	1.0000
SYMIG	42	-1512.6008	3284.1491	0.0000	1.0000
K2P	36	-1539.0198	3299.9946	0.0000	1.0000
JC69G	36	-1543.9268	3309.8086	0.0000	1.0000
JC69I	36	-1544.4591	3310.8732	0.0000	1.0000
JC69IG	37	-1543.6780	3315.4765	0.0000	1.0000
F81	38	-1544.6357	3323.5573	0.0000	1.0000
SYM	40	-1542.1757	3330.9681	0.0000	1.0000
JC69	35	-1571.4964	3358.7825	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14426.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14426.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14426.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14426.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14426.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_14426.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

