
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:07 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/chr2_24792.phylip" (19 taxa, 952 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-3547.8966	7190.3628	0.6032	0.6032
GTRG	44	-3549.4187	7191.2035	0.3962	0.9993
GTRI	44	-3556.2559	7204.8779	0.0004	0.9998
HKYIG	41	-3560.5830	7206.9506	0.0002	0.9999
HKYG	40	-3562.2433	7208.0871	0.0001	1.0000
HKYI	40	-3568.3857	7220.3717	0.0000	1.0000
K2PIG	38	-3600.5856	7280.4176	0.0000	1.0000
SYMIG	42	-3596.2950	7280.5636	0.0000	1.0000
K2PG	37	-3602.4787	7282.0339	0.0000	1.0000
SYMG	41	-3598.1895	7282.1636	0.0000	1.0000
K2PI	37	-3609.6070	7296.2906	0.0000	1.0000
SYMI	41	-3606.1689	7298.1224	0.0000	1.0000
F81IG	40	-3610.3035	7304.2074	0.0000	1.0000
F81G	39	-3611.8286	7305.0782	0.0000	1.0000
F81I	39	-3617.0805	7315.5820	0.0000	1.0000
GTR	43	-3637.1032	7364.3737	0.0000	1.0000
JC69IG	37	-3646.2349	7369.5464	0.0000	1.0000
JC69G	36	-3647.9110	7370.7334	0.0000	1.0000
JC69I	36	-3654.6435	7384.1985	0.0000	1.0000
HKY	39	-3653.5072	7388.4355	0.0000	1.0000
SYM	40	-3693.3706	7470.3415	0.0000	1.0000
K2P	36	-3701.6621	7478.2358	0.0000	1.0000
F81	38	-3703.2766	7485.7997	0.0000	1.0000
JC69	35	-3743.5822	7559.9154	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-3547.8966	7185.7932	0.6273	0.6273
GTRG	44	-3549.4187	7186.8374	0.3722	0.9994
GTRI	44	-3556.2559	7200.5118	0.0004	0.9998
HKYIG	41	-3560.5830	7203.1660	0.0001	0.9999
HKYG	40	-3562.2433	7204.4867	0.0001	1.0000
HKYI	40	-3568.3857	7216.7713	0.0000	1.0000
SYMIG	42	-3596.2950	7276.5900	0.0000	1.0000
K2PIG	38	-3600.5856	7277.1711	0.0000	1.0000
SYMG	41	-3598.1895	7278.3790	0.0000	1.0000
K2PG	37	-3602.4787	7278.9573	0.0000	1.0000
K2PI	37	-3609.6070	7293.2140	0.0000	1.0000
SYMI	41	-3606.1689	7294.3377	0.0000	1.0000
F81IG	40	-3610.3035	7300.6069	0.0000	1.0000
F81G	39	-3611.8286	7301.6572	0.0000	1.0000
F81I	39	-3617.0805	7312.1609	0.0000	1.0000
GTR	43	-3637.1032	7360.2063	0.0000	1.0000
JC69IG	37	-3646.2349	7366.4698	0.0000	1.0000
JC69G	36	-3647.9110	7367.8220	0.0000	1.0000
JC69I	36	-3654.6435	7381.2871	0.0000	1.0000
HKY	39	-3653.5072	7385.0145	0.0000	1.0000
SYM	40	-3693.3706	7466.7411	0.0000	1.0000
K2P	36	-3701.6621	7475.3243	0.0000	1.0000
F81	38	-3703.2766	7482.5533	0.0000	1.0000
JC69	35	-3743.5822	7557.1643	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3562.2433	7398.8293	0.6084	0.6084
GTRG	44	-3549.4187	7400.6143	0.2492	0.8577
HKYIG	41	-3560.5830	7402.3671	0.1037	0.9614
GTRIG	45	-3547.8966	7404.4287	0.0370	0.9984
HKYI	40	-3568.3857	7411.1139	0.0013	0.9997
GTRI	44	-3556.2559	7414.2887	0.0003	1.0000
K2PG	37	-3602.4787	7458.7242	0.0000	1.0000
K2PIG	38	-3600.5856	7461.7966	0.0000	1.0000
K2PI	37	-3609.6070	7472.9809	0.0000	1.0000
SYMG	41	-3598.1895	7477.5802	0.0000	1.0000
SYMIG	42	-3596.2950	7480.6497	0.0000	1.0000
F81G	39	-3611.8286	7491.1412	0.0000	1.0000
SYMI	41	-3606.1689	7493.5389	0.0000	1.0000
F81IG	40	-3610.3035	7494.9495	0.0000	1.0000
F81I	39	-3617.0805	7501.6450	0.0000	1.0000
JC69G	36	-3647.9110	7542.7303	0.0000	1.0000
JC69IG	37	-3646.2349	7546.2367	0.0000	1.0000
JC69I	36	-3654.6435	7556.1954	0.0000	1.0000
GTR	43	-3637.1032	7569.1246	0.0000	1.0000
HKY	39	-3653.5072	7574.4985	0.0000	1.0000
K2P	36	-3701.6621	7650.2326	0.0000	1.0000
SYM	40	-3693.3706	7661.0837	0.0000	1.0000
F81	38	-3703.2766	7667.1787	0.0000	1.0000
JC69	35	-3743.5822	7727.2141	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24792.phylip = 1 - 952;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24792.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24792.phylip = 1 - 952;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24792.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.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/chr2_24792.phylip = 1 - 952;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24792.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

