I am trying to solve a MILP supply network optimization problem using R\'s OMPR package. I need to define a 3D binary variable (x[i,j,q]) in the model to map Customers (i) to DC