
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:36 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/chr3_9719.phylip" (19 taxa, 510 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	-1344.8041	2776.6019	0.4111	0.4111
HKYG	40	-1344.8961	2776.7858	0.3750	0.7861
HKYIG	41	-1344.8959	2779.1507	0.1149	0.9010
GTRI	44	-1342.2737	2781.0635	0.0442	0.9452
GTRG	44	-1342.3203	2781.1567	0.0422	0.9873
GTRIG	45	-1342.3200	2783.5625	0.0127	1.0000
K2PI	37	-1359.4410	2798.8397	0.0000	1.0000
K2PG	37	-1359.5005	2798.9586	0.0000	1.0000
HKY	39	-1357.7656	2800.1694	0.0000	1.0000
K2PIG	38	-1359.5002	2801.2935	0.0000	1.0000
SYMI	41	-1356.9834	2803.3258	0.0000	1.0000
SYMG	41	-1356.9986	2803.3562	0.0000	1.0000
GTR	43	-1355.3177	2804.7555	0.0000	1.0000
SYMIG	42	-1357.0140	2805.7624	0.0000	1.0000
F81I	39	-1362.2002	2809.0388	0.0000	1.0000
F81G	39	-1362.2131	2809.0645	0.0000	1.0000
F81IG	40	-1362.2182	2811.4301	0.0000	1.0000
K2P	36	-1372.8540	2823.3402	0.0000	1.0000
SYM	40	-1370.4195	2827.8325	0.0000	1.0000
JC69I	36	-1375.7293	2829.0907	0.0000	1.0000
JC69G	36	-1375.7327	2829.0975	0.0000	1.0000
JC69IG	37	-1375.7381	2831.4339	0.0000	1.0000
F81	38	-1375.2176	2832.7283	0.0000	1.0000
JC69	35	-1389.0116	2853.3397	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1344.8041	2769.6083	0.3599	0.3599
HKYG	40	-1344.8961	2769.7922	0.3283	0.6883
HKYIG	41	-1344.8959	2771.7917	0.1208	0.8091
GTRI	44	-1342.2737	2772.5473	0.0828	0.8919
GTRG	44	-1342.3203	2772.6406	0.0790	0.9709
GTRIG	45	-1342.3200	2774.6400	0.0291	1.0000
K2PI	37	-1359.4410	2792.8821	0.0000	1.0000
K2PG	37	-1359.5005	2793.0010	0.0000	1.0000
HKY	39	-1357.7656	2793.5311	0.0000	1.0000
K2PIG	38	-1359.5002	2795.0005	0.0000	1.0000
SYMI	41	-1356.9834	2795.9668	0.0000	1.0000
SYMG	41	-1356.9986	2795.9972	0.0000	1.0000
GTR	43	-1355.3177	2796.6353	0.0000	1.0000
SYMIG	42	-1357.0140	2798.0280	0.0000	1.0000
F81I	39	-1362.2002	2802.4005	0.0000	1.0000
F81G	39	-1362.2131	2802.4262	0.0000	1.0000
F81IG	40	-1362.2182	2804.4365	0.0000	1.0000
K2P	36	-1372.8540	2817.7080	0.0000	1.0000
SYM	40	-1370.4195	2820.8389	0.0000	1.0000
JC69I	36	-1375.7293	2823.4586	0.0000	1.0000
JC69G	36	-1375.7327	2823.4653	0.0000	1.0000
JC69IG	37	-1375.7381	2825.4763	0.0000	1.0000
F81	38	-1375.2176	2826.4353	0.0000	1.0000
JC69	35	-1389.0116	2848.0232	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1344.8041	2938.9847	0.5095	0.5095
HKYG	40	-1344.8961	2939.1686	0.4647	0.9742
HKYIG	41	-1344.8959	2945.4026	0.0206	0.9948
K2PI	37	-1359.4410	2949.5553	0.0026	0.9974
K2PG	37	-1359.5005	2949.6742	0.0024	0.9998
K2PIG	38	-1359.5002	2955.9081	0.0001	0.9999
HKY	39	-1357.7656	2958.6732	0.0000	0.9999
GTRI	44	-1342.2737	2958.8614	0.0000	1.0000
GTRG	44	-1342.3203	2958.9547	0.0000	1.0000
GTRIG	45	-1342.3200	2965.1885	0.0000	1.0000
F81I	39	-1362.2002	2967.5425	0.0000	1.0000
F81G	39	-1362.2131	2967.5682	0.0000	1.0000
SYMI	41	-1356.9834	2969.5776	0.0000	1.0000
SYMG	41	-1356.9986	2969.6081	0.0000	1.0000
K2P	36	-1372.8540	2970.1468	0.0000	1.0000
F81IG	40	-1362.2182	2973.8129	0.0000	1.0000
SYMIG	42	-1357.0140	2975.8732	0.0000	1.0000
JC69I	36	-1375.7293	2975.8973	0.0000	1.0000
JC69G	36	-1375.7327	2975.9041	0.0000	1.0000
GTR	43	-1355.3177	2978.7150	0.0000	1.0000
JC69IG	37	-1375.7381	2982.1495	0.0000	1.0000
F81	38	-1375.2176	2987.3429	0.0000	1.0000
SYM	40	-1370.4195	2990.2153	0.0000	1.0000
JC69	35	-1389.0116	2996.2276	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/chr3_9719.phylip = 1 - 510;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_9719.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/chr3_9719.phylip = 1 - 510;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_9719.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/chr3_9719.phylip = 1 - 510;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_9719.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

