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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-2599.0608	5341.1125	0.8625	0.8625
GTRIG	65	-2599.6556	5344.7985	0.1366	0.9991
SYMG	61	-2609.7690	5355.0937	0.0008	0.9999
SYMIG	62	-2610.7521	5359.5294	0.0001	1.0000
HKYG	60	-2615.1280	5363.3509	0.0000	1.0000
HKYIG	61	-2614.4351	5364.4257	0.0000	1.0000
GTRI	64	-2611.4559	5365.9028	0.0000	1.0000
K2PG	57	-2620.5062	5366.7775	0.0000	1.0000
K2PIG	58	-2622.1593	5372.5183	0.0000	1.0000
SYMI	61	-2621.1732	5377.9019	0.0000	1.0000
HKYI	60	-2623.2121	5379.5190	0.0000	1.0000
K2PI	57	-2631.5869	5388.9389	0.0000	1.0000
GTR	63	-2659.3851	5459.2738	0.0000	1.0000
SYM	60	-2669.1623	5471.4194	0.0000	1.0000
F81G	59	-2670.8111	5472.2651	0.0000	1.0000
JC69G	56	-2674.7868	5472.9128	0.0000	1.0000
F81IG	60	-2670.2072	5473.5092	0.0000	1.0000
JC69IG	57	-2674.2322	5474.2294	0.0000	1.0000
HKY	59	-2673.9077	5478.4583	0.0000	1.0000
K2P	56	-2679.3644	5482.0681	0.0000	1.0000
F81I	59	-2679.6634	5489.9697	0.0000	1.0000
JC69I	56	-2683.7377	5490.8146	0.0000	1.0000
JC69	55	-2729.1804	5579.2828	0.0000	1.0000
F81	58	-2725.5506	5579.3008	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2599.0608	5326.1215	0.8309	0.8309
GTRIG	65	-2599.6556	5329.3111	0.1686	0.9996
SYMG	61	-2609.7690	5341.5381	0.0004	0.9999
SYMIG	62	-2610.7521	5345.5043	0.0001	1.0000
HKYG	60	-2615.1280	5350.2561	0.0000	1.0000
HKYIG	61	-2614.4351	5350.8701	0.0000	1.0000
GTRI	64	-2611.4559	5350.9118	0.0000	1.0000
K2PG	57	-2620.5062	5355.0124	0.0000	1.0000
K2PIG	58	-2622.1593	5360.3186	0.0000	1.0000
SYMI	61	-2621.1732	5364.3464	0.0000	1.0000
HKYI	60	-2623.2121	5366.4242	0.0000	1.0000
K2PI	57	-2631.5869	5377.1737	0.0000	1.0000
GTR	63	-2659.3851	5444.7702	0.0000	1.0000
SYM	60	-2669.1623	5458.3246	0.0000	1.0000
F81G	59	-2670.8111	5459.6223	0.0000	1.0000
F81IG	60	-2670.2072	5460.4143	0.0000	1.0000
JC69G	56	-2674.7868	5461.5735	0.0000	1.0000
JC69IG	57	-2674.2322	5462.4643	0.0000	1.0000
HKY	59	-2673.9077	5465.8154	0.0000	1.0000
K2P	56	-2679.3644	5470.7288	0.0000	1.0000
F81I	59	-2679.6634	5477.3269	0.0000	1.0000
JC69I	56	-2683.7377	5479.4754	0.0000	1.0000
F81	58	-2725.5506	5567.1011	0.0000	1.0000
JC69	55	-2729.1804	5568.3608	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2620.5062	5607.5064	0.6671	0.6671
GTRG	64	-2599.0608	5609.6236	0.2314	0.8985
SYMG	61	-2609.7690	5611.7510	0.0799	0.9784
HKYG	60	-2615.1280	5616.0392	0.0094	0.9878
K2PIG	58	-2622.1593	5617.2424	0.0051	0.9929
GTRIG	65	-2599.6556	5617.2429	0.0051	0.9980
SYMIG	62	-2610.7521	5620.1469	0.0012	0.9992
HKYIG	61	-2614.4351	5621.0830	0.0008	1.0000
K2PI	57	-2631.5869	5629.6678	0.0000	1.0000
HKYI	60	-2623.2121	5632.2074	0.0000	1.0000
GTRI	64	-2611.4559	5634.4139	0.0000	1.0000
SYMI	61	-2621.1732	5634.5592	0.0000	1.0000
JC69G	56	-2674.7868	5709.6378	0.0000	1.0000
JC69IG	57	-2674.2322	5714.9583	0.0000	1.0000
K2P	56	-2679.3644	5718.7931	0.0000	1.0000
F81G	59	-2670.8111	5720.9757	0.0000	1.0000
GTR	63	-2659.3851	5723.8425	0.0000	1.0000
SYM	60	-2669.1623	5724.1078	0.0000	1.0000
F81IG	60	-2670.2072	5726.1975	0.0000	1.0000
HKY	59	-2673.9077	5727.1688	0.0000	1.0000
JC69I	56	-2683.7377	5727.5397	0.0000	1.0000
F81I	59	-2679.6634	5738.6803	0.0000	1.0000
JC69	55	-2729.1804	5811.9954	0.0000	1.0000
F81	58	-2725.5506	5824.0249	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/7444_hg18.phylip = 1 - 620;
 Partition Dummy = 1:261-aggressive-phylip/7444_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/7444_hg18.phylip = 1 - 620;
 Partition Dummy = 1:261-aggressive-phylip/7444_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 (K2PG)]

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

