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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1110.6677	2321.5955	0.5977	0.5977
GTRIG	45	-1110.1347	2323.1265	0.2780	0.8757
GTRI	44	-1112.2378	2324.7356	0.1243	1.0000
HKYG	40	-1131.9907	2354.0120	0.0000	1.0000
HKYI	40	-1132.5038	2355.0382	0.0000	1.0000
HKYIG	41	-1131.7390	2356.0424	0.0000	1.0000
SYMG	41	-1133.4841	2359.5327	0.0000	1.0000
SYMI	41	-1134.3716	2361.3076	0.0000	1.0000
SYMIG	42	-1133.5159	2362.1455	0.0000	1.0000
GTR	43	-1137.0411	2371.7612	0.0000	1.0000
K2PG	37	-1148.6629	2379.8471	0.0000	1.0000
K2PI	37	-1148.9962	2380.5135	0.0000	1.0000
K2PIG	38	-1148.4404	2381.8898	0.0000	1.0000
F81G	39	-1149.5962	2386.7047	0.0000	1.0000
F81I	39	-1150.0475	2387.6071	0.0000	1.0000
F81IG	40	-1149.3725	2388.7755	0.0000	1.0000
HKY	39	-1156.5603	2400.6327	0.0000	1.0000
SYM	40	-1158.1101	2406.2507	0.0000	1.0000
JC69G	36	-1164.2960	2408.6403	0.0000	1.0000
JC69I	36	-1164.6981	2409.4445	0.0000	1.0000
JC69IG	37	-1164.1889	2410.8991	0.0000	1.0000
K2P	36	-1173.4904	2427.0292	0.0000	1.0000
F81	38	-1174.0835	2433.1760	0.0000	1.0000
JC69	35	-1188.6413	2454.8729	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1110.6677	2309.3354	0.5450	0.5450
GTRIG	45	-1110.1347	2310.2693	0.3417	0.8866
GTRI	44	-1112.2378	2312.4756	0.1134	1.0000
HKYG	40	-1131.9907	2343.9814	0.0000	1.0000
HKYI	40	-1132.5038	2345.0076	0.0000	1.0000
HKYIG	41	-1131.7390	2345.4780	0.0000	1.0000
SYMG	41	-1133.4841	2348.9682	0.0000	1.0000
SYMI	41	-1134.3716	2350.7432	0.0000	1.0000
SYMIG	42	-1133.5159	2351.0317	0.0000	1.0000
GTR	43	-1137.0411	2360.0822	0.0000	1.0000
K2PG	37	-1148.6629	2371.3259	0.0000	1.0000
K2PI	37	-1148.9962	2371.9923	0.0000	1.0000
K2PIG	38	-1148.4404	2372.8807	0.0000	1.0000
F81G	39	-1149.5962	2377.1925	0.0000	1.0000
F81I	39	-1150.0475	2378.0949	0.0000	1.0000
F81IG	40	-1149.3725	2378.7449	0.0000	1.0000
HKY	39	-1156.5603	2391.1205	0.0000	1.0000
SYM	40	-1158.1101	2396.2201	0.0000	1.0000
JC69G	36	-1164.2960	2400.5919	0.0000	1.0000
JC69I	36	-1164.6981	2401.3962	0.0000	1.0000
JC69IG	37	-1164.1889	2402.3779	0.0000	1.0000
K2P	36	-1173.4904	2418.9808	0.0000	1.0000
F81	38	-1174.0835	2424.1669	0.0000	1.0000
JC69	35	-1188.6413	2447.2825	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1110.6677	2481.2911	0.7710	0.7710
GTRI	44	-1112.2378	2484.4312	0.1604	0.9314
GTRIG	45	-1110.1347	2486.1331	0.0685	0.9999
HKYG	40	-1131.9907	2500.3047	0.0001	1.0000
HKYI	40	-1132.5038	2501.3309	0.0000	1.0000
HKYIG	41	-1131.7390	2505.7094	0.0000	1.0000
SYMG	41	-1133.4841	2509.1996	0.0000	1.0000
SYMI	41	-1134.3716	2510.9746	0.0000	1.0000
SYMIG	42	-1133.5159	2515.1712	0.0000	1.0000
K2PG	37	-1148.6629	2515.9250	0.0000	1.0000
K2PI	37	-1148.9962	2516.5914	0.0000	1.0000
K2PIG	38	-1148.4404	2521.3879	0.0000	1.0000
GTR	43	-1137.0411	2528.1298	0.0000	1.0000
F81G	39	-1149.5962	2529.6077	0.0000	1.0000
F81I	39	-1150.0475	2530.5101	0.0000	1.0000
F81IG	40	-1149.3725	2535.0682	0.0000	1.0000
JC69G	36	-1164.2960	2541.2829	0.0000	1.0000
JC69I	36	-1164.6981	2542.0871	0.0000	1.0000
HKY	39	-1156.5603	2543.5358	0.0000	1.0000
JC69IG	37	-1164.1889	2546.9769	0.0000	1.0000
SYM	40	-1158.1101	2552.5434	0.0000	1.0000
K2P	36	-1173.4904	2559.6718	0.0000	1.0000
F81	38	-1174.0835	2572.6741	0.0000	1.0000
JC69	35	-1188.6413	2584.0654	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_12990.phylip = 1 - 368;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_12990.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_12990.phylip = 1 - 368;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_12990.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 (GTRG)]

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


-------------------------------------------------------------
End of Output

