
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:23 08/01/2011

Input data from file "261-aggressive-phylip/1553_hg18.phylip" (29 taxa, 541 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-2228.1775	4601.8341	0.6247	0.6247
GTRIG	65	-2227.8724	4603.8079	0.2329	0.8576
HKYG	60	-2235.0887	4605.4275	0.1036	0.9612
HKYIG	61	-2234.8822	4607.5557	0.0357	0.9970
GTRI	64	-2233.8488	4613.1765	0.0022	0.9991
HKYI	60	-2239.8540	4614.9580	0.0009	1.0000
K2PG	57	-2260.3284	4648.3463	0.0000	1.0000
K2PIG	58	-2259.9228	4650.0448	0.0000	1.0000
SYMG	61	-2257.3484	4652.4879	0.0000	1.0000
SYMIG	62	-2256.8979	4654.1390	0.0000	1.0000
K2PI	57	-2265.1250	4657.9394	0.0000	1.0000
SYMI	61	-2262.3154	4662.4221	0.0000	1.0000
GTR	63	-2278.5041	4699.9139	0.0000	1.0000
HKY	59	-2286.2573	4705.2340	0.0000	1.0000
K2P	56	-2313.1031	4751.3962	0.0000	1.0000
SYM	60	-2310.1307	4755.5113	0.0000	1.0000
F81G	59	-2323.1454	4779.0102	0.0000	1.0000
F81IG	60	-2322.9174	4781.0849	0.0000	1.0000
F81I	59	-2327.5109	4787.7412	0.0000	1.0000
JC69G	56	-2344.3477	4813.8856	0.0000	1.0000
JC69IG	57	-2344.0588	4815.8070	0.0000	1.0000
JC69I	56	-2348.7543	4822.6987	0.0000	1.0000
F81	58	-2371.9418	4874.0827	0.0000	1.0000
JC69	55	-2394.4698	4911.6406	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2228.1775	4584.3551	0.6321	0.6321
GTRIG	65	-2227.8724	4585.7448	0.3155	0.9476
HKYG	60	-2235.0887	4590.1775	0.0344	0.9820
HKYIG	61	-2234.8822	4591.7645	0.0156	0.9975
GTRI	64	-2233.8488	4595.6975	0.0022	0.9997
HKYI	60	-2239.8540	4599.7080	0.0003	1.0000
K2PG	57	-2260.3284	4634.6569	0.0000	1.0000
K2PIG	58	-2259.9228	4635.8456	0.0000	1.0000
SYMG	61	-2257.3484	4636.6967	0.0000	1.0000
SYMIG	62	-2256.8979	4637.7959	0.0000	1.0000
K2PI	57	-2265.1250	4644.2500	0.0000	1.0000
SYMI	61	-2262.3154	4646.6309	0.0000	1.0000
GTR	63	-2278.5041	4683.0083	0.0000	1.0000
HKY	59	-2286.2573	4690.5147	0.0000	1.0000
K2P	56	-2313.1031	4738.2061	0.0000	1.0000
SYM	60	-2310.1307	4740.2613	0.0000	1.0000
F81G	59	-2323.1454	4764.2908	0.0000	1.0000
F81IG	60	-2322.9174	4765.8349	0.0000	1.0000
F81I	59	-2327.5109	4773.0219	0.0000	1.0000
JC69G	56	-2344.3477	4800.6955	0.0000	1.0000
JC69IG	57	-2344.0588	4802.1176	0.0000	1.0000
JC69I	56	-2348.7543	4809.5087	0.0000	1.0000
F81	58	-2371.9418	4859.8835	0.0000	1.0000
JC69	55	-2394.4698	4898.9395	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2235.0887	4847.7826	0.9390	0.9390
HKYIG	61	-2234.8822	4853.6631	0.0496	0.9886
HKYI	60	-2239.8540	4857.3132	0.0080	0.9966
GTRG	64	-2228.1775	4859.1339	0.0032	0.9998
GTRIG	65	-2227.8724	4864.8170	0.0002	1.0000
GTRI	64	-2233.8488	4870.4763	0.0000	1.0000
K2PG	57	-2260.3284	4879.3818	0.0000	1.0000
K2PIG	58	-2259.9228	4884.8639	0.0000	1.0000
K2PI	57	-2265.1250	4888.9749	0.0000	1.0000
SYMG	61	-2257.3484	4898.5953	0.0000	1.0000
SYMIG	62	-2256.8979	4903.9879	0.0000	1.0000
SYMI	61	-2262.3154	4908.5294	0.0000	1.0000
HKY	59	-2286.2573	4943.8264	0.0000	1.0000
GTR	63	-2278.5041	4953.4937	0.0000	1.0000
K2P	56	-2313.1031	4978.6376	0.0000	1.0000
SYM	60	-2310.1307	4997.8665	0.0000	1.0000
F81G	59	-2323.1454	5017.6026	0.0000	1.0000
F81IG	60	-2322.9174	5023.4401	0.0000	1.0000
F81I	59	-2327.5109	5026.3336	0.0000	1.0000
JC69G	56	-2344.3477	5041.1270	0.0000	1.0000
JC69IG	57	-2344.0588	5046.8425	0.0000	1.0000
JC69I	56	-2348.7543	5049.9401	0.0000	1.0000
F81	58	-2371.9418	5108.9019	0.0000	1.0000
JC69	55	-2394.4698	5135.0776	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 261-aggressive-phylip/1553_hg18.phylip = 1 - 541;
 Partition Dummy = 1:261-aggressive-phylip/1553_hg18.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 261-aggressive-phylip/1553_hg18.phylip = 1 - 541;
 Partition Dummy = 1:261-aggressive-phylip/1553_hg18.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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1553_hg18.phylip = 1 - 541;
 Partition Dummy = 1:261-aggressive-phylip/1553_hg18.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

