
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:25 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/chr4_19912.phylip" (19 taxa, 692 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	44	-1743.9528	3582.0261	0.4861	0.4861
GTRG	44	-1744.1741	3582.4687	0.3896	0.8757
GTRIG	45	-1744.1738	3584.7562	0.1241	0.9998
HKYG	40	-1757.0490	3599.1363	0.0001	0.9999
HKYI	40	-1757.0749	3599.1883	0.0001	1.0000
HKYIG	41	-1757.0659	3601.4302	0.0000	1.0000
GTR	43	-1760.2272	3612.2939	0.0000	1.0000
HKY	39	-1772.3809	3627.5470	0.0000	1.0000
SYMG	41	-1779.2108	3645.7200	0.0000	1.0000
SYMI	41	-1779.3372	3645.9728	0.0000	1.0000
SYMIG	42	-1779.2051	3647.9756	0.0000	1.0000
F81G	39	-1782.6677	3648.1207	0.0000	1.0000
F81I	39	-1782.7981	3648.3814	0.0000	1.0000
K2PG	37	-1785.9510	3650.2016	0.0000	1.0000
F81IG	40	-1782.6580	3650.3544	0.0000	1.0000
K2PI	37	-1786.1945	3650.6886	0.0000	1.0000
K2PIG	38	-1785.8908	3652.3207	0.0000	1.0000
F81	38	-1797.0270	3674.5930	0.0000	1.0000
SYM	40	-1795.0657	3675.1697	0.0000	1.0000
K2P	36	-1801.8830	3679.8331	0.0000	1.0000
JC69G	36	-1808.4373	3692.9417	0.0000	1.0000
JC69I	36	-1808.6647	3693.3966	0.0000	1.0000
JC69IG	37	-1808.3933	3695.0863	0.0000	1.0000
JC69	35	-1823.1499	3720.1412	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	44	-1743.9528	3575.9056	0.4769	0.4769
GTRG	44	-1744.1741	3576.3482	0.3823	0.8592
GTRIG	45	-1744.1738	3578.3475	0.1407	0.9999
HKYG	40	-1757.0490	3594.0979	0.0001	0.9999
HKYI	40	-1757.0749	3594.1499	0.0001	1.0000
HKYIG	41	-1757.0659	3596.1318	0.0000	1.0000
GTR	43	-1760.2272	3606.4544	0.0000	1.0000
HKY	39	-1772.3809	3622.7617	0.0000	1.0000
SYMG	41	-1779.2108	3640.4216	0.0000	1.0000
SYMI	41	-1779.3372	3640.6743	0.0000	1.0000
SYMIG	42	-1779.2051	3642.4102	0.0000	1.0000
F81G	39	-1782.6677	3643.3355	0.0000	1.0000
F81I	39	-1782.7981	3643.5962	0.0000	1.0000
F81IG	40	-1782.6580	3645.3160	0.0000	1.0000
K2PG	37	-1785.9510	3645.9019	0.0000	1.0000
K2PI	37	-1786.1945	3646.3889	0.0000	1.0000
K2PIG	38	-1785.8908	3647.7816	0.0000	1.0000
F81	38	-1797.0270	3670.0540	0.0000	1.0000
SYM	40	-1795.0657	3670.1313	0.0000	1.0000
K2P	36	-1801.8830	3675.7659	0.0000	1.0000
JC69G	36	-1808.4373	3688.8745	0.0000	1.0000
JC69I	36	-1808.6647	3689.3294	0.0000	1.0000
JC69IG	37	-1808.3933	3690.7866	0.0000	1.0000
JC69	35	-1823.1499	3716.2998	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	44	-1743.9528	3775.6474	0.2625	0.2625
HKYG	40	-1757.0490	3775.6813	0.2581	0.5205
HKYI	40	-1757.0749	3775.7333	0.2514	0.7720
GTRG	44	-1744.1741	3776.0900	0.2104	0.9824
HKYIG	41	-1757.0659	3782.2548	0.0096	0.9920
GTRIG	45	-1744.1738	3782.6289	0.0080	1.0000
HKY	39	-1772.3809	3799.8056	0.0000	1.0000
GTR	43	-1760.2272	3801.6566	0.0000	1.0000
K2PG	37	-1785.9510	3813.8666	0.0000	1.0000
K2PI	37	-1786.1945	3814.3536	0.0000	1.0000
K2PIG	38	-1785.8908	3820.2859	0.0000	1.0000
F81G	39	-1782.6677	3820.3793	0.0000	1.0000
F81I	39	-1782.7981	3820.6400	0.0000	1.0000
SYMG	41	-1779.2108	3826.5446	0.0000	1.0000
SYMI	41	-1779.3372	3826.7973	0.0000	1.0000
F81IG	40	-1782.6580	3826.8994	0.0000	1.0000
SYMIG	42	-1779.2051	3833.0728	0.0000	1.0000
K2P	36	-1801.8830	3839.1910	0.0000	1.0000
F81	38	-1797.0270	3842.5582	0.0000	1.0000
SYM	40	-1795.0657	3851.7147	0.0000	1.0000
JC69G	36	-1808.4373	3852.2996	0.0000	1.0000
JC69I	36	-1808.6647	3852.7545	0.0000	1.0000
JC69IG	37	-1808.3933	3858.7513	0.0000	1.0000
JC69	35	-1823.1499	3875.1853	0.0000	1.0000



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

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

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


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

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


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

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


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

