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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: F81G


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	-1394.8159	2830.7872	0.3541	0.3541
GTRI	20	-1394.8167	2830.7889	0.3538	0.7079
GTRIG	21	-1394.8165	2832.9057	0.1228	0.8307
F81G	15	-1402.3200	2835.2956	0.0372	0.8679
F81I	15	-1402.3248	2835.3054	0.0370	0.9049
HKYG	16	-1401.7154	2836.1751	0.0239	0.9288
HKYI	16	-1401.7207	2836.1857	0.0238	0.9526
GTR	19	-1398.6257	2836.2954	0.0225	0.9752
F81IG	16	-1402.3030	2837.3503	0.0133	0.9885
HKYIG	17	-1401.6997	2838.2379	0.0085	0.9970
F81	14	-1406.3810	2841.3350	0.0018	0.9988
HKY	15	-1405.7740	2842.2038	0.0012	1.0000
SYMI	17	-1419.5350	2873.9084	0.0000	1.0000
SYMG	17	-1419.5391	2873.9165	0.0000	1.0000
SYMIG	18	-1419.5488	2876.0358	0.0000	1.0000
SYM	16	-1423.3459	2879.4359	0.0000	1.0000
JC69G	12	-1428.3615	2881.1474	0.0000	1.0000
JC69I	12	-1428.3654	2881.1552	0.0000	1.0000
K2PG	13	-1427.8920	2882.2799	0.0000	1.0000
K2PI	13	-1427.8963	2882.2886	0.0000	1.0000
JC69IG	13	-1428.3541	2883.2042	0.0000	1.0000
K2PIG	14	-1427.8834	2884.3398	0.0000	1.0000
JC69	11	-1432.3174	2886.9934	0.0000	1.0000
K2P	12	-1431.8673	2888.1590	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1394.8159	2829.6318	0.3626	0.3626
GTRI	20	-1394.8167	2829.6334	0.3623	0.7250
GTRIG	21	-1394.8165	2831.6330	0.1333	0.8583
F81G	15	-1402.3200	2834.6399	0.0296	0.8880
F81I	15	-1402.3248	2834.6497	0.0295	0.9175
GTR	19	-1398.6257	2835.2515	0.0218	0.9393
HKYG	16	-1401.7154	2835.4309	0.0200	0.9593
HKYI	16	-1401.7207	2835.4415	0.0199	0.9791
F81IG	16	-1402.3030	2836.6061	0.0111	0.9902
HKYIG	17	-1401.6997	2837.3995	0.0075	0.9977
F81	14	-1406.3810	2840.7621	0.0014	0.9991
HKY	15	-1405.7740	2841.5480	0.0009	1.0000
SYMI	17	-1419.5350	2873.0701	0.0000	1.0000
SYMG	17	-1419.5391	2873.0781	0.0000	1.0000
SYMIG	18	-1419.5488	2875.0976	0.0000	1.0000
SYM	16	-1423.3459	2878.6917	0.0000	1.0000
JC69G	12	-1428.3615	2880.7229	0.0000	1.0000
JC69I	12	-1428.3654	2880.7307	0.0000	1.0000
K2PG	13	-1427.8920	2881.7840	0.0000	1.0000
K2PI	13	-1427.8963	2881.7927	0.0000	1.0000
JC69IG	13	-1428.3541	2882.7083	0.0000	1.0000
K2PIG	14	-1427.8834	2883.7668	0.0000	1.0000
JC69	11	-1432.3174	2886.6347	0.0000	1.0000
K2P	12	-1431.8673	2887.7345	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
F81G	15	-1402.3200	2903.9009	0.3743	0.3743
F81I	15	-1402.3248	2903.9107	0.3725	0.7468
F81	14	-1406.3810	2905.4057	0.1764	0.9232
HKYG	16	-1401.7154	2909.3093	0.0251	0.9483
HKYI	16	-1401.7207	2909.3199	0.0249	0.9732
F81IG	16	-1402.3030	2910.4845	0.0139	0.9871
HKY	15	-1405.7740	2910.8091	0.0118	0.9990
HKYIG	17	-1401.6997	2915.8954	0.0009	0.9999
GTRG	20	-1394.8159	2921.9798	0.0000	0.9999
GTRI	20	-1394.8167	2921.9815	0.0000	1.0000
GTR	19	-1398.6257	2922.9821	0.0000	1.0000
GTRIG	21	-1394.8165	2928.5985	0.0000	1.0000
JC69G	12	-1428.3615	2936.1318	0.0000	1.0000
JC69I	12	-1428.3654	2936.1396	0.0000	1.0000
JC69	11	-1432.3174	2937.4262	0.0000	1.0000
K2PG	13	-1427.8920	2941.8102	0.0000	1.0000
K2PI	13	-1427.8963	2941.8189	0.0000	1.0000
JC69IG	13	-1428.3541	2942.7345	0.0000	1.0000
K2P	12	-1431.8673	2943.1434	0.0000	1.0000
K2PIG	14	-1427.8834	2948.4105	0.0000	1.0000
SYMI	17	-1419.5350	2951.5659	0.0000	1.0000
SYMG	17	-1419.5391	2951.5740	0.0000	1.0000
SYM	16	-1423.3459	2952.5702	0.0000	1.0000
SYMIG	18	-1419.5488	2958.2108	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/chr5_11955.phylip = 1 - 748;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11955.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/chr5_11955.phylip = 1 - 748;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11955.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 (F81G)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11955.phylip = 1 - 748;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11955.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

