How to fasten going through all independent variable combinations?

Question

I want to write function combination_rsquare(y, data, factor_number) where

y - A vector - dependent variable

data - A da