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

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

Minimum AIC  model: GTRIG
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	-2443.2953	5034.4950	0.5721	0.5721
GTRIG	65	-2442.2500	5035.0756	0.4279	1.0000
SYMIG	62	-2462.2171	5067.0342	0.0000	1.0000
SYMG	61	-2463.5666	5067.0999	0.0000	1.0000
GTRI	64	-2464.0369	5075.9780	0.0000	1.0000
HKYIG	61	-2473.3841	5086.7349	0.0000	1.0000
HKYG	60	-2475.0068	5087.3596	0.0000	1.0000
K2PIG	58	-2480.0267	5092.1949	0.0000	1.0000
K2PG	57	-2481.7825	5093.1227	0.0000	1.0000
SYMI	61	-2480.8761	5101.7189	0.0000	1.0000
HKYI	60	-2493.1480	5123.6419	0.0000	1.0000
K2PI	57	-2499.5348	5128.6273	0.0000	1.0000
JC69IG	57	-2527.0722	5183.7020	0.0000	1.0000
F81G	59	-2525.2486	5185.2348	0.0000	1.0000
F81IG	60	-2523.9684	5185.2827	0.0000	1.0000
JC69G	56	-2530.9613	5188.9086	0.0000	1.0000
GTR	63	-2533.0661	5211.3780	0.0000	1.0000
F81I	59	-2539.6783	5214.0941	0.0000	1.0000
JC69I	56	-2545.8486	5218.6830	0.0000	1.0000
SYM	60	-2552.2239	5241.7939	0.0000	1.0000
HKY	59	-2565.4499	5265.6373	0.0000	1.0000
K2P	56	-2574.1866	5275.3591	0.0000	1.0000
F81	58	-2613.2051	5358.5517	0.0000	1.0000
JC69	55	-2619.3556	5363.1374	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2442.2500	5014.5001	0.5113	0.5113
GTRG	64	-2443.2953	5014.5907	0.4887	1.0000
SYMIG	62	-2462.2171	5048.4342	0.0000	1.0000
SYMG	61	-2463.5666	5049.1331	0.0000	1.0000
GTRI	64	-2464.0369	5056.0737	0.0000	1.0000
HKYIG	61	-2473.3841	5068.7681	0.0000	1.0000
HKYG	60	-2475.0068	5070.0136	0.0000	1.0000
K2PIG	58	-2480.0267	5076.0534	0.0000	1.0000
K2PG	57	-2481.7825	5077.5651	0.0000	1.0000
SYMI	61	-2480.8761	5083.7521	0.0000	1.0000
HKYI	60	-2493.1480	5106.2959	0.0000	1.0000
K2PI	57	-2499.5348	5113.0696	0.0000	1.0000
F81IG	60	-2523.9684	5167.9367	0.0000	1.0000
JC69IG	57	-2527.0722	5168.1443	0.0000	1.0000
F81G	59	-2525.2486	5168.4972	0.0000	1.0000
JC69G	56	-2530.9613	5173.9227	0.0000	1.0000
GTR	63	-2533.0661	5192.1322	0.0000	1.0000
F81I	59	-2539.6783	5197.3565	0.0000	1.0000
JC69I	56	-2545.8486	5203.6971	0.0000	1.0000
SYM	60	-2552.2239	5224.4479	0.0000	1.0000
HKY	59	-2565.4499	5248.8997	0.0000	1.0000
K2P	56	-2574.1866	5260.3732	0.0000	1.0000
F81	58	-2613.2051	5342.4102	0.0000	1.0000
JC69	55	-2619.3556	5348.7111	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2443.2953	5282.1117	0.8854	0.8854
GTRIG	65	-2442.2500	5286.2011	0.1146	1.0000
SYMG	61	-2463.5666	5304.1141	0.0000	1.0000
SYMIG	62	-2462.2171	5307.5952	0.0000	1.0000
K2PG	57	-2481.7825	5315.8260	0.0000	1.0000
K2PIG	58	-2480.0267	5318.4943	0.0000	1.0000
HKYG	60	-2475.0068	5320.8146	0.0000	1.0000
GTRI	64	-2464.0369	5323.5948	0.0000	1.0000
HKYIG	61	-2473.3841	5323.7491	0.0000	1.0000
SYMI	61	-2480.8761	5338.7332	0.0000	1.0000
K2PI	57	-2499.5348	5351.3306	0.0000	1.0000
HKYI	60	-2493.1480	5357.0969	0.0000	1.0000
JC69IG	57	-2527.0722	5406.4053	0.0000	1.0000
JC69G	56	-2530.9613	5408.0036	0.0000	1.0000
F81G	59	-2525.2486	5415.1182	0.0000	1.0000
F81IG	60	-2523.9684	5418.7377	0.0000	1.0000
JC69I	56	-2545.8486	5437.7781	0.0000	1.0000
F81I	59	-2539.6783	5443.9775	0.0000	1.0000
GTR	63	-2533.0661	5455.4732	0.0000	1.0000
SYM	60	-2552.2239	5475.2489	0.0000	1.0000
K2P	56	-2574.1866	5494.4541	0.0000	1.0000
HKY	59	-2565.4499	5495.5207	0.0000	1.0000
JC69	55	-2619.3556	5578.6121	0.0000	1.0000
F81	58	-2613.2051	5584.8511	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_17003.phylip = 1 - 483;
 Partition Dummy = 1:183-phylip/chr2_17003.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

