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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	61	-1493.4019	3130.4771	0.5031	0.5031
HKYG	60	-1495.0204	3130.9551	0.3961	0.8992
GTRIG	65	-1490.0070	3134.8837	0.0556	0.9547
GTRG	64	-1491.6962	3135.4387	0.0421	0.9968
HKYI	60	-1499.9341	3140.7826	0.0029	0.9997
GTRI	64	-1497.1742	3146.3946	0.0002	0.9999
K2PIG	58	-1506.6798	3148.8028	0.0001	1.0000
K2PG	57	-1508.3634	3149.4577	0.0000	1.0000
SYMIG	62	-1505.4483	3157.3450	0.0000	1.0000
K2PI	57	-1512.5128	3157.7566	0.0000	1.0000
SYMG	61	-1507.1291	3157.9316	0.0000	1.0000
SYMI	61	-1511.3247	3166.3228	0.0000	1.0000
HKY	59	-1532.4496	3203.0702	0.0000	1.0000
F81G	59	-1532.7720	3203.7150	0.0000	1.0000
F81IG	60	-1531.4104	3203.7350	0.0000	1.0000
GTR	63	-1529.0526	3207.3444	0.0000	1.0000
F81I	59	-1536.8049	3211.7807	0.0000	1.0000
JC69IG	57	-1543.0921	3218.9152	0.0000	1.0000
JC69G	56	-1544.5038	3219.0415	0.0000	1.0000
K2P	56	-1547.6631	3225.3601	0.0000	1.0000
JC69I	56	-1548.2702	3226.5743	0.0000	1.0000
SYM	60	-1546.0874	3233.0891	0.0000	1.0000
F81	58	-1568.9686	3273.3805	0.0000	1.0000
JC69	55	-1582.3334	3292.0189	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1493.4019	3108.8037	0.4230	0.4230
GTRIG	65	-1490.0070	3110.0141	0.2310	0.6540
HKYG	60	-1495.0204	3110.0408	0.2279	0.8819
GTRG	64	-1491.6962	3111.3925	0.1159	0.9978
HKYI	60	-1499.9341	3119.8683	0.0017	0.9995
GTRI	64	-1497.1742	3122.3484	0.0005	1.0000
K2PIG	58	-1506.6798	3129.3597	0.0000	1.0000
K2PG	57	-1508.3634	3130.7268	0.0000	1.0000
SYMIG	62	-1505.4483	3134.8967	0.0000	1.0000
SYMG	61	-1507.1291	3136.2582	0.0000	1.0000
K2PI	57	-1512.5128	3139.0257	0.0000	1.0000
SYMI	61	-1511.3247	3144.6494	0.0000	1.0000
F81IG	60	-1531.4104	3182.8207	0.0000	1.0000
HKY	59	-1532.4496	3182.8993	0.0000	1.0000
F81G	59	-1532.7720	3183.5441	0.0000	1.0000
GTR	63	-1529.0526	3184.1052	0.0000	1.0000
F81I	59	-1536.8049	3191.6097	0.0000	1.0000
JC69IG	57	-1543.0921	3200.1843	0.0000	1.0000
JC69G	56	-1544.5038	3201.0076	0.0000	1.0000
K2P	56	-1547.6631	3207.3262	0.0000	1.0000
JC69I	56	-1548.2702	3208.5404	0.0000	1.0000
SYM	60	-1546.0874	3212.1749	0.0000	1.0000
F81	58	-1568.9686	3253.9373	0.0000	1.0000
JC69	55	-1582.3334	3274.6668	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1495.0204	3351.1564	0.7852	0.7852
HKYIG	61	-1493.4019	3353.9379	0.1954	0.9806
K2PG	57	-1508.3634	3359.7866	0.0105	0.9911
HKYI	60	-1499.9341	3360.9839	0.0058	0.9969
K2PIG	58	-1506.6798	3362.4381	0.0028	0.9997
K2PI	57	-1512.5128	3368.0855	0.0002	0.9998
GTRG	64	-1491.6962	3368.5824	0.0001	1.0000
GTRIG	65	-1490.0070	3371.2227	0.0000	1.0000
GTRI	64	-1497.1742	3379.5383	0.0000	1.0000
SYMG	61	-1507.1291	3381.3924	0.0000	1.0000
SYMIG	62	-1505.4483	3384.0495	0.0000	1.0000
SYMI	61	-1511.3247	3389.7836	0.0000	1.0000
HKY	59	-1532.4496	3419.9963	0.0000	1.0000
F81G	59	-1532.7720	3420.6411	0.0000	1.0000
F81IG	60	-1531.4104	3423.9363	0.0000	1.0000
JC69G	56	-1544.5038	3426.0489	0.0000	1.0000
F81I	59	-1536.8049	3428.7067	0.0000	1.0000
JC69IG	57	-1543.0921	3429.2441	0.0000	1.0000
K2P	56	-1547.6631	3432.3674	0.0000	1.0000
JC69I	56	-1548.2702	3433.5816	0.0000	1.0000
GTR	63	-1529.0526	3437.2766	0.0000	1.0000
SYM	60	-1546.0874	3453.2905	0.0000	1.0000
F81	58	-1568.9686	3487.0157	0.0000	1.0000
JC69	55	-1582.3334	3495.6894	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_6036.phylip = 1 - 411;
 Partition Dummy = 1:183-phylip/chrZ_6036.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_6036.phylip = 1 - 411;
 Partition Dummy = 1:183-phylip/chrZ_6036.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

