Is there an R function to create a cross-sell matrix for x products and y-length combinations of them?

Question

I made a similar post here: link

I\'m now considering if and how this can be done in R.

What I am trying to achieve is a matrix with columns representing standalon