Fix sigma estimate bug.
Return vectors, not matrices for fit methods.
Add simple test data.