
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:50 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/chr7_10361.phylip" (19 taxa, 371 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	-1198.4754	2486.8902	0.4259	0.4259
HKYI	40	-1198.9706	2487.8807	0.2596	0.6855
HKYIG	41	-1198.4708	2489.4096	0.1209	0.8064
GTRG	44	-1194.6995	2489.5463	0.1129	0.9193
GTRI	44	-1195.5242	2491.1956	0.0495	0.9687
GTRIG	45	-1194.6933	2492.1251	0.0311	0.9998
K2PG	37	-1211.0483	2504.5411	0.0001	0.9999
K2PI	37	-1211.4623	2505.3690	0.0000	0.9999
HKY	39	-1209.1351	2505.6962	0.0000	1.0000
K2PIG	38	-1211.0425	2507.0128	0.0000	1.0000
GTR	43	-1205.3280	2508.2278	0.0000	1.0000
SYMG	41	-1209.5578	2511.5836	0.0000	1.0000
SYMI	41	-1209.9925	2512.4531	0.0000	1.0000
SYMIG	42	-1209.5552	2514.1226	0.0000	1.0000
K2P	36	-1221.5165	2523.0090	0.0000	1.0000
SYM	40	-1219.8206	2529.5805	0.0000	1.0000
F81G	39	-1225.0174	2537.4608	0.0000	1.0000
F81I	39	-1225.4877	2538.4014	0.0000	1.0000
F81IG	40	-1225.0146	2539.9686	0.0000	1.0000
JC69G	36	-1235.0932	2550.1625	0.0000	1.0000
JC69I	36	-1235.4894	2550.9548	0.0000	1.0000
JC69IG	37	-1235.0924	2552.6292	0.0000	1.0000
F81	38	-1234.9268	2554.7813	0.0000	1.0000
JC69	35	-1245.0947	2567.7117	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1198.4754	2476.9508	0.2920	0.2920
GTRG	44	-1194.6995	2477.3990	0.2334	0.5254
HKYI	40	-1198.9706	2477.9413	0.1780	0.7033
HKYIG	41	-1198.4708	2478.9415	0.1079	0.8112
GTRI	44	-1195.5242	2479.0484	0.1023	0.9135
GTRIG	45	-1194.6933	2479.3866	0.0864	0.9999
K2PG	37	-1211.0483	2496.0967	0.0000	0.9999
HKY	39	-1209.1351	2496.2702	0.0000	1.0000
GTR	43	-1205.3280	2496.6559	0.0000	1.0000
K2PI	37	-1211.4623	2496.9246	0.0000	1.0000
K2PIG	38	-1211.0425	2498.0850	0.0000	1.0000
SYMG	41	-1209.5578	2501.1156	0.0000	1.0000
SYMI	41	-1209.9925	2501.9851	0.0000	1.0000
SYMIG	42	-1209.5552	2503.1104	0.0000	1.0000
K2P	36	-1221.5165	2515.0330	0.0000	1.0000
SYM	40	-1219.8206	2519.6411	0.0000	1.0000
F81G	39	-1225.0174	2528.0349	0.0000	1.0000
F81I	39	-1225.4877	2528.9754	0.0000	1.0000
F81IG	40	-1225.0146	2530.0292	0.0000	1.0000
JC69G	36	-1235.0932	2542.1865	0.0000	1.0000
JC69I	36	-1235.4894	2542.9788	0.0000	1.0000
JC69IG	37	-1235.0924	2544.1848	0.0000	1.0000
F81	38	-1234.9268	2545.8536	0.0000	1.0000
JC69	35	-1245.0947	2560.1893	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1198.4754	2633.5989	0.5865	0.5865
HKYI	40	-1198.9706	2634.5893	0.3575	0.9440
HKYIG	41	-1198.4708	2639.5058	0.0306	0.9746
K2PG	37	-1211.0483	2640.9962	0.0145	0.9891
K2PI	37	-1211.4623	2641.8241	0.0096	0.9987
K2PIG	38	-1211.0425	2646.9007	0.0008	0.9994
HKY	39	-1209.1351	2649.0021	0.0003	0.9997
GTRG	44	-1194.6995	2649.7119	0.0002	0.9999
GTRI	44	-1195.5242	2651.3613	0.0001	1.0000
GTRIG	45	-1194.6933	2655.6157	0.0000	1.0000
K2P	36	-1221.5165	2656.0162	0.0000	1.0000
SYMG	41	-1209.5578	2661.6798	0.0000	1.0000
SYMI	41	-1209.9925	2662.5493	0.0000	1.0000
GTR	43	-1205.3280	2665.0526	0.0000	1.0000
SYMIG	42	-1209.5552	2667.5909	0.0000	1.0000
SYM	40	-1219.8206	2676.2892	0.0000	1.0000
F81G	39	-1225.0174	2680.7667	0.0000	1.0000
F81I	39	-1225.4877	2681.7073	0.0000	1.0000
JC69G	36	-1235.0932	2683.1698	0.0000	1.0000
JC69I	36	-1235.4894	2683.9620	0.0000	1.0000
F81IG	40	-1225.0146	2686.6773	0.0000	1.0000
JC69IG	37	-1235.0924	2689.0843	0.0000	1.0000
F81	38	-1234.9268	2694.6693	0.0000	1.0000
JC69	35	-1245.0947	2697.2564	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/chr7_10361.phylip = 1 - 371;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10361.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/chr7_10361.phylip = 1 - 371;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10361.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/chr7_10361.phylip = 1 - 371;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10361.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

