from sklearn.model_selection import train_test_split
from sklearn.linear_model import LogisticRegression
X, y = data[features), data["is_male"]
X_tr, X_te, y_tr, y_te = train_test_split(
X, y, test size=0.3, random_state=0, stratify=y)
medians = X_tr.median()
X_tr = X_tr.fillna(medians)
X_te = X_te.fillna (medians)
model = LogisticRegression(max_iter=1000)
model.fit(X_tr, y_tr)
score model.score(X_te, y_te)