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

Input data from file "261-aggressive-phylip/8711_hg18.phylip" (29 taxa, 682 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	-4168.9443	8479.3733	0.7724	0.7724
GTRIG	65	-4168.9445	8481.8176	0.2275	1.0000
HKYG	60	-4184.9286	8501.6446	0.0000	1.0000
HKYIG	61	-4184.8146	8503.8293	0.0000	1.0000
GTRI	64	-4200.2028	8541.8901	0.0000	1.0000
HKYI	60	-4212.6093	8557.0060	0.0000	1.0000
K2PG	57	-4219.9095	8564.4152	0.0000	1.0000
K2PIG	58	-4219.5486	8566.0828	0.0000	1.0000
SYMG	61	-4217.4174	8569.0349	0.0000	1.0000
SYMIG	62	-4217.0616	8570.7435	0.0000	1.0000
K2PI	57	-4246.5712	8617.7385	0.0000	1.0000
SYMI	61	-4243.1094	8620.4189	0.0000	1.0000
GTR	63	-4281.9052	8702.8590	0.0000	1.0000
HKY	59	-4301.2569	8731.8964	0.0000	1.0000
K2P	56	-4346.8837	8815.9817	0.0000	1.0000
SYM	60	-4342.4898	8816.7671	0.0000	1.0000
F81G	59	-4352.8967	8835.1760	0.0000	1.0000
F81IG	60	-4352.6537	8837.0949	0.0000	1.0000
JC69G	56	-4375.0794	8872.3732	0.0000	1.0000
JC69IG	57	-4374.7392	8874.0746	0.0000	1.0000
F81I	59	-4376.6965	8882.7757	0.0000	1.0000
JC69I	56	-4399.0432	8920.3007	0.0000	1.0000
F81	58	-4466.9814	9060.9484	0.0000	1.0000
JC69	55	-4490.6736	9101.1876	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-4168.9443	8465.8887	0.7311	0.7311
GTRIG	65	-4168.9445	8467.8891	0.2689	1.0000
HKYG	60	-4184.9286	8489.8571	0.0000	1.0000
HKYIG	61	-4184.8146	8491.6293	0.0000	1.0000
GTRI	64	-4200.2028	8528.4055	0.0000	1.0000
HKYI	60	-4212.6093	8545.2186	0.0000	1.0000
K2PG	57	-4219.9095	8553.8191	0.0000	1.0000
K2PIG	58	-4219.5486	8555.0973	0.0000	1.0000
SYMG	61	-4217.4174	8556.8349	0.0000	1.0000
SYMIG	62	-4217.0616	8558.1231	0.0000	1.0000
K2PI	57	-4246.5712	8607.1424	0.0000	1.0000
SYMI	61	-4243.1094	8608.2189	0.0000	1.0000
GTR	63	-4281.9052	8689.8104	0.0000	1.0000
HKY	59	-4301.2569	8720.5137	0.0000	1.0000
SYM	60	-4342.4898	8804.9796	0.0000	1.0000
K2P	56	-4346.8837	8805.7673	0.0000	1.0000
F81G	59	-4352.8967	8823.7933	0.0000	1.0000
F81IG	60	-4352.6537	8825.3075	0.0000	1.0000
JC69G	56	-4375.0794	8862.1588	0.0000	1.0000
JC69IG	57	-4374.7392	8863.4784	0.0000	1.0000
F81I	59	-4376.6965	8871.3930	0.0000	1.0000
JC69I	56	-4399.0432	8910.0863	0.0000	1.0000
F81	58	-4466.9814	9049.9628	0.0000	1.0000
JC69	55	-4490.6736	9091.3473	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-4168.9443	8755.4906	0.9143	0.9143
HKYG	60	-4184.9286	8761.3589	0.0486	0.9629
GTRIG	65	-4168.9445	8762.0160	0.0350	0.9979
HKYIG	61	-4184.8146	8767.6561	0.0021	1.0000
K2PG	57	-4219.9095	8811.7458	0.0000	1.0000
HKYI	60	-4212.6093	8816.7203	0.0000	1.0000
K2PIG	58	-4219.5486	8817.5490	0.0000	1.0000
GTRI	64	-4200.2028	8818.0074	0.0000	1.0000
SYMG	61	-4217.4174	8832.8617	0.0000	1.0000
SYMIG	62	-4217.0616	8838.6749	0.0000	1.0000
K2PI	57	-4246.5712	8865.0691	0.0000	1.0000
SYMI	61	-4243.1094	8884.2457	0.0000	1.0000
GTR	63	-4281.9052	8974.8873	0.0000	1.0000
HKY	59	-4301.2569	8987.4905	0.0000	1.0000
K2P	56	-4346.8837	9059.1690	0.0000	1.0000
SYM	60	-4342.4898	9076.4814	0.0000	1.0000
F81G	59	-4352.8967	9090.7701	0.0000	1.0000
F81IG	60	-4352.6537	9096.8093	0.0000	1.0000
JC69G	56	-4375.0794	9115.5604	0.0000	1.0000
JC69IG	57	-4374.7392	9121.4051	0.0000	1.0000
F81I	59	-4376.6965	9138.3698	0.0000	1.0000
JC69I	56	-4399.0432	9163.4880	0.0000	1.0000
F81	58	-4466.9814	9312.4145	0.0000	1.0000
JC69	55	-4490.6736	9340.2239	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/8711_hg18.phylip = 1 - 682;
 Partition Dummy = 1:261-aggressive-phylip/8711_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/8711_hg18.phylip = 1 - 682;
 Partition Dummy = 1:261-aggressive-phylip/8711_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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8711_hg18.phylip = 1 - 682;
 Partition Dummy = 1:261-aggressive-phylip/8711_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;


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

