
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:28 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/chr4_9746.phylip" (19 taxa, 1017 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	-3737.8051	7569.8739	0.5455	0.5455
GTRG	44	-3739.0987	7570.2714	0.4472	0.9927
HKYIG	41	-3747.1430	7579.8183	0.0038	0.9965
HKYG	40	-3748.3561	7580.0728	0.0033	0.9998
GTRI	44	-3747.1454	7586.3649	0.0001	0.9999
SYMIG	42	-3751.0391	7589.7866	0.0000	1.0000
SYMG	41	-3752.2805	7590.0933	0.0000	1.0000
HKYI	40	-3755.1530	7593.6667	0.0000	1.0000
K2PIG	38	-3759.0303	7597.0913	0.0000	1.0000
K2PG	37	-3760.2591	7597.3905	0.0000	1.0000
SYMI	41	-3759.0337	7603.5997	0.0000	1.0000
K2PI	37	-3766.8321	7610.5364	0.0000	1.0000
GTR	43	-3824.9142	7739.7174	0.0000	1.0000
HKY	39	-3832.3348	7745.8630	0.0000	1.0000
SYM	40	-3836.3823	7756.1252	0.0000	1.0000
K2P	36	-3844.6898	7764.0980	0.0000	1.0000
F81IG	40	-3848.8447	7781.0501	0.0000	1.0000
F81G	39	-3849.9456	7781.0847	0.0000	1.0000
F81I	39	-3856.4077	7794.0088	0.0000	1.0000
JC69IG	37	-3861.8528	7800.5779	0.0000	1.0000
JC69G	36	-3862.9558	7800.6301	0.0000	1.0000
JC69I	36	-3869.1323	7812.9830	0.0000	1.0000
F81	38	-3930.1531	7939.3369	0.0000	1.0000
JC69	35	-3942.0764	7956.7217	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-3737.8051	7565.6103	0.5699	0.5699
GTRG	44	-3739.0987	7566.1973	0.4249	0.9949
HKYIG	41	-3747.1430	7576.2860	0.0027	0.9976
HKYG	40	-3748.3561	7576.7121	0.0022	0.9998
GTRI	44	-3747.1454	7582.2908	0.0001	1.0000
SYMIG	42	-3751.0391	7586.0781	0.0000	1.0000
SYMG	41	-3752.2805	7586.5610	0.0000	1.0000
HKYI	40	-3755.1530	7590.3060	0.0000	1.0000
K2PIG	38	-3759.0303	7594.0607	0.0000	1.0000
K2PG	37	-3760.2591	7594.5182	0.0000	1.0000
SYMI	41	-3759.0337	7600.0674	0.0000	1.0000
K2PI	37	-3766.8321	7607.6641	0.0000	1.0000
GTR	43	-3824.9142	7735.8284	0.0000	1.0000
HKY	39	-3832.3348	7742.6696	0.0000	1.0000
SYM	40	-3836.3823	7752.7646	0.0000	1.0000
K2P	36	-3844.6898	7761.3797	0.0000	1.0000
F81IG	40	-3848.8447	7777.6894	0.0000	1.0000
F81G	39	-3849.9456	7777.8912	0.0000	1.0000
F81I	39	-3856.4077	7790.8154	0.0000	1.0000
JC69IG	37	-3861.8528	7797.7056	0.0000	1.0000
JC69G	36	-3862.9558	7797.9117	0.0000	1.0000
JC69I	36	-3869.1323	7810.2646	0.0000	1.0000
F81	38	-3930.1531	7936.3063	0.0000	1.0000
JC69	35	-3942.0764	7954.1529	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3748.3561	7773.6966	0.7342	0.7342
K2PG	37	-3760.2591	7776.7288	0.1612	0.8954
HKYIG	41	-3747.1430	7778.1951	0.0774	0.9729
K2PIG	38	-3759.0303	7781.1959	0.0173	0.9902
GTRG	44	-3739.0987	7782.8803	0.0074	0.9976
GTRIG	45	-3737.8051	7787.2178	0.0009	0.9984
HKYI	40	-3755.1530	7787.2905	0.0008	0.9993
SYMG	41	-3752.2805	7788.4701	0.0005	0.9997
K2PI	37	-3766.8321	7789.8748	0.0002	0.9999
SYMIG	42	-3751.0391	7792.9119	0.0000	1.0000
GTRI	44	-3747.1454	7798.9738	0.0000	1.0000
SYMI	41	-3759.0337	7801.9765	0.0000	1.0000
HKY	39	-3832.3348	7934.7295	0.0000	1.0000
K2P	36	-3844.6898	7938.6657	0.0000	1.0000
GTR	43	-3824.9142	7947.5867	0.0000	1.0000
SYM	40	-3836.3823	7949.7491	0.0000	1.0000
F81G	39	-3849.9456	7969.9511	0.0000	1.0000
F81IG	40	-3848.8447	7974.6739	0.0000	1.0000
JC69G	36	-3862.9558	7975.1977	0.0000	1.0000
JC69IG	37	-3861.8528	7979.9162	0.0000	1.0000
F81I	39	-3856.4077	7982.8753	0.0000	1.0000
JC69I	36	-3869.1323	7987.5507	0.0000	1.0000
F81	38	-3930.1531	8123.4415	0.0000	1.0000
JC69	35	-3942.0764	8126.5143	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/chr4_9746.phylip = 1 - 1017;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_9746.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/chr4_9746.phylip = 1 - 1017;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_9746.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/chr4_9746.phylip = 1 - 1017;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_9746.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

