• mcp(fact = K) with K = contrMat(...) didn't work spotted by Yves Brostaux