Subsetting data with specific values in glm function

Question

I need to examine expenditure for specific incomes of $20K and 40K. My glm() works, but when I add subset, I get an error:

glm(district21$expend ~ 1 + income, f