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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-3613.1179	7321.0517	0.7285	0.7285
GTRIG	45	-3613.0530	7323.2439	0.2434	0.9719
HKYG	40	-3621.2568	7328.1205	0.0213	0.9931
HKYIG	41	-3621.2439	7330.3851	0.0068	1.0000
GTRI	44	-3623.9112	7342.6382	0.0000	1.0000
K2PG	37	-3634.9510	7348.6843	0.0000	1.0000
HKYI	40	-3632.3661	7350.3389	0.0000	1.0000
K2PIG	38	-3634.9482	7350.9459	0.0000	1.0000
SYMG	41	-3634.0754	7356.0480	0.0000	1.0000
SYMIG	42	-3634.0548	7358.3052	0.0000	1.0000
K2PI	37	-3647.1533	7373.0890	0.0000	1.0000
SYMI	41	-3645.5582	7379.0136	0.0000	1.0000
GTR	43	-3654.1079	7400.7176	0.0000	1.0000
HKY	39	-3666.1679	7415.6601	0.0000	1.0000
K2P	36	-3684.6178	7445.7585	0.0000	1.0000
SYM	40	-3681.0116	7447.6300	0.0000	1.0000
F81G	39	-3688.8854	7461.0950	0.0000	1.0000
F81IG	40	-3688.7886	7463.1841	0.0000	1.0000
JC69G	36	-3697.9815	7472.4860	0.0000	1.0000
JC69IG	37	-3697.9121	7474.6065	0.0000	1.0000
F81I	39	-3698.1235	7479.5711	0.0000	1.0000
JC69I	36	-3707.9313	7492.3856	0.0000	1.0000
F81	38	-3731.3800	7543.8094	0.0000	1.0000
JC69	35	-3743.0249	7560.3211	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3613.1179	7314.2359	0.7070	0.7070
GTRIG	45	-3613.0530	7316.1060	0.2775	0.9845
HKYG	40	-3621.2568	7322.5137	0.0113	0.9958
HKYIG	41	-3621.2439	7324.4879	0.0042	1.0000
GTRI	44	-3623.9112	7335.8224	0.0000	1.0000
K2PG	37	-3634.9510	7343.9019	0.0000	1.0000
HKYI	40	-3632.3661	7344.7321	0.0000	1.0000
K2PIG	38	-3634.9482	7345.8965	0.0000	1.0000
SYMG	41	-3634.0754	7350.1507	0.0000	1.0000
SYMIG	42	-3634.0548	7352.1096	0.0000	1.0000
K2PI	37	-3647.1533	7368.3066	0.0000	1.0000
SYMI	41	-3645.5582	7373.1163	0.0000	1.0000
GTR	43	-3654.1079	7394.2158	0.0000	1.0000
HKY	39	-3666.1679	7410.3359	0.0000	1.0000
K2P	36	-3684.6178	7441.2355	0.0000	1.0000
SYM	40	-3681.0116	7442.0232	0.0000	1.0000
F81G	39	-3688.8854	7455.7708	0.0000	1.0000
F81IG	40	-3688.7886	7457.5772	0.0000	1.0000
JC69G	36	-3697.9815	7467.9631	0.0000	1.0000
JC69IG	37	-3697.9121	7469.8242	0.0000	1.0000
F81I	39	-3698.1235	7474.2469	0.0000	1.0000
JC69I	36	-3707.9313	7487.8627	0.0000	1.0000
F81	38	-3731.3800	7538.7600	0.0000	1.0000
JC69	35	-3743.0249	7556.0499	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3621.2568	7500.0877	0.9363	0.9363
HKYIG	41	-3621.2439	7506.5012	0.0379	0.9742
K2PG	37	-3634.9510	7508.1579	0.0166	0.9908
GTRG	44	-3613.1179	7509.5673	0.0082	0.9990
K2PIG	38	-3634.9482	7514.5918	0.0007	0.9996
GTRIG	45	-3613.0530	7515.8768	0.0003	1.0000
HKYI	40	-3632.3661	7522.3061	0.0000	1.0000
GTRI	44	-3623.9112	7531.1538	0.0000	1.0000
SYMG	41	-3634.0754	7532.1641	0.0000	1.0000
K2PI	37	-3647.1533	7532.5626	0.0000	1.0000
SYMIG	42	-3634.0548	7538.5624	0.0000	1.0000
SYMI	41	-3645.5582	7555.1297	0.0000	1.0000
HKY	39	-3666.1679	7583.4705	0.0000	1.0000
GTR	43	-3654.1079	7585.1079	0.0000	1.0000
K2P	36	-3684.6178	7601.0522	0.0000	1.0000
SYM	40	-3681.0116	7619.5972	0.0000	1.0000
JC69G	36	-3697.9815	7627.7797	0.0000	1.0000
F81G	39	-3688.8854	7628.9055	0.0000	1.0000
JC69IG	37	-3697.9121	7634.0801	0.0000	1.0000
F81IG	40	-3688.7886	7635.1513	0.0000	1.0000
F81I	39	-3698.1235	7647.3816	0.0000	1.0000
JC69I	36	-3707.9313	7647.6793	0.0000	1.0000
F81	38	-3731.3800	7707.4553	0.0000	1.0000
JC69	35	-3743.0249	7711.4271	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/chr2_22628.phylip = 1 - 626;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22628.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/chr2_22628.phylip = 1 - 626;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22628.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 (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22628.phylip = 1 - 626;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22628.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

