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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1761.2882	3564.0763	0.3974	0.3974
GTRI	20	-1761.6211	3564.7422	0.2849	0.6822
GTRIG	21	-1761.2938	3566.2405	0.1347	0.8169
HKYG	16	-1767.1716	3567.3077	0.0790	0.8959
HKYI	16	-1767.4006	3567.7657	0.0628	0.9587
HKYIG	17	-1767.1735	3569.4341	0.0273	0.9860
GTR	19	-1766.1297	3571.6140	0.0092	0.9951
HKY	15	-1771.0187	3572.8870	0.0049	1.0000
K2PG	13	-1792.6330	3611.9079	0.0000	1.0000
K2PI	13	-1792.7449	3612.1319	0.0000	1.0000
K2PIG	14	-1792.6157	3613.9734	0.0000	1.0000
F81G	15	-1792.2004	3615.2503	0.0000	1.0000
F81I	15	-1792.3396	3615.5287	0.0000	1.0000
SYMG	17	-1790.3139	3615.7149	0.0000	1.0000
SYMI	17	-1790.4792	3616.0455	0.0000	1.0000
K2P	12	-1796.2947	3617.1386	0.0000	1.0000
F81IG	16	-1792.1930	3617.3505	0.0000	1.0000
SYMIG	18	-1790.3025	3617.8221	0.0000	1.0000
F81	14	-1795.9653	3620.6726	0.0000	1.0000
SYM	16	-1794.1255	3621.2155	0.0000	1.0000
JC69G	12	-1813.8716	3652.2924	0.0000	1.0000
JC69I	12	-1813.9535	3652.4563	0.0000	1.0000
JC69IG	13	-1813.8532	3654.3484	0.0000	1.0000
JC69	11	-1817.5613	3657.5866	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1761.2882	3562.5763	0.4095	0.4095
GTRI	20	-1761.6211	3563.2422	0.2936	0.7031
GTRIG	21	-1761.2938	3564.5875	0.1498	0.8529
HKYG	16	-1767.1716	3566.3431	0.0623	0.9152
HKYI	16	-1767.4006	3566.8011	0.0495	0.9647
HKYIG	17	-1767.1735	3568.3471	0.0229	0.9876
GTR	19	-1766.1297	3570.2593	0.0088	0.9964
HKY	15	-1771.0187	3572.0374	0.0036	1.0000
K2PG	13	-1792.6330	3611.2660	0.0000	1.0000
K2PI	13	-1792.7449	3611.4899	0.0000	1.0000
K2PIG	14	-1792.6157	3613.2313	0.0000	1.0000
F81G	15	-1792.2004	3614.4008	0.0000	1.0000
SYMG	17	-1790.3139	3614.6279	0.0000	1.0000
F81I	15	-1792.3396	3614.6792	0.0000	1.0000
SYMI	17	-1790.4792	3614.9584	0.0000	1.0000
F81IG	16	-1792.1930	3616.3860	0.0000	1.0000
K2P	12	-1796.2947	3616.5893	0.0000	1.0000
SYMIG	18	-1790.3025	3616.6050	0.0000	1.0000
F81	14	-1795.9653	3619.9306	0.0000	1.0000
SYM	16	-1794.1255	3620.2509	0.0000	1.0000
JC69G	12	-1813.8716	3651.7431	0.0000	1.0000
JC69I	12	-1813.9535	3651.9070	0.0000	1.0000
JC69IG	13	-1813.8532	3653.7064	0.0000	1.0000
JC69	11	-1817.5613	3657.1226	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1767.1716	3636.1791	0.4250	0.4250
HKYI	16	-1767.4006	3636.6371	0.3380	0.7629
HKY	15	-1771.0187	3637.5087	0.2186	0.9815
HKYIG	17	-1767.1735	3642.5478	0.0176	0.9991
GTRG	20	-1761.2882	3649.8713	0.0005	0.9996
GTRI	20	-1761.6211	3650.5372	0.0003	0.9999
GTR	19	-1766.1297	3653.1896	0.0001	1.0000
GTRIG	21	-1761.2938	3656.2473	0.0000	1.0000
K2PG	13	-1792.6330	3668.0077	0.0000	1.0000
K2PI	13	-1792.7449	3668.2316	0.0000	1.0000
K2P	12	-1796.2947	3668.9663	0.0000	1.0000
K2PIG	14	-1792.6157	3674.3378	0.0000	1.0000
F81G	15	-1792.2004	3679.8720	0.0000	1.0000
F81I	15	-1792.3396	3680.1504	0.0000	1.0000
F81	14	-1795.9653	3681.0371	0.0000	1.0000
F81IG	16	-1792.1930	3686.2220	0.0000	1.0000
SYMG	17	-1790.3139	3688.8286	0.0000	1.0000
SYMI	17	-1790.4792	3689.1592	0.0000	1.0000
SYM	16	-1794.1255	3690.0870	0.0000	1.0000
SYMIG	18	-1790.3025	3695.1706	0.0000	1.0000
JC69G	12	-1813.8716	3704.1201	0.0000	1.0000
JC69I	12	-1813.9535	3704.2840	0.0000	1.0000
JC69	11	-1817.5613	3705.1348	0.0000	1.0000
JC69IG	13	-1813.8532	3710.4482	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32248.phylip = 1 - 581;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32248.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32248.phylip = 1 - 581;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32248.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

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

