# score the model
log.score(X, y)