Scikit-learn implementation

Now using sci-kit Learn implementation. Way faster (C libraries I
believe), but still similar problems to statsmodels implementation: no
more singular matrices, but parameters do run out of control and
moreover from the initial (random) estimation the error actually
increases on the first run, which should not be possible. Maybe the
scale difference between learning rate and the other parameters causes
problems?