
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:50 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_10370.phylip" (19 taxa, 579 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	-1394.2006	2874.4978	0.3950	0.3950
HKYI	40	-1394.2574	2874.6114	0.3732	0.7681
HKYIG	41	-1394.1155	2876.6443	0.1350	0.9031
GTRG	44	-1391.7805	2878.9768	0.0421	0.9452
GTRI	44	-1391.8806	2879.1770	0.0381	0.9833
GTRIG	45	-1391.5265	2880.8204	0.0167	1.0000
K2PG	37	-1413.4051	2906.0080	0.0000	1.0000
K2PI	37	-1413.4879	2906.1736	0.0000	1.0000
HKY	39	-1412.0438	2907.8761	0.0000	1.0000
K2PIG	38	-1413.2847	2908.0582	0.0000	1.0000
F81G	39	-1413.6712	2911.1310	0.0000	1.0000
SYMG	41	-1411.4177	2911.2487	0.0000	1.0000
F81I	39	-1413.7791	2911.3467	0.0000	1.0000
SYMI	41	-1411.5449	2911.5033	0.0000	1.0000
GTR	43	-1410.0355	2913.1439	0.0000	1.0000
F81IG	40	-1413.5253	2913.1473	0.0000	1.0000
SYMIG	42	-1411.2099	2913.1587	0.0000	1.0000
K2P	36	-1430.1841	2937.2834	0.0000	1.0000
JC69G	36	-1431.0018	2938.9186	0.0000	1.0000
JC69I	36	-1431.1197	2939.1546	0.0000	1.0000
F81	38	-1429.4627	2940.4143	0.0000	1.0000
JC69IG	37	-1430.8558	2940.9094	0.0000	1.0000
SYM	40	-1428.1061	2942.3089	0.0000	1.0000
JC69	35	-1446.4851	2967.6110	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1394.2006	2868.4011	0.3527	0.3527
HKYI	40	-1394.2574	2868.5147	0.3332	0.6859
HKYIG	41	-1394.1155	2870.2309	0.1413	0.8272
GTRG	44	-1391.7805	2871.5611	0.0726	0.8998
GTRI	44	-1391.8806	2871.7613	0.0657	0.9655
GTRIG	45	-1391.5265	2873.0531	0.0345	1.0000
K2PG	37	-1413.4051	2900.8102	0.0000	1.0000
K2PI	37	-1413.4879	2900.9758	0.0000	1.0000
HKY	39	-1412.0438	2902.0876	0.0000	1.0000
K2PIG	38	-1413.2847	2902.5693	0.0000	1.0000
SYMG	41	-1411.4177	2904.8353	0.0000	1.0000
SYMI	41	-1411.5449	2905.0899	0.0000	1.0000
F81G	39	-1413.6712	2905.3425	0.0000	1.0000
F81I	39	-1413.7791	2905.5582	0.0000	1.0000
GTR	43	-1410.0355	2906.0710	0.0000	1.0000
SYMIG	42	-1411.2099	2906.4199	0.0000	1.0000
F81IG	40	-1413.5253	2907.0506	0.0000	1.0000
K2P	36	-1430.1841	2932.3682	0.0000	1.0000
JC69G	36	-1431.0018	2934.0035	0.0000	1.0000
JC69I	36	-1431.1197	2934.2394	0.0000	1.0000
F81	38	-1429.4627	2934.9254	0.0000	1.0000
JC69IG	37	-1430.8558	2935.7117	0.0000	1.0000
SYM	40	-1428.1061	2936.2122	0.0000	1.0000
JC69	35	-1446.4851	2962.9701	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1394.2006	3042.8532	0.5025	0.5025
HKYI	40	-1394.2574	3042.9668	0.4747	0.9772
HKYIG	41	-1394.1155	3049.0443	0.0227	0.9999
K2PG	37	-1413.4051	3062.1784	0.0000	0.9999
K2PI	37	-1413.4879	3062.3440	0.0000	1.0000
GTRG	44	-1391.7805	3063.4584	0.0000	1.0000
GTRI	44	-1391.8806	3063.6586	0.0000	1.0000
K2PIG	38	-1413.2847	3068.2988	0.0000	1.0000
GTRIG	45	-1391.5265	3069.3117	0.0000	1.0000
HKY	39	-1412.0438	3072.1784	0.0000	1.0000
F81G	39	-1413.6712	3075.4333	0.0000	1.0000
F81I	39	-1413.7791	3075.6490	0.0000	1.0000
F81IG	40	-1413.5253	3081.5027	0.0000	1.0000
SYMG	41	-1411.4177	3083.6487	0.0000	1.0000
SYMI	41	-1411.5449	3083.9033	0.0000	1.0000
K2P	36	-1430.1841	3089.3751	0.0000	1.0000
SYMIG	42	-1411.2099	3089.5946	0.0000	1.0000
JC69G	36	-1431.0018	3091.0104	0.0000	1.0000
JC69I	36	-1431.1197	3091.2463	0.0000	1.0000
GTR	43	-1410.0355	3093.6070	0.0000	1.0000
JC69IG	37	-1430.8558	3097.0799	0.0000	1.0000
F81	38	-1429.4627	3100.6549	0.0000	1.0000
SYM	40	-1428.1061	3110.6643	0.0000	1.0000
JC69	35	-1446.4851	3115.6157	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_10370.phylip = 1 - 579;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10370.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_10370.phylip = 1 - 579;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10370.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_10370.phylip = 1 - 579;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10370.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

