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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	40	-1361.9771	2813.1676	0.3492	0.3492
HKYI	40	-1361.9980	2813.2095	0.3420	0.6912
GTRG	44	-1358.2857	2815.8214	0.0926	0.7838
HKYIG	41	-1362.0747	2815.8508	0.0913	0.8751
GTRI	44	-1358.3246	2815.8992	0.0891	0.9642
GTRIG	45	-1357.9639	2817.7226	0.0358	1.0000
K2PG	37	-1383.4717	2848.7763	0.0000	1.0000
K2PI	37	-1383.6828	2849.1984	0.0000	1.0000
K2PIG	38	-1383.4716	2851.2225	0.0000	1.0000
SYMG	41	-1381.2599	2854.2212	0.0000	1.0000
SYMI	41	-1381.3741	2854.4496	0.0000	1.0000
SYMIG	42	-1381.3704	2856.9441	0.0000	1.0000
HKY	39	-1389.7927	2866.3249	0.0000	1.0000
GTR	43	-1385.3987	2867.5170	0.0000	1.0000
F81I	39	-1393.8037	2874.3470	0.0000	1.0000
F81G	39	-1393.8478	2874.4351	0.0000	1.0000
F81IG	40	-1393.8477	2876.9089	0.0000	1.0000
JC69G	36	-1412.2978	2903.9956	0.0000	1.0000
JC69I	36	-1412.3638	2904.1277	0.0000	1.0000
JC69IG	37	-1412.2977	2906.4282	0.0000	1.0000
K2P	36	-1414.6014	2908.6028	0.0000	1.0000
SYM	40	-1412.7451	2914.7037	0.0000	1.0000
F81	38	-1421.5176	2927.3145	0.0000	1.0000
JC69	35	-1441.6785	2960.3375	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1361.9771	2803.9542	0.2423	0.2423
HKYI	40	-1361.9980	2803.9960	0.2373	0.4796
GTRG	44	-1358.2857	2804.5714	0.1780	0.6576
GTRI	44	-1358.3246	2804.6492	0.1712	0.8288
GTRIG	45	-1357.9639	2805.9277	0.0903	0.9191
HKYIG	41	-1362.0747	2806.1494	0.0809	1.0000
K2PG	37	-1383.4717	2840.9435	0.0000	1.0000
K2PI	37	-1383.6828	2841.3656	0.0000	1.0000
K2PIG	38	-1383.4716	2842.9432	0.0000	1.0000
SYMG	41	-1381.2599	2844.5198	0.0000	1.0000
SYMI	41	-1381.3741	2844.7482	0.0000	1.0000
SYMIG	42	-1381.3704	2846.7407	0.0000	1.0000
GTR	43	-1385.3987	2856.7974	0.0000	1.0000
HKY	39	-1389.7927	2857.5854	0.0000	1.0000
F81I	39	-1393.8037	2865.6075	0.0000	1.0000
F81G	39	-1393.8478	2865.6956	0.0000	1.0000
F81IG	40	-1393.8477	2867.6954	0.0000	1.0000
JC69G	36	-1412.2978	2896.5956	0.0000	1.0000
JC69I	36	-1412.3638	2896.7277	0.0000	1.0000
JC69IG	37	-1412.2977	2898.5953	0.0000	1.0000
K2P	36	-1414.6014	2901.2028	0.0000	1.0000
SYM	40	-1412.7451	2905.4902	0.0000	1.0000
F81	38	-1421.5176	2919.0352	0.0000	1.0000
JC69	35	-1441.6785	2953.3569	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1361.9771	2963.3116	0.4937	0.4937
HKYI	40	-1361.9980	2963.3534	0.4835	0.9773
HKYIG	41	-1362.0747	2969.4907	0.0225	0.9997
GTRG	44	-1358.2857	2979.8646	0.0001	0.9999
GTRI	44	-1358.3246	2979.9424	0.0001	1.0000
GTRIG	45	-1357.9639	2985.2049	0.0000	1.0000
K2PG	37	-1383.4717	2988.3491	0.0000	1.0000
K2PI	37	-1383.6828	2988.7712	0.0000	1.0000
K2PIG	38	-1383.4716	2994.3327	0.0000	1.0000
SYMG	41	-1381.2599	3007.8612	0.0000	1.0000
SYMI	41	-1381.3741	3008.0896	0.0000	1.0000
HKY	39	-1389.7927	3012.9589	0.0000	1.0000
SYMIG	42	-1381.3704	3014.0660	0.0000	1.0000
F81I	39	-1393.8037	3020.9810	0.0000	1.0000
F81G	39	-1393.8478	3021.0691	0.0000	1.0000
F81IG	40	-1393.8477	3027.0529	0.0000	1.0000
GTR	43	-1385.3987	3028.1067	0.0000	1.0000
JC69G	36	-1412.2978	3040.0173	0.0000	1.0000
JC69I	36	-1412.3638	3040.1494	0.0000	1.0000
K2P	36	-1414.6014	3044.6245	0.0000	1.0000
JC69IG	37	-1412.2977	3046.0010	0.0000	1.0000
SYM	40	-1412.7451	3064.8477	0.0000	1.0000
F81	38	-1421.5176	3070.4248	0.0000	1.0000
JC69	35	-1441.6785	3092.7947	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10052.phylip = 1 - 397;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10052.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10052.phylip = 1 - 397;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10052.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10052.phylip = 1 - 397;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_10052.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

