my design matrix Φ has a column of one's to handle the intercept. so anytime I Z-score transform, I get NaNs after dividing by a zero variance. would be nice if transform had an option to handle this, ie. not transform when sigma=0.0.
feature_transform = fit(ZScoreTransform, Φ, dims=1)
Φ̂ = StatsBase.transform(feature_transform, Φ)