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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	44	-5974.4842	12040.5296	0.7185	0.7185
GTRIG	45	-5974.3385	12042.4033	0.2815	1.0000
HKYG	40	-5999.5366	12082.0122	0.0000	1.0000
HKYIG	41	-5999.1523	12083.3934	0.0000	1.0000
GTRI	44	-6004.7073	12100.9758	0.0000	1.0000
HKYI	40	-6023.0632	12129.0654	0.0000	1.0000
SYMG	41	-6035.2567	12155.6021	0.0000	1.0000
SYMIG	42	-6034.9125	12157.0673	0.0000	1.0000
K2PG	37	-6045.8507	12168.2143	0.0000	1.0000
K2PIG	38	-6045.4185	12169.4882	0.0000	1.0000
SYMI	41	-6056.8039	12198.6967	0.0000	1.0000
K2PI	37	-6066.9002	12210.3133	0.0000	1.0000
GTR	43	-6124.5815	12338.5628	0.0000	1.0000
F81G	39	-6143.9225	12368.6383	0.0000	1.0000
F81IG	40	-6143.5768	12370.0927	0.0000	1.0000
HKY	39	-6146.3052	12373.4036	0.0000	1.0000
F81I	39	-6161.3707	12403.5346	0.0000	1.0000
JC69G	36	-6175.2023	12424.7833	0.0000	1.0000
JC69IG	37	-6174.8647	12426.2423	0.0000	1.0000
SYM	40	-6184.9453	12452.8297	0.0000	1.0000
JC69I	36	-6190.7001	12455.7787	0.0000	1.0000
K2P	36	-6194.5681	12463.5148	0.0000	1.0000
F81	38	-6281.3445	12641.3402	0.0000	1.0000
JC69	35	-6313.4799	12699.2079	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-5974.4842	12036.9684	0.7015	0.7015
GTRIG	45	-5974.3385	12038.6770	0.2985	1.0000
HKYG	40	-5999.5366	12079.0732	0.0000	1.0000
HKYIG	41	-5999.1523	12080.3047	0.0000	1.0000
GTRI	44	-6004.7073	12097.4146	0.0000	1.0000
HKYI	40	-6023.0632	12126.1264	0.0000	1.0000
SYMG	41	-6035.2567	12152.5133	0.0000	1.0000
SYMIG	42	-6034.9125	12153.8249	0.0000	1.0000
K2PG	37	-6045.8507	12165.7014	0.0000	1.0000
K2PIG	38	-6045.4185	12166.8370	0.0000	1.0000
SYMI	41	-6056.8039	12195.6079	0.0000	1.0000
K2PI	37	-6066.9002	12207.8003	0.0000	1.0000
GTR	43	-6124.5815	12335.1630	0.0000	1.0000
F81G	39	-6143.9225	12365.8451	0.0000	1.0000
F81IG	40	-6143.5768	12367.1536	0.0000	1.0000
HKY	39	-6146.3052	12370.6104	0.0000	1.0000
F81I	39	-6161.3707	12400.7414	0.0000	1.0000
JC69G	36	-6175.2023	12422.4047	0.0000	1.0000
JC69IG	37	-6174.8647	12423.7294	0.0000	1.0000
SYM	40	-6184.9453	12449.8906	0.0000	1.0000
JC69I	36	-6190.7001	12453.4001	0.0000	1.0000
K2P	36	-6194.5681	12461.1362	0.0000	1.0000
F81	38	-6281.3445	12638.6890	0.0000	1.0000
JC69	35	-6313.4799	12696.9599	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-5974.4842	12259.3262	0.9671	0.9671
GTRIG	45	-5974.3385	12266.0883	0.0329	1.0000
HKYG	40	-5999.5366	12281.2166	0.0000	1.0000
HKYIG	41	-5999.1523	12287.5017	0.0000	1.0000
GTRI	44	-6004.7073	12319.7724	0.0000	1.0000
HKYI	40	-6023.0632	12328.2698	0.0000	1.0000
K2PG	37	-6045.8507	12352.6841	0.0000	1.0000
K2PIG	38	-6045.4185	12358.8733	0.0000	1.0000
SYMG	41	-6035.2567	12359.7103	0.0000	1.0000
SYMIG	42	-6034.9125	12366.0755	0.0000	1.0000
K2PI	37	-6066.9002	12394.7830	0.0000	1.0000
SYMI	41	-6056.8039	12402.8049	0.0000	1.0000
GTR	43	-6124.5815	12552.4672	0.0000	1.0000
F81G	39	-6143.9225	12562.9349	0.0000	1.0000
HKY	39	-6146.3052	12567.7002	0.0000	1.0000
F81IG	40	-6143.5768	12569.2970	0.0000	1.0000
F81I	39	-6161.3707	12597.8312	0.0000	1.0000
JC69G	36	-6175.2023	12604.3338	0.0000	1.0000
JC69IG	37	-6174.8647	12610.7121	0.0000	1.0000
JC69I	36	-6190.7001	12635.3292	0.0000	1.0000
K2P	36	-6194.5681	12643.0653	0.0000	1.0000
SYM	40	-6184.9453	12652.0340	0.0000	1.0000
F81	38	-6281.3445	12830.7253	0.0000	1.0000
JC69	35	-6313.4799	12873.8354	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3231.phylip = 1 - 1157;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3231.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3231.phylip = 1 - 1157;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3231.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3231.phylip = 1 - 1157;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3231.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

