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

Input data from file "591-tomodel-phylip/chr3_17860.phylip" (7 taxa, 586 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	-1574.8731	3182.7023	0.2897	0.2897
HKYG	16	-1575.1270	3183.2100	0.2248	0.5145
GTRI	20	-1571.0823	3183.6513	0.1803	0.6948
GTRG	20	-1571.2517	3183.9900	0.1522	0.8470
HKYIG	17	-1574.8797	3184.8368	0.0997	0.9467
GTRIG	21	-1571.2582	3186.1547	0.0516	0.9982
HKY	15	-1581.8320	3194.5062	0.0008	0.9990
GTR	19	-1577.6028	3194.5483	0.0008	0.9998
K2PI	13	-1586.1066	3198.8496	0.0001	0.9999
K2PG	13	-1586.4122	3199.4608	0.0001	0.9999
K2PIG	14	-1586.1148	3200.9651	0.0000	1.0000
SYMI	17	-1583.9475	3202.9725	0.0000	1.0000
SYMG	17	-1584.2168	3203.5111	0.0000	1.0000
SYMIG	18	-1583.9546	3205.1155	0.0000	1.0000
K2P	12	-1593.0040	3210.5524	0.0000	1.0000
SYM	16	-1590.3380	3213.6321	0.0000	1.0000
F81I	15	-1594.2238	3219.2896	0.0000	1.0000
F81G	15	-1594.5412	3219.9244	0.0000	1.0000
F81IG	16	-1594.2320	3221.4201	0.0000	1.0000
F81	14	-1600.8554	3230.4463	0.0000	1.0000
JC69I	12	-1604.0615	3232.6675	0.0000	1.0000
JC69G	12	-1604.4080	3233.3605	0.0000	1.0000
JC69IG	13	-1604.0707	3234.7778	0.0000	1.0000
JC69	11	-1610.6860	3243.8319	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1574.8731	3181.7463	0.2568	0.2568
GTRI	20	-1571.0823	3182.1646	0.2083	0.4651
HKYG	16	-1575.1270	3182.2540	0.1992	0.6643
GTRG	20	-1571.2517	3182.5033	0.1759	0.8402
HKYIG	17	-1574.8797	3183.7594	0.0939	0.9341
GTRIG	21	-1571.2582	3184.5164	0.0643	0.9983
GTR	19	-1577.6028	3193.2056	0.0008	0.9992
HKY	15	-1581.8320	3193.6641	0.0007	0.9998
K2PI	13	-1586.1066	3198.2132	0.0001	0.9999
K2PG	13	-1586.4122	3198.8245	0.0001	1.0000
K2PIG	14	-1586.1148	3200.2295	0.0000	1.0000
SYMI	17	-1583.9475	3201.8950	0.0000	1.0000
SYMG	17	-1584.2168	3202.4336	0.0000	1.0000
SYMIG	18	-1583.9546	3203.9092	0.0000	1.0000
K2P	12	-1593.0040	3210.0079	0.0000	1.0000
SYM	16	-1590.3380	3212.6761	0.0000	1.0000
F81I	15	-1594.2238	3218.4475	0.0000	1.0000
F81G	15	-1594.5412	3219.0823	0.0000	1.0000
F81IG	16	-1594.2320	3220.4640	0.0000	1.0000
F81	14	-1600.8554	3229.7107	0.0000	1.0000
JC69I	12	-1604.0615	3232.1230	0.0000	1.0000
JC69G	12	-1604.4080	3232.8160	0.0000	1.0000
JC69IG	13	-1604.0707	3234.1414	0.0000	1.0000
JC69	11	-1610.6860	3243.3720	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1574.8731	3251.7194	0.4591	0.4591
HKYG	16	-1575.1270	3252.2271	0.3562	0.8153
K2PI	13	-1586.1066	3255.0664	0.0861	0.9014
K2PG	13	-1586.4122	3255.6776	0.0634	0.9648
HKYIG	17	-1574.8797	3258.1058	0.0188	0.9837
HKY	15	-1581.8320	3259.2639	0.0106	0.9942
K2PIG	14	-1586.1148	3261.4560	0.0035	0.9978
K2P	12	-1593.0040	3262.4878	0.0021	0.9999
GTRI	20	-1571.0823	3269.6310	0.0001	0.9999
GTRG	20	-1571.2517	3269.9697	0.0000	1.0000
SYMI	17	-1583.9475	3276.2414	0.0000	1.0000
GTR	19	-1577.6028	3276.2987	0.0000	1.0000
GTRIG	21	-1571.2582	3276.3561	0.0000	1.0000
SYMG	17	-1584.2168	3276.7801	0.0000	1.0000
SYMIG	18	-1583.9546	3282.6290	0.0000	1.0000
SYM	16	-1590.3380	3282.6492	0.0000	1.0000
F81I	15	-1594.2238	3284.0473	0.0000	1.0000
JC69I	12	-1604.0615	3284.6028	0.0000	1.0000
F81G	15	-1594.5412	3284.6821	0.0000	1.0000
JC69G	12	-1604.4080	3285.2958	0.0000	1.0000
F81IG	16	-1594.2320	3290.4371	0.0000	1.0000
F81	14	-1600.8554	3290.9372	0.0000	1.0000
JC69IG	13	-1604.0707	3290.9946	0.0000	1.0000
JC69	11	-1610.6860	3291.4785	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/chr3_17860.phylip = 1 - 586;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17860.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/chr3_17860.phylip = 1 - 586;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17860.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/chr3_17860.phylip = 1 - 586;
 Partition Dummy = 1:591-tomodel-phylip/chr3_17860.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

