
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:27 08/01/2011

Input data from file "261-aggressive-phylip/4956_hg18.phylip" (29 taxa, 461 characters)

Minimum AIC  model: GTRIG
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	-1621.6021	3392.2144	0.2457	0.2457
HKYIG	61	-1625.6883	3392.3341	0.2315	0.4772
GTRIG	65	-1620.4017	3392.5250	0.2104	0.6876
HKYG	60	-1627.1712	3392.6423	0.1984	0.8860
HKYI	60	-1628.3547	3395.0095	0.0607	0.9467
GTRI	64	-1623.1601	3395.3304	0.0517	0.9984
K2PIG	58	-1634.9267	3402.8783	0.0012	0.9996
K2PG	57	-1637.7060	3405.8189	0.0003	0.9999
K2PI	57	-1639.6868	3409.7805	0.0000	0.9999
SYMIG	62	-1633.3505	3410.3292	0.0000	1.0000
SYMG	61	-1634.8095	3410.5764	0.0000	1.0000
SYMI	61	-1635.7042	3412.3658	0.0000	1.0000
F81IG	60	-1668.1980	3474.6959	0.0000	1.0000
F81G	59	-1669.7139	3475.0837	0.0000	1.0000
F81I	59	-1671.0297	3477.7153	0.0000	1.0000
JC69G	56	-1677.6408	3483.0837	0.0000	1.0000
JC69I	56	-1679.0448	3485.8915	0.0000	1.0000
JC69IG	57	-1678.4404	3487.2878	0.0000	1.0000
GTR	63	-1670.9495	3488.2113	0.0000	1.0000
HKY	59	-1678.2393	3492.1344	0.0000	1.0000
K2P	56	-1691.7124	3511.2267	0.0000	1.0000
SYM	60	-1687.0065	3512.3130	0.0000	1.0000
F81	58	-1719.5056	3572.0362	0.0000	1.0000
JC69	55	-1730.6540	3586.5180	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1620.4017	3370.8035	0.4016	0.4016
GTRG	64	-1621.6021	3371.2043	0.3286	0.7302
HKYIG	61	-1625.6883	3373.3767	0.1109	0.8411
GTRI	64	-1623.1601	3374.3203	0.0692	0.9103
HKYG	60	-1627.1712	3374.3423	0.0684	0.9788
HKYI	60	-1628.3547	3376.7095	0.0210	0.9997
K2PIG	58	-1634.9267	3385.8535	0.0002	0.9999
K2PG	57	-1637.7060	3389.4120	0.0000	1.0000
SYMIG	62	-1633.3505	3390.7010	0.0000	1.0000
SYMG	61	-1634.8095	3391.6190	0.0000	1.0000
K2PI	57	-1639.6868	3393.3735	0.0000	1.0000
SYMI	61	-1635.7042	3393.4084	0.0000	1.0000
F81IG	60	-1668.1980	3456.3959	0.0000	1.0000
F81G	59	-1669.7139	3457.4278	0.0000	1.0000
F81I	59	-1671.0297	3460.0595	0.0000	1.0000
JC69G	56	-1677.6408	3467.2817	0.0000	1.0000
GTR	63	-1670.9495	3467.8990	0.0000	1.0000
JC69I	56	-1679.0448	3470.0895	0.0000	1.0000
JC69IG	57	-1678.4404	3470.8808	0.0000	1.0000
HKY	59	-1678.2393	3474.4785	0.0000	1.0000
SYM	60	-1687.0065	3494.0130	0.0000	1.0000
K2P	56	-1691.7124	3495.4247	0.0000	1.0000
F81	58	-1719.5056	3555.0113	0.0000	1.0000
JC69	55	-1730.6540	3571.3081	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1627.1712	3622.3462	0.4975	0.4975
HKYI	60	-1628.3547	3624.7134	0.1523	0.6498
K2PG	57	-1637.7060	3625.0157	0.1310	0.7808
HKYIG	61	-1625.6883	3625.5140	0.1021	0.8828
K2PIG	58	-1634.9267	3625.5905	0.0982	0.9811
K2PI	57	-1639.6868	3628.9772	0.0181	0.9991
GTRG	64	-1621.6021	3635.7418	0.0006	0.9998
GTRI	64	-1623.1601	3638.8577	0.0001	0.9999
GTRIG	65	-1620.4017	3639.4743	0.0001	1.0000
SYMG	61	-1634.8095	3643.7563	0.0000	1.0000
SYMI	61	-1635.7042	3645.5457	0.0000	1.0000
SYMIG	62	-1633.3505	3646.9717	0.0000	1.0000
JC69G	56	-1677.6408	3698.7520	0.0000	1.0000
F81G	59	-1669.7139	3701.2983	0.0000	1.0000
JC69I	56	-1679.0448	3701.5598	0.0000	1.0000
F81I	59	-1671.0297	3703.9299	0.0000	1.0000
F81IG	60	-1668.1980	3704.3998	0.0000	1.0000
JC69IG	57	-1678.4404	3706.4845	0.0000	1.0000
HKY	59	-1678.2393	3718.3490	0.0000	1.0000
K2P	56	-1691.7124	3726.8950	0.0000	1.0000
GTR	63	-1670.9495	3728.3031	0.0000	1.0000
SYM	60	-1687.0065	3742.0169	0.0000	1.0000
F81	58	-1719.5056	3794.7484	0.0000	1.0000
JC69	55	-1730.6540	3798.6450	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4956_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/4956_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4956_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/4956_hg18.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 261-aggressive-phylip/4956_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/4956_hg18.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

