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

Input data from file "183-phylip/chr7_4747.phylip" (29 taxa, 421 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	-1375.8924	2892.1182	0.5004	0.5004
HKYIG	61	-1375.3606	2893.7909	0.2168	0.7172
GTRG	64	-1371.7200	2894.8107	0.1302	0.8474
HKYI	60	-1377.6414	2895.6162	0.0870	0.9345
GTRIG	65	-1371.2916	2896.7522	0.0493	0.9838
GTRI	64	-1373.8049	2898.9806	0.0162	1.0000
K2PG	57	-1396.8825	2925.9798	0.0000	1.0000
HKY	59	-1394.4109	2926.4340	0.0000	1.0000
K2PIG	58	-1396.2799	2927.4658	0.0000	1.0000
GTR	63	-1390.1779	2928.9440	0.0000	1.0000
K2PI	57	-1398.7045	2929.6238	0.0000	1.0000
SYMG	61	-1396.0778	2935.2253	0.0000	1.0000
SYMIG	62	-1395.5156	2936.8524	0.0000	1.0000
SYMI	61	-1397.8837	2938.8370	0.0000	1.0000
F81G	59	-1401.4572	2940.5266	0.0000	1.0000
F81IG	60	-1400.9354	2942.2042	0.0000	1.0000
F81I	59	-1403.1384	2943.8890	0.0000	1.0000
K2P	56	-1415.9948	2961.5281	0.0000	1.0000
JC69G	56	-1419.4897	2968.5179	0.0000	1.0000
JC69IG	57	-1418.9169	2970.0486	0.0000	1.0000
SYM	60	-1414.8660	2970.0652	0.0000	1.0000
JC69I	56	-1421.2683	2972.0750	0.0000	1.0000
F81	58	-1419.6603	2974.2267	0.0000	1.0000
JC69	55	-1438.3114	3003.4995	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1371.7200	2871.4400	0.3121	0.3121
HKYG	60	-1375.8924	2871.7848	0.2627	0.5748
GTRIG	65	-1371.2916	2872.5832	0.1762	0.7510
HKYIG	61	-1375.3606	2872.7213	0.1645	0.9155
HKYI	60	-1377.6414	2875.2828	0.0457	0.9612
GTRI	64	-1373.8049	2875.6098	0.0388	1.0000
GTR	63	-1390.1779	2906.3558	0.0000	1.0000
HKY	59	-1394.4109	2906.8219	0.0000	1.0000
K2PG	57	-1396.8825	2907.7649	0.0000	1.0000
K2PIG	58	-1396.2799	2908.5597	0.0000	1.0000
K2PI	57	-1398.7045	2911.4089	0.0000	1.0000
SYMG	61	-1396.0778	2914.1556	0.0000	1.0000
SYMIG	62	-1395.5156	2915.0312	0.0000	1.0000
SYMI	61	-1397.8837	2917.7674	0.0000	1.0000
F81G	59	-1401.4572	2920.9144	0.0000	1.0000
F81IG	60	-1400.9354	2921.8708	0.0000	1.0000
F81I	59	-1403.1384	2924.2768	0.0000	1.0000
K2P	56	-1415.9948	2943.9896	0.0000	1.0000
SYM	60	-1414.8660	2949.7319	0.0000	1.0000
JC69G	56	-1419.4897	2950.9794	0.0000	1.0000
JC69IG	57	-1418.9169	2951.8337	0.0000	1.0000
JC69I	56	-1421.2683	2954.5365	0.0000	1.0000
F81	58	-1419.6603	2955.3207	0.0000	1.0000
JC69	55	-1438.3114	2986.6228	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1375.8924	3114.3428	0.7953	0.7953
HKYI	60	-1377.6414	3117.8408	0.1383	0.9337
HKYIG	61	-1375.3606	3119.3219	0.0660	0.9996
GTRG	64	-1371.7200	3130.1685	0.0003	0.9999
GTRI	64	-1373.8049	3134.3383	0.0000	1.0000
GTRIG	65	-1371.2916	3135.3544	0.0000	1.0000
K2PG	57	-1396.8825	3138.1950	0.0000	1.0000
K2PI	57	-1398.7045	3141.8390	0.0000	1.0000
K2PIG	58	-1396.2799	3143.0324	0.0000	1.0000
HKY	59	-1394.4109	3145.3372	0.0000	1.0000
F81G	59	-1401.4572	3159.4298	0.0000	1.0000
SYMG	61	-1396.0778	3160.7562	0.0000	1.0000
GTR	63	-1390.1779	3161.0417	0.0000	1.0000
F81I	59	-1403.1384	3162.7921	0.0000	1.0000
SYMI	61	-1397.8837	3164.3680	0.0000	1.0000
F81IG	60	-1400.9354	3164.4288	0.0000	1.0000
SYMIG	62	-1395.5156	3165.6744	0.0000	1.0000
K2P	56	-1415.9948	3170.3771	0.0000	1.0000
JC69G	56	-1419.4897	3177.3669	0.0000	1.0000
JC69I	56	-1421.2683	3180.9240	0.0000	1.0000
JC69IG	57	-1418.9169	3182.2638	0.0000	1.0000
F81	58	-1419.6603	3189.7934	0.0000	1.0000
SYM	60	-1414.8660	3192.2899	0.0000	1.0000
JC69	55	-1438.3114	3208.9676	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/chr7_4747.phylip = 1 - 421;
 Partition Dummy = 1:183-phylip/chr7_4747.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/chr7_4747.phylip = 1 - 421;
 Partition Dummy = 1:183-phylip/chr7_4747.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/chr7_4747.phylip = 1 - 421;
 Partition Dummy = 1:183-phylip/chr7_4747.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

