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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	60	-1083.5117	2310.6363	0.5669	0.5669
HKYG	60	-1084.3573	2312.3275	0.2433	0.8102
HKYIG	61	-1083.5593	2313.5976	0.1290	0.9392
GTRI	64	-1080.4515	2316.0925	0.0370	0.9762
GTRG	64	-1081.1245	2317.4386	0.0189	0.9951
GTRIG	65	-1081.2428	2320.6168	0.0039	0.9989
HKY	59	-1091.2731	2323.3116	0.0010	1.0000
GTR	63	-1088.7731	2329.8134	0.0000	1.0000
F81I	59	-1096.4812	2333.7277	0.0000	1.0000
F81G	59	-1097.3178	2335.4009	0.0000	1.0000
F81IG	60	-1096.5047	2336.6224	0.0000	1.0000
F81	58	-1103.2351	2344.4062	0.0000	1.0000
SYMI	61	-1103.1814	2352.8417	0.0000	1.0000
K2PI	57	-1109.1456	2353.4158	0.0000	1.0000
SYMG	61	-1103.7136	2353.9062	0.0000	1.0000
K2PG	57	-1109.8705	2354.8657	0.0000	1.0000
K2PIG	58	-1109.1674	2356.2708	0.0000	1.0000
SYMIG	62	-1103.7856	2356.9349	0.0000	1.0000
K2P	56	-1116.2450	2364.8212	0.0000	1.0000
SYM	60	-1110.8131	2365.2391	0.0000	1.0000
JC69I	56	-1119.6104	2371.5521	0.0000	1.0000
JC69G	56	-1120.3424	2373.0160	0.0000	1.0000
JC69IG	57	-1119.6320	2374.3886	0.0000	1.0000
JC69	55	-1126.1953	2381.9463	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-1083.5117	2287.0234	0.4104	0.4104
HKYG	60	-1084.3573	2288.7146	0.1762	0.5866
GTRI	64	-1080.4515	2288.9030	0.1603	0.7469
HKYIG	61	-1083.5593	2289.1186	0.1440	0.8909
GTRG	64	-1081.1245	2290.2491	0.0818	0.9727
GTRIG	65	-1081.2428	2292.4856	0.0267	0.9994
HKY	59	-1091.2731	2300.5463	0.0005	0.9999
GTR	63	-1088.7731	2303.5463	0.0001	1.0000
F81I	59	-1096.4812	2310.9624	0.0000	1.0000
F81G	59	-1097.3178	2312.6357	0.0000	1.0000
F81IG	60	-1096.5047	2313.0095	0.0000	1.0000
F81	58	-1103.2351	2322.4703	0.0000	1.0000
SYMI	61	-1103.1814	2328.3627	0.0000	1.0000
SYMG	61	-1103.7136	2329.4273	0.0000	1.0000
SYMIG	62	-1103.7856	2331.5712	0.0000	1.0000
K2PI	57	-1109.1456	2332.2912	0.0000	1.0000
K2PG	57	-1109.8705	2333.7411	0.0000	1.0000
K2PIG	58	-1109.1674	2334.3349	0.0000	1.0000
SYM	60	-1110.8131	2341.6262	0.0000	1.0000
K2P	56	-1116.2450	2344.4900	0.0000	1.0000
JC69I	56	-1119.6104	2351.2209	0.0000	1.0000
JC69G	56	-1120.3424	2352.6848	0.0000	1.0000
JC69IG	57	-1119.6320	2353.2640	0.0000	1.0000
JC69	55	-1126.1953	2362.3907	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1083.5117	2521.9955	0.6724	0.6724
HKYG	60	-1084.3573	2523.6867	0.2886	0.9610
HKYIG	61	-1083.5593	2528.0070	0.0333	0.9943
HKY	59	-1091.2731	2531.6022	0.0055	0.9998
GTRI	64	-1080.4515	2539.5399	0.0001	0.9999
GTRG	64	-1081.1245	2540.8860	0.0001	1.0000
F81I	59	-1096.4812	2542.0183	0.0000	1.0000
F81G	59	-1097.3178	2543.6916	0.0000	1.0000
GTRIG	65	-1081.2428	2547.0388	0.0000	1.0000
F81IG	60	-1096.5047	2547.9816	0.0000	1.0000
F81	58	-1103.2351	2549.6100	0.0000	1.0000
GTR	63	-1088.7731	2550.2670	0.0000	1.0000
K2PI	57	-1109.1456	2555.5148	0.0000	1.0000
K2PG	57	-1109.8705	2556.9646	0.0000	1.0000
K2PIG	58	-1109.1674	2561.4746	0.0000	1.0000
K2P	56	-1116.2450	2563.7973	0.0000	1.0000
SYMI	61	-1103.1814	2567.2511	0.0000	1.0000
SYMG	61	-1103.7136	2568.3156	0.0000	1.0000
JC69I	56	-1119.6104	2570.5282	0.0000	1.0000
JC69G	56	-1120.3424	2571.9921	0.0000	1.0000
SYMIG	62	-1103.7856	2574.3757	0.0000	1.0000
JC69IG	57	-1119.6320	2576.4875	0.0000	1.0000
SYM	60	-1110.8131	2576.5983	0.0000	1.0000
JC69	55	-1126.1953	2577.7818	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10841.phylip = 1 - 371;
 Partition Dummy = 1:183-phylip/chr5_10841.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10841.phylip = 1 - 371;
 Partition Dummy = 1:183-phylip/chr5_10841.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10841.phylip = 1 - 371;
 Partition Dummy = 1:183-phylip/chr5_10841.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

