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

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

Minimum AIC  model: HKYG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	57	-1624.2629	3378.4200	0.3648	0.3648
HKYG	60	-1620.9259	3379.5757	0.2047	0.5694
K2PI	57	-1625.3089	3380.5121	0.1281	0.6976
K2PIG	58	-1624.1963	3380.8842	0.1064	0.8040
HKYI	60	-1622.1520	3382.0280	0.0601	0.8640
HKYIG	61	-1620.8442	3382.0476	0.0595	0.9235
SYMG	61	-1621.4447	3383.2487	0.0326	0.9561
GTRG	64	-1618.3078	3384.9580	0.0139	0.9700
SYMI	61	-1622.3932	3385.1456	0.0126	0.9826
SYMIG	62	-1621.3912	3385.7896	0.0092	0.9918
GTRI	64	-1619.4748	3387.2919	0.0043	0.9961
GTRIG	65	-1618.2368	3387.5030	0.0039	1.0000
K2P	56	-1640.8816	3409.0726	0.0000	1.0000
HKY	59	-1637.7154	3410.5323	0.0000	1.0000
SYM	60	-1638.1336	3413.9912	0.0000	1.0000
GTR	63	-1635.1821	3416.0325	0.0000	1.0000
JC69G	56	-1655.8324	3438.9741	0.0000	1.0000
JC69I	56	-1656.8269	3440.9632	0.0000	1.0000
JC69IG	57	-1655.7727	3441.4397	0.0000	1.0000
F81G	59	-1653.1978	3441.4971	0.0000	1.0000
F81I	59	-1654.2636	3443.6286	0.0000	1.0000
F81IG	60	-1653.1324	3443.9888	0.0000	1.0000
JC69	55	-1671.7586	3468.2539	0.0000	1.0000
F81	58	-1669.1626	3470.8167	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1620.9259	3361.8517	0.2663	0.2663
K2PG	57	-1624.2629	3362.5257	0.1901	0.4564
HKYIG	61	-1620.8442	3363.6884	0.1063	0.5627
HKYI	60	-1622.1520	3364.3041	0.0781	0.6408
K2PIG	58	-1624.1963	3364.3926	0.0748	0.7156
GTRG	64	-1618.3078	3364.6157	0.0669	0.7825
K2PI	57	-1625.3089	3364.6179	0.0668	0.8493
SYMG	61	-1621.4447	3364.8894	0.0583	0.9076
GTRIG	65	-1618.2368	3366.4736	0.0264	0.9340
SYMIG	62	-1621.3912	3366.7823	0.0226	0.9566
SYMI	61	-1622.3932	3366.7863	0.0226	0.9792
GTRI	64	-1619.4748	3366.9496	0.0208	1.0000
HKY	59	-1637.7154	3393.4309	0.0000	1.0000
K2P	56	-1640.8816	3393.7632	0.0000	1.0000
SYM	60	-1638.1336	3396.2672	0.0000	1.0000
GTR	63	-1635.1821	3396.3642	0.0000	1.0000
JC69G	56	-1655.8324	3423.6647	0.0000	1.0000
F81G	59	-1653.1978	3424.3956	0.0000	1.0000
JC69IG	57	-1655.7727	3425.5454	0.0000	1.0000
JC69I	56	-1656.8269	3425.6539	0.0000	1.0000
F81IG	60	-1653.1324	3426.2648	0.0000	1.0000
F81I	59	-1654.2636	3426.5271	0.0000	1.0000
JC69	55	-1671.7586	3453.5171	0.0000	1.0000
F81	58	-1669.1626	3454.3252	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1624.2629	3599.7146	0.7122	0.7122
K2PI	57	-1625.3089	3601.8067	0.2502	0.9624
K2PIG	58	-1624.1963	3605.7427	0.0350	0.9973
HKYG	60	-1620.9259	3611.5241	0.0019	0.9993
HKYI	60	-1622.1520	3613.9765	0.0006	0.9998
HKYIG	61	-1620.8442	3617.5221	0.0001	0.9999
SYMG	61	-1621.4447	3618.7231	0.0001	1.0000
SYMI	61	-1622.3932	3620.6200	0.0000	1.0000
SYMIG	62	-1621.3912	3624.7772	0.0000	1.0000
K2P	56	-1640.8816	3626.7908	0.0000	1.0000
GTRG	64	-1618.3078	3630.9330	0.0000	1.0000
GTRI	64	-1619.4748	3633.2669	0.0000	1.0000
GTRIG	65	-1618.2368	3636.9520	0.0000	1.0000
HKY	59	-1637.7154	3638.9421	0.0000	1.0000
SYM	60	-1638.1336	3645.9397	0.0000	1.0000
JC69G	56	-1655.8324	3656.6924	0.0000	1.0000
GTR	63	-1635.1821	3658.5202	0.0000	1.0000
JC69I	56	-1656.8269	3658.6815	0.0000	1.0000
JC69IG	57	-1655.7727	3662.7342	0.0000	1.0000
F81G	59	-1653.1978	3669.9069	0.0000	1.0000
F81I	59	-1654.2636	3672.0384	0.0000	1.0000
F81IG	60	-1653.1324	3675.9373	0.0000	1.0000
JC69	55	-1671.7586	3682.3835	0.0000	1.0000
F81	58	-1669.1626	3695.6752	0.0000	1.0000



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

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

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


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1813_hg18.phylip = 1 - 474;
 Partition Dummy = 1:261-aggressive-phylip/1813_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1813_hg18.phylip = 1 - 474;
 Partition Dummy = 1:261-aggressive-phylip/1813_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

