Fixed a "bug" in which I predicted probabilities in LMLR using sigmoid instead of softmax. It doesn't really change anything. It would matter if we could somehow find a loss function which guaranteed summing probabilities to one during training, but I don't see any obvious candidate which is still convex.