
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:54 04/04/2011

Input data from file "591-tomodel-phylip/chr8_446.phylip" (7 taxa, 610 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	20	-1170.6477	2382.7215	0.4150	0.4150
GTRI	20	-1170.6837	2382.7936	0.4003	0.8153
GTRIG	21	-1170.6253	2384.8220	0.1452	0.9605
GTR	19	-1174.8824	2389.0529	0.0175	0.9780
HKYG	16	-1178.7057	2390.3287	0.0093	0.9873
HKYI	16	-1178.7537	2390.4247	0.0088	0.9961
HKYIG	17	-1178.6846	2392.4030	0.0033	0.9994
HKY	15	-1182.4425	2395.6931	0.0006	1.0000
F81G	15	-1198.0648	2426.9378	0.0000	1.0000
F81I	15	-1198.1357	2427.0796	0.0000	1.0000
F81IG	16	-1198.0156	2428.9486	0.0000	1.0000
SYMG	17	-1198.5109	2432.0556	0.0000	1.0000
SYMI	17	-1198.5796	2432.1929	0.0000	1.0000
F81	14	-1201.7889	2432.2837	0.0000	1.0000
SYMIG	18	-1198.4860	2434.1294	0.0000	1.0000
K2PG	13	-1203.8852	2434.3811	0.0000	1.0000
K2PI	13	-1203.9549	2434.5206	0.0000	1.0000
K2PIG	14	-1203.8502	2436.4063	0.0000	1.0000
SYM	16	-1201.9266	2436.7705	0.0000	1.0000
K2P	12	-1207.3742	2439.2710	0.0000	1.0000
JC69G	12	-1221.1386	2466.7998	0.0000	1.0000
JC69I	12	-1221.2168	2466.9563	0.0000	1.0000
JC69IG	13	-1221.0977	2468.8061	0.0000	1.0000
JC69	11	-1224.7510	2471.9435	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1170.6477	2381.2953	0.4130	0.4130
GTRI	20	-1170.6837	2381.3675	0.3983	0.8113
GTRIG	21	-1170.6253	2383.2506	0.1554	0.9667
GTR	19	-1174.8824	2387.7648	0.0163	0.9829
HKYG	16	-1178.7057	2389.4113	0.0071	0.9901
HKYI	16	-1178.7537	2389.5073	0.0068	0.9969
HKYIG	17	-1178.6846	2391.3693	0.0027	0.9995
HKY	15	-1182.4425	2394.8850	0.0005	1.0000
F81G	15	-1198.0648	2426.1297	0.0000	1.0000
F81I	15	-1198.1357	2426.2715	0.0000	1.0000
F81IG	16	-1198.0156	2428.0313	0.0000	1.0000
SYMG	17	-1198.5109	2431.0218	0.0000	1.0000
SYMI	17	-1198.5796	2431.1591	0.0000	1.0000
F81	14	-1201.7889	2431.5778	0.0000	1.0000
SYMIG	18	-1198.4860	2432.9720	0.0000	1.0000
K2PG	13	-1203.8852	2433.7703	0.0000	1.0000
K2PI	13	-1203.9549	2433.9099	0.0000	1.0000
K2PIG	14	-1203.8502	2435.7004	0.0000	1.0000
SYM	16	-1201.9266	2435.8532	0.0000	1.0000
K2P	12	-1207.3742	2438.7484	0.0000	1.0000
JC69G	12	-1221.1386	2466.2772	0.0000	1.0000
JC69I	12	-1221.2168	2466.4337	0.0000	1.0000
JC69IG	13	-1221.0977	2468.1953	0.0000	1.0000
JC69	11	-1224.7510	2471.5020	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1178.7057	2460.0267	0.3842	0.3842
HKYI	16	-1178.7537	2460.1227	0.3662	0.7503
HKY	15	-1182.4425	2461.0869	0.2261	0.9764
HKYIG	17	-1178.6846	2466.3981	0.0159	0.9923
GTRG	20	-1170.6477	2469.5645	0.0033	0.9956
GTRI	20	-1170.6837	2469.6367	0.0031	0.9987
GTR	19	-1174.8824	2471.6205	0.0012	0.9999
GTRIG	21	-1170.6253	2475.9332	0.0001	1.0000
K2PG	13	-1203.8852	2491.1453	0.0000	1.0000
K2PI	13	-1203.9549	2491.2848	0.0000	1.0000
K2P	12	-1207.3742	2491.7099	0.0000	1.0000
F81G	15	-1198.0648	2492.3316	0.0000	1.0000
F81I	15	-1198.1357	2492.4734	0.0000	1.0000
F81	14	-1201.7889	2493.3662	0.0000	1.0000
K2PIG	14	-1203.8502	2497.4888	0.0000	1.0000
F81IG	16	-1198.0156	2498.6466	0.0000	1.0000
SYMG	17	-1198.5109	2506.0506	0.0000	1.0000
SYMI	17	-1198.5796	2506.1879	0.0000	1.0000
SYM	16	-1201.9266	2506.4685	0.0000	1.0000
SYMIG	18	-1198.4860	2512.4143	0.0000	1.0000
JC69G	12	-1221.1386	2519.2387	0.0000	1.0000
JC69I	12	-1221.2168	2519.3952	0.0000	1.0000
JC69	11	-1224.7510	2520.0500	0.0000	1.0000
JC69IG	13	-1221.0977	2525.5703	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 591-tomodel-phylip/chr8_446.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr8_446.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 591-tomodel-phylip/chr8_446.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr8_446.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 591-tomodel-phylip/chr8_446.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr8_446.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

