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

Input data from file "183-phylip/chr22_640.phylip" (29 taxa, 378 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	-1517.6123	3189.8060	0.5192	0.5192
GTRIG	65	-1517.0343	3191.5687	0.2151	0.7343
HKYG	60	-1524.8085	3192.7085	0.1216	0.8559
GTRI	64	-1519.6694	3193.9203	0.0664	0.9223
HKYIG	61	-1524.1935	3194.3237	0.0542	0.9765
HKYI	60	-1526.4543	3196.0000	0.0235	1.0000
K2PG	57	-1538.6984	3212.0592	0.0000	1.0000
K2PIG	58	-1537.9951	3213.4447	0.0000	1.0000
K2PI	57	-1540.4099	3215.4824	0.0000	1.0000
SYMG	61	-1535.3298	3216.5964	0.0000	1.0000
SYMIG	62	-1534.6369	3218.0738	0.0000	1.0000
SYMI	61	-1537.0268	3219.9904	0.0000	1.0000
GTR	63	-1544.9370	3241.5555	0.0000	1.0000
HKY	59	-1553.3529	3246.9699	0.0000	1.0000
SYM	60	-1562.0786	3267.2486	0.0000	1.0000
K2P	56	-1569.2293	3270.3465	0.0000	1.0000
F81G	59	-1568.4026	3277.0693	0.0000	1.0000
F81IG	60	-1567.8853	3278.8621	0.0000	1.0000
F81I	59	-1569.6782	3279.6206	0.0000	1.0000
JC69G	56	-1580.4916	3292.8711	0.0000	1.0000
JC69IG	57	-1579.9394	3294.5413	0.0000	1.0000
JC69I	56	-1581.8087	3295.5052	0.0000	1.0000
F81	58	-1596.7331	3330.9208	0.0000	1.0000
JC69	55	-1609.5970	3348.3243	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1517.6123	3163.2246	0.5376	0.5376
GTRIG	65	-1517.0343	3164.0687	0.3525	0.8901
GTRI	64	-1519.6694	3167.3388	0.0687	0.9588
HKYG	60	-1524.8085	3169.6170	0.0220	0.9808
HKYIG	61	-1524.1935	3170.3870	0.0150	0.9958
HKYI	60	-1526.4543	3172.9085	0.0042	1.0000
K2PG	57	-1538.6984	3191.3967	0.0000	1.0000
K2PIG	58	-1537.9951	3191.9901	0.0000	1.0000
SYMG	61	-1535.3298	3192.6597	0.0000	1.0000
SYMIG	62	-1534.6369	3193.2738	0.0000	1.0000
K2PI	57	-1540.4099	3194.8199	0.0000	1.0000
SYMI	61	-1537.0268	3196.0537	0.0000	1.0000
GTR	63	-1544.9370	3215.8740	0.0000	1.0000
HKY	59	-1553.3529	3224.7057	0.0000	1.0000
SYM	60	-1562.0786	3244.1571	0.0000	1.0000
K2P	56	-1569.2293	3250.4586	0.0000	1.0000
F81G	59	-1568.4026	3254.8051	0.0000	1.0000
F81IG	60	-1567.8853	3255.7706	0.0000	1.0000
F81I	59	-1569.6782	3257.3564	0.0000	1.0000
JC69G	56	-1580.4916	3272.9832	0.0000	1.0000
JC69IG	57	-1579.9394	3273.8788	0.0000	1.0000
JC69I	56	-1581.8087	3275.6174	0.0000	1.0000
F81	58	-1596.7331	3309.4662	0.0000	1.0000
JC69	55	-1609.5970	3329.1939	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1524.8085	3405.7107	0.7644	0.7644
HKYI	60	-1526.4543	3409.0022	0.1474	0.9119
HKYIG	61	-1524.1935	3410.4156	0.0727	0.9846
GTRG	64	-1517.6123	3415.0578	0.0071	0.9917
K2PG	57	-1538.6984	3415.6857	0.0052	0.9969
K2PI	57	-1540.4099	3419.1089	0.0009	0.9979
GTRI	64	-1519.6694	3419.1720	0.0009	0.9988
GTRIG	65	-1517.0343	3419.8368	0.0007	0.9995
K2PIG	58	-1537.9951	3420.2140	0.0005	1.0000
SYMG	61	-1535.3298	3432.6882	0.0000	1.0000
SYMI	61	-1537.0268	3436.0822	0.0000	1.0000
SYMIG	62	-1534.6369	3437.2372	0.0000	1.0000
HKY	59	-1553.3529	3456.8645	0.0000	1.0000
GTR	63	-1544.9370	3463.7724	0.0000	1.0000
K2P	56	-1569.2293	3470.8127	0.0000	1.0000
SYM	60	-1562.0786	3480.2508	0.0000	1.0000
F81G	59	-1568.4026	3486.9639	0.0000	1.0000
F81I	59	-1569.6782	3489.5152	0.0000	1.0000
F81IG	60	-1567.8853	3491.8642	0.0000	1.0000
JC69G	56	-1580.4916	3493.3373	0.0000	1.0000
JC69I	56	-1581.8087	3495.9715	0.0000	1.0000
JC69IG	57	-1579.9394	3498.1678	0.0000	1.0000
F81	58	-1596.7331	3537.6901	0.0000	1.0000
JC69	55	-1609.5970	3545.6131	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/chr22_640.phylip = 1 - 378;
 Partition Dummy = 1:183-phylip/chr22_640.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/chr22_640.phylip = 1 - 378;
 Partition Dummy = 1:183-phylip/chr22_640.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/chr22_640.phylip = 1 - 378;
 Partition Dummy = 1:183-phylip/chr22_640.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

