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

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

Minimum AIC  model: GTRG
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	-1517.6712	3179.5008	0.4342	0.4342
GTRG	64	-1512.4576	3180.7412	0.2335	0.6677
HKYIG	61	-1517.1310	3181.3083	0.1759	0.8435
GTRIG	65	-1511.6042	3182.0004	0.1244	0.9679
HKYI	60	-1520.7548	3185.6680	0.0199	0.9878
GTRI	64	-1515.4938	3186.8137	0.0112	0.9990
K2PG	57	-1528.4404	3192.4887	0.0007	0.9997
K2PIG	58	-1527.8721	3194.1836	0.0003	1.0000
K2PI	57	-1531.7120	3199.0318	0.0000	1.0000
SYMG	61	-1526.8336	3200.7136	0.0000	1.0000
SYMIG	62	-1526.2939	3202.5413	0.0000	1.0000
SYMI	61	-1529.9867	3207.0198	0.0000	1.0000
F81G	59	-1553.6559	3248.6012	0.0000	1.0000
F81IG	60	-1553.1406	3250.4397	0.0000	1.0000
HKY	59	-1555.2344	3251.7583	0.0000	1.0000
F81I	59	-1556.3319	3253.9533	0.0000	1.0000
JC69G	56	-1560.9396	3254.6739	0.0000	1.0000
GTR	63	-1551.6923	3256.2646	0.0000	1.0000
JC69IG	57	-1560.4109	3256.4297	0.0000	1.0000
JC69I	56	-1563.7238	3260.2423	0.0000	1.0000
K2P	56	-1564.6835	3262.1618	0.0000	1.0000
SYM	60	-1563.1264	3270.4113	0.0000	1.0000
F81	58	-1588.9396	3316.3184	0.0000	1.0000
JC69	55	-1595.6645	3321.3290	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1512.4576	3152.9151	0.4149	0.4149
GTRIG	65	-1511.6042	3153.2084	0.3583	0.7732
HKYG	60	-1517.6712	3155.3423	0.1233	0.8965
HKYIG	61	-1517.1310	3156.2619	0.0778	0.9743
GTRI	64	-1515.4938	3158.9876	0.0199	0.9943
HKYI	60	-1520.7548	3161.5096	0.0056	0.9999
K2PG	57	-1528.4404	3170.8808	0.0001	1.0000
K2PIG	58	-1527.8721	3171.7443	0.0000	1.0000
SYMG	61	-1526.8336	3175.6672	0.0000	1.0000
SYMIG	62	-1526.2939	3176.5878	0.0000	1.0000
K2PI	57	-1531.7120	3177.4239	0.0000	1.0000
SYMI	61	-1529.9867	3181.9734	0.0000	1.0000
F81G	59	-1553.6559	3225.3117	0.0000	1.0000
F81IG	60	-1553.1406	3226.2813	0.0000	1.0000
HKY	59	-1555.2344	3228.4688	0.0000	1.0000
GTR	63	-1551.6923	3229.3846	0.0000	1.0000
F81I	59	-1556.3319	3230.6639	0.0000	1.0000
JC69G	56	-1560.9396	3233.8791	0.0000	1.0000
JC69IG	57	-1560.4109	3234.8219	0.0000	1.0000
JC69I	56	-1563.7238	3239.4476	0.0000	1.0000
K2P	56	-1564.6835	3241.3670	0.0000	1.0000
SYM	60	-1563.1264	3246.2529	0.0000	1.0000
F81	58	-1588.9396	3293.8791	0.0000	1.0000
JC69	55	-1595.6645	3301.3290	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1517.6712	3389.1716	0.7677	0.7677
K2PG	57	-1528.4404	3393.0186	0.1122	0.8799
HKYIG	61	-1517.1310	3393.9883	0.0691	0.9490
HKYI	60	-1520.7548	3395.3388	0.0352	0.9841
K2PIG	58	-1527.8721	3397.7792	0.0104	0.9945
K2PI	57	-1531.7120	3399.5617	0.0043	0.9987
GTRG	64	-1512.4576	3402.3330	0.0011	0.9998
GTRIG	65	-1511.6042	3406.5234	0.0001	0.9999
GTRI	64	-1515.4938	3408.4054	0.0001	1.0000
SYMG	61	-1526.8336	3413.3936	0.0000	1.0000
SYMIG	62	-1526.2939	3418.2114	0.0000	1.0000
SYMI	61	-1529.9867	3419.6998	0.0000	1.0000
JC69G	56	-1560.9396	3452.1197	0.0000	1.0000
F81G	59	-1553.6559	3455.2438	0.0000	1.0000
JC69IG	57	-1560.4109	3456.9596	0.0000	1.0000
JC69I	56	-1563.7238	3457.6882	0.0000	1.0000
HKY	59	-1555.2344	3458.4009	0.0000	1.0000
K2P	56	-1564.6835	3459.6076	0.0000	1.0000
F81IG	60	-1553.1406	3460.1105	0.0000	1.0000
F81I	59	-1556.3319	3460.5959	0.0000	1.0000
GTR	63	-1551.6923	3474.9053	0.0000	1.0000
SYM	60	-1563.1264	3480.0821	0.0000	1.0000
JC69	55	-1595.6645	3515.6724	0.0000	1.0000
F81	58	-1588.9396	3519.9140	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/chr11_4777.phylip = 1 - 364;
 Partition Dummy = 1:183-phylip/chr11_4777.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr11_4777.phylip = 1 - 364;
 Partition Dummy = 1:183-phylip/chr11_4777.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 183-phylip/chr11_4777.phylip = 1 - 364;
 Partition Dummy = 1:183-phylip/chr11_4777.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

