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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-1749.1186	3531.1319	0.3644	0.3644
HKYG	16	-1749.1669	3531.2286	0.3472	0.7117
HKYIG	17	-1749.0324	3533.0730	0.1381	0.8497
GTRI	20	-1746.5765	3534.5436	0.0662	0.9159
GTRG	20	-1746.6210	3534.6327	0.0633	0.9792
GTRIG	21	-1746.6634	3536.8591	0.0208	1.0000
HKY	15	-1770.3934	3571.5749	0.0000	1.0000
GTR	19	-1770.0239	3579.3039	0.0000	1.0000
K2PG	13	-1790.4012	3607.3982	0.0000	1.0000
K2PI	13	-1790.4156	3607.4269	0.0000	1.0000
K2PIG	14	-1790.4331	3609.5547	0.0000	1.0000
SYMG	17	-1790.0595	3615.1272	0.0000	1.0000
SYMI	17	-1790.1070	3615.2222	0.0000	1.0000
SYMIG	18	-1790.0860	3617.3008	0.0000	1.0000
F81I	15	-1804.6336	3640.0553	0.0000	1.0000
F81G	15	-1804.7987	3640.3856	0.0000	1.0000
F81IG	16	-1804.8440	3642.5828	0.0000	1.0000
K2P	12	-1811.4213	3647.3524	0.0000	1.0000
SYM	16	-1811.0354	3654.9655	0.0000	1.0000
F81	14	-1822.4603	3673.6091	0.0000	1.0000
JC69I	12	-1838.3708	3701.2515	0.0000	1.0000
JC69G	12	-1838.4633	3701.4363	0.0000	1.0000
JC69IG	13	-1838.5175	3703.6308	0.0000	1.0000
JC69	11	-1856.3245	3735.0797	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1749.1186	3530.2371	0.3463	0.3463
HKYG	16	-1749.1669	3530.3338	0.3300	0.6763
HKYIG	17	-1749.0324	3532.0648	0.1389	0.8151
GTRI	20	-1746.5765	3533.1529	0.0806	0.8957
GTRG	20	-1746.6210	3533.2420	0.0771	0.9728
GTRIG	21	-1746.6634	3535.3268	0.0272	1.0000
HKY	15	-1770.3934	3570.7867	0.0000	1.0000
GTR	19	-1770.0239	3578.0477	0.0000	1.0000
K2PG	13	-1790.4012	3606.8025	0.0000	1.0000
K2PI	13	-1790.4156	3606.8311	0.0000	1.0000
K2PIG	14	-1790.4331	3608.8662	0.0000	1.0000
SYMG	17	-1790.0595	3614.1190	0.0000	1.0000
SYMI	17	-1790.1070	3614.2140	0.0000	1.0000
SYMIG	18	-1790.0860	3616.1721	0.0000	1.0000
F81I	15	-1804.6336	3639.2671	0.0000	1.0000
F81G	15	-1804.7987	3639.5974	0.0000	1.0000
F81IG	16	-1804.8440	3641.6881	0.0000	1.0000
K2P	12	-1811.4213	3646.8426	0.0000	1.0000
SYM	16	-1811.0354	3654.0708	0.0000	1.0000
F81	14	-1822.4603	3672.9205	0.0000	1.0000
JC69I	12	-1838.3708	3700.7417	0.0000	1.0000
JC69G	12	-1838.4633	3700.9265	0.0000	1.0000
JC69IG	13	-1838.5175	3703.0350	0.0000	1.0000
JC69	11	-1856.3245	3734.6491	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1749.1186	3601.2411	0.5009	0.5009
HKYG	16	-1749.1669	3601.3378	0.4772	0.9781
HKYIG	17	-1749.0324	3607.5065	0.0218	1.0000
GTRI	20	-1746.5765	3621.9080	0.0000	1.0000
GTRG	20	-1746.6210	3621.9970	0.0000	1.0000
GTRIG	21	-1746.6634	3628.5195	0.0000	1.0000
HKY	15	-1770.3934	3637.3530	0.0000	1.0000
GTR	19	-1770.0239	3662.3650	0.0000	1.0000
K2PG	13	-1790.4012	3664.4933	0.0000	1.0000
K2PI	13	-1790.4156	3664.5219	0.0000	1.0000
K2PIG	14	-1790.4331	3670.9947	0.0000	1.0000
SYMG	17	-1790.0595	3689.5607	0.0000	1.0000
SYMI	17	-1790.1070	3689.6558	0.0000	1.0000
SYMIG	18	-1790.0860	3696.0516	0.0000	1.0000
K2P	12	-1811.4213	3700.0956	0.0000	1.0000
F81I	15	-1804.6336	3705.8334	0.0000	1.0000
F81G	15	-1804.7987	3706.1637	0.0000	1.0000
F81IG	16	-1804.8440	3712.6921	0.0000	1.0000
SYM	16	-1811.0354	3725.0748	0.0000	1.0000
F81	14	-1822.4603	3735.0491	0.0000	1.0000
JC69I	12	-1838.3708	3753.9947	0.0000	1.0000
JC69G	12	-1838.4633	3754.1795	0.0000	1.0000
JC69IG	13	-1838.5175	3760.7258	0.0000	1.0000
JC69	11	-1856.3245	3783.4643	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10694.phylip = 1 - 625;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10694.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10694.phylip = 1 - 625;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10694.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10694.phylip = 1 - 625;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10694.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

