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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1644.3878	3321.6704	0.4683	0.4683
HKYI	16	-1644.9237	3322.7422	0.2740	0.7423
HKYIG	17	-1644.2734	3323.5550	0.1825	0.9248
GTRG	20	-1642.6097	3326.6101	0.0396	0.9644
GTRI	20	-1643.2930	3327.9767	0.0200	0.9844
GTRIG	21	-1642.4724	3328.4771	0.0156	1.0000
HKY	15	-1656.9898	3344.7679	0.0000	1.0000
GTR	19	-1654.4548	3348.1658	0.0000	1.0000
F81G	15	-1678.5563	3387.9008	0.0000	1.0000
F81I	15	-1679.2551	3389.2984	0.0000	1.0000
F81IG	16	-1678.4751	3389.8448	0.0000	1.0000
F81	14	-1689.6183	3407.9251	0.0000	1.0000
SYMG	17	-1686.6431	3408.2945	0.0000	1.0000
SYMIG	18	-1686.5158	3410.1603	0.0000	1.0000
SYMI	17	-1687.6506	3410.3095	0.0000	1.0000
K2PG	13	-1693.3434	3413.2826	0.0000	1.0000
K2PI	13	-1693.9770	3414.5498	0.0000	1.0000
K2PIG	14	-1693.2485	3415.1855	0.0000	1.0000
SYM	16	-1697.5698	3428.0343	0.0000	1.0000
K2P	12	-1704.7937	3434.0973	0.0000	1.0000
JC69G	12	-1721.1708	3466.8513	0.0000	1.0000
JC69I	12	-1721.9088	3468.3274	0.0000	1.0000
JC69IG	13	-1721.0911	3468.7779	0.0000	1.0000
JC69	11	-1731.7879	3486.0064	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1644.3878	3320.7757	0.4532	0.4532
HKYI	16	-1644.9237	3321.8475	0.2652	0.7184
HKYIG	17	-1644.2734	3322.5468	0.1869	0.9053
GTRG	20	-1642.6097	3325.2193	0.0491	0.9545
GTRI	20	-1643.2930	3326.5860	0.0248	0.9793
GTRIG	21	-1642.4724	3326.9447	0.0207	1.0000
HKY	15	-1656.9898	3343.9797	0.0000	1.0000
GTR	19	-1654.4548	3346.9096	0.0000	1.0000
F81G	15	-1678.5563	3387.1126	0.0000	1.0000
F81I	15	-1679.2551	3388.5102	0.0000	1.0000
F81IG	16	-1678.4751	3388.9501	0.0000	1.0000
F81	14	-1689.6183	3407.2366	0.0000	1.0000
SYMG	17	-1686.6431	3407.2863	0.0000	1.0000
SYMIG	18	-1686.5158	3409.0315	0.0000	1.0000
SYMI	17	-1687.6506	3409.3013	0.0000	1.0000
K2PG	13	-1693.3434	3412.6869	0.0000	1.0000
K2PI	13	-1693.9770	3413.9540	0.0000	1.0000
K2PIG	14	-1693.2485	3414.4970	0.0000	1.0000
SYM	16	-1697.5698	3427.1396	0.0000	1.0000
K2P	12	-1704.7937	3433.5875	0.0000	1.0000
JC69G	12	-1721.1708	3466.3415	0.0000	1.0000
JC69I	12	-1721.9088	3467.8176	0.0000	1.0000
JC69IG	13	-1721.0911	3468.1822	0.0000	1.0000
JC69	11	-1731.7879	3485.5758	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1644.3878	3391.7797	0.6135	0.6135
HKYI	16	-1644.9237	3392.8515	0.3590	0.9724
HKYIG	17	-1644.2734	3397.9885	0.0275	0.9999
HKY	15	-1656.9898	3410.5460	0.0001	1.0000
GTRG	20	-1642.6097	3413.9744	0.0000	1.0000
GTRI	20	-1643.2930	3415.3410	0.0000	1.0000
GTRIG	21	-1642.4724	3420.1375	0.0000	1.0000
GTR	19	-1654.4548	3431.2269	0.0000	1.0000
F81G	15	-1678.5563	3453.6789	0.0000	1.0000
F81I	15	-1679.2551	3455.0765	0.0000	1.0000
F81IG	16	-1678.4751	3459.9541	0.0000	1.0000
F81	14	-1689.6183	3469.3651	0.0000	1.0000
K2PG	13	-1693.3434	3470.3776	0.0000	1.0000
K2PI	13	-1693.9770	3471.6448	0.0000	1.0000
K2PIG	14	-1693.2485	3476.6255	0.0000	1.0000
SYMG	17	-1686.6431	3482.7281	0.0000	1.0000
SYMI	17	-1687.6506	3484.7430	0.0000	1.0000
K2P	12	-1704.7937	3486.8405	0.0000	1.0000
SYMIG	18	-1686.5158	3488.9111	0.0000	1.0000
SYM	16	-1697.5698	3498.1436	0.0000	1.0000
JC69G	12	-1721.1708	3519.5945	0.0000	1.0000
JC69I	12	-1721.9088	3521.0706	0.0000	1.0000
JC69IG	13	-1721.0911	3525.8730	0.0000	1.0000
JC69	11	-1731.7879	3534.3910	0.0000	1.0000



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

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

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


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

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


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

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

