For n predictor variables, there are 2n -1 possible combinations. The best combination of variables was selected by comparing all model subsets using AIC with the function stepAIC in R package "MASS", but i don't know how to calculate the Akaike weight for each variable based on all possible combinations of predictor variables