n_cols = 1 + len(ESTIMATORS)
plt.figure(figsize=(2. * n_cols, 2.26 * n_faces))
plt.suptitle("Face completion with multi-output estimators", size=16)
true_face = np.hstack((X_test[i], y_test[i]))
sub = plt.subplot(n_faces, n_cols, i * n_cols + 1)
sub = plt.subplot(n_faces, n_cols, i * n_cols + 1,
sub.imshow(true_face.reshape(image_shape),
for j, est in enumerate(sorted(ESTIMATORS)):
completed_face = np.hstack((X_test[i], y_test_predict[est][i]))
sub = plt.subplot(n_faces, n_cols, i * n_cols + 2 + j)
sub = plt.subplot(n_faces, n_cols, i * n_cols + 2 + j,
sub.imshow(completed_face.reshape(image_shape),