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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRIG


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
GTRIG	65	-1554.7066	3258.1878	0.9988	0.9988
GTRG	64	-1562.8020	3271.7699	0.0011	0.9999
GTRI	64	-1565.4281	3277.0222	0.0001	1.0000
HKYIG	61	-1569.9616	3278.3311	0.0000	1.0000
HKYG	60	-1576.2719	3288.3880	0.0000	1.0000
HKYI	60	-1576.6021	3289.0485	0.0000	1.0000
SYMIG	62	-1587.4511	3315.8848	0.0000	1.0000
SYMG	61	-1591.3990	3321.2057	0.0000	1.0000
SYMI	61	-1594.0204	3326.4486	0.0000	1.0000
K2PIG	58	-1603.7010	3338.1519	0.0000	1.0000
K2PI	57	-1610.9099	3350.0392	0.0000	1.0000
K2PG	57	-1611.3287	3350.8768	0.0000	1.0000
GTR	63	-1635.3306	3414.2298	0.0000	1.0000
F81IG	60	-1639.2666	3414.3773	0.0000	1.0000
F81I	59	-1645.3157	3423.9230	0.0000	1.0000
F81G	59	-1645.6408	3424.5731	0.0000	1.0000
HKY	59	-1645.7790	3424.8496	0.0000	1.0000
SYM	60	-1665.2889	3466.4219	0.0000	1.0000
JC69IG	57	-1673.5264	3475.2722	0.0000	1.0000
JC69I	56	-1678.5029	3482.7053	0.0000	1.0000
JC69G	56	-1679.1544	3484.0084	0.0000	1.0000
K2P	56	-1683.1371	3491.9737	0.0000	1.0000
F81	58	-1715.1052	3560.9604	0.0000	1.0000
JC69	55	-1747.5871	3618.3648	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1554.7066	3239.4132	0.9991	0.9991
GTRG	64	-1562.8020	3253.6040	0.0008	0.9999
GTRI	64	-1565.4281	3258.8562	0.0001	1.0000
HKYIG	61	-1569.9616	3261.9233	0.0000	1.0000
HKYG	60	-1576.2719	3272.5439	0.0000	1.0000
HKYI	60	-1576.6021	3273.2043	0.0000	1.0000
SYMIG	62	-1587.4511	3298.9022	0.0000	1.0000
SYMG	61	-1591.3990	3304.7979	0.0000	1.0000
SYMI	61	-1594.0204	3310.0408	0.0000	1.0000
K2PIG	58	-1603.7010	3323.4019	0.0000	1.0000
K2PI	57	-1610.9099	3335.8198	0.0000	1.0000
K2PG	57	-1611.3287	3336.6574	0.0000	1.0000
GTR	63	-1635.3306	3396.6612	0.0000	1.0000
F81IG	60	-1639.2666	3398.5332	0.0000	1.0000
F81I	59	-1645.3157	3408.6314	0.0000	1.0000
F81G	59	-1645.6408	3409.2815	0.0000	1.0000
HKY	59	-1645.7790	3409.5580	0.0000	1.0000
SYM	60	-1665.2889	3450.5778	0.0000	1.0000
JC69IG	57	-1673.5264	3461.0529	0.0000	1.0000
JC69I	56	-1678.5029	3469.0057	0.0000	1.0000
JC69G	56	-1679.1544	3470.3088	0.0000	1.0000
K2P	56	-1683.1371	3478.2742	0.0000	1.0000
F81	58	-1715.1052	3546.2104	0.0000	1.0000
JC69	55	-1747.5871	3605.1742	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRIG	65	-1554.7066	3516.2860	0.9286	0.9286
HKYIG	61	-1569.9616	3521.7578	0.0602	0.9888
GTRG	64	-1562.8020	3526.2172	0.0065	0.9952
HKYG	60	-1576.2719	3528.1187	0.0025	0.9977
HKYI	60	-1576.6021	3528.7792	0.0018	0.9995
GTRI	64	-1565.4281	3531.4695	0.0005	1.0000
SYMIG	62	-1587.4511	3562.9963	0.0000	1.0000
SYMG	61	-1591.3990	3564.6324	0.0000	1.0000
SYMI	61	-1594.0204	3569.8752	0.0000	1.0000
K2PIG	58	-1603.7010	3570.4577	0.0000	1.0000
K2PI	57	-1610.9099	3578.6160	0.0000	1.0000
K2PG	57	-1611.3287	3579.4536	0.0000	1.0000
F81IG	60	-1639.2666	3654.1081	0.0000	1.0000
F81I	59	-1645.3157	3659.9467	0.0000	1.0000
F81G	59	-1645.6408	3660.5968	0.0000	1.0000
HKY	59	-1645.7790	3660.8733	0.0000	1.0000
GTR	63	-1635.3306	3665.0149	0.0000	1.0000
JC69IG	57	-1673.5264	3703.8490	0.0000	1.0000
SYM	60	-1665.2889	3706.1526	0.0000	1.0000
JC69I	56	-1678.5029	3707.5423	0.0000	1.0000
JC69G	56	-1679.1544	3708.8453	0.0000	1.0000
K2P	56	-1683.1371	3716.8107	0.0000	1.0000
F81	58	-1715.1052	3793.2662	0.0000	1.0000
JC69	55	-1747.5871	3839.4512	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/3860_hg18.phylip = 1 - 523;
 Partition Dummy = 1:261-aggressive-phylip/3860_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 (GTRIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/3860_hg18.phylip = 1 - 523;
 Partition Dummy = 1:261-aggressive-phylip/3860_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 BIC model (GTRIG)]

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


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

