
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:43 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/chr6_3430.phylip" (19 taxa, 1245 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	-4101.1385	8293.5770	0.6832	0.6832
GTRIG	45	-4100.8329	8295.1187	0.3160	0.9992
GTRI	44	-4107.9275	8307.1549	0.0008	1.0000
HKYG	40	-4116.7692	8316.2627	0.0000	1.0000
HKYIG	41	-4116.5640	8317.9908	0.0000	1.0000
HKYI	40	-4123.5890	8329.9023	0.0000	1.0000
GTR	43	-4165.1424	8419.4354	0.0000	1.0000
HKY	39	-4176.3865	8433.3622	0.0000	1.0000
F81G	39	-4188.6397	8457.8686	0.0000	1.0000
F81IG	40	-4188.5253	8459.7748	0.0000	1.0000
F81I	39	-4194.2886	8469.1664	0.0000	1.0000
SYMG	41	-4195.2654	8475.3937	0.0000	1.0000
SYMIG	42	-4195.1381	8477.2812	0.0000	1.0000
K2PG	37	-4205.5301	8487.3900	0.0000	1.0000
SYMI	41	-4201.5883	8488.0395	0.0000	1.0000
K2PIG	38	-4205.4003	8489.2583	0.0000	1.0000
K2PI	37	-4212.1273	8500.5843	0.0000	1.0000
F81	38	-4240.9495	8560.3567	0.0000	1.0000
SYM	40	-4252.1150	8586.9542	0.0000	1.0000
K2P	36	-4260.9852	8596.1757	0.0000	1.0000
JC69G	36	-4264.1849	8602.5752	0.0000	1.0000
JC69IG	37	-4264.1042	8604.5382	0.0000	1.0000
JC69I	36	-4269.3460	8612.8973	0.0000	1.0000
JC69	35	-4314.4378	8700.9600	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-4101.1385	8290.2770	0.6664	0.6664
GTRIG	45	-4100.8329	8291.6658	0.3328	0.9992
GTRI	44	-4107.9275	8303.8549	0.0008	1.0000
HKYG	40	-4116.7692	8313.5385	0.0000	1.0000
HKYIG	41	-4116.5640	8315.1279	0.0000	1.0000
HKYI	40	-4123.5890	8327.1780	0.0000	1.0000
GTR	43	-4165.1424	8416.2847	0.0000	1.0000
HKY	39	-4176.3865	8430.7729	0.0000	1.0000
F81G	39	-4188.6397	8455.2794	0.0000	1.0000
F81IG	40	-4188.5253	8457.0505	0.0000	1.0000
F81I	39	-4194.2886	8466.5772	0.0000	1.0000
SYMG	41	-4195.2654	8472.5308	0.0000	1.0000
SYMIG	42	-4195.1381	8474.2762	0.0000	1.0000
K2PG	37	-4205.5301	8485.0602	0.0000	1.0000
SYMI	41	-4201.5883	8485.1766	0.0000	1.0000
K2PIG	38	-4205.4003	8486.8006	0.0000	1.0000
K2PI	37	-4212.1273	8498.2546	0.0000	1.0000
F81	38	-4240.9495	8557.8990	0.0000	1.0000
SYM	40	-4252.1150	8584.2300	0.0000	1.0000
K2P	36	-4260.9852	8593.9704	0.0000	1.0000
JC69G	36	-4264.1849	8600.3699	0.0000	1.0000
JC69IG	37	-4264.1042	8602.2084	0.0000	1.0000
JC69I	36	-4269.3460	8610.6920	0.0000	1.0000
JC69	35	-4314.4378	8698.8756	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-4101.1385	8515.8602	0.7686	0.7686
HKYG	40	-4116.7692	8518.6141	0.1940	0.9626
GTRIG	45	-4100.8329	8522.3759	0.0296	0.9922
HKYIG	41	-4116.5640	8525.3305	0.0067	0.9989
GTRI	44	-4107.9275	8529.4381	0.0009	0.9998
HKYI	40	-4123.5890	8532.2537	0.0002	1.0000
HKY	39	-4176.3865	8630.7217	0.0000	1.0000
GTR	43	-4165.1424	8636.7410	0.0000	1.0000
F81G	39	-4188.6397	8655.2282	0.0000	1.0000
F81IG	40	-4188.5253	8662.1262	0.0000	1.0000
F81I	39	-4194.2886	8666.5259	0.0000	1.0000
K2PG	37	-4205.5301	8674.7552	0.0000	1.0000
K2PIG	38	-4205.4003	8681.6225	0.0000	1.0000
SYMG	41	-4195.2654	8682.7333	0.0000	1.0000
K2PI	37	-4212.1273	8687.9496	0.0000	1.0000
SYMIG	42	-4195.1381	8689.6056	0.0000	1.0000
SYMI	41	-4201.5883	8695.3791	0.0000	1.0000
F81	38	-4240.9495	8752.7209	0.0000	1.0000
K2P	36	-4260.9852	8778.5384	0.0000	1.0000
JC69G	36	-4264.1849	8784.9379	0.0000	1.0000
SYM	40	-4252.1150	8789.3056	0.0000	1.0000
JC69IG	37	-4264.1042	8791.9034	0.0000	1.0000
JC69I	36	-4269.3460	8795.2600	0.0000	1.0000
JC69	35	-4314.4378	8878.3168	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/chr6_3430.phylip = 1 - 1245;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3430.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/chr6_3430.phylip = 1 - 1245;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3430.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/chr6_3430.phylip = 1 - 1245;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3430.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

