fitting raw linear regression to probabilities yields illegal values outside
[0,1] and is outlier-sensitive; LR guarantees valid probabilities, and the CE gradient
p−y decays only linearly as
z→±∞ (never saturates), whereas squared loss
21(p−y)2 gives vanishing gradients near
p→y — slow to train; adding
2λ∥w∥2 yields regularized LR for high-dimensional collinearity.