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

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

Minimum AIC  model: GTRIG
Minimum AICc model: HKYG
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
HKYG	60	-1570.0014	3284.3216	0.3427	0.3427
GTRIG	65	-1562.7071	3284.4007	0.3294	0.6721
GTRG	64	-1564.3636	3284.7406	0.2779	0.9500
HKYIG	61	-1570.7238	3288.6609	0.0391	0.9892
K2PG	57	-1577.8633	3291.4766	0.0096	0.9988
K2PIG	58	-1579.1238	3296.8351	0.0007	0.9994
SYMG	61	-1575.4129	3298.0392	0.0004	0.9998
HKYI	60	-1577.7327	3299.7844	0.0002	0.9999
GTRI	64	-1573.0852	3302.1838	0.0000	1.0000
SYMIG	62	-1576.8099	3303.7469	0.0000	1.0000
K2PI	57	-1585.7484	3307.2467	0.0000	1.0000
SYMI	61	-1583.4102	3314.0338	0.0000	1.0000
F81G	59	-1592.7147	3326.8731	0.0000	1.0000
JC69G	56	-1598.6787	3330.2884	0.0000	1.0000
F81IG	60	-1594.7468	3333.8126	0.0000	1.0000
JC69IG	57	-1601.0788	3337.9077	0.0000	1.0000
F81I	59	-1601.3377	3344.1190	0.0000	1.0000
JC69I	56	-1607.5157	3347.9625	0.0000	1.0000
HKY	59	-1615.3688	3372.1813	0.0000	1.0000
GTR	63	-1610.2326	3373.5257	0.0000	1.0000
K2P	56	-1623.6245	3380.1802	0.0000	1.0000
SYM	60	-1621.0374	3386.3938	0.0000	1.0000
F81	58	-1637.3668	3413.3210	0.0000	1.0000
JC69	55	-1643.6320	3417.3948	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1562.7071	3255.4143	0.6102	0.6102
GTRG	64	-1564.3636	3256.7271	0.3165	0.9268
HKYG	60	-1570.0014	3260.0027	0.0615	0.9883
HKYIG	61	-1570.7238	3263.4476	0.0110	0.9993
K2PG	57	-1577.8633	3269.7266	0.0005	0.9998
SYMG	61	-1575.4129	3272.8258	0.0001	0.9999
GTRI	64	-1573.0852	3274.1703	0.0001	0.9999
K2PIG	58	-1579.1238	3274.2477	0.0000	1.0000
HKYI	60	-1577.7327	3275.4654	0.0000	1.0000
SYMIG	62	-1576.8099	3277.6198	0.0000	1.0000
K2PI	57	-1585.7484	3285.4967	0.0000	1.0000
SYMI	61	-1583.4102	3288.8204	0.0000	1.0000
F81G	59	-1592.7147	3303.4294	0.0000	1.0000
JC69G	56	-1598.6787	3309.3573	0.0000	1.0000
F81IG	60	-1594.7468	3309.4937	0.0000	1.0000
JC69IG	57	-1601.0788	3316.1577	0.0000	1.0000
F81I	59	-1601.3377	3320.6753	0.0000	1.0000
JC69I	56	-1607.5157	3327.0313	0.0000	1.0000
GTR	63	-1610.2326	3346.4653	0.0000	1.0000
HKY	59	-1615.3688	3348.7376	0.0000	1.0000
K2P	56	-1623.6245	3359.2491	0.0000	1.0000
SYM	60	-1621.0374	3362.0748	0.0000	1.0000
F81	58	-1637.3668	3390.7335	0.0000	1.0000
JC69	55	-1643.6320	3397.2640	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1577.8633	3491.5503	0.7127	0.7127
HKYG	60	-1570.0014	3493.5014	0.2687	0.9813
K2PIG	58	-1579.1238	3499.9630	0.0106	0.9920
HKYIG	61	-1570.7238	3500.8379	0.0069	0.9988
GTRG	64	-1564.3636	3505.7923	0.0006	0.9994
K2PI	57	-1585.7484	3507.3205	0.0003	0.9997
GTRIG	65	-1562.7071	3508.3711	0.0002	0.9998
HKYI	60	-1577.7327	3508.9641	0.0001	0.9999
SYMG	61	-1575.4129	3510.2161	0.0001	1.0000
SYMIG	62	-1576.8099	3518.9017	0.0000	1.0000
GTRI	64	-1573.0852	3523.2355	0.0000	1.0000
SYMI	61	-1583.4102	3526.2107	0.0000	1.0000
JC69G	56	-1598.6787	3527.2894	0.0000	1.0000
F81G	59	-1592.7147	3533.0364	0.0000	1.0000
JC69IG	57	-1601.0788	3537.9814	0.0000	1.0000
F81IG	60	-1594.7468	3542.9923	0.0000	1.0000
JC69I	56	-1607.5157	3544.9634	0.0000	1.0000
F81I	59	-1601.3377	3550.2823	0.0000	1.0000
K2P	56	-1623.6245	3577.1812	0.0000	1.0000
HKY	59	-1615.3688	3578.3446	0.0000	1.0000
GTR	63	-1610.2326	3591.6389	0.0000	1.0000
SYM	60	-1621.0374	3595.5735	0.0000	1.0000
JC69	55	-1643.6320	3611.3045	0.0000	1.0000
F81	58	-1637.3668	3616.4489	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 183-phylip/chr2_11740.phylip = 1 - 362;
 Partition Dummy = 1:183-phylip/chr2_11740.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_11740.phylip = 1 - 362;
 Partition Dummy = 1:183-phylip/chr2_11740.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 BIC model (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_11740.phylip = 1 - 362;
 Partition Dummy = 1:183-phylip/chr2_11740.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

