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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: K2PG


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
HKYIG	41	-1131.2639	2355.3922	0.4144	0.4144
HKYG	40	-1133.0597	2356.4340	0.2461	0.6605
HKYI	40	-1133.2604	2356.8353	0.2014	0.8619
GTRG	44	-1129.1773	2358.9660	0.0694	0.9313
GTRI	44	-1129.7156	2360.0427	0.0405	0.9719
GTRIG	45	-1129.0053	2361.2375	0.0223	0.9941
K2PG	37	-1141.1060	2364.9721	0.0034	0.9976
K2PIG	38	-1140.9601	2367.1828	0.0011	0.9987
K2PI	37	-1142.3141	2367.3884	0.0010	0.9998
SYMG	41	-1139.1718	2371.2081	0.0002	0.9999
SYMIG	42	-1139.0258	2373.4820	0.0000	1.0000
SYMI	41	-1140.4502	2373.7647	0.0000	1.0000
HKY	39	-1147.8506	2383.4818	0.0000	1.0000
GTR	43	-1144.2847	2386.5821	0.0000	1.0000
K2P	36	-1156.3376	2392.9485	0.0000	1.0000
SYM	40	-1154.3349	2398.9842	0.0000	1.0000
F81I	39	-1171.7961	2431.3727	0.0000	1.0000
F81G	39	-1171.7996	2431.3799	0.0000	1.0000
F81IG	40	-1171.9555	2434.2254	0.0000	1.0000
JC69G	36	-1178.1229	2436.5191	0.0000	1.0000
JC69IG	37	-1178.0192	2438.7985	0.0000	1.0000
JC69I	36	-1179.3552	2438.9837	0.0000	1.0000
F81	38	-1185.0209	2455.3042	0.0000	1.0000
JC69	35	-1191.5809	2460.9637	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-1131.2639	2344.5278	0.3793	0.3793
HKYG	40	-1133.0597	2346.1195	0.1711	0.5504
GTRG	44	-1129.1773	2346.3545	0.1522	0.7026
HKYI	40	-1133.2604	2346.5209	0.1400	0.8426
GTRI	44	-1129.7156	2347.4312	0.0888	0.9314
GTRIG	45	-1129.0053	2348.0107	0.0665	0.9979
K2PG	37	-1141.1060	2356.2120	0.0011	0.9990
K2PIG	38	-1140.9601	2357.9203	0.0005	0.9994
K2PI	37	-1142.3141	2358.6283	0.0003	0.9998
SYMG	41	-1139.1718	2360.3437	0.0001	0.9999
SYMIG	42	-1139.0258	2362.0516	0.0001	1.0000
SYMI	41	-1140.4502	2362.9004	0.0000	1.0000
HKY	39	-1147.8506	2373.7013	0.0000	1.0000
GTR	43	-1144.2847	2374.5694	0.0000	1.0000
K2P	36	-1156.3376	2384.6752	0.0000	1.0000
SYM	40	-1154.3349	2388.6697	0.0000	1.0000
F81I	39	-1171.7961	2421.5921	0.0000	1.0000
F81G	39	-1171.7996	2421.5993	0.0000	1.0000
F81IG	40	-1171.9555	2423.9110	0.0000	1.0000
JC69G	36	-1178.1229	2428.2458	0.0000	1.0000
JC69IG	37	-1178.0192	2430.0384	0.0000	1.0000
JC69I	36	-1179.3552	2430.7104	0.0000	1.0000
F81	38	-1185.0209	2446.0417	0.0000	1.0000
JC69	35	-1191.5809	2453.1619	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1141.1060	2499.8949	0.4272	0.4272
HKYG	40	-1133.0597	2501.4524	0.1961	0.6233
HKYI	40	-1133.2604	2501.8538	0.1604	0.7838
K2PI	37	-1142.3141	2502.3112	0.1276	0.9114
HKYIG	41	-1131.2639	2503.7440	0.0624	0.9738
K2PIG	38	-1140.9601	2505.4865	0.0261	0.9998
GTRG	44	-1129.1773	2517.2207	0.0001	0.9999
GTRI	44	-1129.7156	2518.2974	0.0000	1.0000
SYMG	41	-1139.1718	2519.5599	0.0000	1.0000
SYMI	41	-1140.4502	2522.1166	0.0000	1.0000
GTRIG	45	-1129.0053	2522.7602	0.0000	1.0000
K2P	36	-1156.3376	2524.4748	0.0000	1.0000
HKY	39	-1147.8506	2525.1509	0.0000	1.0000
SYMIG	42	-1139.0258	2525.1512	0.0000	1.0000
GTR	43	-1144.2847	2541.5523	0.0000	1.0000
SYM	40	-1154.3349	2544.0026	0.0000	1.0000
JC69G	36	-1178.1229	2568.0454	0.0000	1.0000
JC69I	36	-1179.3552	2570.5100	0.0000	1.0000
F81I	39	-1171.7961	2573.0417	0.0000	1.0000
F81G	39	-1171.7996	2573.0489	0.0000	1.0000
JC69IG	37	-1178.0192	2573.7213	0.0000	1.0000
F81IG	40	-1171.9555	2579.2439	0.0000	1.0000
JC69	35	-1191.5809	2589.0781	0.0000	1.0000
F81	38	-1185.0209	2593.6080	0.0000	1.0000



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

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

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


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

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


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

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


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

