# NOT RUN {
data("mortgage")
fm <- deny ~ black + p_irat + hse_inc + ccred + mcred + pubrec + ltv_med +
ltv_high + denpmi + selfemp + single + hischl
test <- SPE(fm = fm, data = mortgage, var.T = "black", method = "logit",
us = c(1:9)/10)
SPEplot(output = test, main="APE and SPE of Being Black on the prob of Mortgage Denial",
sub="Logit Model", ylab="Change in Probability")
# }
Run the code above in your browser using DataLab