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

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

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


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	-2970.4631	6081.6480	0.5936	0.5936
GTRIG	65	-2969.6338	6082.4069	0.4061	0.9997
HKYG	60	-2983.5069	6098.1385	0.0002	0.9999
HKYIG	61	-2982.7483	6099.0095	0.0001	1.0000
K2PG	57	-3005.4423	6134.8875	0.0000	1.0000
K2PIG	58	-3005.4503	6137.2702	0.0000	1.0000
SYMG	61	-3003.2438	6140.0005	0.0000	1.0000
SYMIG	62	-3003.2517	6142.4120	0.0000	1.0000
GTRI	64	-3036.4956	6213.7129	0.0000	1.0000
HKYI	60	-3050.2468	6231.6182	0.0000	1.0000
F81IG	60	-3055.0714	6241.2675	0.0000	1.0000
F81G	59	-3058.4783	6245.7001	0.0000	1.0000
GTR	63	-3059.5193	6257.3501	0.0000	1.0000
JC69G	56	-3069.9952	6261.6338	0.0000	1.0000
JC69IG	57	-3070.0030	6264.0090	0.0000	1.0000
K2PI	57	-3073.7386	6271.4802	0.0000	1.0000
HKY	59	-3071.9997	6272.7430	0.0000	1.0000
SYMI	61	-3071.4343	6276.3816	0.0000	1.0000
K2P	56	-3095.8602	6313.3638	0.0000	1.0000
SYM	60	-3092.9257	6316.9760	0.0000	1.0000
F81I	59	-3124.5816	6377.9068	0.0000	1.0000
JC69I	56	-3137.2193	6396.0821	0.0000	1.0000
F81	58	-3147.2940	6420.9577	0.0000	1.0000
JC69	55	-3158.2530	6435.7972	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2970.4631	6068.9262	0.5425	0.5425
GTRIG	65	-2969.6338	6069.2676	0.4574	0.9999
HKYG	60	-2983.5069	6087.0139	0.0001	0.9999
HKYIG	61	-2982.7483	6087.4965	0.0001	1.0000
K2PG	57	-3005.4423	6124.8845	0.0000	1.0000
K2PIG	58	-3005.4503	6126.9005	0.0000	1.0000
SYMG	61	-3003.2438	6128.4876	0.0000	1.0000
SYMIG	62	-3003.2517	6130.5034	0.0000	1.0000
GTRI	64	-3036.4956	6200.9911	0.0000	1.0000
HKYI	60	-3050.2468	6220.4936	0.0000	1.0000
F81IG	60	-3055.0714	6230.1429	0.0000	1.0000
F81G	59	-3058.4783	6234.9565	0.0000	1.0000
GTR	63	-3059.5193	6245.0387	0.0000	1.0000
JC69G	56	-3069.9952	6251.9903	0.0000	1.0000
JC69IG	57	-3070.0030	6254.0060	0.0000	1.0000
K2PI	57	-3073.7386	6261.4771	0.0000	1.0000
HKY	59	-3071.9997	6261.9995	0.0000	1.0000
SYMI	61	-3071.4343	6264.8686	0.0000	1.0000
K2P	56	-3095.8602	6303.7203	0.0000	1.0000
SYM	60	-3092.9257	6305.8513	0.0000	1.0000
F81I	59	-3124.5816	6367.1633	0.0000	1.0000
JC69I	56	-3137.2193	6386.4386	0.0000	1.0000
F81	58	-3147.2940	6410.5880	0.0000	1.0000
JC69	55	-3158.2530	6426.5061	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2983.5069	6361.6855	0.4878	0.4878
GTRG	64	-2970.4631	6361.9094	0.4361	0.9239
HKYIG	61	-2982.7483	6366.7461	0.0388	0.9627
GTRIG	65	-2969.6338	6366.8285	0.0373	1.0000
K2PG	57	-3005.4423	6385.8226	0.0000	1.0000
K2PIG	58	-3005.4503	6392.4165	0.0000	1.0000
SYMG	61	-3003.2438	6407.7371	0.0000	1.0000
SYMIG	62	-3003.2517	6414.3308	0.0000	1.0000
GTRI	64	-3036.4956	6493.9743	0.0000	1.0000
HKYI	60	-3050.2468	6495.1653	0.0000	1.0000
F81IG	60	-3055.0714	6504.8145	0.0000	1.0000
F81G	59	-3058.4783	6505.0503	0.0000	1.0000
JC69G	56	-3069.9952	6508.3505	0.0000	1.0000
JC69IG	57	-3070.0030	6514.9441	0.0000	1.0000
K2PI	57	-3073.7386	6522.4152	0.0000	1.0000
HKY	59	-3071.9997	6532.0933	0.0000	1.0000
GTR	63	-3059.5193	6533.4439	0.0000	1.0000
SYMI	61	-3071.4343	6544.1182	0.0000	1.0000
K2P	56	-3095.8602	6560.0806	0.0000	1.0000
SYM	60	-3092.9257	6580.5230	0.0000	1.0000
F81I	59	-3124.5816	6637.2571	0.0000	1.0000
JC69I	56	-3137.2193	6642.7988	0.0000	1.0000
F81	58	-3147.2940	6676.1040	0.0000	1.0000
JC69	55	-3158.2530	6678.2884	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/11272_hg18.phylip = 1 - 719;
 Partition Dummy = 1:261-aggressive-phylip/11272_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/11272_hg18.phylip = 1 - 719;
 Partition Dummy = 1:261-aggressive-phylip/11272_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 (HKYG)]

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


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

