
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:47 04/04/2011

Input data from file "591-tomodel-phylip/chr7_10997.phylip" (7 taxa, 285 characters)

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: HKY


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
HKY	15	-672.3005	1376.3853	0.4543	0.4543
HKYI	16	-671.9016	1377.8331	0.2203	0.6745
HKYG	16	-671.9675	1377.9649	0.2062	0.8807
HKYIG	17	-671.9039	1380.0999	0.0709	0.9516
GTR	19	-670.7746	1382.4170	0.0223	0.9739
GTRI	20	-670.2592	1383.7002	0.0117	0.9856
GTRG	20	-670.3452	1383.8722	0.0108	0.9964
GTRIG	21	-670.2622	1386.0378	0.0036	1.0000
F81	14	-688.3620	1406.2796	0.0000	1.0000
F81I	15	-688.0519	1407.8882	0.0000	1.0000
F81G	15	-688.1318	1408.0479	0.0000	1.0000
F81IG	16	-688.0536	1410.1370	0.0000	1.0000
K2P	12	-692.6649	1410.4769	0.0000	1.0000
K2PI	13	-692.2170	1411.7771	0.0000	1.0000
K2PG	13	-692.2732	1411.8896	0.0000	1.0000
K2PIG	14	-692.2187	1413.9930	0.0000	1.0000
SYM	16	-691.9860	1418.0018	0.0000	1.0000
SYMI	17	-691.4989	1419.2900	0.0000	1.0000
SYMG	17	-691.5594	1419.4110	0.0000	1.0000
SYMIG	18	-691.5414	1421.6542	0.0000	1.0000
JC69	11	-706.0281	1435.0233	0.0000	1.0000
JC69I	12	-705.6656	1436.4782	0.0000	1.0000
JC69G	12	-705.7421	1436.6312	0.0000	1.0000
JC69IG	13	-705.6680	1438.6791	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-672.3005	1374.6009	0.4056	0.4056
HKYI	16	-671.9016	1375.8032	0.2224	0.6280
HKYG	16	-671.9675	1375.9351	0.2082	0.8361
HKYIG	17	-671.9039	1377.8077	0.0816	0.9178
GTR	19	-670.7746	1379.5491	0.0342	0.9519
GTRI	20	-670.2592	1380.5184	0.0210	0.9730
GTRG	20	-670.3452	1380.6904	0.0193	0.9923
GTRIG	21	-670.2622	1382.5245	0.0077	1.0000
F81	14	-688.3620	1404.7240	0.0000	1.0000
F81I	15	-688.0519	1406.1038	0.0000	1.0000
F81G	15	-688.1318	1406.2635	0.0000	1.0000
F81IG	16	-688.0536	1408.1071	0.0000	1.0000
K2P	12	-692.6649	1409.3298	0.0000	1.0000
K2PI	13	-692.2170	1410.4339	0.0000	1.0000
K2PG	13	-692.2732	1410.5464	0.0000	1.0000
K2PIG	14	-692.2187	1412.4375	0.0000	1.0000
SYM	16	-691.9860	1415.9720	0.0000	1.0000
SYMI	17	-691.4989	1416.9979	0.0000	1.0000
SYMG	17	-691.5594	1417.1188	0.0000	1.0000
SYMIG	18	-691.5414	1419.0828	0.0000	1.0000
JC69	11	-706.0281	1434.0563	0.0000	1.0000
JC69I	12	-705.6656	1435.3311	0.0000	1.0000
JC69G	12	-705.7421	1435.4841	0.0000	1.0000
JC69IG	13	-705.6680	1437.3359	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-672.3005	1429.3883	0.8502	0.8502
HKYI	16	-671.9016	1434.2430	0.0750	0.9252
HKYG	16	-671.9675	1434.3749	0.0703	0.9955
HKYIG	17	-671.9039	1439.9000	0.0044	0.9999
GTR	19	-670.7746	1448.9464	0.0000	1.0000
K2P	12	-692.6649	1453.1597	0.0000	1.0000
GTRI	20	-670.2592	1453.5682	0.0000	1.0000
GTRG	20	-670.3452	1453.7402	0.0000	1.0000
F81	14	-688.3620	1455.8588	0.0000	1.0000
K2PI	13	-692.2170	1457.9163	0.0000	1.0000
K2PG	13	-692.2732	1458.0288	0.0000	1.0000
GTRIG	21	-670.2622	1459.2267	0.0000	1.0000
F81I	15	-688.0519	1460.8911	0.0000	1.0000
F81G	15	-688.1318	1461.0508	0.0000	1.0000
K2PIG	14	-692.2187	1463.5723	0.0000	1.0000
F81IG	16	-688.0536	1466.5470	0.0000	1.0000
JC69	11	-706.0281	1474.2337	0.0000	1.0000
SYM	16	-691.9860	1474.4118	0.0000	1.0000
SYMI	17	-691.4989	1479.0902	0.0000	1.0000
JC69I	12	-705.6656	1479.1610	0.0000	1.0000
SYMG	17	-691.5594	1479.2112	0.0000	1.0000
JC69G	12	-705.7421	1479.3140	0.0000	1.0000
JC69IG	13	-705.6680	1484.8183	0.0000	1.0000
SYMIG	18	-691.5414	1484.8276	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10997.phylip = 1 - 285;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10997.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10997.phylip = 1 - 285;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10997.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10997.phylip = 1 - 285;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10997.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

