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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1581.0983	3322.4447	0.7599	0.7599
GTRIG	65	-1580.6853	3324.7558	0.2393	0.9992
HKYG	60	-1594.6942	3337.3274	0.0004	0.9997
GTRI	64	-1589.2763	3338.8007	0.0002	0.9999
HKYIG	61	-1594.5789	3340.1387	0.0001	1.0000
HKYI	60	-1601.8618	3351.6626	0.0000	1.0000
SYMG	61	-1611.9710	3374.9229	0.0000	1.0000
SYMIG	62	-1611.9397	3377.9256	0.0000	1.0000
SYMI	61	-1623.5663	3398.1133	0.0000	1.0000
K2PG	57	-1632.9966	3404.9441	0.0000	1.0000
K2PIG	58	-1634.1296	3410.1834	0.0000	1.0000
GTR	63	-1628.2833	3413.7018	0.0000	1.0000
HKY	59	-1635.0333	3414.9868	0.0000	1.0000
K2PI	57	-1643.4922	3425.9354	0.0000	1.0000
SYM	60	-1658.7364	3465.4117	0.0000	1.0000
F81G	59	-1662.1540	3469.2282	0.0000	1.0000
F81IG	60	-1661.8615	3471.6618	0.0000	1.0000
F81I	59	-1670.0611	3485.0423	0.0000	1.0000
K2P	56	-1678.6037	3493.2074	0.0000	1.0000
F81	58	-1708.7218	3559.3678	0.0000	1.0000
JC69G	56	-1719.3057	3574.6115	0.0000	1.0000
JC69IG	57	-1722.3290	3583.6088	0.0000	1.0000
JC69I	56	-1734.7005	3605.4009	0.0000	1.0000
JC69	55	-1761.5782	3656.2275	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1581.0983	3290.1967	0.6425	0.6425
GTRIG	65	-1580.6853	3291.3706	0.3572	0.9998
GTRI	64	-1589.2763	3306.5526	0.0002	0.9999
HKYG	60	-1594.6942	3309.3884	0.0000	1.0000
HKYIG	61	-1594.5789	3311.1578	0.0000	1.0000
HKYI	60	-1601.8618	3323.7236	0.0000	1.0000
SYMG	61	-1611.9710	3345.9420	0.0000	1.0000
SYMIG	62	-1611.9397	3347.8794	0.0000	1.0000
SYMI	61	-1623.5663	3369.1325	0.0000	1.0000
K2PG	57	-1632.9966	3379.9932	0.0000	1.0000
GTR	63	-1628.2833	3382.5667	0.0000	1.0000
K2PIG	58	-1634.1296	3384.2591	0.0000	1.0000
HKY	59	-1635.0333	3388.0666	0.0000	1.0000
K2PI	57	-1643.4922	3400.9844	0.0000	1.0000
SYM	60	-1658.7364	3437.4728	0.0000	1.0000
F81G	59	-1662.1540	3442.3080	0.0000	1.0000
F81IG	60	-1661.8615	3443.7229	0.0000	1.0000
F81I	59	-1670.0611	3458.1221	0.0000	1.0000
K2P	56	-1678.6037	3469.2074	0.0000	1.0000
F81	58	-1708.7218	3533.4436	0.0000	1.0000
JC69G	56	-1719.3057	3550.6115	0.0000	1.0000
JC69IG	57	-1722.3290	3558.6579	0.0000	1.0000
JC69I	56	-1734.7005	3581.4009	0.0000	1.0000
JC69	55	-1761.5782	3633.1564	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1581.0983	3531.9664	0.8180	0.8180
HKYG	60	-1594.6942	3536.0476	0.1063	0.9243
GTRIG	65	-1580.6853	3536.9180	0.0688	0.9931
HKYIG	61	-1594.5789	3541.5946	0.0066	0.9997
GTRI	64	-1589.2763	3548.3223	0.0002	0.9999
HKYI	60	-1601.8618	3550.3828	0.0001	1.0000
SYMG	61	-1611.9710	3576.3788	0.0000	1.0000
SYMIG	62	-1611.9397	3582.0939	0.0000	1.0000
K2PG	57	-1632.9966	3595.3194	0.0000	1.0000
SYMI	61	-1623.5663	3599.5693	0.0000	1.0000
K2PIG	58	-1634.1296	3603.3630	0.0000	1.0000
HKY	59	-1635.0333	3610.9481	0.0000	1.0000
K2PI	57	-1643.4922	3616.3106	0.0000	1.0000
GTR	63	-1628.2833	3620.5588	0.0000	1.0000
SYM	60	-1658.7364	3664.1319	0.0000	1.0000
F81G	59	-1662.1540	3665.1895	0.0000	1.0000
F81IG	60	-1661.8615	3670.3820	0.0000	1.0000
K2P	56	-1678.6037	3680.7559	0.0000	1.0000
F81I	59	-1670.0611	3681.0036	0.0000	1.0000
F81	58	-1708.7218	3752.5474	0.0000	1.0000
JC69G	56	-1719.3057	3762.1600	0.0000	1.0000
JC69IG	57	-1722.3290	3773.9841	0.0000	1.0000
JC69I	56	-1734.7005	3792.9495	0.0000	1.0000
JC69	55	-1761.5782	3840.9273	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/chr6_3147.phylip = 1 - 323;
 Partition Dummy = 1:183-phylip/chr6_3147.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/chr6_3147.phylip = 1 - 323;
 Partition Dummy = 1:183-phylip/chr6_3147.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 (GTRG)]

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


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

