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

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

Minimum AIC  model: HKYG
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	-1554.1689	3242.5790	0.7588	0.7588
HKYIG	61	-1554.1765	3245.0977	0.2154	0.9742
GTRG	64	-1552.7584	3249.8305	0.0202	0.9944
GTRIG	65	-1552.7660	3252.3886	0.0056	1.0000
HKY	59	-1572.5785	3276.9047	0.0000	1.0000
HKYI	60	-1572.5830	3279.4073	0.0000	1.0000
GTR	63	-1571.3502	3284.4812	0.0000	1.0000
GTRI	64	-1571.3547	3287.0231	0.0000	1.0000
F81G	59	-1587.0988	3305.9451	0.0000	1.0000
F81IG	60	-1587.1064	3308.4540	0.0000	1.0000
K2PG	57	-1599.6873	3326.1638	0.0000	1.0000
K2PIG	58	-1599.6950	3328.6536	0.0000	1.0000
SYMG	61	-1598.5104	3333.7654	0.0000	1.0000
F81	58	-1603.5006	3336.2647	0.0000	1.0000
SYMIG	62	-1598.5181	3336.2940	0.0000	1.0000
F81I	59	-1603.5054	3338.7584	0.0000	1.0000
K2P	56	-1617.5472	3359.4187	0.0000	1.0000
K2PI	57	-1617.5518	3361.8927	0.0000	1.0000
SYM	60	-1616.3010	3366.8432	0.0000	1.0000
SYMI	61	-1616.3056	3369.3558	0.0000	1.0000
JC69G	56	-1625.0634	3374.4511	0.0000	1.0000
JC69IG	57	-1625.0711	3376.9314	0.0000	1.0000
JC69	55	-1642.0640	3405.9969	0.0000	1.0000
JC69I	56	-1642.0687	3408.4617	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1554.1689	3228.3378	0.6814	0.6814
HKYIG	61	-1554.1765	3230.3530	0.2488	0.9302
GTRG	64	-1552.7584	3233.5168	0.0511	0.9813
GTRIG	65	-1552.7660	3235.5320	0.0187	1.0000
HKY	59	-1572.5785	3263.1571	0.0000	1.0000
HKYI	60	-1572.5830	3265.1661	0.0000	1.0000
GTR	63	-1571.3502	3268.7004	0.0000	1.0000
GTRI	64	-1571.3547	3270.7094	0.0000	1.0000
F81G	59	-1587.0988	3292.1975	0.0000	1.0000
F81IG	60	-1587.1064	3294.2128	0.0000	1.0000
K2PG	57	-1599.6873	3313.3746	0.0000	1.0000
K2PIG	58	-1599.6950	3315.3901	0.0000	1.0000
SYMG	61	-1598.5104	3319.0207	0.0000	1.0000
SYMIG	62	-1598.5181	3321.0362	0.0000	1.0000
F81	58	-1603.5006	3323.0011	0.0000	1.0000
F81I	59	-1603.5054	3325.0108	0.0000	1.0000
K2P	56	-1617.5472	3347.0944	0.0000	1.0000
K2PI	57	-1617.5518	3349.1035	0.0000	1.0000
SYM	60	-1616.3010	3352.6020	0.0000	1.0000
SYMI	61	-1616.3056	3354.6111	0.0000	1.0000
JC69G	56	-1625.0634	3362.1268	0.0000	1.0000
JC69IG	57	-1625.0711	3364.1423	0.0000	1.0000
JC69	55	-1642.0640	3394.1279	0.0000	1.0000
JC69I	56	-1642.0687	3396.1374	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1554.1689	3489.6000	0.9602	0.9602
HKYIG	61	-1554.1765	3495.9696	0.0397	1.0000
GTRG	64	-1552.7584	3512.1965	0.0000	1.0000
GTRIG	65	-1552.7660	3518.5661	0.0000	1.0000
HKY	59	-1572.5785	3520.0649	0.0000	1.0000
HKYI	60	-1572.5830	3526.4283	0.0000	1.0000
GTR	63	-1571.3502	3543.0257	0.0000	1.0000
F81G	59	-1587.0988	3549.1053	0.0000	1.0000
GTRI	64	-1571.3547	3549.3891	0.0000	1.0000
F81IG	60	-1587.1064	3555.4750	0.0000	1.0000
K2PG	57	-1599.6873	3561.5737	0.0000	1.0000
K2PIG	58	-1599.6950	3567.9435	0.0000	1.0000
F81	58	-1603.5006	3575.5546	0.0000	1.0000
F81I	59	-1603.5054	3581.9186	0.0000	1.0000
SYMG	61	-1598.5104	3584.6373	0.0000	1.0000
K2P	56	-1617.5472	3590.9391	0.0000	1.0000
SYMIG	62	-1598.5181	3591.0071	0.0000	1.0000
K2PI	57	-1617.5518	3597.3026	0.0000	1.0000
JC69G	56	-1625.0634	3605.9715	0.0000	1.0000
JC69IG	57	-1625.0711	3612.3414	0.0000	1.0000
SYM	60	-1616.3010	3613.8642	0.0000	1.0000
SYMI	61	-1616.3056	3620.2277	0.0000	1.0000
JC69	55	-1642.0640	3633.6183	0.0000	1.0000
JC69I	56	-1642.0687	3639.9821	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/3835_hg18.phylip = 1 - 575;
 Partition Dummy = 1:261-aggressive-phylip/3835_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 AICc model (HKYG)]

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

