
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:44 08/01/2011

Input data from file "261-aggressive-phylip/5723_hg18.phylip" (29 taxa, 665 characters)

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


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	-2592.2784	5326.4234	0.7243	0.7243
GTRIG	65	-2592.0264	5328.3767	0.2728	0.9971
SYMIG	62	-2600.7382	5338.4531	0.0018	0.9989
SYMG	61	-2602.4900	5339.5239	0.0010	0.9999
GTRI	64	-2602.1474	5346.1614	0.0000	0.9999
K2PG	57	-2610.8061	5346.5052	0.0000	1.0000
HKYG	60	-2608.4516	5349.0225	0.0000	1.0000
HKYIG	61	-2607.2672	5349.0784	0.0000	1.0000
K2PIG	58	-2611.9429	5351.1795	0.0000	1.0000
SYMI	61	-2610.6904	5355.9247	0.0000	1.0000
HKYI	60	-2616.8086	5365.7364	0.0000	1.0000
K2PI	57	-2621.5810	5368.0549	0.0000	1.0000
GTR	63	-2651.3128	5442.0433	0.0000	1.0000
SYM	60	-2656.3243	5444.7677	0.0000	1.0000
K2P	56	-2663.4619	5449.4237	0.0000	1.0000
JC69IG	57	-2663.5475	5451.9879	0.0000	1.0000
JC69G	56	-2665.0009	5452.5018	0.0000	1.0000
F81IG	60	-2660.2040	5452.5271	0.0000	1.0000
F81G	59	-2661.6773	5453.0571	0.0000	1.0000
HKY	59	-2662.9328	5455.5681	0.0000	1.0000
JC69I	56	-2673.1163	5468.7327	0.0000	1.0000
F81I	59	-2669.7835	5469.2695	0.0000	1.0000
JC69	55	-2716.4479	5553.0106	0.0000	1.0000
F81	58	-2713.7152	5554.7242	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2592.2784	5312.5567	0.6777	0.6777
GTRIG	65	-2592.0264	5314.0528	0.3207	0.9984
SYMIG	62	-2600.7382	5325.4763	0.0011	0.9994
SYMG	61	-2602.4900	5326.9800	0.0005	0.9999
GTRI	64	-2602.1474	5332.2948	0.0000	1.0000
K2PG	57	-2610.8061	5335.6122	0.0000	1.0000
HKYIG	61	-2607.2672	5336.5344	0.0000	1.0000
HKYG	60	-2608.4516	5336.9033	0.0000	1.0000
K2PIG	58	-2611.9429	5339.8858	0.0000	1.0000
SYMI	61	-2610.6904	5343.3808	0.0000	1.0000
HKYI	60	-2616.8086	5353.6172	0.0000	1.0000
K2PI	57	-2621.5810	5357.1619	0.0000	1.0000
GTR	63	-2651.3128	5428.6257	0.0000	1.0000
SYM	60	-2656.3243	5432.6485	0.0000	1.0000
K2P	56	-2663.4619	5438.9237	0.0000	1.0000
F81IG	60	-2660.2040	5440.4079	0.0000	1.0000
JC69IG	57	-2663.5475	5441.0950	0.0000	1.0000
F81G	59	-2661.6773	5441.3546	0.0000	1.0000
JC69G	56	-2665.0009	5442.0018	0.0000	1.0000
HKY	59	-2662.9328	5443.8657	0.0000	1.0000
F81I	59	-2669.7835	5457.5670	0.0000	1.0000
JC69I	56	-2673.1163	5458.2327	0.0000	1.0000
JC69	55	-2716.4479	5542.8957	0.0000	1.0000
F81	58	-2713.7152	5543.4305	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2610.8061	5592.1001	0.9616	0.9616
GTRG	64	-2592.2784	5600.5431	0.0141	0.9758
K2PIG	58	-2611.9429	5600.8734	0.0120	0.9877
SYMG	61	-2602.4900	5601.4670	0.0089	0.9966
SYMIG	62	-2600.7382	5604.4631	0.0020	0.9986
GTRIG	65	-2592.0264	5606.5389	0.0007	0.9993
HKYG	60	-2608.4516	5606.8905	0.0006	0.9999
HKYIG	61	-2607.2672	5611.0214	0.0001	1.0000
K2PI	57	-2621.5810	5613.6498	0.0000	1.0000
SYMI	61	-2610.6904	5617.8678	0.0000	1.0000
GTRI	64	-2602.1474	5620.2811	0.0000	1.0000
HKYI	60	-2616.8086	5623.6044	0.0000	1.0000
K2P	56	-2663.4619	5690.9118	0.0000	1.0000
JC69G	56	-2665.0009	5693.9899	0.0000	1.0000
JC69IG	57	-2663.5475	5697.5829	0.0000	1.0000
SYM	60	-2656.3243	5702.6357	0.0000	1.0000
F81G	59	-2661.6773	5706.8421	0.0000	1.0000
HKY	59	-2662.9328	5709.3531	0.0000	1.0000
JC69I	56	-2673.1163	5710.2207	0.0000	1.0000
F81IG	60	-2660.2040	5710.3952	0.0000	1.0000
GTR	63	-2651.3128	5712.1123	0.0000	1.0000
F81I	59	-2669.7835	5723.0545	0.0000	1.0000
JC69	55	-2716.4479	5790.3840	0.0000	1.0000
F81	58	-2713.7152	5804.4181	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 261-aggressive-phylip/5723_hg18.phylip = 1 - 665;
 Partition Dummy = 1:261-aggressive-phylip/5723_hg18.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 261-aggressive-phylip/5723_hg18.phylip = 1 - 665;
 Partition Dummy = 1:261-aggressive-phylip/5723_hg18.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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5723_hg18.phylip = 1 - 665;
 Partition Dummy = 1:261-aggressive-phylip/5723_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

