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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1588.8403	3218.8926	0.3663	0.3663
GTRG	20	-1588.9074	3219.0268	0.3425	0.7089
GTRIG	21	-1588.9204	3221.1760	0.1170	0.8258
HKYI	16	-1594.6294	3222.0393	0.0760	0.9018
HKYG	16	-1594.6714	3222.1233	0.0728	0.9746
HKYIG	17	-1594.6773	3224.2339	0.0254	1.0000
GTR	19	-1599.1498	3237.3947	0.0000	1.0000
HKY	15	-1604.7948	3240.2772	0.0000	1.0000
F81G	15	-1615.4513	3261.5903	0.0000	1.0000
F81I	15	-1615.4796	3261.6469	0.0000	1.0000
K2PG	13	-1618.0498	3262.6195	0.0000	1.0000
K2PI	13	-1618.0650	3262.6501	0.0000	1.0000
F81IG	16	-1615.4566	3263.6936	0.0000	1.0000
K2PIG	14	-1618.0609	3264.7226	0.0000	1.0000
SYMG	17	-1615.7152	3266.3097	0.0000	1.0000
SYMI	17	-1615.7354	3266.3500	0.0000	1.0000
SYMIG	18	-1615.7205	3268.4252	0.0000	1.0000
F81	14	-1625.5911	3279.7831	0.0000	1.0000
K2P	12	-1628.4966	3281.4383	0.0000	1.0000
SYM	16	-1625.5694	3283.9193	0.0000	1.0000
JC69G	12	-1636.7303	3297.9057	0.0000	1.0000
JC69I	12	-1636.8007	3298.0465	0.0000	1.0000
JC69IG	13	-1636.6473	3299.8146	0.0000	1.0000
JC69	11	-1647.1536	3316.6833	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1588.8403	3217.6805	0.3758	0.3758
GTRG	20	-1588.9074	3217.8147	0.3514	0.7273
GTRIG	21	-1588.9204	3219.8407	0.1276	0.8549
HKYI	16	-1594.6294	3221.2588	0.0628	0.9177
HKYG	16	-1594.6714	3221.3428	0.0602	0.9779
HKYIG	17	-1594.6773	3223.3546	0.0220	1.0000
GTR	19	-1599.1498	3236.2996	0.0000	1.0000
HKY	15	-1604.7948	3239.5895	0.0000	1.0000
F81G	15	-1615.4513	3260.9026	0.0000	1.0000
F81I	15	-1615.4796	3260.9593	0.0000	1.0000
K2PG	13	-1618.0498	3262.0995	0.0000	1.0000
K2PI	13	-1618.0650	3262.1301	0.0000	1.0000
F81IG	16	-1615.4566	3262.9131	0.0000	1.0000
K2PIG	14	-1618.0609	3264.1217	0.0000	1.0000
SYMG	17	-1615.7152	3265.4303	0.0000	1.0000
SYMI	17	-1615.7354	3265.4707	0.0000	1.0000
SYMIG	18	-1615.7205	3267.4410	0.0000	1.0000
F81	14	-1625.5911	3279.1822	0.0000	1.0000
K2P	12	-1628.4966	3280.9933	0.0000	1.0000
SYM	16	-1625.5694	3283.1388	0.0000	1.0000
JC69G	12	-1636.7303	3297.4607	0.0000	1.0000
JC69I	12	-1636.8007	3297.6015	0.0000	1.0000
JC69IG	13	-1636.6473	3299.2946	0.0000	1.0000
JC69	11	-1647.1536	3316.3072	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1594.6294	3294.3930	0.5008	0.5008
HKYG	16	-1594.6714	3294.4769	0.4802	0.9810
HKYIG	17	-1594.6773	3301.0596	0.0179	0.9989
HKY	15	-1604.7948	3308.1528	0.0005	0.9994
GTRI	20	-1588.8403	3309.0982	0.0003	0.9997
GTRG	20	-1588.9074	3309.2324	0.0003	1.0000
GTRIG	21	-1588.9204	3315.8292	0.0000	1.0000
K2PG	13	-1618.0498	3321.5210	0.0000	1.0000
K2PI	13	-1618.0650	3321.5516	0.0000	1.0000
GTR	19	-1599.1498	3323.1464	0.0000	1.0000
K2PIG	14	-1618.0609	3328.1141	0.0000	1.0000
F81G	15	-1615.4513	3329.4658	0.0000	1.0000
F81I	15	-1615.4796	3329.5225	0.0000	1.0000
K2P	12	-1628.4966	3335.8439	0.0000	1.0000
F81IG	16	-1615.4566	3336.0472	0.0000	1.0000
SYMG	17	-1615.7152	3343.1354	0.0000	1.0000
F81	14	-1625.5911	3343.1746	0.0000	1.0000
SYMI	17	-1615.7354	3343.1757	0.0000	1.0000
SYMIG	18	-1615.7205	3349.7169	0.0000	1.0000
JC69G	12	-1636.7303	3352.3113	0.0000	1.0000
JC69I	12	-1636.8007	3352.4521	0.0000	1.0000
SYM	16	-1625.5694	3356.2730	0.0000	1.0000
JC69IG	13	-1636.6473	3358.7160	0.0000	1.0000
JC69	11	-1647.1536	3366.5869	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_7577.phylip = 1 - 714;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7577.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_7577.phylip = 1 - 714;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7577.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

