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

Minimum AIC  model: K2PI
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	-507.5403	1111.7581	0.8972	0.8972
K2PG	37	-509.7954	1116.2683	0.0941	0.9913
HKYI	40	-507.2280	1121.5635	0.0067	0.9980
K2PIG	38	-512.5323	1125.1621	0.0011	0.9991
HKYG	40	-509.4578	1126.0231	0.0007	0.9998
SYMG	41	-509.9490	1130.5980	0.0001	0.9998
SYMI	41	-510.0133	1130.7266	0.0001	0.9999
JC69G	36	-519.5301	1132.3723	0.0000	0.9999
JC69I	36	-519.5684	1132.4489	0.0000	1.0000
HKYIG	41	-512.1464	1134.9929	0.0000	1.0000
SYMIG	42	-510.3816	1135.1161	0.0000	1.0000
JC69IG	37	-519.4342	1135.5458	0.0000	1.0000
GTRG	44	-509.3369	1140.5200	0.0000	1.0000
GTRI	44	-509.3916	1140.6293	0.0000	1.0000
F81G	39	-519.3754	1142.3246	0.0000	1.0000
F81I	39	-519.4115	1142.3968	0.0000	1.0000
GTRIG	45	-509.7733	1145.2362	0.0000	1.0000
F81IG	40	-519.2763	1145.6600	0.0000	1.0000
K2P	36	-534.7148	1162.7417	0.0000	1.0000
HKY	39	-534.6236	1172.8210	0.0000	1.0000
SYM	40	-535.5550	1178.2174	0.0000	1.0000
GTR	43	-534.9031	1187.8739	0.0000	1.0000
JC69	35	-551.1858	1192.3716	0.0000	1.0000
F81	38	-551.0680	1202.2336	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	37	-507.5403	1089.0807	0.8426	0.8426
K2PG	37	-509.7954	1093.5909	0.0884	0.9309
HKYI	40	-507.2280	1094.4560	0.0573	0.9882
HKYG	40	-509.4578	1098.9157	0.0062	0.9944
K2PIG	38	-512.5323	1101.0645	0.0021	0.9965
SYMG	41	-509.9490	1101.8980	0.0014	0.9979
SYMI	41	-510.0133	1102.0266	0.0013	0.9992
SYMIG	42	-510.3816	1104.7632	0.0003	0.9995
HKYIG	41	-512.1464	1106.2929	0.0002	0.9997
GTRG	44	-509.3369	1106.6739	0.0001	0.9998
GTRI	44	-509.3916	1106.7832	0.0001	0.9999
GTRIG	45	-509.7733	1109.5465	0.0000	1.0000
JC69G	36	-519.5301	1111.0603	0.0000	1.0000
JC69I	36	-519.5684	1111.1369	0.0000	1.0000
JC69IG	37	-519.4342	1112.8684	0.0000	1.0000
F81G	39	-519.3754	1116.7508	0.0000	1.0000
F81I	39	-519.4115	1116.8230	0.0000	1.0000
F81IG	40	-519.2763	1118.5526	0.0000	1.0000
K2P	36	-534.7148	1141.4297	0.0000	1.0000
HKY	39	-534.6236	1147.2472	0.0000	1.0000
SYM	40	-535.5550	1151.1100	0.0000	1.0000
GTR	43	-534.9031	1155.8061	0.0000	1.0000
JC69	35	-551.1858	1172.3716	0.0000	1.0000
F81	38	-551.0680	1178.1361	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-507.5403	1203.3217	0.9039	0.9039
K2PG	37	-509.7954	1207.8319	0.0948	0.9987
HKYI	40	-507.2280	1217.9599	0.0006	0.9993
K2PIG	38	-512.5323	1218.3932	0.0005	0.9998
JC69G	36	-519.5301	1222.2137	0.0001	0.9999
JC69I	36	-519.5684	1222.2903	0.0001	0.9999
HKYG	40	-509.4578	1222.4195	0.0001	1.0000
JC69IG	37	-519.4342	1227.1094	0.0000	1.0000
SYMG	41	-509.9490	1228.4895	0.0000	1.0000
SYMI	41	-510.0133	1228.6180	0.0000	1.0000
HKYIG	41	-512.1464	1232.8843	0.0000	1.0000
SYMIG	42	-510.3816	1234.4422	0.0000	1.0000
F81G	39	-519.3754	1237.1670	0.0000	1.0000
F81I	39	-519.4115	1237.2393	0.0000	1.0000
F81IG	40	-519.2763	1242.0564	0.0000	1.0000
GTRG	44	-509.3369	1242.5281	0.0000	1.0000
GTRI	44	-509.3916	1242.6374	0.0000	1.0000
GTRIG	45	-509.7733	1248.4883	0.0000	1.0000
K2P	36	-534.7148	1252.5831	0.0000	1.0000
HKY	39	-534.6236	1267.6634	0.0000	1.0000
SYM	40	-535.5550	1274.6138	0.0000	1.0000
JC69	35	-551.1858	1280.4375	0.0000	1.0000
GTR	43	-534.9031	1288.5728	0.0000	1.0000
F81	38	-551.0680	1295.4647	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_30414.phylip = 1 - 162;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_30414.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 AICc model (K2PI)]

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

