
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:35 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_9895.phylip" (19 taxa, 397 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-1148.4374	2386.0882	0.4125	0.4125
HKYG	40	-1148.5283	2386.2700	0.3767	0.7892
HKYIG	41	-1148.5521	2388.8057	0.1060	0.8952
GTRI	44	-1145.5413	2390.3327	0.0494	0.9446
GTRG	44	-1145.6936	2390.6372	0.0424	0.9870
GTRIG	45	-1145.7178	2393.2304	0.0116	0.9986
HKY	39	-1155.5264	2397.7924	0.0012	0.9998
GTR	43	-1152.4369	2401.5934	0.0002	1.0000
F81I	39	-1160.7624	2408.2644	0.0000	1.0000
F81G	39	-1160.8058	2408.3511	0.0000	1.0000
F81IG	40	-1160.8328	2410.8791	0.0000	1.0000
K2PI	37	-1167.3901	2416.6130	0.0000	1.0000
K2PG	37	-1167.4876	2416.8081	0.0000	1.0000
SYMI	41	-1163.4166	2418.5346	0.0000	1.0000
SYMG	41	-1163.4936	2418.6887	0.0000	1.0000
K2PIG	38	-1167.5065	2419.2922	0.0000	1.0000
F81	38	-1168.1542	2420.5877	0.0000	1.0000
SYMIG	42	-1163.5002	2421.2037	0.0000	1.0000
K2P	36	-1174.2072	2427.8143	0.0000	1.0000
SYM	40	-1169.5771	2428.3676	0.0000	1.0000
JC69I	36	-1181.6697	2442.7393	0.0000	1.0000
JC69G	36	-1181.7390	2442.8781	0.0000	1.0000
JC69IG	37	-1181.7560	2445.3449	0.0000	1.0000
JC69	35	-1188.2961	2453.5728	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1148.4374	2376.8748	0.3375	0.3375
HKYG	40	-1148.5283	2377.0565	0.3082	0.6457
GTRI	44	-1145.5413	2379.0827	0.1119	0.7576
HKYIG	41	-1148.5521	2379.1043	0.1107	0.8683
GTRG	44	-1145.6936	2379.3872	0.0961	0.9644
GTRIG	45	-1145.7178	2381.4355	0.0345	0.9989
HKY	39	-1155.5264	2389.0529	0.0008	0.9997
GTR	43	-1152.4369	2390.8739	0.0003	1.0000
F81I	39	-1160.7624	2399.5249	0.0000	1.0000
F81G	39	-1160.8058	2399.6116	0.0000	1.0000
F81IG	40	-1160.8328	2401.6656	0.0000	1.0000
K2PI	37	-1167.3901	2408.7801	0.0000	1.0000
SYMI	41	-1163.4166	2408.8332	0.0000	1.0000
K2PG	37	-1167.4876	2408.9752	0.0000	1.0000
SYMG	41	-1163.4936	2408.9873	0.0000	1.0000
SYMIG	42	-1163.5002	2411.0003	0.0000	1.0000
K2PIG	38	-1167.5065	2411.0129	0.0000	1.0000
F81	38	-1168.1542	2412.3084	0.0000	1.0000
SYM	40	-1169.5771	2419.1541	0.0000	1.0000
K2P	36	-1174.2072	2420.4143	0.0000	1.0000
JC69I	36	-1181.6697	2435.3393	0.0000	1.0000
JC69G	36	-1181.7390	2435.4781	0.0000	1.0000
JC69IG	37	-1181.7560	2437.5120	0.0000	1.0000
JC69	35	-1188.2961	2446.5922	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1148.4374	2536.2322	0.5063	0.5063
HKYG	40	-1148.5283	2536.4140	0.4623	0.9687
HKYIG	41	-1148.5521	2542.4457	0.0227	0.9913
HKY	39	-1155.5264	2544.4264	0.0084	0.9998
GTRI	44	-1145.5413	2554.3759	0.0001	0.9998
GTRG	44	-1145.6936	2554.6804	0.0000	0.9999
F81I	39	-1160.7624	2554.8984	0.0000	0.9999
F81G	39	-1160.8058	2554.9851	0.0000	0.9999
K2PI	37	-1167.3901	2556.1858	0.0000	1.0000
K2PG	37	-1167.4876	2556.3809	0.0000	1.0000
GTRIG	45	-1145.7178	2560.7127	0.0000	1.0000
F81IG	40	-1160.8328	2561.0230	0.0000	1.0000
GTR	43	-1152.4369	2562.1831	0.0000	1.0000
K2PIG	38	-1167.5065	2562.4025	0.0000	1.0000
F81	38	-1168.1542	2563.6979	0.0000	1.0000
K2P	36	-1174.2072	2563.8360	0.0000	1.0000
SYMI	41	-1163.4166	2572.1746	0.0000	1.0000
SYMG	41	-1163.4936	2572.3287	0.0000	1.0000
SYMIG	42	-1163.5002	2578.3256	0.0000	1.0000
SYM	40	-1169.5771	2578.5116	0.0000	1.0000
JC69I	36	-1181.6697	2578.7610	0.0000	1.0000
JC69G	36	-1181.7390	2578.8998	0.0000	1.0000
JC69IG	37	-1181.7560	2584.9177	0.0000	1.0000
JC69	35	-1188.2961	2586.0299	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9895.phylip = 1 - 397;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9895.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9895.phylip = 1 - 397;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9895.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9895.phylip = 1 - 397;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_9895.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

