def score_simulation(preds, endog_expected):
    """Score a simulation based on mean squared error."""
    MSE = 0
    for i in range(len(preds)):
        MSE += (preds[i] - endog_expected[i]) ** 2
    return MSE