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

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

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


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	-1668.3791	3369.8334	0.5422	0.5422
HKYIG	17	-1668.4137	3372.0392	0.1800	0.7222
HKYG	16	-1669.9596	3372.9944	0.1116	0.8338
GTRI	20	-1665.8700	3373.4133	0.0905	0.9244
GTRIG	21	-1665.9050	3375.6543	0.0295	0.9539
GTRG	20	-1667.4925	3376.6583	0.0179	0.9718
K2PI	13	-1674.9955	3376.7061	0.0174	0.9892
K2PIG	14	-1675.0312	3378.8892	0.0059	0.9951
K2PG	13	-1676.5466	3379.8082	0.0037	0.9988
SYMI	17	-1673.8580	3382.9278	0.0008	0.9996
SYMIG	18	-1673.8936	3385.1443	0.0003	0.9998
SYMG	17	-1675.3948	3386.0014	0.0002	1.0000
HKY	15	-1680.2210	3391.3888	0.0000	1.0000
GTR	19	-1677.7098	3394.9306	0.0000	1.0000
K2P	12	-1686.3070	3397.2257	0.0000	1.0000
SYM	16	-1684.9161	3402.9074	0.0000	1.0000
F81I	15	-1694.7447	3420.4360	0.0000	1.0000
F81IG	16	-1694.7775	3422.6300	0.0000	1.0000
F81G	15	-1696.2085	3423.3638	0.0000	1.0000
JC69I	12	-1699.8871	3424.3859	0.0000	1.0000
JC69IG	13	-1699.9229	3426.5609	0.0000	1.0000
JC69G	12	-1701.3941	3427.3999	0.0000	1.0000
F81	14	-1705.5387	3439.9042	0.0000	1.0000
JC69	11	-1710.4238	3443.3643	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1668.3791	3368.7583	0.5114	0.5114
HKYIG	17	-1668.4137	3370.8273	0.1817	0.6931
GTRI	20	-1665.8700	3371.7400	0.1152	0.8083
HKYG	16	-1669.9596	3371.9193	0.1053	0.9136
GTRIG	21	-1665.9050	3373.8100	0.0409	0.9545
GTRG	20	-1667.4925	3374.9850	0.0227	0.9772
K2PI	13	-1674.9955	3375.9910	0.0137	0.9910
K2PIG	14	-1675.0312	3378.0624	0.0049	0.9958
K2PG	13	-1676.5466	3379.0931	0.0029	0.9988
SYMI	17	-1673.8580	3381.7159	0.0008	0.9995
SYMIG	18	-1673.8936	3383.7872	0.0003	0.9998
SYMG	17	-1675.3948	3384.7895	0.0002	1.0000
HKY	15	-1680.2210	3390.4420	0.0000	1.0000
GTR	19	-1677.7098	3393.4196	0.0000	1.0000
K2P	12	-1686.3070	3396.6139	0.0000	1.0000
SYM	16	-1684.9161	3401.8323	0.0000	1.0000
F81I	15	-1694.7447	3419.4893	0.0000	1.0000
F81IG	16	-1694.7775	3421.5549	0.0000	1.0000
F81G	15	-1696.2085	3422.4170	0.0000	1.0000
JC69I	12	-1699.8871	3423.7741	0.0000	1.0000
JC69IG	13	-1699.9229	3425.8458	0.0000	1.0000
JC69G	12	-1701.3941	3426.7881	0.0000	1.0000
F81	14	-1705.5387	3439.0774	0.0000	1.0000
JC69	11	-1710.4238	3442.8476	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1674.9955	3431.3656	0.7505	0.7505
K2PG	13	-1676.5466	3434.4677	0.1591	0.9096
HKYI	16	-1668.3791	3436.9116	0.0469	0.9565
K2PIG	14	-1675.0312	3437.6965	0.0317	0.9881
HKYG	16	-1669.9596	3440.0726	0.0097	0.9978
HKYIG	17	-1668.4137	3443.2402	0.0020	0.9998
K2P	12	-1686.3070	3447.7289	0.0002	1.0000
SYMI	17	-1673.8580	3454.1288	0.0000	1.0000
HKY	15	-1680.2210	3454.3357	0.0000	1.0000
GTRI	20	-1665.8700	3456.9317	0.0000	1.0000
SYMG	17	-1675.3948	3457.2024	0.0000	1.0000
GTRG	20	-1667.4925	3460.1767	0.0000	1.0000
SYMIG	18	-1673.8936	3460.4596	0.0000	1.0000
GTRIG	21	-1665.9050	3463.2612	0.0000	1.0000
SYM	16	-1684.9161	3469.9856	0.0000	1.0000
GTR	19	-1677.7098	3474.3517	0.0000	1.0000
JC69I	12	-1699.8871	3474.8891	0.0000	1.0000
JC69G	12	-1701.3941	3477.9031	0.0000	1.0000
JC69IG	13	-1699.9229	3481.2203	0.0000	1.0000
F81I	15	-1694.7447	3483.3830	0.0000	1.0000
F81G	15	-1696.2085	3486.3107	0.0000	1.0000
JC69	11	-1710.4238	3489.7030	0.0000	1.0000
F81IG	16	-1694.7775	3489.7082	0.0000	1.0000
F81	14	-1705.5387	3498.7115	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/chr2_6336.phylip = 1 - 523;
 Partition Dummy = 1:591-tomodel-phylip/chr2_6336.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/chr2_6336.phylip = 1 - 523;
 Partition Dummy = 1:591-tomodel-phylip/chr2_6336.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 (K2PI)]

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


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

