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

Input data from file "591-tomodel-phylip/chr2_11804.phylip" (7 taxa, 770 characters)

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


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	-2264.3255	4569.7724	0.3074	0.3074
HKYG	16	-2268.6617	4570.0459	0.2681	0.5756
HKYI	16	-2269.5226	4571.7677	0.1134	0.6889
GTRI	20	-2265.3690	4571.8594	0.1083	0.7972
GTRIG	21	-2264.3146	4571.8645	0.1080	0.9052
HKYIG	17	-2268.6562	4572.1263	0.0948	1.0000
HKY	15	-2279.8471	4590.3308	0.0000	1.0000
GTR	19	-2276.0891	4591.1915	0.0000	1.0000
K2PG	13	-2283.9410	4594.3635	0.0000	1.0000
K2PI	13	-2284.8281	4596.1377	0.0000	1.0000
K2PIG	14	-2283.9356	4596.4274	0.0000	1.0000
SYMG	17	-2282.4431	4599.7001	0.0000	1.0000
SYMI	17	-2283.3521	4601.5181	0.0000	1.0000
SYMIG	18	-2282.4367	4601.7842	0.0000	1.0000
K2P	12	-2294.8625	4614.1372	0.0000	1.0000
SYM	16	-2293.5285	4619.7795	0.0000	1.0000
F81G	15	-2299.0403	4628.7172	0.0000	1.0000
F81I	15	-2299.9078	4630.4522	0.0000	1.0000
F81IG	16	-2299.0390	4630.8005	0.0000	1.0000
F81	14	-2308.7920	4646.1403	0.0000	1.0000
JC69G	12	-2313.0396	4650.4914	0.0000	1.0000
JC69I	12	-2313.8982	4652.2086	0.0000	1.0000
JC69IG	13	-2313.0394	4652.5603	0.0000	1.0000
JC69	11	-2322.7532	4667.8547	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2264.3255	4568.6509	0.3328	0.3328
HKYG	16	-2268.6617	4569.3235	0.2378	0.5705
GTRIG	21	-2264.3146	4570.6292	0.1238	0.6943
GTRI	20	-2265.3690	4570.7379	0.1172	0.8115
HKYI	16	-2269.5226	4571.0452	0.1005	0.9120
HKYIG	17	-2268.6562	4571.3124	0.0879	1.0000
HKY	15	-2279.8471	4589.6942	0.0000	1.0000
GTR	19	-2276.0891	4590.1781	0.0000	1.0000
K2PG	13	-2283.9410	4593.8820	0.0000	1.0000
K2PI	13	-2284.8281	4595.6562	0.0000	1.0000
K2PIG	14	-2283.9356	4595.8711	0.0000	1.0000
SYMG	17	-2282.4431	4598.8862	0.0000	1.0000
SYMI	17	-2283.3521	4600.7043	0.0000	1.0000
SYMIG	18	-2282.4367	4600.8734	0.0000	1.0000
K2P	12	-2294.8625	4613.7250	0.0000	1.0000
SYM	16	-2293.5285	4619.0571	0.0000	1.0000
F81G	15	-2299.0403	4628.0806	0.0000	1.0000
F81I	15	-2299.9078	4629.8156	0.0000	1.0000
F81IG	16	-2299.0390	4630.0781	0.0000	1.0000
F81	14	-2308.7920	4645.5840	0.0000	1.0000
JC69G	12	-2313.0396	4650.0793	0.0000	1.0000
JC69I	12	-2313.8982	4651.7965	0.0000	1.0000
JC69IG	13	-2313.0394	4652.0789	0.0000	1.0000
JC69	11	-2322.7532	4667.5064	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2268.6617	4643.6657	0.6818	0.6818
HKYI	16	-2269.5226	4645.3875	0.2882	0.9700
HKYIG	17	-2268.6562	4650.3011	0.0247	0.9947
K2PG	13	-2283.9410	4654.2851	0.0034	0.9981
K2PI	13	-2284.8281	4656.0593	0.0014	0.9995
HKY	15	-2279.8471	4659.3901	0.0003	0.9998
K2PIG	14	-2283.9356	4660.9206	0.0001	0.9999
GTRG	20	-2264.3255	4661.5787	0.0001	1.0000
GTRI	20	-2265.3690	4663.6658	0.0000	1.0000
GTRIG	21	-2264.3146	4668.2034	0.0000	1.0000
K2P	12	-2294.8625	4669.4817	0.0000	1.0000
SYMG	17	-2282.4431	4677.8749	0.0000	1.0000
GTR	19	-2276.0891	4678.4595	0.0000	1.0000
SYMI	17	-2283.3521	4679.6929	0.0000	1.0000
SYMIG	18	-2282.4367	4684.5084	0.0000	1.0000
SYM	16	-2293.5285	4693.3993	0.0000	1.0000
F81G	15	-2299.0403	4697.7764	0.0000	1.0000
F81I	15	-2299.9078	4699.5115	0.0000	1.0000
F81IG	16	-2299.0390	4704.4203	0.0000	1.0000
JC69G	12	-2313.0396	4705.8359	0.0000	1.0000
JC69I	12	-2313.8982	4707.5531	0.0000	1.0000
F81	14	-2308.7920	4710.6335	0.0000	1.0000
JC69IG	13	-2313.0394	4712.4819	0.0000	1.0000
JC69	11	-2322.7532	4718.6167	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 591-tomodel-phylip/chr2_11804.phylip = 1 - 770;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11804.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 591-tomodel-phylip/chr2_11804.phylip = 1 - 770;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11804.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_11804.phylip = 1 - 770;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11804.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;


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

