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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	44	-1156.0782	2410.3626	0.3003	0.3003
GTRG	44	-1156.3769	2410.9600	0.2227	0.5230
HKYI	40	-1161.4721	2411.3115	0.1868	0.7098
HKYG	40	-1161.7716	2411.9105	0.1385	0.8483
GTRIG	45	-1156.0856	2412.8688	0.0858	0.9340
HKYIG	41	-1161.4800	2413.7681	0.0547	0.9887
GTR	43	-1161.0477	2417.8229	0.0072	0.9959
HKY	39	-1166.5171	2418.9731	0.0041	1.0000
K2PI	37	-1188.7298	2458.5786	0.0000	1.0000
K2PG	37	-1188.9225	2458.9640	0.0000	1.0000
K2PIG	38	-1188.7350	2460.9928	0.0000	1.0000
SYMI	41	-1185.3876	2461.5834	0.0000	1.0000
SYMG	41	-1185.5681	2461.9443	0.0000	1.0000
SYMIG	42	-1185.3921	2464.0457	0.0000	1.0000
K2P	36	-1193.5538	2465.8348	0.0000	1.0000
SYM	40	-1190.0339	2468.4351	0.0000	1.0000
F81I	39	-1193.5359	2473.0107	0.0000	1.0000
F81G	39	-1193.7389	2473.4167	0.0000	1.0000
F81IG	40	-1193.5412	2475.4497	0.0000	1.0000
F81	38	-1198.2280	2479.9788	0.0000	1.0000
JC69I	36	-1216.0385	2510.8043	0.0000	1.0000
JC69G	36	-1216.1897	2511.1067	0.0000	1.0000
JC69IG	37	-1216.1973	2513.5136	0.0000	1.0000
JC69	35	-1220.7781	2517.9038	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	44	-1156.0782	2400.1564	0.3770	0.3770
GTRG	44	-1156.3769	2400.7538	0.2796	0.6566
GTRIG	45	-1156.0856	2402.1711	0.1377	0.7943
HKYI	40	-1161.4721	2402.9442	0.0935	0.8878
HKYG	40	-1161.7716	2403.5432	0.0693	0.9571
HKYIG	41	-1161.4800	2404.9599	0.0341	0.9912
GTR	43	-1161.0477	2408.0954	0.0071	0.9984
HKY	39	-1166.5171	2411.0341	0.0016	1.0000
K2PI	37	-1188.7298	2451.4597	0.0000	1.0000
K2PG	37	-1188.9225	2451.8450	0.0000	1.0000
SYMI	41	-1185.3876	2452.7752	0.0000	1.0000
SYMG	41	-1185.5681	2453.1361	0.0000	1.0000
K2PIG	38	-1188.7350	2453.4699	0.0000	1.0000
SYMIG	42	-1185.3921	2454.7841	0.0000	1.0000
K2P	36	-1193.5538	2459.1075	0.0000	1.0000
SYM	40	-1190.0339	2460.0678	0.0000	1.0000
F81I	39	-1193.5359	2465.0717	0.0000	1.0000
F81G	39	-1193.7389	2465.4778	0.0000	1.0000
F81IG	40	-1193.5412	2467.0824	0.0000	1.0000
F81	38	-1198.2280	2472.4560	0.0000	1.0000
JC69I	36	-1216.0385	2504.0770	0.0000	1.0000
JC69G	36	-1216.1897	2504.3795	0.0000	1.0000
JC69IG	37	-1216.1973	2506.3946	0.0000	1.0000
JC69	35	-1220.7781	2511.5562	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1161.4721	2565.7737	0.5194	0.5194
HKYG	40	-1161.7716	2566.3727	0.3850	0.9044
HKY	39	-1166.5171	2569.7929	0.0696	0.9741
HKYIG	41	-1161.4800	2571.8602	0.0248	0.9988
GTRI	44	-1156.0782	2579.2689	0.0006	0.9994
GTRG	44	-1156.3769	2579.8662	0.0005	0.9999
GTR	43	-1161.0477	2583.1371	0.0001	1.0000
GTRIG	45	-1156.0856	2585.3543	0.0000	1.0000
K2PI	37	-1188.7298	2602.0770	0.0000	1.0000
K2PG	37	-1188.9225	2602.4623	0.0000	1.0000
K2P	36	-1193.5538	2605.6541	0.0000	1.0000
K2PIG	38	-1188.7350	2608.1580	0.0000	1.0000
SYMI	41	-1185.3876	2619.6755	0.0000	1.0000
SYMG	41	-1185.5681	2620.0364	0.0000	1.0000
SYM	40	-1190.0339	2622.8973	0.0000	1.0000
F81I	39	-1193.5359	2623.8305	0.0000	1.0000
F81G	39	-1193.7389	2624.2366	0.0000	1.0000
SYMIG	42	-1185.3921	2625.7551	0.0000	1.0000
F81	38	-1198.2280	2627.1440	0.0000	1.0000
F81IG	40	-1193.5412	2629.9119	0.0000	1.0000
JC69I	36	-1216.0385	2650.6235	0.0000	1.0000
JC69G	36	-1216.1897	2650.9260	0.0000	1.0000
JC69	35	-1220.7781	2654.0320	0.0000	1.0000
JC69IG	37	-1216.1973	2657.0119	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/chr1_31743.phylip = 1 - 433;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_31743.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 (GTRI)]

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

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

