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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-1426.0077	3016.8296	0.8082	0.8082
K2PG	57	-1439.8372	3022.0520	0.0594	0.8676
HKYG	60	-1435.2312	3022.2885	0.0527	0.9203
GTRIG	65	-1427.2655	3022.6643	0.0437	0.9640
HKYIG	61	-1434.4680	3023.9666	0.0228	0.9868
K2PIG	58	-1440.8136	3027.1272	0.0047	0.9915
GTRI	64	-1431.4593	3027.7328	0.0035	0.9950
SYMG	61	-1436.6092	3028.2491	0.0027	0.9976
SYMIG	62	-1435.7991	3029.8614	0.0012	0.9988
HKYI	60	-1439.4770	3030.7801	0.0008	0.9996
K2PI	57	-1445.7811	3033.9398	0.0002	0.9997
SYMI	61	-1439.6331	3034.2967	0.0001	0.9999
F81G	59	-1443.8441	3036.3376	0.0000	0.9999
JC69G	56	-1448.7348	3036.7517	0.0000	1.0000
F81IG	60	-1443.0726	3037.9713	0.0000	1.0000
JC69IG	57	-1447.9958	3038.3692	0.0000	1.0000
F81I	59	-1447.8660	3044.3813	0.0000	1.0000
JC69I	56	-1452.9475	3045.1770	0.0000	1.0000
HKY	59	-1475.2217	3099.0928	0.0000	1.0000
GTR	63	-1469.0280	3099.5801	0.0000	1.0000
K2P	56	-1480.2059	3099.6938	0.0000	1.0000
SYM	60	-1475.2018	3102.2298	0.0000	1.0000
JC69	55	-1487.8006	3111.8140	0.0000	1.0000
F81	58	-1483.5431	3112.5863	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1426.0077	2980.0154	0.8922	0.8922
GTRIG	65	-1427.2655	2984.5309	0.0933	0.9856
HKYG	60	-1435.2312	2990.4624	0.0048	0.9904
GTRI	64	-1431.4593	2990.9186	0.0038	0.9942
HKYIG	61	-1434.4680	2990.9360	0.0038	0.9980
K2PG	57	-1439.8372	2993.6743	0.0010	0.9990
SYMG	61	-1436.6092	2995.2185	0.0004	0.9994
SYMIG	62	-1435.7991	2995.5983	0.0004	0.9998
K2PIG	58	-1440.8136	2997.6272	0.0001	0.9999
HKYI	60	-1439.4770	2998.9540	0.0001	1.0000
SYMI	61	-1439.6331	3001.2661	0.0000	1.0000
K2PI	57	-1445.7811	3005.5621	0.0000	1.0000
F81G	59	-1443.8441	3005.6882	0.0000	1.0000
F81IG	60	-1443.0726	3006.1453	0.0000	1.0000
JC69G	56	-1448.7348	3009.4696	0.0000	1.0000
JC69IG	57	-1447.9958	3009.9915	0.0000	1.0000
F81I	59	-1447.8660	3013.7319	0.0000	1.0000
JC69I	56	-1452.9475	3017.8950	0.0000	1.0000
GTR	63	-1469.0280	3064.0559	0.0000	1.0000
HKY	59	-1475.2217	3068.4435	0.0000	1.0000
SYM	60	-1475.2018	3070.4037	0.0000	1.0000
K2P	56	-1480.2059	3072.4117	0.0000	1.0000
F81	58	-1483.5431	3083.0863	0.0000	1.0000
JC69	55	-1487.8006	3085.6012	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1439.8372	3203.0537	0.9495	0.9495
K2PIG	58	-1440.8136	3210.6799	0.0210	0.9704
HKYG	60	-1435.2312	3210.8618	0.0191	0.9896
K2PI	57	-1445.7811	3214.9415	0.0025	0.9921
HKYIG	61	-1434.4680	3215.0088	0.0024	0.9945
GTRG	64	-1426.0077	3215.1081	0.0023	0.9968
JC69G	56	-1448.7348	3215.1757	0.0022	0.9990
SYMG	61	-1436.6092	3219.2912	0.0003	0.9993
HKYI	60	-1439.4770	3219.3534	0.0003	0.9995
JC69IG	57	-1447.9958	3219.3709	0.0003	0.9998
F81G	59	-1443.8441	3222.4143	0.0001	0.9999
GTRIG	65	-1427.2655	3223.2969	0.0000	0.9999
SYMIG	62	-1435.7991	3223.3443	0.0000	0.9999
JC69I	56	-1452.9475	3223.6011	0.0000	1.0000
SYMI	61	-1439.6331	3225.3389	0.0000	1.0000
GTRI	64	-1431.4593	3226.0113	0.0000	1.0000
F81IG	60	-1443.0726	3226.5447	0.0000	1.0000
F81I	59	-1447.8660	3230.4580	0.0000	1.0000
K2P	56	-1480.2059	3278.1178	0.0000	1.0000
HKY	59	-1475.2217	3285.1695	0.0000	1.0000
JC69	55	-1487.8006	3287.6340	0.0000	1.0000
SYM	60	-1475.2018	3290.8031	0.0000	1.0000
GTR	63	-1469.0280	3295.4753	0.0000	1.0000
F81	58	-1483.5431	3296.1390	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/chr12_5908.phylip = 1 - 291;
 Partition Dummy = 1:183-phylip/chr12_5908.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/chr12_5908.phylip = 1 - 291;
 Partition Dummy = 1:183-phylip/chr12_5908.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5908.phylip = 1 - 291;
 Partition Dummy = 1:183-phylip/chr12_5908.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

