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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKY
Minimum BIC  model: F81


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	59	-627.2329	1399.5922	0.3143	0.3143
HKYI	60	-626.0881	1400.3301	0.2174	0.5317
HKYG	60	-626.1218	1400.3973	0.2102	0.7419
F81	58	-630.1096	1402.3413	0.0795	0.8214
F81I	59	-629.0451	1403.2166	0.0513	0.8727
F81G	59	-629.0829	1403.2922	0.0494	0.9221
HKYIG	61	-626.1213	1403.4472	0.0457	0.9679
F81IG	60	-629.0913	1406.3364	0.0108	0.9786
GTR	63	-625.0899	1407.5572	0.0059	0.9845
GTRI	64	-623.8423	1408.1845	0.0043	0.9888
GTRG	64	-623.8796	1408.2593	0.0041	0.9929
K2P	56	-636.8437	1409.8692	0.0018	0.9947
K2PI	57	-635.7757	1410.6920	0.0012	0.9960
K2PG	57	-635.8131	1410.7669	0.0012	0.9971
GTRIG	65	-623.8794	1411.4059	0.0009	0.9980
JC69	55	-639.3207	1411.8867	0.0007	0.9987
JC69I	56	-638.2896	1412.7611	0.0004	0.9991
JC69G	56	-638.3290	1412.8398	0.0004	0.9995
K2PIG	58	-635.8217	1413.7655	0.0003	0.9998
JC69IG	57	-638.3372	1415.8151	0.0001	0.9999
SYM	60	-634.5175	1417.1888	0.0000	0.9999
SYMI	61	-633.4022	1418.0091	0.0000	1.0000
SYMG	61	-633.4471	1418.0989	0.0000	1.0000
SYMIG	62	-633.4469	1421.1729	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-626.0881	1372.1762	0.2387	0.2387
HKYG	60	-626.1218	1372.2435	0.2308	0.4695
HKY	59	-627.2329	1372.4657	0.2065	0.6760
HKYIG	61	-626.1213	1374.2426	0.0849	0.7610
GTRI	64	-623.8423	1375.6845	0.0413	0.8023
GTRG	64	-623.8796	1375.7593	0.0398	0.8421
F81I	59	-629.0451	1376.0902	0.0337	0.8758
F81G	59	-629.0829	1376.1658	0.0325	0.9083
GTR	63	-625.0899	1376.1797	0.0322	0.9405
F81	58	-630.1096	1376.2192	0.0316	0.9721
GTRIG	65	-623.8794	1377.7588	0.0146	0.9868
F81IG	60	-629.0913	1378.1825	0.0118	0.9986
K2PI	57	-635.7757	1385.5513	0.0003	0.9989
K2PG	57	-635.8131	1385.6262	0.0003	0.9992
K2P	56	-636.8437	1385.6874	0.0003	0.9995
K2PIG	58	-635.8217	1387.6433	0.0001	0.9996
JC69I	56	-638.2896	1388.5793	0.0001	0.9997
JC69	55	-639.3207	1388.6415	0.0001	0.9997
JC69G	56	-638.3290	1388.6579	0.0001	0.9998
SYMI	61	-633.4022	1388.8045	0.0001	0.9998
SYMG	61	-633.4471	1388.8942	0.0001	0.9999
SYM	60	-634.5175	1389.0349	0.0001	1.0000
JC69IG	57	-638.3372	1390.6745	0.0000	1.0000
SYMIG	62	-633.4469	1390.8938	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
F81	58	-630.1096	1594.9628	0.2544	0.2544
HKY	59	-627.2329	1594.9808	0.2521	0.5064
JC69	55	-639.3207	1596.0707	0.1462	0.6526
K2P	56	-636.8437	1596.8881	0.0971	0.7497
HKYI	60	-626.0881	1598.4627	0.0442	0.7939
HKYG	60	-626.1218	1598.5300	0.0427	0.8367
F81I	59	-629.0451	1598.6052	0.0412	0.8778
F81G	59	-629.0829	1598.6808	0.0396	0.9175
JC69I	56	-638.2896	1599.7800	0.0229	0.9404
JC69G	56	-638.3290	1599.8586	0.0220	0.9623
K2PI	57	-635.7757	1600.5235	0.0158	0.9781
K2PG	57	-635.8131	1600.5983	0.0152	0.9933
HKYIG	61	-626.1213	1604.3005	0.0024	0.9957
F81IG	60	-629.0913	1604.4690	0.0022	0.9979
JC69IG	57	-638.3372	1605.6466	0.0012	0.9991
K2PIG	58	-635.8217	1606.3869	0.0008	1.0000
GTR	63	-625.0899	1613.7805	0.0000	1.0000
SYM	60	-634.5175	1615.3214	0.0000	1.0000
GTRI	64	-623.8423	1617.0567	0.0000	1.0000
GTRG	64	-623.8796	1617.1315	0.0000	1.0000
SYMI	61	-633.4022	1618.8624	0.0000	1.0000
SYMG	61	-633.4471	1618.9521	0.0000	1.0000
GTRIG	65	-623.8794	1622.9025	0.0000	1.0000
SYMIG	62	-633.4469	1624.7232	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/chr1_32215.phylip = 1 - 321;
 Partition Dummy = 1:183-phylip/chr1_32215.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 (HKY)]

BEGIN MRBAYES;
 Charset 183-phylip/chr1_32215.phylip = 1 - 321;
 Partition Dummy = 1:183-phylip/chr1_32215.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 (F81)]

BEGIN MRBAYES;
 Charset 183-phylip/chr1_32215.phylip = 1 - 321;
 Partition Dummy = 1:183-phylip/chr1_32215.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0);
END;


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

