
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:16 04/04/2011

Input data from file "683-tomodel-phylip/chr3_3067.phylip" (7 taxa, 561 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	20	-1644.3526	3330.2607	0.4925	0.4925
GTRI	20	-1645.1085	3331.7726	0.2313	0.7238
GTRIG	21	-1644.3349	3332.3841	0.1703	0.8942
HKYG	16	-1651.4983	3335.9966	0.0280	0.9221
SYMG	17	-1650.4504	3336.0278	0.0276	0.9497
SYMI	17	-1651.4189	3337.9649	0.0105	0.9602
HKYIG	17	-1651.4892	3338.1054	0.0097	0.9699
SYMIG	18	-1650.4463	3338.1546	0.0095	0.9794
HKYI	16	-1652.6171	3338.2343	0.0091	0.9886
K2PG	13	-1656.1161	3338.8977	0.0066	0.9951
K2PI	13	-1657.0493	3340.7640	0.0026	0.9977
K2PIG	14	-1656.1125	3340.9942	0.0023	1.0000
GTR	19	-1657.3728	3354.1504	0.0000	1.0000
SYM	16	-1661.8767	3356.7535	0.0000	1.0000
HKY	15	-1663.6033	3358.0873	0.0000	1.0000
K2P	12	-1667.4712	3359.5117	0.0000	1.0000
F81G	15	-1700.2762	3431.4331	0.0000	1.0000
JC69G	12	-1703.4581	3431.4856	0.0000	1.0000
JC69I	12	-1703.8145	3432.1983	0.0000	1.0000
F81I	15	-1700.6616	3432.2039	0.0000	1.0000
F81IG	16	-1700.2827	3433.5654	0.0000	1.0000
JC69IG	13	-1703.4609	3433.5873	0.0000	1.0000
JC69	11	-1712.9528	3448.3864	0.0000	1.0000
F81	14	-1710.0233	3448.8159	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1644.3526	3328.7051	0.4972	0.4972
GTRI	20	-1645.1085	3330.2170	0.2335	0.7307
GTRIG	21	-1644.3349	3330.6698	0.1862	0.9169
SYMG	17	-1650.4504	3334.9007	0.0224	0.9394
HKYG	16	-1651.4983	3334.9966	0.0214	0.9608
SYMI	17	-1651.4189	3336.8378	0.0085	0.9693
SYMIG	18	-1650.4463	3336.8926	0.0083	0.9776
HKYIG	17	-1651.4892	3336.9783	0.0079	0.9855
HKYI	16	-1652.6171	3337.2343	0.0070	0.9925
K2PG	13	-1656.1161	3338.2322	0.0042	0.9968
K2PI	13	-1657.0493	3340.0985	0.0017	0.9984
K2PIG	14	-1656.1125	3340.2250	0.0016	1.0000
GTR	19	-1657.3728	3352.7456	0.0000	1.0000
SYM	16	-1661.8767	3355.7535	0.0000	1.0000
HKY	15	-1663.6033	3357.2066	0.0000	1.0000
K2P	12	-1667.4712	3358.9424	0.0000	1.0000
F81G	15	-1700.2762	3430.5524	0.0000	1.0000
JC69G	12	-1703.4581	3430.9163	0.0000	1.0000
F81I	15	-1700.6616	3431.3231	0.0000	1.0000
JC69I	12	-1703.8145	3431.6289	0.0000	1.0000
F81IG	16	-1700.2827	3432.5654	0.0000	1.0000
JC69IG	13	-1703.4609	3432.9219	0.0000	1.0000
JC69	11	-1712.9528	3447.9056	0.0000	1.0000
F81	14	-1710.0233	3448.0466	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1656.1161	3394.5186	0.6907	0.6907
K2PI	13	-1657.0493	3396.3849	0.2717	0.9624
K2PIG	14	-1656.1125	3400.8411	0.0293	0.9917
HKYG	16	-1651.4983	3404.2722	0.0053	0.9969
HKYI	16	-1652.6171	3406.5098	0.0017	0.9986
SYMG	17	-1650.4504	3408.5060	0.0006	0.9993
SYMI	17	-1651.4189	3410.4430	0.0002	0.9995
HKYIG	17	-1651.4892	3410.5836	0.0002	0.9997
K2P	12	-1667.4712	3410.8990	0.0002	0.9999
SYMIG	18	-1650.4463	3414.8275	0.0000	1.0000
GTRG	20	-1644.3526	3415.2995	0.0000	1.0000
GTRI	20	-1645.1085	3416.8115	0.0000	1.0000
GTRIG	21	-1644.3349	3421.5940	0.0000	1.0000
HKY	15	-1663.6033	3422.1524	0.0000	1.0000
SYM	16	-1661.8767	3425.0290	0.0000	1.0000
GTR	19	-1657.3728	3435.0103	0.0000	1.0000
JC69G	12	-1703.4581	3482.8729	0.0000	1.0000
JC69I	12	-1703.8145	3483.5856	0.0000	1.0000
JC69IG	13	-1703.4609	3489.2083	0.0000	1.0000
F81G	15	-1700.2762	3495.4982	0.0000	1.0000
JC69	11	-1712.9528	3495.5325	0.0000	1.0000
F81I	15	-1700.6616	3496.2690	0.0000	1.0000
F81IG	16	-1700.2827	3501.8409	0.0000	1.0000
F81	14	-1710.0233	3508.6627	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 683-tomodel-phylip/chr3_3067.phylip = 1 - 561;
 Partition Dummy = 1:683-tomodel-phylip/chr3_3067.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 683-tomodel-phylip/chr3_3067.phylip = 1 - 561;
 Partition Dummy = 1:683-tomodel-phylip/chr3_3067.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 683-tomodel-phylip/chr3_3067.phylip = 1 - 561;
 Partition Dummy = 1:683-tomodel-phylip/chr3_3067.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

