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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1622.7379	3391.5236	0.5265	0.5265
GTRIG	65	-1622.7370	3394.1261	0.1433	0.6698
SYMG	61	-1627.9998	3394.3013	0.1313	0.8011
GTRI	64	-1624.1792	3394.4062	0.1246	0.9257
SYMIG	62	-1628.0056	3396.8838	0.0361	0.9617
SYMI	61	-1629.3944	3397.0906	0.0325	0.9943
HKYG	60	-1632.8849	3401.5118	0.0036	0.9979
HKYI	60	-1634.1477	3404.0373	0.0010	0.9989
HKYIG	61	-1632.8876	3404.0769	0.0010	0.9999
K2PG	57	-1640.4269	3408.9820	0.0001	1.0000
K2PIG	58	-1640.4286	3411.5125	0.0000	1.0000
K2PI	57	-1641.7142	3411.5566	0.0000	1.0000
GTR	63	-1636.9010	3417.2566	0.0000	1.0000
SYM	60	-1641.4111	3418.5642	0.0000	1.0000
HKY	59	-1651.7674	3436.7280	0.0000	1.0000
K2P	56	-1659.7581	3445.1281	0.0000	1.0000
F81G	59	-1659.0860	3451.3652	0.0000	1.0000
F81IG	60	-1659.0881	3453.9182	0.0000	1.0000
F81I	59	-1660.4098	3454.0127	0.0000	1.0000
JC69G	56	-1667.5761	3460.7642	0.0000	1.0000
JC69IG	57	-1667.5761	3463.2803	0.0000	1.0000
JC69I	56	-1668.8522	3463.3163	0.0000	1.0000
F81	58	-1679.1508	3488.9568	0.0000	1.0000
JC69	55	-1686.3511	3495.8086	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1622.7379	3373.4759	0.5629	0.5629
GTRIG	65	-1622.7370	3375.4739	0.2073	0.7702
GTRI	64	-1624.1792	3376.3585	0.1332	0.9034
SYMG	61	-1627.9998	3377.9996	0.0586	0.9620
SYMIG	62	-1628.0056	3380.0112	0.0214	0.9835
SYMI	61	-1629.3944	3380.7889	0.0145	0.9980
HKYG	60	-1632.8849	3385.7698	0.0012	0.9992
HKYIG	61	-1632.8876	3387.7752	0.0004	0.9996
HKYI	60	-1634.1477	3388.2954	0.0003	1.0000
K2PG	57	-1640.4269	3394.8538	0.0000	1.0000
K2PIG	58	-1640.4286	3396.8573	0.0000	1.0000
K2PI	57	-1641.7142	3397.4284	0.0000	1.0000
GTR	63	-1636.9010	3399.8020	0.0000	1.0000
SYM	60	-1641.4111	3402.8223	0.0000	1.0000
HKY	59	-1651.7674	3421.5348	0.0000	1.0000
K2P	56	-1659.7581	3431.5162	0.0000	1.0000
F81G	59	-1659.0860	3436.1720	0.0000	1.0000
F81IG	60	-1659.0881	3438.1763	0.0000	1.0000
F81I	59	-1660.4098	3438.8196	0.0000	1.0000
JC69G	56	-1667.5761	3447.1523	0.0000	1.0000
JC69IG	57	-1667.5761	3449.1521	0.0000	1.0000
JC69I	56	-1668.8522	3449.7043	0.0000	1.0000
F81	58	-1679.1508	3474.3016	0.0000	1.0000
JC69	55	-1686.3511	3482.7022	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1640.4269	3637.9760	0.3690	0.3690
SYMG	61	-1627.9998	3638.1830	0.3328	0.7018
K2PI	57	-1641.7142	3640.5506	0.1019	0.8037
SYMI	61	-1629.3944	3640.9723	0.0825	0.8862
HKYG	60	-1632.8849	3641.6879	0.0577	0.9439
HKYI	60	-1634.1477	3644.2134	0.0163	0.9602
K2PIG	58	-1640.4286	3644.2448	0.0161	0.9762
SYMIG	62	-1628.0056	3644.4599	0.0144	0.9907
GTRG	64	-1622.7379	3646.4552	0.0053	0.9960
HKYIG	61	-1632.8876	3647.9586	0.0025	0.9985
GTRI	64	-1624.1792	3649.3377	0.0013	0.9998
GTRIG	65	-1622.7370	3652.7185	0.0002	1.0000
SYM	60	-1641.4111	3658.7404	0.0000	1.0000
GTR	63	-1636.9010	3668.5160	0.0000	1.0000
K2P	56	-1659.7581	3670.3731	0.0000	1.0000
HKY	59	-1651.7674	3673.1876	0.0000	1.0000
JC69G	56	-1667.5761	3686.0091	0.0000	1.0000
F81G	59	-1659.0860	3687.8248	0.0000	1.0000
JC69I	56	-1668.8522	3688.5612	0.0000	1.0000
F81I	59	-1660.4098	3690.4724	0.0000	1.0000
JC69IG	57	-1667.5761	3692.2743	0.0000	1.0000
F81IG	60	-1659.0881	3694.0943	0.0000	1.0000
JC69	55	-1686.3511	3717.2938	0.0000	1.0000
F81	58	-1679.1508	3721.6890	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/chr2_18144.phylip = 1 - 526;
 Partition Dummy = 1:183-phylip/chr2_18144.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/chr2_18144.phylip = 1 - 526;
 Partition Dummy = 1:183-phylip/chr2_18144.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 (K2PG)]

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

