
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:07 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/chr8_5177.phylip" (19 taxa, 380 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	-1476.2965	3042.2686	0.6655	0.6655
HKYIG	41	-1476.2963	3044.7820	0.1894	0.8549
GTRG	44	-1473.4665	3046.7538	0.0707	0.9256
HKYI	40	-1478.8878	3047.4510	0.0499	0.9754
GTRIG	45	-1473.4646	3049.3245	0.0195	0.9950
GTRI	44	-1476.1169	3052.0548	0.0050	1.0000
HKY	39	-1487.5479	3062.2723	0.0000	1.0000
GTR	43	-1484.7303	3066.7225	0.0000	1.0000
K2PG	37	-1495.1960	3072.6142	0.0000	1.0000
K2PIG	38	-1495.1942	3075.0804	0.0000	1.0000
K2PI	37	-1497.8285	3077.8792	0.0000	1.0000
SYMG	41	-1493.6757	3079.5407	0.0000	1.0000
SYMIG	42	-1493.6697	3082.0574	0.0000	1.0000
SYMI	41	-1496.1349	3084.4591	0.0000	1.0000
F81G	39	-1502.7447	3092.6659	0.0000	1.0000
K2P	36	-1507.1510	3094.0687	0.0000	1.0000
F81IG	40	-1502.7450	3095.1655	0.0000	1.0000
F81I	39	-1505.1619	3097.5002	0.0000	1.0000
SYM	40	-1505.4520	3100.5795	0.0000	1.0000
F81	38	-1514.2631	3113.2184	0.0000	1.0000
JC69G	36	-1520.3362	3120.4391	0.0000	1.0000
JC69IG	37	-1520.3344	3122.8909	0.0000	1.0000
JC69I	36	-1522.6845	3125.1357	0.0000	1.0000
JC69	35	-1532.0754	3141.4764	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1476.2965	3032.5931	0.5292	0.5292
HKYIG	41	-1476.2963	3034.5926	0.1947	0.7239
GTRG	44	-1473.4665	3034.9329	0.1643	0.8882
GTRIG	45	-1473.4646	3036.9293	0.0605	0.9487
HKYI	40	-1478.8878	3037.7755	0.0397	0.9884
GTRI	44	-1476.1169	3040.2339	0.0116	1.0000
HKY	39	-1487.5479	3053.0959	0.0000	1.0000
GTR	43	-1484.7303	3055.4606	0.0000	1.0000
K2PG	37	-1495.1960	3064.3920	0.0000	1.0000
K2PIG	38	-1495.1942	3066.3883	0.0000	1.0000
SYMG	41	-1493.6757	3069.3514	0.0000	1.0000
K2PI	37	-1497.8285	3069.6570	0.0000	1.0000
SYMIG	42	-1493.6697	3071.3393	0.0000	1.0000
SYMI	41	-1496.1349	3074.2698	0.0000	1.0000
F81G	39	-1502.7447	3083.4894	0.0000	1.0000
F81IG	40	-1502.7450	3085.4899	0.0000	1.0000
K2P	36	-1507.1510	3086.3020	0.0000	1.0000
F81I	39	-1505.1619	3088.3237	0.0000	1.0000
SYM	40	-1505.4520	3090.9040	0.0000	1.0000
F81	38	-1514.2631	3104.5263	0.0000	1.0000
JC69G	36	-1520.3362	3112.6723	0.0000	1.0000
JC69IG	37	-1520.3344	3114.6687	0.0000	1.0000
JC69I	36	-1522.6845	3117.3690	0.0000	1.0000
JC69	35	-1532.0754	3134.1508	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1476.2965	3190.1999	0.8876	0.8876
HKYI	40	-1478.8878	3195.3824	0.0665	0.9541
HKYIG	41	-1476.2963	3196.1396	0.0455	0.9996
HKY	39	-1487.5479	3206.7625	0.0002	0.9998
GTRG	44	-1473.4665	3208.3005	0.0001	0.9999
K2PG	37	-1495.1960	3210.1783	0.0000	1.0000
GTRI	44	-1476.1169	3213.6014	0.0000	1.0000
GTRIG	45	-1473.4646	3214.2370	0.0000	1.0000
K2PI	37	-1497.8285	3215.4433	0.0000	1.0000
K2PIG	38	-1495.1942	3216.1148	0.0000	1.0000
GTR	43	-1484.7303	3224.8880	0.0000	1.0000
K2P	36	-1507.1510	3228.1481	0.0000	1.0000
SYMG	41	-1493.6757	3230.8984	0.0000	1.0000
SYMI	41	-1496.1349	3235.8168	0.0000	1.0000
SYMIG	42	-1493.6697	3236.8265	0.0000	1.0000
F81G	39	-1502.7447	3237.1561	0.0000	1.0000
F81I	39	-1505.1619	3241.9904	0.0000	1.0000
F81IG	40	-1502.7450	3243.0968	0.0000	1.0000
SYM	40	-1505.4520	3248.5108	0.0000	1.0000
F81	38	-1514.2631	3254.2528	0.0000	1.0000
JC69G	36	-1520.3362	3254.5185	0.0000	1.0000
JC69I	36	-1522.6845	3259.2151	0.0000	1.0000
JC69IG	37	-1520.3344	3260.4550	0.0000	1.0000
JC69	35	-1532.0754	3272.0568	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/chr8_5177.phylip = 1 - 380;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_5177.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/chr8_5177.phylip = 1 - 380;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_5177.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/chr8_5177.phylip = 1 - 380;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_5177.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

