
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:09 04/04/2011

Input data from file "683-tomodel-phylip/chr12_6056.phylip" (7 taxa, 807 characters)

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: HKY


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
HKY	15	-1871.7737	3774.1543	0.3497	0.3497
HKYG	16	-1871.0722	3774.8329	0.2491	0.5988
HKYI	16	-1871.2192	3775.1269	0.2150	0.8138
HKYIG	17	-1871.0716	3776.9188	0.0878	0.9016
GTRG	20	-1868.8666	3778.8020	0.0342	0.9358
GTRI	20	-1869.0134	3779.0955	0.0296	0.9654
GTR	19	-1870.3271	3779.6199	0.0227	0.9881
GTRIG	21	-1868.8732	3780.9234	0.0119	1.0000
F81	14	-1897.2566	3823.0436	0.0000	1.0000
F81G	15	-1896.7208	3824.0484	0.0000	1.0000
F81I	15	-1896.7414	3824.0897	0.0000	1.0000
F81IG	16	-1896.7245	3826.1376	0.0000	1.0000
K2P	12	-1905.4154	3835.2238	0.0000	1.0000
K2PG	13	-1904.8081	3836.0753	0.0000	1.0000
K2PI	13	-1904.9503	3836.3597	0.0000	1.0000
K2PIG	14	-1904.8073	3838.1449	0.0000	1.0000
SYM	16	-1903.6798	3840.0481	0.0000	1.0000
SYMG	17	-1902.6702	3840.1160	0.0000	1.0000
SYMI	17	-1902.7955	3840.3667	0.0000	1.0000
SYMIG	18	-1902.6843	3842.2365	0.0000	1.0000
JC69	11	-1927.0831	3876.4983	0.0000	1.0000
JC69G	12	-1926.5426	3877.4782	0.0000	1.0000
JC69I	12	-1926.5919	3877.5767	0.0000	1.0000
JC69IG	13	-1926.5542	3879.5675	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-1871.7737	3773.5475	0.3324	0.3324
HKYG	16	-1871.0722	3774.1443	0.2466	0.5790
HKYI	16	-1871.2192	3774.4383	0.2129	0.7920
HKYIG	17	-1871.0716	3776.1432	0.0908	0.8828
GTRG	20	-1868.8666	3777.7333	0.0410	0.9238
GTRI	20	-1869.0134	3778.0268	0.0354	0.9591
GTR	19	-1870.3271	3778.6542	0.0259	0.9850
GTRIG	21	-1868.8732	3779.7463	0.0150	1.0000
F81	14	-1897.2566	3822.5133	0.0000	1.0000
F81G	15	-1896.7208	3823.4416	0.0000	1.0000
F81I	15	-1896.7414	3823.4828	0.0000	1.0000
F81IG	16	-1896.7245	3825.4490	0.0000	1.0000
K2P	12	-1905.4154	3834.8308	0.0000	1.0000
K2PG	13	-1904.8081	3835.6163	0.0000	1.0000
K2PI	13	-1904.9503	3835.9007	0.0000	1.0000
K2PIG	14	-1904.8073	3837.6146	0.0000	1.0000
SYMG	17	-1902.6702	3839.3403	0.0000	1.0000
SYM	16	-1903.6798	3839.3595	0.0000	1.0000
SYMI	17	-1902.7955	3839.5911	0.0000	1.0000
SYMIG	18	-1902.6843	3841.3685	0.0000	1.0000
JC69	11	-1927.0831	3876.1663	0.0000	1.0000
JC69G	12	-1926.5426	3877.0852	0.0000	1.0000
JC69I	12	-1926.5919	3877.1838	0.0000	1.0000
JC69IG	13	-1926.5542	3879.1084	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1871.7737	3843.9473	0.8812	0.8812
HKYG	16	-1871.0722	3849.2375	0.0626	0.9438
HKYI	16	-1871.2192	3849.5315	0.0540	0.9978
HKYIG	17	-1871.0716	3855.9297	0.0022	1.0000
GTR	19	-1870.3271	3867.8273	0.0000	1.0000
GTRG	20	-1868.8666	3871.5998	0.0000	1.0000
GTRI	20	-1869.0134	3871.8933	0.0000	1.0000
GTRIG	21	-1868.8732	3878.3061	0.0000	1.0000
F81	14	-1897.2566	3888.2198	0.0000	1.0000
K2P	12	-1905.4154	3891.1507	0.0000	1.0000
F81G	15	-1896.7208	3893.8414	0.0000	1.0000
F81I	15	-1896.7414	3893.8827	0.0000	1.0000
K2PG	13	-1904.8081	3896.6295	0.0000	1.0000
K2PI	13	-1904.9503	3896.9139	0.0000	1.0000
F81IG	16	-1896.7245	3900.5422	0.0000	1.0000
K2PIG	14	-1904.8073	3903.3211	0.0000	1.0000
SYM	16	-1903.6798	3914.4527	0.0000	1.0000
SYMG	17	-1902.6702	3919.1268	0.0000	1.0000
SYMI	17	-1902.7955	3919.3776	0.0000	1.0000
SYMIG	18	-1902.6843	3925.8483	0.0000	1.0000
JC69	11	-1927.0831	3927.7928	0.0000	1.0000
JC69G	12	-1926.5426	3933.4051	0.0000	1.0000
JC69I	12	-1926.5919	3933.5036	0.0000	1.0000
JC69IG	13	-1926.5542	3940.1216	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr12_6056.phylip = 1 - 807;
 Partition Dummy = 1:683-tomodel-phylip/chr12_6056.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr12_6056.phylip = 1 - 807;
 Partition Dummy = 1:683-tomodel-phylip/chr12_6056.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr12_6056.phylip = 1 - 807;
 Partition Dummy = 1:683-tomodel-phylip/chr12_6056.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

