
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:35 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_9464.phylip" (19 taxa, 440 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	-1547.4893	3193.0039	0.6794	0.6794
GTRIG	45	-1547.3215	3195.1506	0.2323	0.9116
GTRI	44	-1549.5320	3197.0894	0.0881	0.9997
SYMG	41	-1559.4457	3209.5447	0.0002	0.9999
SYMIG	42	-1559.2106	3211.5194	0.0001	1.0000
SYMI	41	-1561.4202	3213.4936	0.0000	1.0000
K2PG	37	-1574.5479	3230.0909	0.0000	1.0000
HKYG	40	-1571.2885	3230.7975	0.0000	1.0000
K2PIG	38	-1574.1560	3231.7035	0.0000	1.0000
HKYIG	41	-1570.8813	3232.4158	0.0000	1.0000
GTR	43	-1568.5986	3232.7527	0.0000	1.0000
K2PI	37	-1577.1108	3235.2167	0.0000	1.0000
HKYI	40	-1574.2217	3236.6640	0.0000	1.0000
SYM	40	-1579.9063	3248.0331	0.0000	1.0000
F81G	39	-1582.5321	3250.8643	0.0000	1.0000
JC69G	36	-1586.4196	3251.4496	0.0000	1.0000
F81IG	40	-1582.2933	3252.8072	0.0000	1.0000
JC69IG	37	-1585.9075	3252.8101	0.0000	1.0000
F81I	39	-1584.6510	3255.1020	0.0000	1.0000
JC69I	36	-1588.3914	3255.3933	0.0000	1.0000
K2P	36	-1597.4298	3273.4700	0.0000	1.0000
HKY	39	-1594.0247	3273.8493	0.0000	1.0000
F81	38	-1602.5020	3288.3954	0.0000	1.0000
JC69	35	-1606.6939	3289.6254	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1547.4893	3182.9786	0.6390	0.6390
GTRIG	45	-1547.3215	3184.6430	0.2780	0.9170
GTRI	44	-1549.5320	3187.0641	0.0829	0.9999
SYMG	41	-1559.4457	3200.8915	0.0001	1.0000
SYMIG	42	-1559.2106	3202.4211	0.0000	1.0000
SYMI	41	-1561.4202	3204.8403	0.0000	1.0000
HKYG	40	-1571.2885	3222.5770	0.0000	1.0000
K2PG	37	-1574.5479	3223.0958	0.0000	1.0000
GTR	43	-1568.5986	3223.1972	0.0000	1.0000
HKYIG	41	-1570.8813	3223.7626	0.0000	1.0000
K2PIG	38	-1574.1560	3224.3120	0.0000	1.0000
K2PI	37	-1577.1108	3228.2217	0.0000	1.0000
HKYI	40	-1574.2217	3228.4435	0.0000	1.0000
SYM	40	-1579.9063	3239.8125	0.0000	1.0000
F81G	39	-1582.5321	3243.0643	0.0000	1.0000
F81IG	40	-1582.2933	3244.5866	0.0000	1.0000
JC69G	36	-1586.4196	3244.8392	0.0000	1.0000
JC69IG	37	-1585.9075	3245.8150	0.0000	1.0000
F81I	39	-1584.6510	3247.3020	0.0000	1.0000
JC69I	36	-1588.3914	3248.7829	0.0000	1.0000
HKY	39	-1594.0247	3266.0493	0.0000	1.0000
K2P	36	-1597.4298	3266.8595	0.0000	1.0000
F81	38	-1602.5020	3281.0039	0.0000	1.0000
JC69	35	-1606.6939	3283.3878	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1547.4893	3362.7967	0.7932	0.7932
GTRI	44	-1549.5320	3366.8821	0.1029	0.8960
SYMG	41	-1559.4457	3368.4492	0.0470	0.9430
GTRIG	45	-1547.3215	3368.5479	0.0447	0.9878
SYMI	41	-1561.4202	3372.3981	0.0065	0.9943
SYMIG	42	-1559.2106	3374.0657	0.0028	0.9971
K2PG	37	-1574.5479	3374.3065	0.0025	0.9996
K2PI	37	-1577.1108	3379.4323	0.0002	0.9998
K2PIG	38	-1574.1560	3379.6094	0.0002	1.0000
HKYG	40	-1571.2885	3386.0479	0.0000	1.0000
HKYIG	41	-1570.8813	3391.3203	0.0000	1.0000
HKYI	40	-1574.2217	3391.9144	0.0000	1.0000
JC69G	36	-1586.4196	3391.9631	0.0000	1.0000
JC69I	36	-1588.3914	3395.9068	0.0000	1.0000
JC69IG	37	-1585.9075	3397.0257	0.0000	1.0000
GTR	43	-1568.5986	3398.9285	0.0000	1.0000
F81G	39	-1582.5321	3402.4485	0.0000	1.0000
SYM	40	-1579.9063	3403.2835	0.0000	1.0000
F81I	39	-1584.6510	3406.6862	0.0000	1.0000
F81IG	40	-1582.2933	3408.0576	0.0000	1.0000
K2P	36	-1597.4298	3413.9834	0.0000	1.0000
HKY	39	-1594.0247	3425.4336	0.0000	1.0000
JC69	35	-1606.6939	3426.4249	0.0000	1.0000
F81	38	-1602.5020	3436.3014	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_9464.phylip = 1 - 440;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_9464.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_9464.phylip = 1 - 440;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_9464.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_9464.phylip = 1 - 440;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_9464.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

