def METHOD_NAME():
    family_status = np.asarray(
        ["single", "single", "single", "single", "couple", "couple", "couple", "couple"]
    )
    housing_occupancy_status = np.asarray(
        ["owner", "owner", "tenant", "tenant", "owner", "owner", "tenant", "tenant"]
    )
    zone = np.asarray(["z1", "z2", "z1", "z2", "z1", "z2", "z1", "z2"])
    assert_near(
        P[family_status][housing_occupancy_status][zone],
        [100, 200, 300, 400, 500, 600, 700, 800],
    )