Pull requests to linearmodel or statsmodel

[nick-zhy]

Thanks for your hard work. This function works prettey well. I suggest pulling it to linearmodel or statsmodel for better using experience.