
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:45 04/04/2011

Input data from file "591-tomodel-phylip/chr6_1350.phylip" (7 taxa, 1148 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-3619.5427	7279.0853	0.5822	0.5822
GTRIG	21	-3619.5498	7281.0996	0.2127	0.7948
GTRG	20	-3620.5950	7281.1900	0.2033	0.9981
HKYI	16	-3629.8118	7291.6236	0.0011	0.9992
HKYIG	17	-3629.8221	7293.6443	0.0004	0.9996
HKYG	16	-3630.8352	7293.6703	0.0004	1.0000
GTR	19	-3645.3734	7328.7469	0.0000	1.0000
HKY	15	-3654.6618	7339.3237	0.0000	1.0000
SYMI	17	-3669.4270	7372.8541	0.0000	1.0000
K2PI	13	-3673.6060	7373.2119	0.0000	1.0000
SYMIG	18	-3669.4419	7374.8838	0.0000	1.0000
SYMG	17	-3670.6097	7375.2194	0.0000	1.0000
K2PIG	14	-3673.6195	7375.2391	0.0000	1.0000
K2PG	13	-3674.6871	7375.3743	0.0000	1.0000
SYM	16	-3693.7167	7419.4334	0.0000	1.0000
K2P	12	-3698.2727	7420.5455	0.0000	1.0000
F81I	15	-3746.3066	7522.6132	0.0000	1.0000
F81G	15	-3747.2863	7524.5727	0.0000	1.0000
F81IG	16	-3746.3213	7524.6426	0.0000	1.0000
F81	14	-3767.4266	7562.8531	0.0000	1.0000
JC69I	12	-3779.5724	7583.1447	0.0000	1.0000
JC69G	12	-3780.5100	7585.0199	0.0000	1.0000
JC69IG	13	-3779.5857	7585.1713	0.0000	1.0000
JC69	11	-3800.8602	7623.7205	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3629.8118	7372.3560	0.7048	0.7048
HKYG	16	-3630.8352	7374.4027	0.2533	0.9582
HKYIG	17	-3629.8221	7379.4225	0.0206	0.9787
GTRI	20	-3619.5427	7380.0009	0.0154	0.9942
GTRG	20	-3620.5950	7382.1055	0.0054	0.9995
GTRIG	21	-3619.5498	7387.0609	0.0005	1.0000
HKY	15	-3654.6618	7415.0103	0.0000	1.0000
GTR	19	-3645.3734	7424.6166	0.0000	1.0000
K2PI	13	-3673.6060	7438.8070	0.0000	1.0000
K2PG	13	-3674.6871	7440.9694	0.0000	1.0000
K2PIG	14	-3673.6195	7445.8800	0.0000	1.0000
SYMI	17	-3669.4270	7458.6323	0.0000	1.0000
SYMG	17	-3670.6097	7460.9976	0.0000	1.0000
SYMIG	18	-3669.4419	7465.7078	0.0000	1.0000
K2P	12	-3698.2727	7481.0948	0.0000	1.0000
SYM	16	-3693.7167	7500.1658	0.0000	1.0000
F81I	15	-3746.3066	7598.2998	0.0000	1.0000
F81G	15	-3747.2863	7600.2593	0.0000	1.0000
F81IG	16	-3746.3213	7605.3750	0.0000	1.0000
F81	14	-3767.4266	7633.4940	0.0000	1.0000
JC69I	12	-3779.5724	7643.6940	0.0000	1.0000
JC69G	12	-3780.5100	7645.5692	0.0000	1.0000
JC69IG	13	-3779.5857	7650.7664	0.0000	1.0000
JC69	11	-3800.8602	7679.2240	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_1350.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr6_1350.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_1350.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr6_1350.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_1350.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr6_1350.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

