Como posso adaptar o ANOVA para dados binários?

Tenho quatro modelos concorrentes que utilizo para prever uma variável de resultado binário (por exemplo, status de emprego após a graduação, 1 = empregado, 0 = não empregado) para n indivíduos. Uma métrica natural do desempenho do modelo é a taxa de acertos, que é a porcentagem de previsões corretas para cada um dos modelos.

Parece-me que não posso usar ANOVA nessa configuração, pois os dados violam as suposições subjacentes à ANOVA. Existe um procedimento equivalente que eu poderia usar em vez da ANOVA na configuração acima para testar a hipótese de que todos os quatro modelos são igualmente eficazes?

Tabela de contingência (qui-quadrado). A regressão logística também é sua amiga - use variáveis ​​fictícias.