RESUMO

OBJETIVO

Predizer o risco de ausência laboral decorrente de morbidades dos docentes que atuam na educação infantil na rede pública municipal, com o uso de algoritmos de machine learning.

MÉTODOS

Trata-se de um estudo transversal utilizando dados secundários, públicos e anônimos da Relação Anual de Informações Sociais, selecionando professores da educação infantil que atuaram na rede pública municipal do estado de São Paulo entre 2014 e 2018 (n = 174.294). Foram também vinculados dados da média de alunos por turma e número de habitantes no município. Os dados foram separados em treinamento e teste, utilizando os registros de 2014 a 2016 (n = 103.357) para treinar cinco modelos preditivos e os dados de 2017 a 2018 (n = 70.937) para testar seus desempenhos em dados novos. A performance preditiva dos algoritmos foi avaliada por meio do valor da área abaixo da curva ROC (AUROC).

RESULTADOS

Todos os cinco algoritmos testados apresentaram área abaixo da curva acima de 0,76. O algoritmo com melhor performance preditiva (redes neurais artificiais) obteve 0,79 de área abaixo da curva, com acurácia de 71,52%, sensibilidade de 72,86%, especificidade de 70,52% e kappa de 0,427 nos dados de teste.

CONCLUSÃO

É possível predizer casos de afastamentos por morbidade em docentes da rede pública com machine learning usando dados públicos. O melhor algoritmo apresentou melhor resultado da área abaixo da curva quando comparado ao modelo de referência (regressão logística). Os algoritmos podem contribuir para predições mais assertivas na área da saúde pública e da saúde do trabalhador, permitindo acompanhar e ajudar a prevenir afastamentos por morbidade desses trabalhadores.

Absenteísmo; Fatores de Risco; Aprendizado de Máquina Supervisionado; Professores Escolares; Educação Infantil

INTRODUÇÃO

Profissionais da área da educação básica apresentam altos índices de afastamento por doença relacionada ao trabalho, situando-se entre as primeiras posições no Brasil¹1. DIEESE. Anuário do Sistema Público de Emprego, Trabalho e Renda: mercado de trabalho 2016. São Paulo; 2016.. Entre as suas principais doenças, estão os distúrbios de voz, doenças do aparelho respiratório e do sistema osteomuscular, além de transtornos mentais e comportamentais²2. Gasparini SM, Barreto SM, Assunção AA. O professor, as condições de trabalho e os efeitos sobre sua saúde. Educ Pesqui. 2005;31(2):189-99. https://doi.org/10.1590/S1517-97022005000200003
https://doi.org/10.1590/S1517-9702200500... ,³3. Medeiros AM, Vieira MT. Ausência ao trabalho por distúrbio vocal de professores da Educação Básica no Brasil. Cad Saude Publica. 2019;35 Supl 1:e00171717. https://doi.org/10.1590/0102-311x00171717
https://doi.org/10.1590/0102-311x0017171... . Houve um crescimento recente de notificações relacionadas à saúde mental, como no caso do estado de São Paulo⁴4. Arcoverde, L, Franco E, Galvão D, Prado G. Número de professores afastados por transtornos em SP quase dobra em 2016 e vai a 50 mil. G1 Globo News (São Paulo Ed.). 21 nov 2017 [cited 2020 Feb 7]. Available from: https://g1.globo.com/sp/sao-paulo/noticia/numero-de-professores-afastados-por-transtornos-em-sp-quase-dobra-em-2016-e-vai-a-50-mil.ghtml
https://g1.globo.com/sp/sao-paulo/notici... . As condições de trabalho a que esses profissionais estão expostos, com longas jornadas, turmas numerosas⁵5. Rodríguez-Loureiro L, Artazcoz L, López-Ruiz M, Assunção AA, Benavides FG. Joint effect of paid working hours and multiple job holding on work absence due to health problems among basic education teachers in Brazil: the Educatel Study. Cad Saude Publica. 2019;35 Supl 1:e00081118. https://doi.org/10.1590/0102-311x00081118
https://doi.org/10.1590/0102-311x0008111... e falta de reconhecimento, acabam agravando o problema, com consequências diretas na vida particular dos docentes⁶6. Silva J, Fischer FM. Invasão multiforme da vida pelo trabalho entre professores de educação básica e repercussões sobre a saúde. Rev Saude Publica. 2020;54:03. https://doi.org/10.11606/s1518-8787.2020054001547
https://doi.org/10.11606/s1518-8787.2020... .

No entanto, ainda há poucos indicadores que analisam as condições dos docentes em cada etapa do ensino básico, desde a educação infantil até o ensino médio, áreas que apresentam distinções e especificidades no contexto de trabalho⁷7. Assunção AA, Oliveira DA. Intensificação do trabalho e saúde dos professores. Educ Soc. 2009;30(107):349-72. https://doi.org/10.1590/S0101-73302009000200003
https://doi.org/10.1590/S0101-7330200900... . Estudos que buscam gerar indicadores normalmente utilizam valores agregados das diferentes etapas do ensino básico¹1. DIEESE. Anuário do Sistema Público de Emprego, Trabalho e Renda: mercado de trabalho 2016. São Paulo; 2016. ou são baseados em inquéritos em regiões específicas que dependem da ação do respondente, podendo gerar altos custos operacionais⁸8. Porto LA, Oliveira NF, Carvalho FM, Araújo TM. Construção de um índice de morbidade para professoras da educação básica. Rev Baiana Saude Publica. 2008;32(2):282-96. https://doi.org/10.22278/2318-2660.2008.v32.n2.a1449
https://doi.org/10.22278/2318-2660.2008.... ,⁹9. Maia EG, Claro RM, Assunção AA. Múltiplas exposições ao risco de faltar ao trabalho nas escolas da Educação Básica no Brasil. Cad Saude Publica. 2019;35 Supl 1:e00166517. https://doi.org/10.1590/0102-311x00166517
https://doi.org/10.1590/0102-311x0016651... .

O uso de indicadores agregados para a análise de doenças pode omitir situações críticas de morbidade e não corresponder à realidade das diferentes condições de trabalho às quais os docentes estão expostos. O problema se agrava ao analisar professores da rede pública municipal, cujos dados de afastamento somente são acessíveis por meio de secretarias de saúde ou departamentos de perícias médicas de cada município²2. Gasparini SM, Barreto SM, Assunção AA. O professor, as condições de trabalho e os efeitos sobre sua saúde. Educ Pesqui. 2005;31(2):189-99. https://doi.org/10.1590/S1517-97022005000200003
https://doi.org/10.1590/S1517-9702200500... . Uma forma de analisar a situação de morbidade por etapa de ensino e a relação do trabalho com o adoecimento de professores é buscar fontes de dados oficiais públicas, que possuam registros individualizados e atualizados com ampla cobertura geográfica.

Nesse cenário, a partir de 2011, o acesso aos dados provenientes do governo federal foi ampliado, mediante a publicação da Lei nº 12.527 de 2011, conhecida como Lei de Acesso à Informação (LAI). Essa legislação permite que sejam solicitadas informações a qualquer órgão público e prevê que as entidades e órgãos públicos se antecipem e publiquem seus dados e informações na internet em formatos abertos e não proprietários, conhecidos como dados abertos¹⁰10. Portal Brasileiro de Dados Abertos. Primeira Lei de Acesso no mundo que prevê dados abertos. Brasília, DF; 2018 [cited 2018 Oct 1]. Available from: http://dados.gov.br/noticia/primeira-lei-de-acesso-no-mundo-que-preve-dados-abertos
http://dados.gov.br/noticia/primeira-lei... . Um exemplo é a Relação Anual de Informações Sociais (RAIS), instituída pelo Decreto nº 76.900, de 23 de dezembro de 1975, atualmente sob responsabilidade do Ministério da Economia. Os dados são divulgados anualmente de forma pública e anônima, contendo registros individualizados com cobertura em todo o território nacional¹¹11. Ministério do Trabalho (BR). PDET – Progama de Disseminação das Estatísticas do Trabalho. Microdados RAIS e CAGED. Brasília, DF; 2019 [cited 2020 Feb 7]. Available from: http://pdet.mte.gov.br/microdados-rais-e-caged
http://pdet.mte.gov.br/microdados-rais-e... .

A RAIS tem como principal finalidade subsidiar a geração de estatísticas do trabalho e do mercado formal (privado e estatutário) para as entidades governamentais. No entanto, também são disponibilizados os tipos de afastamento, como doença ou doença relacionada ao trabalho, sem informações da doença específica. O preenchimento é obrigatório a todas as empresas públicas e privadas desde 1977¹²12. Brasil. Decreto Nº 76.900, de 23 de dezembro de 1975. Institui a Relação Anual de Informações Sociais – RAIS e dá outras providências. Brasília, DF; 1975 [cited 2020 Feb 28]. Available from: http://www.planalto.gov.br/ccivil_03/decreto/Antigos/D76900.htm
http://www.planalto.gov.br/ccivil_03/dec... . Nos dados de afastamento de professores na RAIS, verifica-se que em média mais de 40% dos professores da educação infantil se afastaram do trabalho ao menos uma vez ao ano devido a morbidades no período entre 2014 e 2018, superando as demais etapas do ensino básico e apontando para a necessidade de analisar esse grupo específico.

O presente estudo propõe construir modelos preditivos capazes de estimar o risco de ausência laboral dos docentes da educação infantil que atuam na rede pública municipal, considerando todos os municípios do estado de São Paulo, utilizando inteligência artificial e algoritmos de aprendizagem supervisionada (machine learning) conhecidos pelo alto desempenho, principalmente na área da saúde¹³13. Fernandes FT, Chiavegatto Filho ADP. Perspectivas do uso de mineração de dados e aprendizado de máquina em saúde e segurança no trabalho. Rev Bras Saude Ocup. 2019;44:e13. https://doi.org/10.1590/2317-6369000019418
https://doi.org/10.1590/2317-63690000194... . Espera-se que o melhor modelo proposto possa ser aplicado para estimar o número de afastamentos por doença em professores, utilizando dados projetados dos trabalhadores (por exemplo, tempo de ocupação, carga horária e renda) e do ambiente de trabalho (por exemplo, alunos por turma), de forma a subsidiar políticas públicas direcionadas de prevenção de ausência laboral, evitando novos adoecimentos, melhorando a qualidade de vida desses profissionais e, consequentemente, a qualidade do ensino na rede pública.

MÉTODOS

Trata-se de um estudo transversal para predizer risco de afastamento por morbidade (doença ou doença relacionada ao trabalho) em professores que atuam na educação infantil (período pré-escolar) de toda a rede pública municipal do estado de São Paulo, pelo uso de dados secundários, públicos e anônimos.

Fontes de Dados

Foram utilizados dados individualizados em seu menor nível de desagregação, conhecidos também como microdados¹¹11. Ministério do Trabalho (BR). PDET – Progama de Disseminação das Estatísticas do Trabalho. Microdados RAIS e CAGED. Brasília, DF; 2019 [cited 2020 Feb 7]. Available from: http://pdet.mte.gov.br/microdados-rais-e-caged
http://pdet.mte.gov.br/microdados-rais-e... , de vínculos empregatícios formais fornecidos pela RAIS no período entre 2014 e 2018. A RAIS mantém informações fornecidas anualmente pelas empresas públicas e privadas sobre os vínculos empregatícios formais de seus funcionários, como contratações e demissões realizadas no ano corrente, além de dados sobre jornada de trabalho e afastamentos por doença ou doença relacionada ao trabalho. Foram incluídos professores da rede pública municipal de todos os municípios do estado de São Paulo que informaram dados à RAIS, nos quais os docentes possuíam nível superior e vínculo empregatício ativo até 31 de dezembro do ano corrente atuando no ensino pré-escolar, classificados na RAIS com a CBO 231105 (Professores de Ensino Superior no Ensino Infantil, quatro a seis anos).

Como forma de identificação do desfecho, foram utilizados os campos de causa de afastamento presentes na RAIS, filtrando-se pelos códigos 30 (doença de trabalho) e 40 (doença). Não foram consideradas aposentadorias decorrentes de doenças. Para a identificação do tipo de estabelecimento (federal, estadual, municipal ou privado), foi utilizada a informação da natureza jurídica do estabelecimento (códigos 1031, 1066, 1120, 1155, 1180 e 1244). Para identificação da atividade do estabelecimento, foram utilizados o código 841160 (administração pública em geral) e a divisão 85 (educação) da Classificação Nacional de Atividades Econômicas versão 2.0 (CNAE 2.0) fornecidos na RAIS.

Cerca de 99,7% (n = 174.346) dos registros estavam cadastrados no CNAE 841160, sendo que o restante, 0,03% (n = 52), pertencia a um CNAE não relacionado ao ensino infantil pré-escolar, sendo excluído da seleção. Também foram excluídos registros em que o município de exercício do trabalhador não correspondia à unidade federativa de São Paulo, resultando no total de 174.294 registros. Nas situações em que o município de exercício do trabalhador não foi informado, foi adotado o município do estabelecimento.

Para capturar parte das condições de trabalho a que esses profissionais estavam expostos, foi vinculada a informação da média de alunos por turma na etapa pré-escolar, fornecida pelo Instituto Nacional de Estudos e Pesquisas Educacionais Anísio Teixeira (Inep)¹⁴14. Instituto Nacional de Estudos e Pesquisas Educacionais Anísio Teixeira. Indicadores Educacionais. Brasília, DF: INEP; 2019 [cited 2019 Sep 18]. Available from: http://portal.inep.gov.br/web/guest/indicadores-educacionais
http://portal.inep.gov.br/web/guest/indi... , considerando que algumas doenças comuns aos docentes estão associadas à acústica do ambiente e à complexidade na atuação do docente³3. Medeiros AM, Vieira MT. Ausência ao trabalho por distúrbio vocal de professores da Educação Básica no Brasil. Cad Saude Publica. 2019;35 Supl 1:e00171717. https://doi.org/10.1590/0102-311x00171717
https://doi.org/10.1590/0102-311x0017171... , e adicionado o número de habitantes por município, fornecido pela Fundação Seade¹⁵15. Fundação Sistema Estadual de Análise de Dados. Informações dos Municípios Paulistas – IMP. São Paulo: SEADE; 2019 [cited 2019 Feb 18]. Available from: http://www.imp.seade.gov.br/frontend/#/
http://www.imp.seade.gov.br/frontend/#/... . Os dados do Inep e da Fundação Seade foram vinculados à amostra selecionada da RAIS por meio do código do município no Instituto Brasileiro de Geografia e Estatística (IBGE).

Ao todo, foram selecionadas 11 variáveis preditoras: sexo, idade, pós-graduação, tamanho do estabelecimento, tipo de vínculo empregatício, primeiro emprego, tempo no emprego, horas contratadas, média anual da quantidade de salários mínimos, média de alunos por turma e número de habitantes no município de atuação do trabalhador.

Preparação dos Dados

Os dados passaram por uma etapa inicial de pré-processamento para tratar as variáveis ausentes e transformá-los para que pudessem ser utilizados na construção dos diferentes tipos de modelos preditivos de machine learning¹⁶16. Kuhn M, Johnson K. Applied predictive modeling. New York: Springer Science & Business Media; 2013.. Todas as variáveis selecionadas estavam totalmente preenchidas, sem valores ausentes (missing). Variáveis com mais de duas categorias foram representadas por um conjunto de variáveis chamadas dummy, em que para cada categoria é gerada uma nova variável com valores de 0 ou 1. As variáveis contínuas foram padronizadas por meio de escore-z (z-score). Foi testada a correlação entre as variáveis numéricas e constatada alta correlação entre o número de habitantes no município e a média de alunos (0,93), portanto optou-se por dicotomizar a variável de número de habitantes utilizando a definição de 500 mil habitantes para identificar atuação em cidade grande¹⁷17. OECD Data. Urban population by city size. Paris; 2018 [cited 2019 Feb 18]. Available from: https://data.oecd.org/popregion/urban-population-by-city-size.htm
https://data.oecd.org/popregion/urban-po... .

Para a fase de pré-processamento e carga de dados, foi utilizado o sistema gerenciador de banco de dados MS SQL Server. Para a análise de dados e construção dos modelos preditivos, foi utilizado o software R.

Construção dos Modelos Preditivos

Ao todo, foram desenvolvidos cinco algoritmos de aprendizado supervisionado de machine learning: regressão logística, árvores de decisão¹⁶16. Kuhn M, Johnson K. Applied predictive modeling. New York: Springer Science & Business Media; 2013., random forest¹⁸18. Breiman L. Random forests. Mach Learn. 2001;45:5-32. https://doi.org/10.1023/A:1010933404324
https://doi.org/10.1023/A:1010933404324... , XGBoost¹⁹19. Chen T, Guestrin C. XGBoost: a scalable tree boosting system. In: Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining; 2016 Aug; San Francisco, CA. New York: Association for Computing Machinery; 2016. p.785-94. e redes neurais artificiais²⁰20. Bishop C. Neural networks for pattern recognition. Oxford (UK): Oxford University Press; 1995.. Algoritmos de machine learning estão sujeitos a problemas como o subajuste (underfitting) – quando o modelo não consegue se ajustar à variabilidade inerente dos dados, gerando estimativas ruins – e sobreajuste (overfitting) – quando o modelo se ajusta muito bem para a amostra utilizada no treinamento, mas não obtém bons resultados para novos dados. Uma forma de melhorar o desempenho dos algoritmos e evitar problemas como overfitting ou undefitting é aplicar técnicas de reamostragem (resampling) e validação cruzada (cross-validation) para a seleção de hiperparâmetros¹⁶16. Kuhn M, Johnson K. Applied predictive modeling. New York: Springer Science & Business Media; 2013..

Na validação cruzada, o conjunto de dados é separado em duas partes. A primeira parte é destinada ao treinamento do algoritmo e a segunda parte é utilizada para ajustes dos hiperparâmetros do modelo, simulando novos dados e selecionando hiperparâmetros que otimizem a métrica de performance escolhida. Essa fase de treinamento permite que o algoritmo obtenha um melhor desempenho preditivo quando são apresentados novos conjuntos de dados. Uma das técnicas de validação cruzada mais conhecidas é a k-fold, na qual k significa o número de subconjuntos dos dados de testes que serão utilizados para ajustes dos hiperparâmetros do modelo durante a fase de treinamento²¹21. Hastie T, Tibshirani R, Friedman J. The elements of statistical learning. 2. ed. New York: Springer Science & Business Media; 2016..

Para a construção dos modelos, separou-se uma parte dos dados (60%) referente ao período entre 2014 e 2016 (n = 103.357) para treinamento, incluindo validação cruzada com a técnica k-fold utilizando 10 partições para definição dos hiperparâmetros. A outra parte dos registros (40%), referente ao período entre 2017 e 2018 (n = 70.937), foi utilizada para a avaliação de desempenho dos modelos, omitindo-se o resultado sobre afastamento por doença e aplicando o modelo treinado para estimar o desfecho em dados futuros. Para avaliar o desempenho dos modelos, foram analisadas medidas como acurácia, sensibilidade, especificidade, valor preditivo positivo (PPV) e valor preditivo negativo (NPV). Para selecionar o melhor modelo, foi utilizado o valor da área abaixo da curva receiver operating characteristic (AUROC)²²22. Bradley AP. The use of the area under the ROC curve in the evaluation of machine learning algorithms. Pattern Recognit. 1997;30(7):1145-59. https://doi.org/10.1016/S0031-3203(96)00142-2
https://doi.org/10.1016/S0031-3203(96)00... . Foi utilizado o coeficiente kappa de Cohen para avaliar a concordância entre os valores previstos e observados.

Este estudo segue as orientações de descrição de modelos preditivos multivariados para diagnóstico ou prognóstico (Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis – TRIPOD)²³23. Moons KGM, Altman DG, Reitsma JB, Ioannidis JPA, Macaskill P, Steyerberg EW, et al. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): explanation and elaboration. Ann Intern Med. 2015;162(1):W1-73. https://doi.org/10.7326/M14-0698
https://doi.org/10.7326/M14-0698... .

Critérios de Confidencialidade

O projeto foi aprovado pelo Comitê de Ética em Pesquisa da Faculdade de Saúde Pública da Universidade de São Paulo, sob o nº 4.031.362, CAAE 30786620.5.0000.5421.

RESULTADOS

População de Estudo

Foi realizada inicialmente a análise descritiva dos dados de treinamento (2014 a 2016) e teste (2017 a 2018). Os dados de treinamento compreendem 103.357 registros de vínculos empregatícios de professores, havendo 45.419 (43,94%) registros de afastamentos por doença ou doença relacionada ao trabalho. Os dados de teste correspondem a 70.937 registros e 30.261 casos de afastamento por doença (42,65%). A Tabela 1 mostra os resultados dessa análise.

Em seguida, foi realizada a análise das variáveis independentes para verificar se há relação significativa com o desfecho, utilizando o teste de associação pelo qui-quadrado de Pearson para as variáveis categóricas. A Tabela 2 mostra os resultados e permite observar que a maioria das variáveis apresenta associação significativa com o desfecho.

Análise dos Modelos Preditivos

A Figura 1 mostra as curvas ROC de cada modelo desenvolvido, nos quais foram incluídas 11 variáveis (tipo de vínculo empregatício, número de horas contratadas, tempo de emprego em meses, média da quantidade de salários mínimos do ano, sexo, tamanho do estabelecimento, possuir pós-graduação, idade, primeiro emprego, média de alunos por turma pré-escolar no município e número de habitantes).

A Tabela 3 mostra os resultados do desempenho de cada algoritmo. O melhor modelo foi o de redes neurais artificiais, que obteve 0,79 de área abaixo da curva, com acurácia de 71,52%, sensibilidade de 72,86%, especificidade de 70,52%, PPV de 64,77% e NPV de 77,74%. O modelo também apresentou a melhor coeficiente kappa, de 0,4266.

Os hiperparâmetros dos modelos de machine learning que obtiveram os maiores valores de AUC foram os seguintes: árvore de decisão (cp-complexity-parameter = 0.000385), random forest (mtry = 7), XGBoost (nrounds = 100, max_depth = 4, gamma = 0.6, colsample_bytree = 0.8, min_child_weight = 2, subsample = 1) e redes neurais artificiais (size = 19, decay = 0.2).

Cada algoritmo utiliza um conjunto de regras diferentes para se ajustar aos dados com o menor erro possível, sendo que algumas variáveis se tornam mais importantes que outras para a capacidade preditiva do algoritmo. A importância das variáveis refere-se a quanto a variável contribuiu para realizar predições acuradas em cada modelo. Para os modelos lineares, utiliza-se de estatística t tradicional para contabilizar a importância das variáveis. Nos demais modelos, a importância depende da implementação de cada algoritmo.

Por exemplo, no caso de árvores de decisão, se uma variável é utilizada com maior frequência nas divisões de folhas contribuindo para maior pureza do nó, ela terá uma importância maior para o modelo final, significando que sua remoção prejudicará a performance geral do modelo. Nos modelos que utilizam múltiplas árvores, como o random forest, a importância das variáveis é calculada a cada divisão das folhas de uma árvore e somada com o resultado das demais árvores²¹21. Hastie T, Tibshirani R, Friedman J. The elements of statistical learning. 2. ed. New York: Springer Science & Business Media; 2016.. Nos modelos de redes neurais, para este estudo foi utilizado o método “weights”²⁴24. Gevrey M, Dimopoulos I, Lek S. Review and comparison of methods to study the contribution of variables in artificial neural network models. Ecol Modell. 2003;160(3):249-64. https://doi.org/10.1016/S0304-3800(02)00257-0
https://doi.org/10.1016/S0304-3800(02)00... , em que se utiliza uma combinação dos valores de saída dos neurônios das camadas ocultas e os valores dos pesos de cada conexão para calcular a maior contribuição e representatividade de cada variável de entrada no resultado da classificação. Neste estudo, foi utilizado o pacote caret²⁵25. Kuhn M, Wing J, Weston S, Williams A, Keefer C, Engelhardt A, et al. A Short Introduction to the caret package. 2020 [cited 2020 Sep 22]. Available from: https://cran.r-project.org/web/packages/caret/vignettes/caret.html
https://cran.r-project.org/web/packages/... do R, que já implementa o cálculo de importância de variáveis de acordo com cada modelo²⁶26. Kuhn M. Variable importance. 2019 [cited 2020 Aug 11]. Available from: https://topepo.github.io/caret/variable-importance.html
https://topepo.github.io/caret/variable-... .

Apesar de não ter interpretação causal²⁷27. Santos HG. Comparação da performance de algoritmos de machine learning para a análise preditiva em saúde pública e medicina [tese]. São Paulo: Faculdade de Saúde Pública da Universidade de São Paulo; 2018., as variáveis podem dar informações relevantes sobre a sua relação com o desfecho. A Figura 2 demonstra como cada variável contribuiu para os modelos desenvolvidos.

As variáveis que mais contribuíram para o melhor modelo (redes neurais artificiais) foram: quantidade de salários mínimos, cidade grande, média de alunos por turma, tempo no emprego e empresa grande.

DISCUSSÃO

A RAIS pode fornecer indícios da situação de morbidade dos trabalhadores por meio do campo de causa de afastamento. A ampla cobertura e o registro individualizado facilitam a seleção e análise de dados de uma população de trabalhadores formais com a mesma ocupação em um determinado estado, o que permite acompanhar ao longo dos anos a evolução de afastamentos por doença e predizer quais profissionais apresentam maior risco de adoecimento nos próximos anos. Para o melhor algoritmo (redes neurais artificiais), obteve-se 0,79 da área abaixo da curva ROC, sensibilidade de 72,86% e 71,52% de acurácia, usando apenas 11 variáveis preditoras.

Os resultados mostram a viabilidade da utilização de fontes anonimizadas para a predição de morbidade de docentes do ensino pré-escolar com algoritmos de machine learning. As variáveis que mais contribuíram para o melhor modelo podem ser exploradas em estudos futuros, para verificar se há relação causal. Porte da cidade em que se trabalha, renda média do trabalhador e média de alunos por turma podem ter impactos diretos na saúde dos professores. No caso da média de alunos, há evidências na literatura que corroboram este raciocínio, como problemas na acústica do ambiente e elevado ruído devido a maior agitação em sala de aula³3. Medeiros AM, Vieira MT. Ausência ao trabalho por distúrbio vocal de professores da Educação Básica no Brasil. Cad Saude Publica. 2019;35 Supl 1:e00171717. https://doi.org/10.1590/0102-311x00171717
https://doi.org/10.1590/0102-311x0017171... ,²⁸28. Rezende BA, Medeiros AM, Silva AM, Assunção AA. Fatores associados à percepção de ruído ocupacional intenso pelos professores da educação básica no Brasil. Rev Bras Epidemiol. 2019;22:e190063. https://doi.org/10.1590/1980-549720190063
https://doi.org/10.1590/1980-54972019006... .

Algoritmos de machine learning têm sido utilizados em diferentes contextos na saúde do trabalhador para auxiliar na predição de morbidades relacionadas às atividades destes profissionais²⁹29. Aliabadi M, Farhadian M, Darvishi E. Prediction of hearing loss among the noise-exposed workers in a steel factory using artificial intelligence approach. Int Arch Occup Environ Health. 2015;88(6):779-87. https://doi.org/10.1007/s00420-014-1004-z
https://doi.org/10.1007/s00420-014-1004-... ,³⁰30. Lee YC, Huang SC, Huang CH, Wu HH. A new approach to identify high burnout medical staffs by kernel K-means cluster analysis in a regional teaching hospital in Taiwan. Inquiry. 2016;53:0046958016679306. https://doi.org/10.1177/0046958016679306
https://doi.org/10.1177/0046958016679306... . Na área da educação, não foram encontrados estudos que utilizaram machine learning ou demais modelos preditivos para risco de absenteísmo docente como principal desfecho, havendo alguns estudos que analisaram a associação de variáveis preditoras com o desfecho³¹31. Ferris G, Bergin TG, Wayne SJ. Personal characteristics, job performance, and absenteeism of public school teachers. J Appl Soc Psychol. 1988;18(7):552-63. https://doi.org/10.1111/j.1559-1816.1988.tb00036.x
https://doi.org/10.1111/j.1559-1816.1988... ,³²32. Rosenblatt Z, Shirom A. Predicting teacher absenteeism by personal background factors. J Educ Adm. 2005;43(2):209-25. https://doi.org/10.1108/09578230510586597
https://doi.org/10.1108/0957823051058659... , não sendo sua predição o principal objetivo. A predição de absenteísmo docente pode identificar professores que precisam de maior atenção, permitindo adotar medidas para evitar afastamentos por doenças relacionadas à ocupação³3. Medeiros AM, Vieira MT. Ausência ao trabalho por distúrbio vocal de professores da Educação Básica no Brasil. Cad Saude Publica. 2019;35 Supl 1:e00171717. https://doi.org/10.1590/0102-311x00171717
https://doi.org/10.1590/0102-311x0017171... ,²⁸28. Rezende BA, Medeiros AM, Silva AM, Assunção AA. Fatores associados à percepção de ruído ocupacional intenso pelos professores da educação básica no Brasil. Rev Bras Epidemiol. 2019;22:e190063. https://doi.org/10.1590/1980-549720190063
https://doi.org/10.1590/1980-54972019006... , diminuindo a ausência ao trabalho e consequentemente aprimorando a qualidade do ensino³³33. Miller RT, Murnane RJ, Willett JB. Do teacher absences impact student achievement? Longitudinal evidence from one urban school district. Educ Eval Policy Anal. 2008;30(2):181-200. https://doi.org/10.3102/0162373708318019
https://doi.org/10.3102/0162373708318019... .

O estudo possui algumas limitações. Por se tratar de dados anônimos, não foi possível verificar se os docentes, que podem ter múltiplos vínculos empregatícios municipais, atuam em mais de uma escola. Também não foi possível identificar se houve um afastamento prévio por doença pré-existente que possa estar relacionado aos afastamentos seguintes. Por fim, foi analisado um número baixo de variáveis preditoras devido à disponibilidade da RAIS. Espera-se que no futuro, com uma coleta mais extensa realizada por meio de questionários, sejam possíveis resultados ainda mais robustos por parte dos algoritmos.

Em conclusão, somente com dados públicos e anônimos, foi possível obter boa sensibilidade (73%) e especificidade (71%) na identificação de risco de absenteísmo por doença, independentemente da causa da morbidade. Apesar de não ter sido criada para fins da saúde, a RAIS mostrou-se viável para análise de absenteísmo por morbidades, principalmente por conter variáveis das relações de trabalho, podendo ser explorada em outras categorias profissionais. Os modelos desenvolvidos podem futuramente ser aplicados nas demais etapas de ensino e em outras categorias ocupacionais. Algoritmos de machine learning podem contribuir com predições mais assertivas na análise de absenteísmo de professores do ensino infantil da rede pública municipal, podendo ser usados na elaboração de indicadores de morbidade ou para subsidiar políticas públicas prevencionistas direcionadas a essa categoria profissional.

Referências bibliográficas

Datas de Publicação

Histórico