
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:32 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_3232.phylip" (19 taxa, 707 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	-3298.9294	6691.8407	0.7623	0.7623
GTRIG	45	-3298.9546	6694.1724	0.2376	0.9998
GTRI	44	-3307.4526	6708.8870	0.0002	1.0000
HKYG	40	-3317.3764	6719.6778	0.0000	1.0000
SYMG	41	-3316.8033	6720.7856	0.0000	1.0000
HKYIG	41	-3317.4025	6721.9839	0.0000	1.0000
SYMIG	42	-3316.8030	6723.0458	0.0000	1.0000
HKYI	40	-3324.3937	6733.7123	0.0000	1.0000
SYMI	41	-3325.3877	6737.9544	0.0000	1.0000
K2PG	37	-3333.6185	6745.4402	0.0000	1.0000
K2PIG	38	-3333.6437	6747.7246	0.0000	1.0000
K2PI	37	-3340.6500	6759.5033	0.0000	1.0000
GTR	43	-3371.8794	6835.4662	0.0000	1.0000
HKY	39	-3385.6609	6853.9994	0.0000	1.0000
SYM	40	-3386.0078	6856.9406	0.0000	1.0000
K2P	36	-3400.4452	6876.8666	0.0000	1.0000
F81G	39	-3415.0568	6912.7912	0.0000	1.0000
F81IG	40	-3417.4751	6919.8751	0.0000	1.0000
F81I	39	-3422.8368	6928.3512	0.0000	1.0000
JC69G	36	-3431.1256	6938.2272	0.0000	1.0000
JC69IG	37	-3433.6021	6945.4074	0.0000	1.0000
JC69I	36	-3438.8636	6953.7034	0.0000	1.0000
F81	38	-3479.6346	7039.7063	0.0000	1.0000
JC69	35	-3493.9284	7061.6125	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3298.9294	6685.8588	0.7359	0.7359
GTRIG	45	-3298.9546	6687.9091	0.2640	0.9999
GTRI	44	-3307.4526	6702.9052	0.0001	1.0000
HKYG	40	-3317.3764	6714.7528	0.0000	1.0000
SYMG	41	-3316.8033	6715.6067	0.0000	1.0000
HKYIG	41	-3317.4025	6716.8049	0.0000	1.0000
SYMIG	42	-3316.8030	6717.6060	0.0000	1.0000
HKYI	40	-3324.3937	6728.7874	0.0000	1.0000
SYMI	41	-3325.3877	6732.7754	0.0000	1.0000
K2PG	37	-3333.6185	6741.2369	0.0000	1.0000
K2PIG	38	-3333.6437	6743.2875	0.0000	1.0000
K2PI	37	-3340.6500	6755.3000	0.0000	1.0000
GTR	43	-3371.8794	6829.7588	0.0000	1.0000
HKY	39	-3385.6609	6849.3218	0.0000	1.0000
SYM	40	-3386.0078	6852.0157	0.0000	1.0000
K2P	36	-3400.4452	6872.8905	0.0000	1.0000
F81G	39	-3415.0568	6908.1136	0.0000	1.0000
F81IG	40	-3417.4751	6914.9502	0.0000	1.0000
F81I	39	-3422.8368	6923.6735	0.0000	1.0000
JC69G	36	-3431.1256	6934.2511	0.0000	1.0000
JC69IG	37	-3433.6021	6941.2042	0.0000	1.0000
JC69I	36	-3438.8636	6949.7273	0.0000	1.0000
F81	38	-3479.6346	7035.2691	0.0000	1.0000
JC69	35	-3493.9284	7057.8569	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3298.9294	6886.5441	0.9594	0.9594
GTRIG	45	-3298.9546	6893.1555	0.0352	0.9946
HKYG	40	-3317.3764	6897.1941	0.0047	0.9993
SYMG	41	-3316.8033	6902.6089	0.0003	0.9996
GTRI	44	-3307.4526	6903.5905	0.0002	0.9998
HKYIG	41	-3317.4025	6903.8072	0.0002	1.0000
SYMIG	42	-3316.8030	6909.1693	0.0000	1.0000
K2PG	37	-3333.6185	6909.9951	0.0000	1.0000
HKYI	40	-3324.3937	6911.2286	0.0000	1.0000
K2PIG	38	-3333.6437	6916.6066	0.0000	1.0000
SYMI	41	-3325.3877	6919.7777	0.0000	1.0000
K2PI	37	-3340.6500	6924.0582	0.0000	1.0000
GTR	43	-3371.8794	7025.8831	0.0000	1.0000
HKY	39	-3385.6609	7027.2020	0.0000	1.0000
SYM	40	-3386.0078	7034.4569	0.0000	1.0000
K2P	36	-3400.4452	7037.0876	0.0000	1.0000
F81G	39	-3415.0568	7085.9938	0.0000	1.0000
F81IG	40	-3417.4751	7097.3914	0.0000	1.0000
JC69G	36	-3431.1256	7098.4482	0.0000	1.0000
F81I	39	-3422.8368	7101.5537	0.0000	1.0000
JC69IG	37	-3433.6021	7109.9623	0.0000	1.0000
JC69I	36	-3438.8636	7113.9244	0.0000	1.0000
F81	38	-3479.6346	7208.5883	0.0000	1.0000
JC69	35	-3493.9284	7217.4930	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_3232.phylip = 1 - 707;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3232.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_3232.phylip = 1 - 707;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3232.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_3232.phylip = 1 - 707;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3232.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

