
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:09 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/chrZ_6679.phylip" (19 taxa, 607 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	-1820.9919	3727.7789	0.4978	0.4978
HKYG	40	-1821.6902	3729.1754	0.2476	0.7454
GTRI	44	-1817.9239	3730.8941	0.1049	0.8503
HKYIG	41	-1821.6901	3731.4758	0.0784	0.9287
GTRG	44	-1818.5804	3732.2070	0.0544	0.9831
GTRIG	45	-1818.5801	3734.5399	0.0169	1.0000
F81I	39	-1847.7214	3778.9454	0.0000	1.0000
F81G	39	-1848.1769	3779.8563	0.0000	1.0000
HKY	39	-1848.5403	3780.5833	0.0000	1.0000
F81IG	40	-1848.1766	3782.1483	0.0000	1.0000
GTR	43	-1845.0236	3782.7683	0.0000	1.0000
K2PI	37	-1856.4589	3791.8598	0.0000	1.0000
SYMI	41	-1852.1134	3792.3223	0.0000	1.0000
K2PG	37	-1856.8562	3792.6545	0.0000	1.0000
SYMG	41	-1852.5245	3793.1446	0.0000	1.0000
K2PIG	38	-1856.8562	3794.9307	0.0000	1.0000
SYMIG	42	-1852.5244	3795.4530	0.0000	1.0000
F81	38	-1875.7460	3832.7104	0.0000	1.0000
JC69I	36	-1879.8277	3836.3290	0.0000	1.0000
JC69G	36	-1880.1864	3837.0465	0.0000	1.0000
JC69IG	37	-1880.1863	3839.3146	0.0000	1.0000
SYM	40	-1882.1131	3850.0212	0.0000	1.0000
K2P	36	-1890.0379	3856.7496	0.0000	1.0000
JC69	35	-1908.8225	3892.0583	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1820.9919	3721.9839	0.4249	0.4249
HKYG	40	-1821.6902	3723.3804	0.2114	0.6362
GTRI	44	-1817.9239	3723.8478	0.1673	0.8035
GTRG	44	-1818.5804	3725.1608	0.0868	0.8903
HKYIG	41	-1821.6901	3725.3802	0.0778	0.9681
GTRIG	45	-1818.5801	3727.1602	0.0319	1.0000
F81I	39	-1847.7214	3773.4427	0.0000	1.0000
F81G	39	-1848.1769	3774.3537	0.0000	1.0000
HKY	39	-1848.5403	3775.0807	0.0000	1.0000
GTR	43	-1845.0236	3776.0471	0.0000	1.0000
F81IG	40	-1848.1766	3776.3533	0.0000	1.0000
SYMI	41	-1852.1134	3786.2268	0.0000	1.0000
K2PI	37	-1856.4589	3786.9178	0.0000	1.0000
SYMG	41	-1852.5245	3787.0490	0.0000	1.0000
K2PG	37	-1856.8562	3787.7125	0.0000	1.0000
SYMIG	42	-1852.5244	3789.0488	0.0000	1.0000
K2PIG	38	-1856.8562	3789.7124	0.0000	1.0000
F81	38	-1875.7460	3827.4921	0.0000	1.0000
JC69I	36	-1879.8277	3831.6553	0.0000	1.0000
JC69G	36	-1880.1864	3832.3729	0.0000	1.0000
JC69IG	37	-1880.1863	3834.3726	0.0000	1.0000
SYM	40	-1882.1131	3844.2262	0.0000	1.0000
K2P	36	-1890.0379	3852.0759	0.0000	1.0000
JC69	35	-1908.8225	3887.6450	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1820.9919	3898.3250	0.6589	0.6589
HKYG	40	-1821.6902	3899.7215	0.3278	0.9866
HKYIG	41	-1821.6901	3906.1299	0.0133	0.9999
GTRI	44	-1817.9239	3917.8231	0.0000	1.0000
GTRG	44	-1818.5804	3919.1360	0.0000	1.0000
GTRIG	45	-1818.5801	3925.5440	0.0000	1.0000
F81I	39	-1847.7214	3945.3753	0.0000	1.0000
F81G	39	-1848.1769	3946.2863	0.0000	1.0000
HKY	39	-1848.5403	3947.0133	0.0000	1.0000
K2PI	37	-1856.4589	3950.0334	0.0000	1.0000
K2PG	37	-1856.8562	3950.8280	0.0000	1.0000
F81IG	40	-1848.1766	3952.6944	0.0000	1.0000
K2PIG	38	-1856.8562	3957.2365	0.0000	1.0000
GTR	43	-1845.0236	3965.6139	0.0000	1.0000
SYMI	41	-1852.1134	3966.9764	0.0000	1.0000
SYMG	41	-1852.5245	3967.7987	0.0000	1.0000
SYMIG	42	-1852.5244	3974.2070	0.0000	1.0000
JC69I	36	-1879.8277	3990.3624	0.0000	1.0000
JC69G	36	-1880.1864	3991.0799	0.0000	1.0000
F81	38	-1875.7460	3995.0162	0.0000	1.0000
JC69IG	37	-1880.1863	3997.4882	0.0000	1.0000
K2P	36	-1890.0379	4010.7829	0.0000	1.0000
SYM	40	-1882.1131	4020.5673	0.0000	1.0000
JC69	35	-1908.8225	4041.9435	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/chrZ_6679.phylip = 1 - 607;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6679.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/chrZ_6679.phylip = 1 - 607;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6679.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/chrZ_6679.phylip = 1 - 607;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6679.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

