
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:10 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/chrUn_random_7197.phylip" (19 taxa, 650 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	-1678.8800	3443.1458	0.4492	0.4492
HKYI	40	-1679.4145	3444.2150	0.2632	0.7125
HKYIG	41	-1678.7134	3445.0913	0.1698	0.8823
GTRG	44	-1676.3193	3447.1841	0.0596	0.9419
GTRI	44	-1676.8136	3448.1726	0.0364	0.9783
GTRIG	45	-1676.1781	3449.2105	0.0217	1.0000
K2PG	37	-1693.0422	3464.6791	0.0000	1.0000
K2PI	37	-1693.5104	3465.6155	0.0000	1.0000
K2PIG	38	-1692.8677	3466.5864	0.0000	1.0000
HKY	39	-1693.1228	3469.3603	0.0000	1.0000
SYMG	41	-1691.6949	3471.0543	0.0000	1.0000
SYMI	41	-1692.1168	3471.8981	0.0000	1.0000
GTR	43	-1690.3389	3472.9220	0.0000	1.0000
SYMIG	42	-1691.5509	3473.0523	0.0000	1.0000
F81G	39	-1704.7042	3492.5231	0.0000	1.0000
K2P	36	-1708.2398	3492.8255	0.0000	1.0000
F81I	39	-1705.0544	3493.2236	0.0000	1.0000
F81IG	40	-1704.6308	3494.6476	0.0000	1.0000
SYM	40	-1706.6044	3498.5947	0.0000	1.0000
JC69G	36	-1718.2892	3512.9243	0.0000	1.0000
JC69I	36	-1718.6538	3513.6535	0.0000	1.0000
JC69IG	37	-1718.1963	3514.9874	0.0000	1.0000
F81	38	-1718.7973	3518.4456	0.0000	1.0000
JC69	35	-1732.9888	3540.0818	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1678.8800	3437.7600	0.3996	0.3996
HKYI	40	-1679.4145	3438.8291	0.2341	0.6337
HKYIG	41	-1678.7134	3439.4268	0.1736	0.8073
GTRG	44	-1676.3193	3440.6387	0.0947	0.9021
GTRI	44	-1676.8136	3441.6272	0.0578	0.9599
GTRIG	45	-1676.1781	3442.3562	0.0401	1.0000
K2PG	37	-1693.0422	3460.0844	0.0000	1.0000
K2PI	37	-1693.5104	3461.0207	0.0000	1.0000
K2PIG	38	-1692.8677	3461.7353	0.0000	1.0000
HKY	39	-1693.1228	3464.2455	0.0000	1.0000
SYMG	41	-1691.6949	3465.3899	0.0000	1.0000
SYMI	41	-1692.1168	3466.2336	0.0000	1.0000
GTR	43	-1690.3389	3466.6778	0.0000	1.0000
SYMIG	42	-1691.5509	3467.1017	0.0000	1.0000
F81G	39	-1704.7042	3487.4083	0.0000	1.0000
F81I	39	-1705.0544	3488.1089	0.0000	1.0000
K2P	36	-1708.2398	3488.4796	0.0000	1.0000
F81IG	40	-1704.6308	3489.2617	0.0000	1.0000
SYM	40	-1706.6044	3493.2089	0.0000	1.0000
JC69G	36	-1718.2892	3508.5785	0.0000	1.0000
JC69I	36	-1718.6538	3509.3077	0.0000	1.0000
JC69IG	37	-1718.1963	3510.3926	0.0000	1.0000
F81	38	-1718.7973	3513.5946	0.0000	1.0000
JC69	35	-1732.9888	3535.9776	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1678.8800	3616.8389	0.6054	0.6054
HKYI	40	-1679.4145	3617.9080	0.3547	0.9600
HKYIG	41	-1678.7134	3622.9826	0.0280	0.9881
K2PG	37	-1693.0422	3625.7323	0.0071	0.9952
K2PI	37	-1693.5104	3626.6687	0.0044	0.9996
K2PIG	38	-1692.8677	3631.8603	0.0003	1.0000
GTRG	44	-1676.3193	3637.6255	0.0000	1.0000
GTRI	44	-1676.8136	3638.6140	0.0000	1.0000
HKY	39	-1693.1228	3638.8474	0.0000	1.0000
GTRIG	45	-1676.1781	3643.8200	0.0000	1.0000
SYMG	41	-1691.6949	3648.9457	0.0000	1.0000
K2P	36	-1708.2398	3649.6506	0.0000	1.0000
SYMI	41	-1692.1168	3649.7895	0.0000	1.0000
SYMIG	42	-1691.5509	3655.1346	0.0000	1.0000
GTR	43	-1690.3389	3659.1876	0.0000	1.0000
F81G	39	-1704.7042	3662.0102	0.0000	1.0000
F81I	39	-1705.0544	3662.7108	0.0000	1.0000
F81IG	40	-1704.6308	3668.3406	0.0000	1.0000
JC69G	36	-1718.2892	3669.7495	0.0000	1.0000
JC69I	36	-1718.6538	3670.4787	0.0000	1.0000
SYM	40	-1706.6044	3672.2878	0.0000	1.0000
JC69IG	37	-1718.1963	3676.0406	0.0000	1.0000
F81	38	-1718.7973	3683.7195	0.0000	1.0000
JC69	35	-1732.9888	3692.6716	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/chrUn_random_7197.phylip = 1 - 650;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrUn_random_7197.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/chrUn_random_7197.phylip = 1 - 650;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrUn_random_7197.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/chrUn_random_7197.phylip = 1 - 650;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrUn_random_7197.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

