
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:58 02/17/2011

Input data from file "183-phylip/chr6_4409.phylip" (29 taxa, 601 characters)

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


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	64	-3184.8054	6513.1332	0.9994	0.9994
GTRIG	65	-3191.0507	6528.1387	0.0006	1.0000
SYMIG	62	-3205.5878	6549.6960	0.0000	1.0000
SYMG	61	-3209.1836	6554.4006	0.0000	1.0000
HKYIG	61	-3209.4105	6554.8543	0.0000	1.0000
HKYG	60	-3211.6472	6556.8500	0.0000	1.0000
K2PIG	58	-3215.6728	6559.9728	0.0000	1.0000
K2PG	57	-3222.5248	6571.2264	0.0000	1.0000
GTRI	64	-3222.3069	6588.1361	0.0000	1.0000
HKYI	60	-3236.4404	6606.4363	0.0000	1.0000
SYMI	61	-3237.2899	6610.6132	0.0000	1.0000
K2PI	57	-3243.3475	6612.8718	0.0000	1.0000
F81IG	60	-3311.0798	6755.7151	0.0000	1.0000
F81G	59	-3314.1203	6759.3274	0.0000	1.0000
JC69G	56	-3321.0564	6765.8482	0.0000	1.0000
JC69IG	57	-3321.4582	6769.0932	0.0000	1.0000
F81I	59	-3336.7878	6804.6626	0.0000	1.0000
JC69I	56	-3344.9406	6813.6165	0.0000	1.0000
GTR	63	-3349.8000	6840.6168	0.0000	1.0000
SYM	60	-3364.7374	6863.0303	0.0000	1.0000
HKY	59	-3367.4433	6865.9736	0.0000	1.0000
K2P	56	-3374.1559	6872.0472	0.0000	1.0000
F81	58	-3462.5379	7053.7030	0.0000	1.0000
JC69	55	-3470.1268	7061.5564	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3184.8054	6497.6108	0.9993	0.9993
GTRIG	65	-3191.0507	6512.1013	0.0007	1.0000
SYMIG	62	-3205.5878	6535.1755	0.0000	1.0000
SYMG	61	-3209.1836	6540.3672	0.0000	1.0000
HKYIG	61	-3209.4105	6540.8209	0.0000	1.0000
HKYG	60	-3211.6472	6543.2945	0.0000	1.0000
K2PIG	58	-3215.6728	6547.3455	0.0000	1.0000
K2PG	57	-3222.5248	6559.0496	0.0000	1.0000
GTRI	64	-3222.3069	6572.6137	0.0000	1.0000
HKYI	60	-3236.4404	6592.8808	0.0000	1.0000
SYMI	61	-3237.2899	6596.5798	0.0000	1.0000
K2PI	57	-3243.3475	6600.6950	0.0000	1.0000
F81IG	60	-3311.0798	6742.1595	0.0000	1.0000
F81G	59	-3314.1203	6746.2405	0.0000	1.0000
JC69G	56	-3321.0564	6754.1129	0.0000	1.0000
JC69IG	57	-3321.4582	6756.9164	0.0000	1.0000
F81I	59	-3336.7878	6791.5757	0.0000	1.0000
JC69I	56	-3344.9406	6801.8812	0.0000	1.0000
GTR	63	-3349.8000	6825.6000	0.0000	1.0000
SYM	60	-3364.7374	6849.4747	0.0000	1.0000
HKY	59	-3367.4433	6852.8867	0.0000	1.0000
K2P	56	-3374.1559	6860.3119	0.0000	1.0000
F81	58	-3462.5379	7041.0757	0.0000	1.0000
JC69	55	-3470.1268	7050.2536	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-3184.8054	6779.1209	0.9999	0.9999
GTRIG	65	-3191.0507	6798.0100	0.0001	1.0000
K2PIG	58	-3215.6728	6802.4640	0.0000	1.0000
HKYG	60	-3211.6472	6807.2102	0.0000	1.0000
SYMIG	62	-3205.5878	6807.8884	0.0000	1.0000
SYMG	61	-3209.1836	6808.6815	0.0000	1.0000
HKYIG	61	-3209.4105	6809.1352	0.0000	1.0000
K2PG	57	-3222.5248	6809.7695	0.0000	1.0000
K2PI	57	-3243.3475	6851.4149	0.0000	1.0000
GTRI	64	-3222.3069	6854.1238	0.0000	1.0000
HKYI	60	-3236.4404	6856.7965	0.0000	1.0000
SYMI	61	-3237.2899	6864.8941	0.0000	1.0000
JC69G	56	-3321.0564	7000.4342	0.0000	1.0000
F81G	59	-3314.1203	7005.7576	0.0000	1.0000
F81IG	60	-3311.0798	7006.0752	0.0000	1.0000
JC69IG	57	-3321.4582	7007.6364	0.0000	1.0000
JC69I	56	-3344.9406	7048.2025	0.0000	1.0000
F81I	59	-3336.7878	7051.0928	0.0000	1.0000
GTR	63	-3349.8000	7102.7115	0.0000	1.0000
K2P	56	-3374.1559	7106.6332	0.0000	1.0000
HKY	59	-3367.4433	7112.4038	0.0000	1.0000
SYM	60	-3364.7374	7113.3904	0.0000	1.0000
JC69	55	-3470.1268	7292.1763	0.0000	1.0000
F81	58	-3462.5379	7296.1942	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 183-phylip/chr6_4409.phylip = 1 - 601;
 Partition Dummy = 1:183-phylip/chr6_4409.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 183-phylip/chr6_4409.phylip = 1 - 601;
 Partition Dummy = 1:183-phylip/chr6_4409.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 (GTRG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr6_4409.phylip = 1 - 601;
 Partition Dummy = 1:183-phylip/chr6_4409.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;


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

