
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:35 02/17/2011

Input data from file "183-phylip/chr1_5279.phylip" (29 taxa, 604 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	60	-1670.5215	3474.5237	0.3897	0.3897
HKYG	60	-1670.5375	3474.5556	0.3835	0.7732
HKYIG	61	-1670.5374	3477.0306	0.1113	0.8844
GTRI	64	-1667.4907	3478.4174	0.0556	0.9400
GTRG	64	-1667.6658	3478.7675	0.0467	0.9867
GTRIG	65	-1667.6658	3481.2795	0.0133	1.0000
HKY	59	-1687.0284	3505.0716	0.0000	1.0000
GTR	63	-1683.3971	3507.7275	0.0000	1.0000
F81G	59	-1706.6875	3544.3898	0.0000	1.0000
K2PG	57	-1709.2016	3544.5130	0.0000	1.0000
F81I	59	-1706.9809	3544.9765	0.0000	1.0000
K2PI	57	-1709.6603	3545.4304	0.0000	1.0000
F81IG	60	-1706.6770	3546.8346	0.0000	1.0000
K2PIG	58	-1709.1572	3546.8721	0.0000	1.0000
SYMG	61	-1706.9052	3549.7660	0.0000	1.0000
SYMI	61	-1707.3021	3550.5598	0.0000	1.0000
SYMIG	62	-1706.8833	3552.2065	0.0000	1.0000
F81	58	-1722.1709	3572.8995	0.0000	1.0000
K2P	56	-1724.9629	3573.5967	0.0000	1.0000
SYM	60	-1721.9173	3577.3153	0.0000	1.0000
JC69G	56	-1738.9372	3601.5454	0.0000	1.0000
JC69I	56	-1739.4310	3602.5329	0.0000	1.0000
JC69IG	57	-1738.8839	3603.8777	0.0000	1.0000
JC69	55	-1754.3933	3630.0275	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-1670.5215	3461.0431	0.3163	0.3163
HKYG	60	-1670.5375	3461.0750	0.3113	0.6277
GTRI	64	-1667.4907	3462.9814	0.1200	0.7477
HKYIG	61	-1670.5374	3463.0749	0.1145	0.8622
GTRG	64	-1667.6658	3463.3315	0.1007	0.9629
GTRIG	65	-1667.6658	3465.3315	0.0371	1.0000
HKY	59	-1687.0284	3492.0569	0.0000	1.0000
GTR	63	-1683.3971	3492.7941	0.0000	1.0000
F81G	59	-1706.6875	3531.3751	0.0000	1.0000
F81I	59	-1706.9809	3531.9618	0.0000	1.0000
K2PG	57	-1709.2016	3532.4031	0.0000	1.0000
K2PI	57	-1709.6603	3533.3205	0.0000	1.0000
F81IG	60	-1706.6770	3533.3539	0.0000	1.0000
K2PIG	58	-1709.1572	3534.3143	0.0000	1.0000
SYMG	61	-1706.9052	3535.8103	0.0000	1.0000
SYMI	61	-1707.3021	3536.6041	0.0000	1.0000
SYMIG	62	-1706.8833	3537.7665	0.0000	1.0000
F81	58	-1722.1709	3560.3417	0.0000	1.0000
K2P	56	-1724.9629	3561.9258	0.0000	1.0000
SYM	60	-1721.9173	3563.8346	0.0000	1.0000
JC69G	56	-1738.9372	3589.8744	0.0000	1.0000
JC69I	56	-1739.4310	3590.8619	0.0000	1.0000
JC69IG	57	-1738.8839	3591.7678	0.0000	1.0000
JC69	55	-1754.3933	3618.7867	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1670.5215	3725.2575	0.4940	0.4940
HKYG	60	-1670.5375	3725.2894	0.4862	0.9802
HKYIG	61	-1670.5374	3731.6929	0.0198	0.9999
GTRI	64	-1667.4907	3744.8102	0.0000	1.0000
GTRG	64	-1667.6658	3745.1603	0.0000	1.0000
GTRIG	65	-1667.6658	3751.5638	0.0000	1.0000
HKY	59	-1687.0284	3751.8677	0.0000	1.0000
GTR	63	-1683.3971	3770.2193	0.0000	1.0000
K2PG	57	-1709.2016	3783.4069	0.0000	1.0000
K2PI	57	-1709.6603	3784.3243	0.0000	1.0000
K2PIG	58	-1709.1572	3789.7216	0.0000	1.0000
F81G	59	-1706.6875	3791.1860	0.0000	1.0000
F81I	59	-1706.9809	3791.7726	0.0000	1.0000
F81IG	60	-1706.6770	3797.5684	0.0000	1.0000
SYMG	61	-1706.9052	3804.4283	0.0000	1.0000
SYMI	61	-1707.3021	3805.2221	0.0000	1.0000
K2P	56	-1724.9629	3808.5259	0.0000	1.0000
SYMIG	62	-1706.8833	3810.7881	0.0000	1.0000
F81	58	-1722.1709	3815.7490	0.0000	1.0000
SYM	60	-1721.9173	3828.0491	0.0000	1.0000
JC69G	56	-1738.9372	3836.4746	0.0000	1.0000
JC69I	56	-1739.4310	3837.4621	0.0000	1.0000
JC69IG	57	-1738.8839	3842.7715	0.0000	1.0000
JC69	55	-1754.3933	3860.9832	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_5279.phylip = 1 - 604;
 Partition Dummy = 1:183-phylip/chr1_5279.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr1_5279.phylip = 1 - 604;
 Partition Dummy = 1:183-phylip/chr1_5279.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_5279.phylip = 1 - 604;
 Partition Dummy = 1:183-phylip/chr1_5279.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

