print("column means:", scores.mean(axis=0))
print("column quartiles:")
print(np.percentile(scores, [25, 50, 75], axis=0))
print("row means with NaN ignored:", np.nanmean(scores_with_missing, axis=1))
print("ordinary mean with NaN:", np.mean(scores_with_missing))