
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:38 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/chr7_10489.phylip" (19 taxa, 461 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	-1369.1263	2826.0622	0.5929	0.5929
HKYIG	41	-1368.9542	2828.1280	0.2111	0.8039
HKYI	40	-1370.6287	2829.0669	0.1320	0.9359
GTRG	44	-1366.9291	2831.3774	0.0416	0.9775
GTRIG	45	-1366.7911	2833.5581	0.0140	0.9915
GTRI	44	-1368.5116	2834.5424	0.0085	1.0000
SYMG	41	-1380.3699	2850.9594	0.0000	1.0000
SYMIG	42	-1380.2742	2853.1895	0.0000	1.0000
SYMI	41	-1382.0064	2854.2324	0.0000	1.0000
K2PG	37	-1387.2862	2855.2202	0.0000	1.0000
K2PIG	38	-1387.1300	2857.2838	0.0000	1.0000
K2PI	37	-1388.7415	2858.1308	0.0000	1.0000
HKY	39	-1387.4658	2860.3426	0.0000	1.0000
GTR	43	-1384.3371	2863.7485	0.0000	1.0000
F81G	39	-1392.1291	2869.6692	0.0000	1.0000
F81IG	40	-1391.9728	2871.7551	0.0000	1.0000
F81I	39	-1393.6609	2872.7328	0.0000	1.0000
SYM	40	-1396.7589	2881.3272	0.0000	1.0000
K2P	36	-1404.7526	2887.7883	0.0000	1.0000
F81	38	-1409.7133	2902.4504	0.0000	1.0000
JC69G	36	-1413.4412	2905.1655	0.0000	1.0000
JC69IG	37	-1413.3503	2907.3485	0.0000	1.0000
JC69I	36	-1414.7476	2907.7782	0.0000	1.0000
JC69	35	-1429.3534	2934.6363	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1369.1263	2818.2526	0.5187	0.5187
HKYIG	41	-1368.9542	2819.9084	0.2267	0.7454
HKYI	40	-1370.6287	2821.2574	0.1155	0.8608
GTRG	44	-1366.9291	2821.8581	0.0855	0.9463
GTRIG	45	-1366.7911	2823.5822	0.0361	0.9824
GTRI	44	-1368.5116	2825.0232	0.0176	1.0000
SYMG	41	-1380.3699	2842.7398	0.0000	1.0000
SYMIG	42	-1380.2742	2844.5484	0.0000	1.0000
SYMI	41	-1382.0064	2846.0128	0.0000	1.0000
K2PG	37	-1387.2862	2848.5724	0.0000	1.0000
K2PIG	38	-1387.1300	2850.2601	0.0000	1.0000
K2PI	37	-1388.7415	2851.4830	0.0000	1.0000
HKY	39	-1387.4658	2852.9316	0.0000	1.0000
GTR	43	-1384.3371	2854.6741	0.0000	1.0000
F81G	39	-1392.1291	2862.2582	0.0000	1.0000
F81IG	40	-1391.9728	2863.9456	0.0000	1.0000
F81I	39	-1393.6609	2865.3219	0.0000	1.0000
SYM	40	-1396.7589	2873.5177	0.0000	1.0000
K2P	36	-1404.7526	2881.5053	0.0000	1.0000
F81	38	-1409.7133	2895.4267	0.0000	1.0000
JC69G	36	-1413.4412	2898.8825	0.0000	1.0000
JC69IG	37	-1413.3503	2900.7007	0.0000	1.0000
JC69I	36	-1414.7476	2901.4952	0.0000	1.0000
JC69	35	-1429.3534	2928.7069	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1369.1263	2983.5886	0.7824	0.7824
HKYI	40	-1370.6287	2986.5933	0.1742	0.9565
HKYIG	41	-1368.9542	2989.3777	0.0433	0.9998
K2PG	37	-1387.2862	3001.5081	0.0001	0.9999
GTRG	44	-1366.9291	3003.7277	0.0000	1.0000
K2PI	37	-1388.7415	3004.4187	0.0000	1.0000
GTRI	44	-1368.5116	3006.8927	0.0000	1.0000
K2PIG	38	-1387.1300	3007.3292	0.0000	1.0000
GTRIG	45	-1366.7911	3009.5851	0.0000	1.0000
SYMG	41	-1380.3699	3012.2091	0.0000	1.0000
HKY	39	-1387.4658	3014.1342	0.0000	1.0000
SYMI	41	-1382.0064	3015.4821	0.0000	1.0000
SYMIG	42	-1380.2742	3018.1511	0.0000	1.0000
F81G	39	-1392.1291	3023.4608	0.0000	1.0000
F81I	39	-1393.6609	3026.5244	0.0000	1.0000
F81IG	40	-1391.9728	3029.2815	0.0000	1.0000
K2P	36	-1404.7526	3030.3076	0.0000	1.0000
GTR	43	-1384.3371	3032.4102	0.0000	1.0000
SYM	40	-1396.7589	3038.8536	0.0000	1.0000
JC69G	36	-1413.4412	3047.6848	0.0000	1.0000
JC69I	36	-1414.7476	3050.2975	0.0000	1.0000
F81	38	-1409.7133	3052.4958	0.0000	1.0000
JC69IG	37	-1413.3503	3053.6364	0.0000	1.0000
JC69	35	-1429.3534	3073.3758	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/chr7_10489.phylip = 1 - 461;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10489.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/chr7_10489.phylip = 1 - 461;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10489.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/chr7_10489.phylip = 1 - 461;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10489.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

