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

Minimum AIC  model: GTRI
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	37	-534.2372	1167.3593	0.4569	0.4569
K2PG	37	-534.8101	1168.5052	0.2576	0.7145
HKYI	40	-530.2063	1170.2307	0.1087	0.8232
HKYG	40	-530.7992	1171.4165	0.0601	0.8833
K2PIG	38	-534.8100	1172.0843	0.0430	0.9263
SYMI	41	-529.8355	1173.2674	0.0238	0.9502
GTRI	44	-524.1972	1173.7529	0.0187	0.9688
SYMG	41	-530.5735	1174.7434	0.0114	0.9802
HKYIG	41	-530.7991	1175.1945	0.0091	0.9893
GTRG	44	-525.0465	1175.4515	0.0080	0.9973
SYMIG	42	-530.5874	1178.6193	0.0016	0.9990
GTRIG	45	-525.0463	1179.5212	0.0010	1.0000
JC69I	36	-548.9785	1193.3254	0.0000	1.0000
JC69G	36	-549.6392	1194.6469	0.0000	1.0000
K2P	36	-550.3481	1196.0645	0.0000	1.0000
F81I	39	-545.7402	1197.5886	0.0000	1.0000
HKY	39	-545.8695	1197.8470	0.0000	1.0000
JC69IG	37	-549.6392	1198.1633	0.0000	1.0000
F81G	39	-546.4282	1198.9645	0.0000	1.0000
F81IG	40	-546.4281	1202.6745	0.0000	1.0000
SYM	40	-546.6445	1203.1071	0.0000	1.0000
GTR	43	-540.9939	1203.3524	0.0000	1.0000
JC69	35	-566.1422	1224.1975	0.0000	1.0000
F81	38	-562.8037	1228.0717	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	44	-524.1972	1136.3944	0.4942	0.4942
GTRG	44	-525.0465	1138.0930	0.2114	0.7056
GTRIG	45	-525.0463	1140.0926	0.0778	0.7834
HKYI	40	-530.2063	1140.4126	0.0663	0.8497
HKYG	40	-530.7992	1141.5984	0.0366	0.8863
SYMI	41	-529.8355	1141.6710	0.0353	0.9216
K2PI	37	-534.2372	1142.4744	0.0236	0.9453
SYMG	41	-530.5735	1143.1471	0.0169	0.9622
HKYIG	41	-530.7991	1143.5982	0.0135	0.9756
K2PG	37	-534.8101	1143.6202	0.0133	0.9890
SYMIG	42	-530.5874	1145.1749	0.0061	0.9951
K2PIG	38	-534.8100	1145.6201	0.0049	1.0000
GTR	43	-540.9939	1167.9879	0.0000	1.0000
F81I	39	-545.7402	1169.4805	0.0000	1.0000
HKY	39	-545.8695	1169.7389	0.0000	1.0000
JC69I	36	-548.9785	1169.9569	0.0000	1.0000
F81G	39	-546.4282	1170.8564	0.0000	1.0000
JC69G	36	-549.6392	1171.2785	0.0000	1.0000
K2P	36	-550.3481	1172.6961	0.0000	1.0000
F81IG	40	-546.4281	1172.8563	0.0000	1.0000
JC69IG	37	-549.6392	1173.2783	0.0000	1.0000
SYM	40	-546.6445	1173.2889	0.0000	1.0000
F81	38	-562.8037	1201.6075	0.0000	1.0000
JC69	35	-566.1422	1202.2845	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-534.2372	1254.1137	0.6008	0.6008
K2PG	37	-534.8101	1255.2596	0.3388	0.9396
K2PIG	38	-534.8100	1260.2767	0.0276	0.9671
HKYI	40	-530.2063	1261.1038	0.0182	0.9854
HKYG	40	-530.7992	1262.2896	0.0101	0.9954
SYMI	41	-529.8355	1265.3795	0.0021	0.9976
SYMG	41	-530.5735	1266.8556	0.0010	0.9986
HKYIG	41	-530.7991	1267.3067	0.0008	0.9994
GTRI	44	-524.1972	1269.1547	0.0003	0.9998
GTRG	44	-525.0465	1270.8533	0.0001	0.9999
SYMIG	42	-530.5874	1271.9006	0.0001	1.0000
GTRIG	45	-525.0463	1275.8702	0.0000	1.0000
JC69I	36	-548.9785	1278.5790	0.0000	1.0000
JC69G	36	-549.6392	1279.9006	0.0000	1.0000
K2P	36	-550.3481	1281.3182	0.0000	1.0000
JC69IG	37	-549.6392	1284.9177	0.0000	1.0000
F81I	39	-545.7402	1287.1544	0.0000	1.0000
HKY	39	-545.8695	1287.4128	0.0000	1.0000
F81G	39	-546.4282	1288.5304	0.0000	1.0000
F81IG	40	-546.4281	1293.5475	0.0000	1.0000
SYM	40	-546.6445	1293.9801	0.0000	1.0000
GTR	43	-540.9939	1297.7309	0.0000	1.0000
JC69	35	-566.1422	1307.8893	0.0000	1.0000
F81	38	-562.8037	1316.2641	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_851.phylip = 1 - 151;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_851.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

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


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

