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

Input data from file "183-phylip/chr8_9158.phylip" (29 taxa, 348 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	-1479.6773	3116.7539	0.4894	0.4894
SYMG	61	-1484.8521	3118.1517	0.2433	0.7327
GTRIG	65	-1479.2879	3119.0013	0.1591	0.8918
SYMIG	62	-1484.5689	3120.5483	0.0734	0.9652
K2PG	57	-1493.6260	3124.0520	0.0127	0.9779
HKYG	60	-1489.5806	3124.6664	0.0094	0.9873
SYMI	61	-1488.9271	3126.3017	0.0041	0.9914
GTRI	64	-1484.4594	3126.3181	0.0041	0.9955
HKYIG	61	-1489.2653	3126.9782	0.0029	0.9985
K2PIG	58	-1494.4035	3128.4886	0.0014	0.9998
K2PI	57	-1498.4556	3133.7113	0.0001	0.9999
HKYI	60	-1494.6073	3134.7198	0.0001	1.0000
SYM	60	-1515.6203	3176.7459	0.0000	1.0000
GTR	63	-1511.5586	3177.5115	0.0000	1.0000
K2P	56	-1523.5922	3181.1225	0.0000	1.0000
HKY	59	-1520.4894	3183.5622	0.0000	1.0000
JC69G	56	-1526.1790	3186.2962	0.0000	1.0000
F81G	59	-1522.7387	3188.0607	0.0000	1.0000
JC69IG	57	-1525.8702	3188.5404	0.0000	1.0000
F81IG	60	-1522.4202	3190.3456	0.0000	1.0000
JC69I	56	-1530.0795	3194.0971	0.0000	1.0000
F81I	59	-1526.7275	3196.0384	0.0000	1.0000
JC69	55	-1555.1404	3241.3768	0.0000	1.0000
F81	58	-1551.7284	3243.1384	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1479.6773	3087.3546	0.5788	0.5788
GTRIG	65	-1479.2879	3088.5757	0.3143	0.8931
SYMG	61	-1484.8521	3091.7042	0.0658	0.9589
SYMIG	62	-1484.5689	3093.1378	0.0321	0.9910
GTRI	64	-1484.4594	3096.9188	0.0048	0.9958
HKYG	60	-1489.5806	3099.1612	0.0016	0.9974
SYMI	61	-1488.9271	3099.8541	0.0011	0.9985
HKYIG	61	-1489.2653	3100.5307	0.0008	0.9993
K2PG	57	-1493.6260	3101.2520	0.0006	0.9999
K2PIG	58	-1494.4035	3104.8069	0.0001	1.0000
HKYI	60	-1494.6073	3109.2146	0.0000	1.0000
K2PI	57	-1498.4556	3110.9113	0.0000	1.0000
GTR	63	-1511.5586	3149.1171	0.0000	1.0000
SYM	60	-1515.6203	3151.2406	0.0000	1.0000
HKY	59	-1520.4894	3158.9789	0.0000	1.0000
K2P	56	-1523.5922	3159.1844	0.0000	1.0000
F81G	59	-1522.7387	3163.4774	0.0000	1.0000
JC69G	56	-1526.1790	3164.3580	0.0000	1.0000
F81IG	60	-1522.4202	3164.8404	0.0000	1.0000
JC69IG	57	-1525.8702	3165.7404	0.0000	1.0000
F81I	59	-1526.7275	3171.4550	0.0000	1.0000
JC69I	56	-1530.0795	3172.1589	0.0000	1.0000
F81	58	-1551.7284	3219.4567	0.0000	1.0000
JC69	55	-1555.1404	3220.2809	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1493.6260	3320.8275	0.9076	0.9076
SYMG	61	-1484.8521	3326.6885	0.0484	0.9561
K2PIG	58	-1494.4035	3328.2347	0.0224	0.9784
HKYG	60	-1489.5806	3330.2933	0.0080	0.9864
K2PI	57	-1498.4556	3330.4868	0.0073	0.9937
SYMIG	62	-1484.5689	3331.9743	0.0034	0.9971
GTRG	64	-1479.6773	3333.8955	0.0013	0.9984
SYMI	61	-1488.9271	3334.8385	0.0008	0.9992
HKYIG	61	-1489.2653	3335.5150	0.0006	0.9998
GTRIG	65	-1479.2879	3338.9689	0.0001	0.9999
HKYI	60	-1494.6073	3340.3468	0.0001	1.0000
GTRI	64	-1484.4594	3343.4597	0.0000	1.0000
K2P	56	-1523.5922	3374.9077	0.0000	1.0000
JC69G	56	-1526.1790	3380.0814	0.0000	1.0000
SYM	60	-1515.6203	3382.3728	0.0000	1.0000
JC69IG	57	-1525.8702	3385.3159	0.0000	1.0000
HKY	59	-1520.4894	3386.2588	0.0000	1.0000
JC69I	56	-1530.0795	3387.8823	0.0000	1.0000
F81G	59	-1522.7387	3390.7573	0.0000	1.0000
GTR	63	-1511.5586	3391.8059	0.0000	1.0000
F81IG	60	-1522.4202	3395.9725	0.0000	1.0000
F81I	59	-1526.7275	3398.7350	0.0000	1.0000
JC69	55	-1555.1404	3432.1520	0.0000	1.0000
F81	58	-1551.7284	3442.8845	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 183-phylip/chr8_9158.phylip = 1 - 348;
 Partition Dummy = 1:183-phylip/chr8_9158.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 183-phylip/chr8_9158.phylip = 1 - 348;
 Partition Dummy = 1:183-phylip/chr8_9158.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 183-phylip/chr8_9158.phylip = 1 - 348;
 Partition Dummy = 1:183-phylip/chr8_9158.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

