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

Input data from file "591-tomodel-phylip/chr2_18926.phylip" (7 taxa, 773 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	-2795.0087	5631.1345	0.6734	0.6734
GTRIG	21	-2795.0397	5633.3098	0.2269	0.9004
GTRG	20	-2796.9540	5635.0251	0.0963	0.9966
HKYI	16	-2804.8959	5642.5114	0.0023	0.9989
HKYIG	17	-2804.9308	5644.6721	0.0008	0.9997
HKYG	16	-2806.8401	5646.3997	0.0003	1.0000
GTR	19	-2814.1126	5667.2344	0.0000	1.0000
HKY	15	-2823.7061	5678.0462	0.0000	1.0000
K2PI	13	-2844.5454	5715.5704	0.0000	1.0000
K2PIG	14	-2844.5784	5717.7108	0.0000	1.0000
K2PG	13	-2846.3286	5719.1367	0.0000	1.0000
SYMI	17	-2842.6897	5720.1901	0.0000	1.0000
SYMIG	18	-2842.7207	5722.3486	0.0000	1.0000
SYMG	17	-2844.3883	5723.5871	0.0000	1.0000
K2P	12	-2864.4217	5753.2538	0.0000	1.0000
SYM	16	-2862.0448	5756.8091	0.0000	1.0000
F81I	15	-2870.6071	5771.8483	0.0000	1.0000
F81IG	16	-2870.6452	5774.0099	0.0000	1.0000
F81G	15	-2872.4475	5775.5292	0.0000	1.0000
F81	14	-2888.5572	5805.6685	0.0000	1.0000
JC69I	12	-2901.2027	5826.8160	0.0000	1.0000
JC69IG	13	-2901.2381	5828.9558	0.0000	1.0000
JC69G	12	-2902.9295	5830.2696	0.0000	1.0000
JC69	11	-2919.7468	5861.8405	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2795.0087	5630.0175	0.6650	0.6650
GTRIG	21	-2795.0397	5632.0795	0.2372	0.9022
GTRG	20	-2796.9540	5633.9081	0.0951	0.9972
HKYI	16	-2804.8959	5641.7918	0.0018	0.9991
HKYIG	17	-2804.9308	5643.8615	0.0007	0.9997
HKYG	16	-2806.8401	5645.6802	0.0003	1.0000
GTR	19	-2814.1126	5666.2251	0.0000	1.0000
HKY	15	-2823.7061	5677.4121	0.0000	1.0000
K2PI	13	-2844.5454	5715.0908	0.0000	1.0000
K2PIG	14	-2844.5784	5717.1567	0.0000	1.0000
K2PG	13	-2846.3286	5718.6571	0.0000	1.0000
SYMI	17	-2842.6897	5719.3795	0.0000	1.0000
SYMIG	18	-2842.7207	5721.4415	0.0000	1.0000
SYMG	17	-2844.3883	5722.7765	0.0000	1.0000
K2P	12	-2864.4217	5752.8433	0.0000	1.0000
SYM	16	-2862.0448	5756.0895	0.0000	1.0000
F81I	15	-2870.6071	5771.2143	0.0000	1.0000
F81IG	16	-2870.6452	5773.2903	0.0000	1.0000
F81G	15	-2872.4475	5774.8951	0.0000	1.0000
F81	14	-2888.5572	5805.1144	0.0000	1.0000
JC69I	12	-2901.2027	5826.4055	0.0000	1.0000
JC69IG	13	-2901.2381	5828.4762	0.0000	1.0000
JC69G	12	-2902.9295	5829.8590	0.0000	1.0000
JC69	11	-2919.7468	5861.4936	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2804.8959	5716.1963	0.8219	0.8219
HKYG	16	-2806.8401	5720.0846	0.1176	0.9396
HKYIG	17	-2804.9308	5722.9163	0.0286	0.9681
GTRI	20	-2795.0087	5723.0231	0.0271	0.9952
GTRG	20	-2796.9540	5726.9137	0.0039	0.9991
GTRIG	21	-2795.0397	5729.7353	0.0009	1.0000
HKY	15	-2823.7061	5747.1663	0.0000	1.0000
GTR	19	-2814.1126	5754.5804	0.0000	1.0000
K2PI	13	-2844.5454	5775.5444	0.0000	1.0000
K2PG	13	-2846.3286	5779.1107	0.0000	1.0000
K2PIG	14	-2844.5784	5782.2606	0.0000	1.0000
SYMI	17	-2842.6897	5798.4342	0.0000	1.0000
SYMG	17	-2844.3883	5801.8313	0.0000	1.0000
SYMIG	18	-2842.7207	5805.1465	0.0000	1.0000
K2P	12	-2864.4217	5808.6466	0.0000	1.0000
SYM	16	-2862.0448	5830.4940	0.0000	1.0000
F81I	15	-2870.6071	5840.9684	0.0000	1.0000
F81G	15	-2872.4475	5844.6493	0.0000	1.0000
F81IG	16	-2870.6452	5847.6948	0.0000	1.0000
F81	14	-2888.5572	5870.2183	0.0000	1.0000
JC69I	12	-2901.2027	5882.2088	0.0000	1.0000
JC69G	12	-2902.9295	5885.6624	0.0000	1.0000
JC69IG	13	-2901.2381	5888.9298	0.0000	1.0000
JC69	11	-2919.7468	5912.6467	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/chr2_18926.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18926.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/chr2_18926.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18926.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/chr2_18926.phylip = 1 - 773;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18926.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

