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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
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
HKYG	60	-3319.6106	6766.6003	0.5615	0.5615
GTRG	64	-3315.4830	6767.3870	0.3789	0.9404
HKYIG	61	-3321.2686	6772.1698	0.0347	0.9751
GTRIG	65	-3317.0855	6772.8639	0.0245	0.9996
HKYI	60	-3327.1788	6781.7367	0.0003	0.9999
GTRI	64	-3323.4882	6783.3974	0.0001	1.0000
K2PG	57	-3345.0621	6810.7694	0.0000	1.0000
SYMG	61	-3342.8479	6815.3286	0.0000	1.0000
K2PIG	58	-3346.6974	6816.2801	0.0000	1.0000
SYMIG	62	-3344.5010	6820.8930	0.0000	1.0000
K2PI	57	-3352.7335	6826.1121	0.0000	1.0000
SYMI	61	-3350.2229	6830.0785	0.0000	1.0000
F81G	59	-3404.1789	6933.4878	0.0000	1.0000
F81IG	60	-3406.4534	6940.2858	0.0000	1.0000
HKY	59	-3410.9264	6946.9827	0.0000	1.0000
GTR	63	-3408.0737	6950.3012	0.0000	1.0000
F81I	59	-3413.0237	6951.1772	0.0000	1.0000
JC69G	56	-3429.1925	6976.7945	0.0000	1.0000
JC69IG	57	-3431.7128	6984.0708	0.0000	1.0000
K2P	56	-3436.2086	6990.8268	0.0000	1.0000
SYM	60	-3433.0057	6993.3904	0.0000	1.0000
JC69I	56	-3437.8481	6994.1058	0.0000	1.0000
F81	58	-3493.8781	7110.6416	0.0000	1.0000
JC69	55	-3517.6210	7151.4205	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3315.4830	6758.9659	0.4959	0.4959
HKYG	60	-3319.6106	6759.2213	0.4364	0.9323
GTRIG	65	-3317.0855	6764.1709	0.0367	0.9690
HKYIG	61	-3321.2686	6764.5371	0.0306	0.9996
HKYI	60	-3327.1788	6774.3577	0.0002	0.9998
GTRI	64	-3323.4882	6774.9763	0.0002	1.0000
K2PG	57	-3345.0621	6804.1242	0.0000	1.0000
SYMG	61	-3342.8479	6807.6959	0.0000	1.0000
K2PIG	58	-3346.6974	6809.3948	0.0000	1.0000
SYMIG	62	-3344.5010	6813.0021	0.0000	1.0000
K2PI	57	-3352.7335	6819.4669	0.0000	1.0000
SYMI	61	-3350.2229	6822.4458	0.0000	1.0000
F81G	59	-3404.1789	6926.3579	0.0000	1.0000
F81IG	60	-3406.4534	6932.9068	0.0000	1.0000
HKY	59	-3410.9264	6939.8528	0.0000	1.0000
GTR	63	-3408.0737	6942.1475	0.0000	1.0000
F81I	59	-3413.0237	6944.0473	0.0000	1.0000
JC69G	56	-3429.1925	6970.3849	0.0000	1.0000
JC69IG	57	-3431.7128	6977.4256	0.0000	1.0000
K2P	56	-3436.2086	6984.4171	0.0000	1.0000
SYM	60	-3433.0057	6986.0113	0.0000	1.0000
JC69I	56	-3437.8481	6987.6962	0.0000	1.0000
F81	58	-3493.8781	7103.7563	0.0000	1.0000
JC69	55	-3517.6210	7145.2419	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3319.6106	7056.7852	0.9936	0.9936
HKYIG	61	-3321.2686	7067.0604	0.0058	0.9994
HKYI	60	-3327.1788	7071.9216	0.0005	0.9999
GTRG	64	-3315.4830	7076.3674	0.0001	1.0000
GTRIG	65	-3317.0855	7086.5318	0.0000	1.0000
K2PG	57	-3345.0621	7086.8099	0.0000	1.0000
GTRI	64	-3323.4882	7092.3778	0.0000	1.0000
K2PIG	58	-3346.6974	7097.0399	0.0000	1.0000
K2PI	57	-3352.7335	7102.1526	0.0000	1.0000
SYMG	61	-3342.8479	7110.2192	0.0000	1.0000
SYMIG	62	-3344.5010	7120.4848	0.0000	1.0000
SYMI	61	-3350.2229	7124.9691	0.0000	1.0000
F81G	59	-3404.1789	7218.9624	0.0000	1.0000
F81IG	60	-3406.4534	7230.4707	0.0000	1.0000
HKY	59	-3410.9264	7232.4573	0.0000	1.0000
F81I	59	-3413.0237	7236.6518	0.0000	1.0000
JC69G	56	-3429.1925	7248.1112	0.0000	1.0000
GTR	63	-3408.0737	7254.5896	0.0000	1.0000
JC69IG	57	-3431.7128	7260.1113	0.0000	1.0000
K2P	56	-3436.2086	7262.1434	0.0000	1.0000
JC69I	56	-3437.8481	7265.4225	0.0000	1.0000
SYM	60	-3433.0057	7283.5752	0.0000	1.0000
F81	58	-3493.8781	7391.4014	0.0000	1.0000
JC69	55	-3517.6210	7418.0088	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/9210_hg18.phylip = 1 - 1053;
 Partition Dummy = 1:261-aggressive-phylip/9210_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 (HKYG)]

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


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

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

