
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:01 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/chr7_10866.phylip" (19 taxa, 593 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	-1655.5852	3397.1124	0.5191	0.5191
GTRI	44	-1652.0456	3399.3174	0.1724	0.6915
HKYIG	41	-1655.6535	3399.5575	0.1529	0.8444
HKYG	40	-1657.5729	3401.0878	0.0711	0.9155
GTRIG	45	-1651.7919	3401.1524	0.0689	0.9844
GTRG	44	-1654.4476	3404.1215	0.0156	1.0000
HKY	39	-1678.0724	3439.7868	0.0000	1.0000
GTR	43	-1673.9995	3440.8915	0.0000	1.0000
K2PI	37	-1688.5545	3456.1756	0.0000	1.0000
K2PG	37	-1689.5873	3458.2413	0.0000	1.0000
SYMG	41	-1685.8279	3459.9062	0.0000	1.0000
SYMI	41	-1685.8477	3459.9458	0.0000	1.0000
K2PIG	38	-1689.6273	3460.6048	0.0000	1.0000
SYMIG	42	-1685.8364	3462.2401	0.0000	1.0000
F81I	39	-1702.9171	3489.4763	0.0000	1.0000
F81IG	40	-1703.1033	3492.1487	0.0000	1.0000
F81G	39	-1704.4409	3492.5237	0.0000	1.0000
SYM	40	-1705.9171	3497.7762	0.0000	1.0000
K2P	36	-1710.5079	3497.8071	0.0000	1.0000
F81	38	-1723.7243	3528.7987	0.0000	1.0000
JC69G	36	-1732.6070	3542.0053	0.0000	1.0000
JC69I	36	-1735.6091	3548.0095	0.0000	1.0000
JC69IG	37	-1735.0280	3549.1228	0.0000	1.0000
JC69	35	-1751.8816	3578.2874	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1655.5852	3391.1704	0.4052	0.4052
GTRI	44	-1652.0456	3392.0911	0.2557	0.6609
HKYIG	41	-1655.6535	3393.3070	0.1392	0.8001
GTRIG	45	-1651.7919	3393.5839	0.1212	0.9213
HKYG	40	-1657.5729	3395.1458	0.0555	0.9769
GTRG	44	-1654.4476	3396.8952	0.0231	1.0000
GTR	43	-1673.9995	3433.9989	0.0000	1.0000
HKY	39	-1678.0724	3434.1448	0.0000	1.0000
K2PI	37	-1688.5545	3451.1089	0.0000	1.0000
K2PG	37	-1689.5873	3453.1746	0.0000	1.0000
SYMG	41	-1685.8279	3453.6557	0.0000	1.0000
SYMI	41	-1685.8477	3453.6953	0.0000	1.0000
K2PIG	38	-1689.6273	3455.2546	0.0000	1.0000
SYMIG	42	-1685.8364	3455.6728	0.0000	1.0000
F81I	39	-1702.9171	3483.8343	0.0000	1.0000
F81IG	40	-1703.1033	3486.2067	0.0000	1.0000
F81G	39	-1704.4409	3486.8817	0.0000	1.0000
SYM	40	-1705.9171	3491.8341	0.0000	1.0000
K2P	36	-1710.5079	3493.0158	0.0000	1.0000
F81	38	-1723.7243	3523.4485	0.0000	1.0000
JC69G	36	-1732.6070	3537.2140	0.0000	1.0000
JC69I	36	-1735.6091	3543.2181	0.0000	1.0000
JC69IG	37	-1735.0280	3544.0561	0.0000	1.0000
JC69	35	-1751.8816	3573.7632	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1655.5852	3566.5781	0.8507	0.8507
HKYG	40	-1657.5729	3570.5536	0.1166	0.9673
HKYIG	41	-1655.6535	3573.1000	0.0326	0.9999
GTRI	44	-1652.0456	3585.0397	0.0001	1.0000
GTRG	44	-1654.4476	3589.8438	0.0000	1.0000
GTRIG	45	-1651.7919	3590.9176	0.0000	1.0000
HKY	39	-1678.0724	3605.1674	0.0000	1.0000
K2PI	37	-1688.5545	3613.3611	0.0000	1.0000
K2PG	37	-1689.5873	3615.4268	0.0000	1.0000
K2PIG	38	-1689.6273	3621.8920	0.0000	1.0000
GTR	43	-1673.9995	3622.5623	0.0000	1.0000
SYMG	41	-1685.8279	3633.4487	0.0000	1.0000
SYMI	41	-1685.8477	3633.4883	0.0000	1.0000
SYMIG	42	-1685.8364	3639.8510	0.0000	1.0000
K2P	36	-1710.5079	3650.8828	0.0000	1.0000
F81I	39	-1702.9171	3654.8569	0.0000	1.0000
F81G	39	-1704.4409	3657.9043	0.0000	1.0000
F81IG	40	-1703.1033	3661.6144	0.0000	1.0000
SYM	40	-1705.9171	3667.2419	0.0000	1.0000
F81	38	-1723.7243	3690.0859	0.0000	1.0000
JC69G	36	-1732.6070	3695.0810	0.0000	1.0000
JC69I	36	-1735.6091	3701.0851	0.0000	1.0000
JC69IG	37	-1735.0280	3706.3083	0.0000	1.0000
JC69	35	-1751.8816	3727.2450	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/chr7_10866.phylip = 1 - 593;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10866.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/chr7_10866.phylip = 1 - 593;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10866.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/chr7_10866.phylip = 1 - 593;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10866.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

