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

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

Minimum AIC  model: GTRIG
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	-1377.9667	2908.8069	0.4652	0.4652
HKYG	60	-1379.7692	2909.2944	0.3646	0.8297
GTRIG	65	-1372.9924	2911.5864	0.1159	0.9456
SYMIG	62	-1379.1485	2914.3134	0.0296	0.9753
GTRG	64	-1376.5289	2915.4379	0.0169	0.9922
SYMG	61	-1382.0592	2916.9918	0.0078	0.9999
K2PIG	58	-1392.4256	2928.4480	0.0000	1.0000
HKYI	60	-1389.4449	2928.6459	0.0000	1.0000
GTRI	64	-1385.5599	2933.4999	0.0000	1.0000
SYMI	61	-1391.0339	2934.9413	0.0000	1.0000
K2PG	57	-1397.7905	2936.1352	0.0000	1.0000
K2PI	57	-1401.7596	2944.0734	0.0000	1.0000
F81IG	60	-1398.1111	2945.9783	0.0000	1.0000
F81G	59	-1403.0799	2952.8237	0.0000	1.0000
F81I	59	-1408.5965	2963.8570	0.0000	1.0000
JC69IG	57	-1420.5377	2981.6296	0.0000	1.0000
JC69G	56	-1425.2591	2988.0542	0.0000	1.0000
JC69I	56	-1434.4741	3006.4842	0.0000	1.0000
SYM	60	-1440.1759	3030.1079	0.0000	1.0000
GTR	63	-1436.2460	3031.6771	0.0000	1.0000
HKY	59	-1446.6940	3040.0520	0.0000	1.0000
K2P	56	-1460.9450	3059.4261	0.0000	1.0000
F81	58	-1470.7791	3085.1549	0.0000	1.0000
JC69	55	-1487.8262	3110.1941	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1372.9924	2875.9847	0.5971	0.5971
HKYIG	61	-1377.9667	2877.9335	0.2254	0.8225
HKYG	60	-1379.7692	2879.5383	0.1010	0.9235
GTRG	64	-1376.5289	2881.0577	0.0473	0.9708
SYMIG	62	-1379.1485	2882.2970	0.0254	0.9962
SYMG	61	-1382.0592	2886.1184	0.0038	1.0000
HKYI	60	-1389.4449	2898.8898	0.0000	1.0000
GTRI	64	-1385.5599	2899.1197	0.0000	1.0000
K2PIG	58	-1392.4256	2900.8512	0.0000	1.0000
SYMI	61	-1391.0339	2904.0679	0.0000	1.0000
K2PG	57	-1397.7905	2909.5810	0.0000	1.0000
F81IG	60	-1398.1111	2916.2222	0.0000	1.0000
K2PI	57	-1401.7596	2917.5191	0.0000	1.0000
F81G	59	-1403.0799	2924.1598	0.0000	1.0000
F81I	59	-1408.5965	2935.1930	0.0000	1.0000
JC69IG	57	-1420.5377	2955.0754	0.0000	1.0000
JC69G	56	-1425.2591	2962.5182	0.0000	1.0000
JC69I	56	-1434.4741	2980.9482	0.0000	1.0000
GTR	63	-1436.2460	2998.4919	0.0000	1.0000
SYM	60	-1440.1759	3000.3518	0.0000	1.0000
HKY	59	-1446.6940	3011.3880	0.0000	1.0000
K2P	56	-1460.9450	3033.8901	0.0000	1.0000
F81	58	-1470.7791	3057.5581	0.0000	1.0000
JC69	55	-1487.8262	3085.6523	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1379.7692	3103.1492	0.7354	0.7354
HKYIG	61	-1377.9667	3105.2712	0.2545	0.9899
SYMIG	62	-1379.1485	3113.3616	0.0045	0.9943
SYMG	61	-1382.0592	3113.4561	0.0043	0.9986
K2PIG	58	-1392.4256	3117.0084	0.0007	0.9993
GTRIG	65	-1372.9924	3118.2298	0.0004	0.9997
GTRG	64	-1376.5289	3119.5760	0.0002	0.9999
K2PG	57	-1397.7905	3122.0113	0.0001	1.0000
HKYI	60	-1389.4449	3122.5006	0.0000	1.0000
K2PI	57	-1401.7596	3129.9495	0.0000	1.0000
SYMI	61	-1391.0339	3131.4056	0.0000	1.0000
GTRI	64	-1385.5599	3137.6380	0.0000	1.0000
F81IG	60	-1398.1111	3139.8330	0.0000	1.0000
F81G	59	-1403.0799	3144.0438	0.0000	1.0000
F81I	59	-1408.5965	3155.0770	0.0000	1.0000
JC69IG	57	-1420.5377	3167.5057	0.0000	1.0000
JC69G	56	-1425.2591	3171.2217	0.0000	1.0000
JC69I	56	-1434.4741	3189.6517	0.0000	1.0000
SYM	60	-1440.1759	3223.9627	0.0000	1.0000
HKY	59	-1446.6940	3231.2720	0.0000	1.0000
GTR	63	-1436.2460	3233.2833	0.0000	1.0000
K2P	56	-1460.9450	3242.5935	0.0000	1.0000
F81	58	-1470.7791	3273.7153	0.0000	1.0000
JC69	55	-1487.8262	3290.6289	0.0000	1.0000



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

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

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


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

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

