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

Input data from file "261-aggressive-phylip/10615_hg18.phylip" (29 taxa, 349 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	-1587.4772	3332.2502	0.6185	0.6185
HKYG	60	-1593.9409	3333.2985	0.3662	0.9847
GTRIG	65	-1590.2220	3340.7619	0.0088	0.9934
HKYIG	61	-1596.7822	3341.9198	0.0049	0.9984
GTRI	64	-1593.5015	3344.2989	0.0015	0.9998
HKYI	60	-1601.7753	3348.9672	0.0001	1.0000
K2PG	57	-1609.4645	3355.6507	0.0000	1.0000
K2PIG	58	-1612.2496	3364.0993	0.0000	1.0000
SYMG	61	-1608.1720	3364.6993	0.0000	1.0000
K2PI	57	-1617.8362	3372.3940	0.0000	1.0000
SYMIG	62	-1610.9785	3373.2718	0.0000	1.0000
SYMI	61	-1616.5071	3381.3696	0.0000	1.0000
GTR	63	-1622.2143	3398.7234	0.0000	1.0000
HKY	59	-1632.4698	3407.4379	0.0000	1.0000
F81G	59	-1645.9890	3434.4762	0.0000	1.0000
K2P	56	-1651.1914	3436.2458	0.0000	1.0000
F81IG	60	-1645.9346	3437.2858	0.0000	1.0000
F81I	59	-1649.0779	3440.6541	0.0000	1.0000
SYM	60	-1648.3208	3442.0583	0.0000	1.0000
JC69G	56	-1654.5623	3442.9877	0.0000	1.0000
JC69IG	57	-1657.1631	3451.0478	0.0000	1.0000
JC69I	56	-1660.6216	3455.1061	0.0000	1.0000
F81	58	-1680.2045	3500.0089	0.0000	1.0000
JC69	55	-1693.0940	3517.2119	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1587.4772	3302.9544	0.8984	0.8984
HKYG	60	-1593.9409	3307.8818	0.0765	0.9749
GTRIG	65	-1590.2220	3310.4439	0.0212	0.9962
GTRI	64	-1593.5015	3315.0031	0.0022	0.9983
HKYIG	61	-1596.7822	3315.5644	0.0016	1.0000
HKYI	60	-1601.7753	3323.5505	0.0000	1.0000
K2PG	57	-1609.4645	3332.9290	0.0000	1.0000
SYMG	61	-1608.1720	3338.3439	0.0000	1.0000
K2PIG	58	-1612.2496	3340.4993	0.0000	1.0000
SYMIG	62	-1610.9785	3345.9571	0.0000	1.0000
K2PI	57	-1617.8362	3349.6724	0.0000	1.0000
SYMI	61	-1616.5071	3355.0142	0.0000	1.0000
GTR	63	-1622.2143	3370.4286	0.0000	1.0000
HKY	59	-1632.4698	3382.9396	0.0000	1.0000
F81G	59	-1645.9890	3409.9780	0.0000	1.0000
F81IG	60	-1645.9346	3411.8692	0.0000	1.0000
K2P	56	-1651.1914	3414.3827	0.0000	1.0000
F81I	59	-1649.0779	3416.1559	0.0000	1.0000
SYM	60	-1648.3208	3416.6417	0.0000	1.0000
JC69G	56	-1654.5623	3421.1247	0.0000	1.0000
JC69IG	57	-1657.1631	3428.3262	0.0000	1.0000
JC69I	56	-1660.6216	3433.2431	0.0000	1.0000
F81	58	-1680.2045	3476.4089	0.0000	1.0000
JC69	55	-1693.0940	3496.1880	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1593.9409	3539.1862	0.9901	0.9901
GTRG	64	-1587.4772	3549.6790	0.0052	0.9953
HKYIG	61	-1596.7822	3550.7238	0.0031	0.9984
K2PG	57	-1609.4645	3552.6681	0.0012	0.9996
HKYI	60	-1601.7753	3554.8548	0.0004	1.0000
GTRIG	65	-1590.2220	3561.0236	0.0000	1.0000
GTRI	64	-1593.5015	3561.7277	0.0000	1.0000
K2PIG	58	-1612.2496	3564.0934	0.0000	1.0000
K2PI	57	-1617.8362	3569.4115	0.0000	1.0000
SYMG	61	-1608.1720	3573.5033	0.0000	1.0000
SYMIG	62	-1610.9785	3584.9715	0.0000	1.0000
SYMI	61	-1616.5071	3590.1736	0.0000	1.0000
HKY	59	-1632.4698	3610.3889	0.0000	1.0000
GTR	63	-1622.2143	3613.2982	0.0000	1.0000
K2P	56	-1651.1914	3630.2668	0.0000	1.0000
JC69G	56	-1654.5623	3637.0087	0.0000	1.0000
F81G	59	-1645.9890	3637.4272	0.0000	1.0000
F81IG	60	-1645.9346	3643.1735	0.0000	1.0000
F81I	59	-1649.0779	3643.6051	0.0000	1.0000
SYM	60	-1648.3208	3647.9460	0.0000	1.0000
JC69IG	57	-1657.1631	3648.0653	0.0000	1.0000
JC69I	56	-1660.6216	3649.1271	0.0000	1.0000
F81	58	-1680.2045	3700.0031	0.0000	1.0000
JC69	55	-1693.0940	3708.2169	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/10615_hg18.phylip = 1 - 349;
 Partition Dummy = 1:261-aggressive-phylip/10615_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/10615_hg18.phylip = 1 - 349;
 Partition Dummy = 1:261-aggressive-phylip/10615_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/10615_hg18.phylip = 1 - 349;
 Partition Dummy = 1:261-aggressive-phylip/10615_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

