
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:34 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_6584.phylip" (19 taxa, 643 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	-3111.7002	6308.8488	0.5619	0.5619
HKYIG	41	-3111.2668	6310.2640	0.2769	0.8388
GTRG	44	-3109.1770	6312.9761	0.0714	0.9102
HKYI	40	-3114.2675	6313.9835	0.0431	0.9533
GTRIG	45	-3108.6570	6314.2487	0.0378	0.9911
GTRI	44	-3111.2586	6317.1393	0.0089	1.0000
K2PG	37	-3127.9010	6334.4500	0.0000	1.0000
K2PIG	38	-3127.7168	6336.3409	0.0000	1.0000
SYMG	41	-3124.9114	6337.5532	0.0000	1.0000
SYMIG	42	-3124.7717	6339.5633	0.0000	1.0000
K2PI	37	-3131.4336	6341.5152	0.0000	1.0000
SYMI	41	-3128.5034	6344.7372	0.0000	1.0000
HKY	39	-3146.7376	6376.6493	0.0000	1.0000
GTR	43	-3145.1693	6382.6559	0.0000	1.0000
K2P	36	-3160.9650	6398.3260	0.0000	1.0000
SYM	40	-3158.7454	6402.9393	0.0000	1.0000
F81G	39	-3221.2719	6525.7179	0.0000	1.0000
F81IG	40	-3220.9745	6527.3975	0.0000	1.0000
F81I	39	-3223.9232	6531.0206	0.0000	1.0000
JC69G	36	-3231.4842	6539.3644	0.0000	1.0000
JC69IG	37	-3231.2385	6541.1249	0.0000	1.0000
JC69I	36	-3234.2827	6544.9614	0.0000	1.0000
F81	38	-3251.5554	6584.0180	0.0000	1.0000
JC69	35	-3261.2664	6596.6844	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-3111.7002	6303.4003	0.4896	0.4896
HKYIG	41	-3111.2668	6304.5336	0.2778	0.7675
GTRG	44	-3109.1770	6306.3540	0.1118	0.8793
GTRIG	45	-3108.6570	6307.3140	0.0692	0.9485
HKYI	40	-3114.2675	6308.5350	0.0376	0.9861
GTRI	44	-3111.2586	6310.5173	0.0139	1.0000
K2PG	37	-3127.9010	6329.8020	0.0000	1.0000
K2PIG	38	-3127.7168	6331.4336	0.0000	1.0000
SYMG	41	-3124.9114	6331.8227	0.0000	1.0000
SYMIG	42	-3124.7717	6333.5433	0.0000	1.0000
K2PI	37	-3131.4336	6336.8673	0.0000	1.0000
SYMI	41	-3128.5034	6339.0068	0.0000	1.0000
HKY	39	-3146.7376	6371.4751	0.0000	1.0000
GTR	43	-3145.1693	6376.3387	0.0000	1.0000
K2P	36	-3160.9650	6393.9300	0.0000	1.0000
SYM	40	-3158.7454	6397.4908	0.0000	1.0000
F81G	39	-3221.2719	6520.5438	0.0000	1.0000
F81IG	40	-3220.9745	6521.9490	0.0000	1.0000
F81I	39	-3223.9232	6525.8464	0.0000	1.0000
JC69G	36	-3231.4842	6534.9684	0.0000	1.0000
JC69IG	37	-3231.2385	6536.4770	0.0000	1.0000
JC69I	36	-3234.2827	6540.5654	0.0000	1.0000
F81	38	-3251.5554	6579.1107	0.0000	1.0000
JC69	35	-3261.2664	6592.5328	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3111.7002	6482.0461	0.8778	0.8778
HKYI	40	-3114.2675	6487.1808	0.0674	0.9452
HKYIG	41	-3111.2668	6487.6455	0.0534	0.9985
K2PG	37	-3127.9010	6495.0494	0.0013	0.9999
K2PIG	38	-3127.7168	6501.1471	0.0001	0.9999
K2PI	37	-3131.4336	6502.1146	0.0000	1.0000
GTRG	44	-3109.1770	6502.8644	0.0000	1.0000
GTRI	44	-3111.2586	6507.0276	0.0000	1.0000
GTRIG	45	-3108.6570	6508.2905	0.0000	1.0000
SYMG	41	-3124.9114	6514.9347	0.0000	1.0000
SYMIG	42	-3124.7717	6521.1214	0.0000	1.0000
SYMI	41	-3128.5034	6522.1187	0.0000	1.0000
HKY	39	-3146.7376	6545.6548	0.0000	1.0000
K2P	36	-3160.9650	6554.7112	0.0000	1.0000
GTR	43	-3145.1693	6568.3829	0.0000	1.0000
SYM	40	-3158.7454	6576.1366	0.0000	1.0000
F81G	39	-3221.2719	6694.7234	0.0000	1.0000
JC69G	36	-3231.4842	6695.7496	0.0000	1.0000
F81I	39	-3223.9232	6700.0261	0.0000	1.0000
F81IG	40	-3220.9745	6700.5948	0.0000	1.0000
JC69I	36	-3234.2827	6701.3466	0.0000	1.0000
JC69IG	37	-3231.2385	6701.7243	0.0000	1.0000
F81	38	-3251.5554	6748.8242	0.0000	1.0000
JC69	35	-3261.2664	6748.8479	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/chrZ_6584.phylip = 1 - 643;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6584.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/chrZ_6584.phylip = 1 - 643;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6584.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/chrZ_6584.phylip = 1 - 643;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6584.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

