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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1315.9815	2779.4935	0.5866	0.5866
GTRIG	65	-1315.4075	2781.0032	0.2757	0.8623
GTRI	64	-1317.4680	2782.4666	0.1327	0.9950
HKYG	60	-1326.5970	2790.2172	0.0028	0.9977
HKYIG	61	-1326.0866	2791.8048	0.0012	0.9990
HKYI	60	-1327.5927	2792.2086	0.0010	1.0000
GTR	63	-1330.4976	2805.8805	0.0000	1.0000
SYMG	61	-1333.5097	2806.6511	0.0000	1.0000
SYMIG	62	-1332.8029	2807.8582	0.0000	1.0000
K2PG	57	-1339.7282	2808.7267	0.0000	1.0000
SYMI	61	-1334.9200	2809.4717	0.0000	1.0000
K2PIG	58	-1339.0203	2809.8833	0.0000	1.0000
K2PI	57	-1340.8731	2811.0165	0.0000	1.0000
HKY	59	-1341.5521	2817.5310	0.0000	1.0000
SYM	60	-1348.6953	2834.4139	0.0000	1.0000
K2P	56	-1356.2927	2839.2951	0.0000	1.0000
F81G	59	-1354.3847	2843.1963	0.0000	1.0000
F81IG	60	-1353.8931	2844.8094	0.0000	1.0000
F81I	59	-1355.3704	2845.1678	0.0000	1.0000
JC69G	56	-1365.8380	2858.3857	0.0000	1.0000
JC69IG	57	-1365.1128	2859.4957	0.0000	1.0000
JC69I	56	-1367.0739	2860.8574	0.0000	1.0000
F81	58	-1368.9838	2869.8102	0.0000	1.0000
JC69	55	-1381.7062	2887.5733	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1315.9815	2759.9630	0.5314	0.5314
GTRIG	65	-1315.4075	2760.8150	0.3470	0.8784
GTRI	64	-1317.4680	2762.9361	0.1202	0.9986
HKYG	60	-1326.5970	2773.1940	0.0007	0.9993
HKYIG	61	-1326.0866	2774.1731	0.0004	0.9997
HKYI	60	-1327.5927	2775.1854	0.0003	1.0000
GTR	63	-1330.4976	2786.9953	0.0000	1.0000
SYMG	61	-1333.5097	2789.0194	0.0000	1.0000
SYMIG	62	-1332.8029	2789.6059	0.0000	1.0000
SYMI	61	-1334.9200	2791.8400	0.0000	1.0000
K2PG	57	-1339.7282	2793.4565	0.0000	1.0000
K2PIG	58	-1339.0203	2794.0407	0.0000	1.0000
K2PI	57	-1340.8731	2795.7463	0.0000	1.0000
HKY	59	-1341.5521	2801.1041	0.0000	1.0000
SYM	60	-1348.6953	2817.3906	0.0000	1.0000
K2P	56	-1356.2927	2824.5854	0.0000	1.0000
F81G	59	-1354.3847	2826.7694	0.0000	1.0000
F81IG	60	-1353.8931	2827.7861	0.0000	1.0000
F81I	59	-1355.3704	2828.7409	0.0000	1.0000
JC69G	56	-1365.8380	2843.6760	0.0000	1.0000
JC69IG	57	-1365.1128	2844.2255	0.0000	1.0000
JC69I	56	-1367.0739	2846.1477	0.0000	1.0000
F81	58	-1368.9838	2853.9676	0.0000	1.0000
JC69	55	-1381.7062	2873.4123	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1326.5970	3024.9806	0.5896	0.5896
HKYI	60	-1327.5927	3026.9720	0.2178	0.8075
GTRG	64	-1315.9815	3028.5354	0.0997	0.9072
HKYIG	61	-1326.0866	3030.1562	0.0443	0.9515
GTRI	64	-1317.4680	3031.5085	0.0225	0.9740
K2PG	57	-1339.7282	3032.6538	0.0127	0.9868
GTRIG	65	-1315.4075	3033.5838	0.0080	0.9947
K2PI	57	-1340.8731	3034.9436	0.0040	0.9988
K2PIG	58	-1339.0203	3037.4344	0.0012	1.0000
SYMG	61	-1333.5097	3045.0025	0.0000	1.0000
SYMI	61	-1334.9200	3047.8231	0.0000	1.0000
HKY	59	-1341.5521	3048.6943	0.0000	1.0000
SYMIG	62	-1332.8029	3049.7854	0.0000	1.0000
GTR	63	-1330.4976	3051.3713	0.0000	1.0000
K2P	56	-1356.2927	3059.5863	0.0000	1.0000
SYM	60	-1348.6953	3069.1773	0.0000	1.0000
F81G	59	-1354.3847	3074.3596	0.0000	1.0000
F81I	59	-1355.3704	3076.3311	0.0000	1.0000
JC69G	56	-1365.8380	3078.6769	0.0000	1.0000
F81IG	60	-1353.8931	3079.5728	0.0000	1.0000
JC69I	56	-1367.0739	3081.1486	0.0000	1.0000
JC69IG	57	-1365.1128	3083.4229	0.0000	1.0000
F81	58	-1368.9838	3097.3614	0.0000	1.0000
JC69	55	-1381.7062	3104.2168	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/chr9_7170.phylip = 1 - 491;
 Partition Dummy = 1:183-phylip/chr9_7170.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/chr9_7170.phylip = 1 - 491;
 Partition Dummy = 1:183-phylip/chr9_7170.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 (HKYG)]

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

