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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-4153.2710	8340.5420	0.3212	0.3212
K2PG	13	-4157.6238	8341.2475	0.2257	0.5470
SYMIG	18	-4153.2769	8342.5538	0.1175	0.6644
GTRG	20	-4151.3620	8342.7240	0.1079	0.7723
K2PIG	14	-4157.6366	8343.2732	0.0820	0.8543
HKYG	16	-4155.9809	8343.9619	0.0581	0.9124
GTRIG	21	-4151.3663	8344.7327	0.0395	0.9519
HKYIG	17	-4155.9951	8345.9902	0.0211	0.9730
SYMI	17	-4156.4772	8346.9545	0.0130	0.9860
K2PI	13	-4160.9766	8347.9531	0.0079	0.9939
GTRI	20	-4154.5954	8349.1909	0.0043	0.9982
HKYI	16	-4159.4406	8350.8812	0.0018	1.0000
K2P	12	-4197.4262	8418.8523	0.0000	1.0000
HKY	15	-4195.4526	8420.9051	0.0000	1.0000
SYM	16	-4194.6258	8421.2517	0.0000	1.0000
GTR	19	-4192.5372	8423.0745	0.0000	1.0000
JC69G	12	-4263.1623	8550.3246	0.0000	1.0000
JC69IG	13	-4263.0616	8552.1232	0.0000	1.0000
F81G	15	-4261.1041	8552.2082	0.0000	1.0000
JC69I	12	-4264.3724	8552.7448	0.0000	1.0000
F81IG	16	-4261.0123	8554.0247	0.0000	1.0000
F81I	15	-4262.4005	8554.8010	0.0000	1.0000
JC69	11	-4297.4579	8616.9159	0.0000	1.0000
F81	14	-4295.0372	8618.0743	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-4157.6238	8406.0489	0.9387	0.9387
K2PI	13	-4160.9766	8412.7544	0.0328	0.9716
K2PIG	14	-4157.6366	8413.0592	0.0282	0.9998
HKYG	16	-4155.9809	8423.7173	0.0001	0.9999
SYMG	17	-4153.2710	8425.2822	0.0001	1.0000
HKYI	16	-4159.4406	8430.6367	0.0000	1.0000
HKYIG	17	-4155.9951	8430.7304	0.0000	1.0000
SYMI	17	-4156.4772	8431.6947	0.0000	1.0000
SYMIG	18	-4153.2769	8432.2787	0.0000	1.0000
GTRG	20	-4151.3620	8442.4184	0.0000	1.0000
GTRI	20	-4154.5954	8448.8852	0.0000	1.0000
GTRIG	21	-4151.3663	8449.4117	0.0000	1.0000
K2P	12	-4197.4262	8478.6689	0.0000	1.0000
HKY	15	-4195.4526	8495.6759	0.0000	1.0000
SYM	16	-4194.6258	8501.0071	0.0000	1.0000
GTR	19	-4192.5372	8517.7841	0.0000	1.0000
JC69G	12	-4263.1623	8610.1412	0.0000	1.0000
JC69I	12	-4264.3724	8612.5614	0.0000	1.0000
JC69IG	13	-4263.0616	8616.9245	0.0000	1.0000
F81G	15	-4261.1041	8626.9790	0.0000	1.0000
F81I	15	-4262.4005	8629.5717	0.0000	1.0000
F81IG	16	-4261.0123	8633.7802	0.0000	1.0000
JC69	11	-4297.4579	8671.7477	0.0000	1.0000
F81	14	-4295.0372	8687.8603	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_16786.phylip = 1 - 1080;
 Partition Dummy = 1:591-tomodel-phylip/chr3_16786.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_16786.phylip = 1 - 1080;
 Partition Dummy = 1:591-tomodel-phylip/chr3_16786.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_16786.phylip = 1 - 1080;
 Partition Dummy = 1:591-tomodel-phylip/chr3_16786.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

