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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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
GTRG	44	-1614.0968	3322.7937	0.4604	0.4604
GTRIG	45	-1613.2260	3323.3634	0.3463	0.8067
GTRI	44	-1615.1410	3324.8819	0.1621	0.9688
SYMG	41	-1621.0071	3329.7256	0.0144	0.9832
SYMIG	42	-1620.1575	3330.3150	0.0107	0.9939
SYMI	41	-1621.9450	3331.6015	0.0056	0.9996
HKYG	40	-1626.6386	3338.7076	0.0002	0.9997
HKYIG	41	-1625.8788	3339.4691	0.0001	0.9998
HKYI	40	-1627.5318	3340.4940	0.0001	0.9999
K2PG	37	-1631.1648	3340.9623	0.0001	0.9999
K2PIG	38	-1630.3863	3341.6638	0.0000	1.0000
K2PI	37	-1632.0629	3342.7585	0.0000	1.0000
F81G	39	-1637.1735	3357.5040	0.0000	1.0000
F81I	39	-1637.6057	3358.3685	0.0000	1.0000
F81IG	40	-1636.8061	3359.0426	0.0000	1.0000
JC69I	36	-1641.5414	3359.4644	0.0000	1.0000
JC69G	36	-1644.1427	3364.6671	0.0000	1.0000
GTR	43	-1636.2234	3364.7430	0.0000	1.0000
JC69IG	37	-1643.4122	3365.4569	0.0000	1.0000
SYM	40	-1642.8963	3371.2230	0.0000	1.0000
HKY	39	-1647.1671	3377.4913	0.0000	1.0000
K2P	36	-1651.6997	3379.7811	0.0000	1.0000
F81	38	-1655.7921	3392.4752	0.0000	1.0000
JC69	35	-1659.8436	3393.8251	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1614.0968	3316.1937	0.4393	0.4393
GTRIG	45	-1613.2260	3316.4519	0.3861	0.8253
GTRI	44	-1615.1410	3318.2819	0.1546	0.9800
SYMG	41	-1621.0071	3324.0142	0.0088	0.9888
SYMIG	42	-1620.1575	3324.3150	0.0076	0.9963
SYMI	41	-1621.9450	3325.8901	0.0034	0.9998
HKYG	40	-1626.6386	3333.2771	0.0001	0.9999
HKYIG	41	-1625.8788	3333.7576	0.0001	0.9999
HKYI	40	-1627.5318	3335.0636	0.0000	1.0000
K2PG	37	-1631.1648	3336.3297	0.0000	1.0000
K2PIG	38	-1630.3863	3336.7727	0.0000	1.0000
K2PI	37	-1632.0629	3338.1258	0.0000	1.0000
F81G	39	-1637.1735	3352.3470	0.0000	1.0000
F81I	39	-1637.6057	3353.2114	0.0000	1.0000
F81IG	40	-1636.8061	3353.6122	0.0000	1.0000
JC69I	36	-1641.5414	3355.0829	0.0000	1.0000
GTR	43	-1636.2234	3358.4468	0.0000	1.0000
JC69G	36	-1644.1427	3360.2855	0.0000	1.0000
JC69IG	37	-1643.4122	3360.8243	0.0000	1.0000
SYM	40	-1642.8963	3365.7926	0.0000	1.0000
HKY	39	-1647.1671	3372.3343	0.0000	1.0000
K2P	36	-1651.6997	3375.3995	0.0000	1.0000
F81	38	-1655.7921	3387.5841	0.0000	1.0000
JC69	35	-1659.8436	3389.6871	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1631.1648	3501.6919	0.6251	0.6251
K2PI	37	-1632.0629	3503.4881	0.2546	0.8797
K2PIG	38	-1630.3863	3506.6042	0.0536	0.9333
SYMG	41	-1621.0071	3507.2534	0.0387	0.9720
SYMI	41	-1621.9450	3509.1293	0.0152	0.9872
SYMIG	42	-1620.1575	3512.0235	0.0036	0.9908
HKYG	40	-1626.6386	3512.0471	0.0035	0.9943
GTRG	44	-1614.0968	3512.8407	0.0024	0.9967
HKYI	40	-1627.5318	3513.8336	0.0014	0.9981
GTRI	44	-1615.1410	3514.9289	0.0008	0.9989
JC69I	36	-1641.5414	3515.9759	0.0005	0.9994
HKYIG	41	-1625.8788	3516.9969	0.0003	0.9997
GTRIG	45	-1613.2260	3517.5682	0.0002	1.0000
JC69G	36	-1644.1427	3521.1785	0.0000	1.0000
JC69IG	37	-1643.4122	3526.1866	0.0000	1.0000
F81G	39	-1637.1735	3526.6478	0.0000	1.0000
F81I	39	-1637.6057	3527.5122	0.0000	1.0000
F81IG	40	-1636.8061	3532.3822	0.0000	1.0000
K2P	36	-1651.6997	3536.2925	0.0000	1.0000
SYM	40	-1642.8963	3544.5626	0.0000	1.0000
JC69	35	-1659.8436	3546.1109	0.0000	1.0000
HKY	39	-1647.1671	3546.6350	0.0000	1.0000
GTR	43	-1636.2234	3550.6246	0.0000	1.0000
F81	38	-1655.7921	3557.4156	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2984.phylip = 1 - 645;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2984.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2984.phylip = 1 - 645;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2984.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2984.phylip = 1 - 645;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr17_2984.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

