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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	41	-3232.0108	6550.2630	0.6064	0.6064
HKYG	40	-3233.7251	6551.4845	0.3292	0.9356
GTRIG	45	-3230.1082	6555.3402	0.0479	0.9835
GTRG	44	-3232.2986	6557.4922	0.0163	0.9998
HKYI	40	-3241.1262	6566.2869	0.0002	1.0000
GTRI	44	-3239.8956	6572.6860	0.0000	1.0000
K2PIG	38	-3262.1027	6603.8421	0.0000	1.0000
K2PG	37	-3263.8626	6605.1714	0.0000	1.0000
SYMIG	42	-3260.0363	6608.5264	0.0000	1.0000
SYMG	41	-3262.2358	6610.7130	0.0000	1.0000
K2PI	37	-3272.7477	6622.9414	0.0000	1.0000
SYMI	41	-3269.0087	6624.2587	0.0000	1.0000
F81G	39	-3329.9239	6741.6808	0.0000	1.0000
F81IG	40	-3328.8497	6741.7337	0.0000	1.0000
F81I	39	-3334.3270	6750.4870	0.0000	1.0000
GTR	43	-3339.2830	6769.2376	0.0000	1.0000
HKY	39	-3344.0982	6770.0293	0.0000	1.0000
JC69IG	37	-3350.7620	6778.9700	0.0000	1.0000
JC69G	36	-3352.0545	6779.3696	0.0000	1.0000
JC69I	36	-3357.3883	6790.0372	0.0000	1.0000
K2P	36	-3367.8501	6810.9609	0.0000	1.0000
SYM	40	-3366.0078	6816.0501	0.0000	1.0000
F81	38	-3428.1306	6935.8980	0.0000	1.0000
JC69	35	-3449.6395	6972.3598	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-3232.0108	6546.0217	0.6061	0.6061
HKYG	40	-3233.7251	6547.4501	0.2967	0.9028
GTRIG	45	-3230.1082	6550.2164	0.0744	0.9772
GTRG	44	-3232.2986	6552.5972	0.0226	0.9998
HKYI	40	-3241.1262	6562.2525	0.0002	1.0000
GTRI	44	-3239.8956	6567.7911	0.0000	1.0000
K2PIG	38	-3262.1027	6600.2053	0.0000	1.0000
K2PG	37	-3263.8626	6601.7253	0.0000	1.0000
SYMIG	42	-3260.0363	6604.0727	0.0000	1.0000
SYMG	41	-3262.2358	6606.4716	0.0000	1.0000
K2PI	37	-3272.7477	6619.4953	0.0000	1.0000
SYMI	41	-3269.0087	6620.0174	0.0000	1.0000
F81IG	40	-3328.8497	6737.6993	0.0000	1.0000
F81G	39	-3329.9239	6737.8479	0.0000	1.0000
F81I	39	-3334.3270	6746.6541	0.0000	1.0000
GTR	43	-3339.2830	6764.5660	0.0000	1.0000
HKY	39	-3344.0982	6766.1964	0.0000	1.0000
JC69IG	37	-3350.7620	6775.5240	0.0000	1.0000
JC69G	36	-3352.0545	6776.1089	0.0000	1.0000
JC69I	36	-3357.3883	6786.7765	0.0000	1.0000
K2P	36	-3367.8501	6807.7002	0.0000	1.0000
SYM	40	-3366.0078	6812.0157	0.0000	1.0000
F81	38	-3428.1306	6932.2612	0.0000	1.0000
JC69	35	-3449.6395	6969.2791	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3233.7251	6737.4473	0.8399	0.8399
HKYIG	41	-3232.0108	6740.7688	0.1596	0.9995
HKYI	40	-3241.1262	6752.2497	0.0005	1.0000
GTRG	44	-3232.2986	6761.5942	0.0000	1.0000
GTRIG	45	-3230.1082	6763.9633	0.0000	1.0000
GTRI	44	-3239.8956	6776.7881	0.0000	1.0000
K2PG	37	-3263.8626	6777.4727	0.0000	1.0000
K2PIG	38	-3262.1027	6780.7027	0.0000	1.0000
K2PI	37	-3272.7477	6795.2428	0.0000	1.0000
SYMG	41	-3262.2358	6801.2188	0.0000	1.0000
SYMIG	42	-3260.0363	6803.5698	0.0000	1.0000
SYMI	41	-3269.0087	6814.7645	0.0000	1.0000
F81G	39	-3329.9239	6923.0952	0.0000	1.0000
F81IG	40	-3328.8497	6927.6965	0.0000	1.0000
F81I	39	-3334.3270	6931.9014	0.0000	1.0000
JC69G	36	-3352.0545	6947.1065	0.0000	1.0000
JC69IG	37	-3350.7620	6951.2714	0.0000	1.0000
HKY	39	-3344.0982	6951.4437	0.0000	1.0000
JC69I	36	-3357.3883	6957.7741	0.0000	1.0000
GTR	43	-3339.2830	6968.8131	0.0000	1.0000
K2P	36	-3367.8501	6978.6977	0.0000	1.0000
SYM	40	-3366.0078	7002.0129	0.0000	1.0000
F81	38	-3428.1306	7112.7586	0.0000	1.0000
JC69	35	-3449.6395	7135.5267	0.0000	1.0000



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

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

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


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

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

