
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:13 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_11196.phylip" (19 taxa, 759 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	-2644.5379	5373.6440	0.3782	0.3782
HKYIG	41	-2643.7775	5374.3583	0.2646	0.6428
GTRG	44	-2640.6507	5374.8476	0.2072	0.8500
GTRIG	45	-2639.9168	5375.6400	0.1394	0.9894
HKYI	40	-2648.5713	5381.7108	0.0067	0.9961
GTRI	44	-2644.6176	5382.7815	0.0039	1.0000
SYMG	41	-2655.2420	5397.2874	0.0000	1.0000
SYMIG	42	-2654.6583	5398.3614	0.0000	1.0000
SYMI	41	-2659.8897	5406.5827	0.0000	1.0000
K2PG	37	-2665.0605	5408.0211	0.0000	1.0000
K2PIG	38	-2664.4422	5409.0011	0.0000	1.0000
K2PI	37	-2669.2718	5416.4437	0.0000	1.0000
F81G	39	-2689.1004	5460.5401	0.0000	1.0000
F81IG	40	-2688.3647	5461.2977	0.0000	1.0000
F81I	39	-2693.7851	5469.9095	0.0000	1.0000
GTR	43	-2690.1942	5471.6807	0.0000	1.0000
HKY	39	-2697.3980	5477.1354	0.0000	1.0000
JC69G	36	-2705.1851	5486.0600	0.0000	1.0000
JC69IG	37	-2704.5518	5487.0038	0.0000	1.0000
SYM	40	-2702.4036	5489.3754	0.0000	1.0000
JC69I	36	-2709.6740	5495.0377	0.0000	1.0000
K2P	36	-2712.7578	5501.2054	0.0000	1.0000
F81	38	-2737.2799	5554.6764	0.0000	1.0000
JC69	35	-2752.1835	5577.8524	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2644.5379	5369.0758	0.2942	0.2942
GTRG	44	-2640.6507	5369.3014	0.2628	0.5569
HKYIG	41	-2643.7775	5369.5550	0.2315	0.7884
GTRIG	45	-2639.9168	5369.8336	0.2014	0.9898
HKYI	40	-2648.5713	5377.1425	0.0052	0.9950
GTRI	44	-2644.6176	5377.2352	0.0050	1.0000
SYMG	41	-2655.2420	5392.4840	0.0000	1.0000
SYMIG	42	-2654.6583	5393.3167	0.0000	1.0000
SYMI	41	-2659.8897	5401.7793	0.0000	1.0000
K2PG	37	-2665.0605	5404.1210	0.0000	1.0000
K2PIG	38	-2664.4422	5404.8845	0.0000	1.0000
K2PI	37	-2669.2718	5412.5436	0.0000	1.0000
F81G	39	-2689.1004	5456.2007	0.0000	1.0000
F81IG	40	-2688.3647	5456.7295	0.0000	1.0000
F81I	39	-2693.7851	5465.5701	0.0000	1.0000
GTR	43	-2690.1942	5466.3884	0.0000	1.0000
HKY	39	-2697.3980	5472.7961	0.0000	1.0000
JC69G	36	-2705.1851	5482.3702	0.0000	1.0000
JC69IG	37	-2704.5518	5483.1037	0.0000	1.0000
SYM	40	-2702.4036	5484.8072	0.0000	1.0000
JC69I	36	-2709.6740	5491.3479	0.0000	1.0000
K2P	36	-2712.7578	5497.5156	0.0000	1.0000
F81	38	-2737.2799	5550.5597	0.0000	1.0000
JC69	35	-2752.1835	5574.3670	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2644.5379	5554.3559	0.9128	0.9128
HKYIG	41	-2643.7775	5559.4671	0.0709	0.9837
HKYI	40	-2648.5713	5562.4226	0.0162	0.9999
GTRG	44	-2640.6507	5573.1095	0.0001	1.0000
K2PG	37	-2665.0605	5575.5051	0.0000	1.0000
GTRIG	45	-2639.9168	5578.2737	0.0000	1.0000
K2PIG	38	-2664.4422	5580.9005	0.0000	1.0000
GTRI	44	-2644.6176	5581.0433	0.0000	1.0000
SYMG	41	-2655.2420	5582.3961	0.0000	1.0000
K2PI	37	-2669.2718	5583.9277	0.0000	1.0000
SYMIG	42	-2654.6583	5587.8608	0.0000	1.0000
SYMI	41	-2659.8897	5591.6914	0.0000	1.0000
F81G	39	-2689.1004	5636.8488	0.0000	1.0000
F81IG	40	-2688.3647	5642.0095	0.0000	1.0000
F81I	39	-2693.7851	5646.2182	0.0000	1.0000
JC69G	36	-2705.1851	5649.1223	0.0000	1.0000
HKY	39	-2697.3980	5653.4441	0.0000	1.0000
JC69IG	37	-2704.5518	5654.4877	0.0000	1.0000
JC69I	36	-2709.6740	5658.1000	0.0000	1.0000
K2P	36	-2712.7578	5664.2677	0.0000	1.0000
GTR	43	-2690.1942	5665.5645	0.0000	1.0000
SYM	40	-2702.4036	5670.0872	0.0000	1.0000
F81	38	-2737.2799	5726.5758	0.0000	1.0000
JC69	35	-2752.1835	5736.4870	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/chr1_11196.phylip = 1 - 759;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_11196.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/chr1_11196.phylip = 1 - 759;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_11196.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/chr1_11196.phylip = 1 - 759;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_11196.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

