
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:15 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/chr1_5438.phylip" (19 taxa, 437 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	-1148.6752	2385.6333	0.8443	0.8443
HKYI	40	-1151.3182	2390.9192	0.0601	0.9044
GTRIG	45	-1145.5602	2391.7087	0.0405	0.9449
HKYIG	41	-1150.5370	2391.7930	0.0388	0.9837
GTRI	44	-1148.4021	2394.9062	0.0082	0.9919
GTRG	44	-1148.4117	2394.9255	0.0081	1.0000
F81I	39	-1166.0088	2417.8766	0.0000	1.0000
F81G	39	-1166.0838	2418.0265	0.0000	1.0000
F81IG	40	-1165.4710	2419.2249	0.0000	1.0000
K2PI	37	-1172.1055	2425.2585	0.0000	1.0000
K2PG	37	-1172.1928	2425.4333	0.0000	1.0000
K2PIG	38	-1171.4648	2426.3769	0.0000	1.0000
SYMI	41	-1171.5081	2433.7353	0.0000	1.0000
SYMG	41	-1171.5990	2433.9169	0.0000	1.0000
SYMIG	42	-1170.9051	2434.9776	0.0000	1.0000
HKY	39	-1180.5976	2447.0542	0.0000	1.0000
JC69I	36	-1184.2365	2447.1329	0.0000	1.0000
JC69G	36	-1184.3251	2447.3103	0.0000	1.0000
JC69IG	37	-1183.7719	2448.5914	0.0000	1.0000
GTR	43	-1178.2314	2452.0914	0.0000	1.0000
F81	38	-1194.5865	2472.6203	0.0000	1.0000
K2P	36	-1201.2528	2481.1655	0.0000	1.0000
SYM	40	-1200.6744	2489.6317	0.0000	1.0000
JC69	35	-1212.9211	2502.1264	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1148.6752	2377.3504	0.7530	0.7530
GTRIG	45	-1145.5602	2381.1204	0.1143	0.8673
HKYI	40	-1151.3182	2382.6364	0.0536	0.9209
HKYIG	41	-1150.5370	2383.0741	0.0430	0.9639
GTRI	44	-1148.4021	2384.8042	0.0181	0.9821
GTRG	44	-1148.4117	2384.8235	0.0179	1.0000
F81I	39	-1166.0088	2410.0176	0.0000	1.0000
F81G	39	-1166.0838	2410.1676	0.0000	1.0000
F81IG	40	-1165.4710	2410.9421	0.0000	1.0000
K2PI	37	-1172.1055	2418.2109	0.0000	1.0000
K2PG	37	-1172.1928	2418.3857	0.0000	1.0000
K2PIG	38	-1171.4648	2418.9296	0.0000	1.0000
SYMI	41	-1171.5081	2425.0163	0.0000	1.0000
SYMG	41	-1171.5990	2425.1980	0.0000	1.0000
SYMIG	42	-1170.9051	2425.8101	0.0000	1.0000
HKY	39	-1180.5976	2439.1953	0.0000	1.0000
JC69I	36	-1184.2365	2440.4729	0.0000	1.0000
JC69G	36	-1184.3251	2440.6503	0.0000	1.0000
JC69IG	37	-1183.7719	2441.5437	0.0000	1.0000
GTR	43	-1178.2314	2442.4629	0.0000	1.0000
F81	38	-1194.5865	2465.1730	0.0000	1.0000
K2P	36	-1201.2528	2474.5055	0.0000	1.0000
SYM	40	-1200.6744	2481.3489	0.0000	1.0000
JC69	35	-1212.9211	2495.8422	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1148.6752	2540.5478	0.9271	0.9271
HKYI	40	-1151.3182	2545.8337	0.0660	0.9931
HKYIG	41	-1150.5370	2550.3513	0.0069	1.0000
GTRI	44	-1148.4021	2564.3212	0.0000	1.0000
GTRG	44	-1148.4117	2564.3405	0.0000	1.0000
GTRIG	45	-1145.5602	2564.7174	0.0000	1.0000
F81I	39	-1166.0088	2569.1350	0.0000	1.0000
K2PI	37	-1172.1055	2569.1684	0.0000	1.0000
F81G	39	-1166.0838	2569.2850	0.0000	1.0000
K2PG	37	-1172.1928	2569.3432	0.0000	1.0000
K2PIG	38	-1171.4648	2573.9671	0.0000	1.0000
F81IG	40	-1165.4710	2574.1394	0.0000	1.0000
JC69I	36	-1184.2365	2587.3505	0.0000	1.0000
JC69G	36	-1184.3251	2587.5279	0.0000	1.0000
SYMI	41	-1171.5081	2592.2936	0.0000	1.0000
SYMG	41	-1171.5990	2592.4752	0.0000	1.0000
JC69IG	37	-1183.7719	2592.5013	0.0000	1.0000
SYMIG	42	-1170.9051	2597.1673	0.0000	1.0000
HKY	39	-1180.5976	2598.3127	0.0000	1.0000
GTR	43	-1178.2314	2617.9000	0.0000	1.0000
F81	38	-1194.5865	2620.2105	0.0000	1.0000
K2P	36	-1201.2528	2621.3831	0.0000	1.0000
JC69	35	-1212.9211	2638.6398	0.0000	1.0000
SYM	40	-1200.6744	2644.5462	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/chr1_5438.phylip = 1 - 437;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_5438.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/chr1_5438.phylip = 1 - 437;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_5438.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/chr1_5438.phylip = 1 - 437;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_5438.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

