R = QQ[x,y]
I = ideal(x^2-1,y-x)
cv = isCellular (I,ReturnCellVars=>true)
cellularBinomialAssociatedPrimes (I,CellVariables=>cv)