How do I write a Python function that is equivalent to the model.matrix function in R?

Question

Rpy2 is not suggested since it does not produce good results.