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

Input data from file "183-phylip/chr3_5834.phylip" (29 taxa, 267 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	-1100.3683	2369.9247	0.6320	0.6320
GTRI	64	-1101.2988	2371.7858	0.2492	0.8812
GTRIG	65	-1100.8706	2374.4278	0.0665	0.9477
HKYG	60	-1110.2914	2376.1168	0.0286	0.9763
HKYI	60	-1110.8413	2377.2166	0.0165	0.9928
HKYIG	61	-1110.0232	2378.9440	0.0070	0.9997
K2PG	57	-1120.8378	2387.3119	0.0001	0.9998
SYMG	61	-1114.7848	2388.4672	0.0001	0.9999
SYMI	61	-1115.2581	2389.4138	0.0000	0.9999
K2PI	57	-1122.0143	2389.6649	0.0000	1.0000
K2PIG	58	-1120.9337	2390.7712	0.0000	1.0000
SYMIG	62	-1114.6018	2391.4976	0.0000	1.0000
F81G	59	-1132.6780	2417.5589	0.0000	1.0000
GTR	63	-1126.3550	2418.4342	0.0000	1.0000
F81I	59	-1133.1636	2418.5300	0.0000	1.0000
HKY	59	-1133.9180	2420.0388	0.0000	1.0000
F81IG	60	-1132.4815	2420.4970	0.0000	1.0000
JC69G	56	-1140.1605	2422.7209	0.0000	1.0000
JC69I	56	-1140.7941	2423.9881	0.0000	1.0000
JC69IG	57	-1140.0012	2425.6387	0.0000	1.0000
K2P	56	-1143.9940	2430.3880	0.0000	1.0000
SYM	60	-1139.1444	2433.8227	0.0000	1.0000
F81	58	-1155.1237	2459.1513	0.0000	1.0000
JC69	55	-1162.6916	2464.5776	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1100.3683	2328.7366	0.6163	0.6163
GTRI	64	-1101.2988	2330.5977	0.2430	0.8594
GTRIG	65	-1100.8706	2331.7412	0.1372	0.9966
HKYG	60	-1110.2914	2340.5828	0.0016	0.9982
HKYI	60	-1110.8413	2341.6826	0.0010	0.9992
HKYIG	61	-1110.0232	2342.0464	0.0008	1.0000
SYMG	61	-1114.7848	2351.5696	0.0000	1.0000
SYMI	61	-1115.2581	2352.5162	0.0000	1.0000
SYMIG	62	-1114.6018	2353.2035	0.0000	1.0000
K2PG	57	-1120.8378	2355.6756	0.0000	1.0000
K2PIG	58	-1120.9337	2357.8674	0.0000	1.0000
K2PI	57	-1122.0143	2358.0286	0.0000	1.0000
GTR	63	-1126.3550	2378.7101	0.0000	1.0000
F81G	59	-1132.6780	2383.3560	0.0000	1.0000
F81I	59	-1133.1636	2384.3271	0.0000	1.0000
F81IG	60	-1132.4815	2384.9630	0.0000	1.0000
HKY	59	-1133.9180	2385.8359	0.0000	1.0000
JC69G	56	-1140.1605	2392.3209	0.0000	1.0000
JC69I	56	-1140.7941	2393.5881	0.0000	1.0000
JC69IG	57	-1140.0012	2394.0023	0.0000	1.0000
SYM	60	-1139.1444	2398.2888	0.0000	1.0000
K2P	56	-1143.9940	2399.9880	0.0000	1.0000
F81	58	-1155.1237	2426.2475	0.0000	1.0000
JC69	55	-1162.6916	2435.3832	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1110.2914	2555.8177	0.4496	0.4496
HKYI	60	-1110.8413	2556.9176	0.2594	0.7090
GTRG	64	-1100.3683	2558.3205	0.1286	0.8377
K2PG	57	-1120.8378	2560.1488	0.0516	0.8892
GTRI	64	-1101.2988	2560.1816	0.0507	0.9400
HKYIG	61	-1110.0232	2560.8686	0.0360	0.9759
K2PI	57	-1122.0143	2562.5017	0.0159	0.9918
GTRIG	65	-1100.8706	2564.9124	0.0048	0.9966
K2PIG	58	-1120.9337	2565.9278	0.0029	0.9995
SYMG	61	-1114.7848	2570.3918	0.0003	0.9998
SYMI	61	-1115.2581	2571.3384	0.0002	1.0000
SYMIG	62	-1114.6018	2575.6129	0.0000	1.0000
JC69G	56	-1140.1605	2593.2069	0.0000	1.0000
JC69I	56	-1140.7941	2594.4741	0.0000	1.0000
F81G	59	-1132.6780	2595.0037	0.0000	1.0000
F81I	59	-1133.1636	2595.9748	0.0000	1.0000
HKY	59	-1133.9180	2597.4836	0.0000	1.0000
JC69IG	57	-1140.0012	2598.4755	0.0000	1.0000
F81IG	60	-1132.4815	2600.1979	0.0000	1.0000
K2P	56	-1143.9940	2600.8739	0.0000	1.0000
GTR	63	-1126.3550	2604.7067	0.0000	1.0000
SYM	60	-1139.1444	2613.5237	0.0000	1.0000
JC69	55	-1162.6916	2632.6819	0.0000	1.0000
F81	58	-1155.1237	2634.3079	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/chr3_5834.phylip = 1 - 267;
 Partition Dummy = 1:183-phylip/chr3_5834.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/chr3_5834.phylip = 1 - 267;
 Partition Dummy = 1:183-phylip/chr3_5834.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/chr3_5834.phylip = 1 - 267;
 Partition Dummy = 1:183-phylip/chr3_5834.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

