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

Input data from file "261-aggressive-phylip/11487_hg18.phylip" (29 taxa, 491 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	-3219.5458	6586.6221	0.8628	0.8628
GTRIG	65	-3220.0622	6590.3126	0.1363	0.9992
SYMG	61	-3230.4966	6600.6248	0.0008	0.9999
SYMIG	62	-3231.9640	6606.1804	0.0000	1.0000
K2PG	57	-3240.7274	6610.7249	0.0000	1.0000
HKYG	60	-3236.8882	6610.7998	0.0000	1.0000
HKYIG	61	-3236.0928	6611.8172	0.0000	1.0000
K2PIG	58	-3241.5528	6614.9482	0.0000	1.0000
GTRI	64	-3237.6547	6622.8399	0.0000	1.0000
SYMI	61	-3246.7662	6633.1641	0.0000	1.0000
K2PI	57	-3257.1640	6643.5982	0.0000	1.0000
HKYI	60	-3253.4538	6643.9309	0.0000	1.0000
JC69G	56	-3342.7312	6812.1720	0.0000	1.0000
JC69IG	57	-3342.3891	6814.0483	0.0000	1.0000
F81G	59	-3340.2278	6814.8824	0.0000	1.0000
F81IG	60	-3340.9293	6818.8818	0.0000	1.0000
GTR	63	-3343.8830	6832.6513	0.0000	1.0000
JC69I	56	-3356.9744	6840.6584	0.0000	1.0000
F81I	59	-3354.5456	6843.5182	0.0000	1.0000
SYM	60	-3357.8839	6852.7910	0.0000	1.0000
HKY	59	-3364.2493	6862.9256	0.0000	1.0000
K2P	56	-3370.8521	6868.4140	0.0000	1.0000
JC69	55	-3468.5202	7061.2013	0.0000	1.0000
F81	58	-3465.2817	7062.4061	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3219.5458	6567.0916	0.8198	0.8198
GTRIG	65	-3220.0622	6570.1244	0.1799	0.9997
SYMG	61	-3230.4966	6582.9931	0.0003	1.0000
SYMIG	62	-3231.9640	6587.9280	0.0000	1.0000
HKYG	60	-3236.8882	6593.7765	0.0000	1.0000
HKYIG	61	-3236.0928	6594.1855	0.0000	1.0000
K2PG	57	-3240.7274	6595.4547	0.0000	1.0000
K2PIG	58	-3241.5528	6599.1056	0.0000	1.0000
GTRI	64	-3237.6547	6603.3094	0.0000	1.0000
SYMI	61	-3246.7662	6615.5324	0.0000	1.0000
HKYI	60	-3253.4538	6626.9077	0.0000	1.0000
K2PI	57	-3257.1640	6628.3280	0.0000	1.0000
JC69G	56	-3342.7312	6797.4623	0.0000	1.0000
F81G	59	-3340.2278	6798.4555	0.0000	1.0000
JC69IG	57	-3342.3891	6798.7781	0.0000	1.0000
F81IG	60	-3340.9293	6801.8585	0.0000	1.0000
GTR	63	-3343.8830	6813.7661	0.0000	1.0000
JC69I	56	-3356.9744	6825.9488	0.0000	1.0000
F81I	59	-3354.5456	6827.0912	0.0000	1.0000
SYM	60	-3357.8839	6835.7677	0.0000	1.0000
HKY	59	-3364.2493	6846.4987	0.0000	1.0000
K2P	56	-3370.8521	6853.7043	0.0000	1.0000
F81	58	-3465.2817	7046.5635	0.0000	1.0000
JC69	55	-3468.5202	7047.0404	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-3240.7274	6834.6520	0.5682	0.5682
GTRG	64	-3219.5458	6835.6640	0.3426	0.9108
SYMG	61	-3230.4966	6838.9762	0.0654	0.9762
K2PIG	58	-3241.5528	6842.4993	0.0112	0.9874
GTRIG	65	-3220.0622	6842.8933	0.0092	0.9966
HKYG	60	-3236.8882	6845.5631	0.0024	0.9991
SYMIG	62	-3231.9640	6848.1076	0.0007	0.9998
HKYIG	61	-3236.0928	6850.1686	0.0002	1.0000
K2PI	57	-3257.1640	6867.5253	0.0000	1.0000
SYMI	61	-3246.7662	6871.5155	0.0000	1.0000
GTRI	64	-3237.6547	6871.8818	0.0000	1.0000
HKYI	60	-3253.4538	6878.6943	0.0000	1.0000
JC69G	56	-3342.7312	7032.4632	0.0000	1.0000
JC69IG	57	-3342.3891	7037.9755	0.0000	1.0000
F81G	59	-3340.2278	7046.0457	0.0000	1.0000
F81IG	60	-3340.9293	7053.6452	0.0000	1.0000
JC69I	56	-3356.9744	7060.9496	0.0000	1.0000
F81I	59	-3354.5456	7074.6814	0.0000	1.0000
GTR	63	-3343.8830	7078.1421	0.0000	1.0000
SYM	60	-3357.8839	7087.5543	0.0000	1.0000
K2P	56	-3370.8521	7088.7052	0.0000	1.0000
HKY	59	-3364.2493	7094.0889	0.0000	1.0000
JC69	55	-3468.5202	7277.8448	0.0000	1.0000
F81	58	-3465.2817	7289.9572	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/11487_hg18.phylip = 1 - 491;
 Partition Dummy = 1:261-aggressive-phylip/11487_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/11487_hg18.phylip = 1 - 491;
 Partition Dummy = 1:261-aggressive-phylip/11487_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/11487_hg18.phylip = 1 - 491;
 Partition Dummy = 1:261-aggressive-phylip/11487_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

