Fatores associados ao nascimento pré-termo: da regressão logística à modelagem com equações estruturais

Factores asociados con el nacimiento prematuro: de la regresión logística al modelo de ecuaciones estructurales

Resumos

O estudo propôs a aplicação da modelagem com equações estruturais (MEE) para estudar variáveis associadas ao nascimento pré-termo com base em um modelo teórico analisado previamente pela regressão logística hierarquizada. Os dados foram provenientes da pesquisa observacional do tipo caso-controle populacional sobre nascidos vivos hospitalares de mães residentes em Londrina, Paraná, Brasil (junho de 2006 a março de 2007). Para a MEE foi considerada a associação de características socioeconômicas e aspectos psicossociais sobre história reprodutiva, trabalho e atividade física, intercorrências durante a gestação e características fetais. Considerou-se, ainda, a relação dessas associações sobre o desfecho nascimento pré-termo mediado pela adequação da assistência pré-natal. Foram utilizados estimadores de mínimos quadrados ajustados pela média e variância (WLSMV), para dados categóricos, e a máxima verossimilhança robusta (MLR), para obter razões de chances. Foram criadas três variáveis latentes: vulnerabilidade socioeconômica, vulnerabilidade familiar e não aceitação da gravidez. O efeito da vulnerabilidade socioeconômica, da família e da não aceitação da gravidez sobre a prematuridade ocorreu de modo indireto por meio da inadequação da assistência pré-natal. A metodologia proposta possibilitou utilizar construtos, verificar o papel de mediação da inadequação da assistência pré-natal e identificar efeitos diretos e indiretos das variáveis sobre o desfecho nascimento pré-termo.

Palavras-chave: Modelos Estatísticos; Modelos Logísticos; Fatores de Risco; Nascimento Prematuro

Este estudio propuso la aplicación de modelos de ecuaciones estructurales (SEM) para investigar las variables asociadas con el parto prematuro basándose en un modelo teórico previamente analizado mediante regresión logística jerárquica. Los datos provienen de un estudio observacional de casos y controles de base poblacional de nacidos vivos en hospitales de madres que residen en Londrina, estado de Paraná, Brasil (junio de 2006 a marzo de 2007). Para el SEM, el estudio consideró la asociación entre las características socioeconómicas y los aspectos psicosociales relacionados con el historial reproductivo, el trabajo y la actividad física, las complicaciones durante el embarazo y las características fetales. También consideró la relación entre estas asociaciones y el parto prematuro mediado por la adecuación de la atención prenatal. Se utilizó el estimador de los mínimos cuadrados ponderados ajustados por la media y variancia (WLSMV) para datos categóricos y la probabilidad máxima robusta (MLR) para los odds ratios. Se crearon tres variables latentes: vulnerabilidad socioeconómica, vulnerabilidad familiar y embarazo no deseado. El efecto de la vulnerabilidad socioeconómica y familiar y el embarazo no deseado en la prematuridad ocurrió indirectamente por la insuficiencia de la atención prenatal. La metodología propuesta permitió usar constructos, verificar el papel de la mediación por la insuficiencia de la atención prenatal e identificar los efectos directos e indirectos de las variables sobre el resultado “parto prematuro”.

Palabras-clave: Modelos Estadísticos; Modelos Logísticos; Factores de Riesgo; Nacimiento Prematuro

This study proposed the application of structural equation modeling (SEM) to investigate variables associated with preterm birth based on a theoretical model analyzed previously by hierarchical logistic regression. The data came from a population-based case-control observational study of hospital births to mothers residing in Londrina, Paraná State, Brazil (June 2006 to March 2007). For the SEM, the study considered the association between socioeconomic characteristics and psychosocial aspects pertaining to reproductive history, work and physical activity, complications during the pregnancy, and fetal characteristics. It also considered the relationship between these associations and the outcome preterm birth mediated by adequacy of prenatal care. The weighted least square mean and variance adjusted estimator (WLSMV) was used for categorical data and robust maximum likelihood (MLR) for odds ratios. Three latent variables were created: socioeconomic vulnerability, family vulnerability, and unwanted pregnancy. The effect of socioeconomic and family vulnerability and unwanted pregnancy on prematurity occurred indirectly through inadequacy of prenatal care. The proposed methodology allowed using constructs, verifying the role of mediation by inadequacy of prenatal care, and identifying the variables’ direct and indirect effects on the outcome preterm birth.

Keywords: Statistical Models; Logistic Models; Risk Factors; Premature Birth

Introdução

Modelos teóricos para estudar nascimento pré-termo mostram a importância da multiplicidade de variáveis associadas ao nascimento pré-termo. A presença de intercorrências, que incluem a presença de sangramento, pressão alta, eclampsia, alteração do volume do líquido amniótico, infecção do trato genital e diabetes, tem efeito direto sobre a idade gestacional ¹. Diversos autores apontaram a gestação múltipla como fator de risco para nascimento pré-termo, inclusive com a presença de intercorrências, fato este que pode estar relacionado à ocorrência de rotura prematura de membranas ou outros agravos materno-fetais ^1,2,3,4. Idade materna, história reprodutiva envolvendo número de gestações e características dos filhos anteriormente gerados ou não (baixo peso e pré-termo) são variáveis apontadas em diversos estudos como risco para o nascimento pré-termo ^1,5,6,7. Aspectos de inadequação do cuidado realizado no pré-natal, como início tardio (após o terceiro trimestre) e frequência insatisfatória (menos de 6 consultas), constituem riscos para o nascimento pré-termo ^5,6,8.

Estudos mostraram ainda a importância da assistência pré-natal relacionada ao nascimento pré-termo. Alguns trabalhos também apontam o papel de variáveis socioeconômicas, intercorrências e história reprodutiva sobre o uso inadequado da assistência pré-natal e ressaltam o papel da assistência pré-natal como um determinante nos indicadores relacionados à mãe e ao bebê, relacionando aspectos como classe econômica, baixa escolaridade, uso de álcool na gestação ^8,9,10,11, baixa idade materna, estado marital, planejamento da gravidez, serviço público de saúde, alta paridade, parto prematuro anterior ^9,10,12 e gravidez não desejada ¹³. Esses determinantes estão na cadeia de fatores associados tanto ao pré-natal inadequado quanto a desfechos negativos como o nascimento pré-termo. Assim, consideramos o papel mediador da assistência pré-natal na modelagem.

A situação de vulnerabilidade das famílias está relacionada à situação de pobreza, exploração e abuso, aspectos psicossociais e culturais. Outro aspecto em paralelo é que as famílias chefiadas por mulheres são em parte decorrentes de uma gravidez precoce ou indesejada, instabilidade familiar e abandono ¹⁴. A vulnerabilidade social que as mulheres enfrentam pode estar também associada a um grau de vulnerabilidade emocional, que envolve sentimentos de abandono, violência e fragilização que não dê condições para fazer um pré-natal adequado e ter acesso a medicamentos ¹⁵. O não planejar a gravidez está associado a fumar durante a gestação, usar menos ácido fólico, iniciar pré-natal depois do terceiro trimestre e interromper a gestação ¹⁶.

Em estudos epidemiológicos sobre nascimento pré-termo é usual utilizar a abordagem considerando-se modelos hierarquizados com regressão logística, a fim de definir os fatores de risco para o nascimento pré-termo que consideram modelos conceituais para representar relações de uma variável sobre as outras em um modelo de causalidade ^17,18.

Nesses modelos, as variáveis são incluídas simultaneamente no modelo e podem ter seus efeitos superestimados ou subestimados. Ao serem excluídas, podem levar consigo parte da informação que faz a ligação entre as variáveis. Apesar de a regressão hierarquizada considerar elementos importantes como a ordenação temporal e a lógica da relação entre as variáveis, tais modelos não consideram adequadamente as relações de variáveis de confusão e de mediação ¹⁹.

Os modelos de equações estruturais (MEE) são mais apropriados, pois permite que múltiplas equações simultâneas incorporem a confusão e a mediação, além de incorporar variáveis latentes para representar medidas mais complexas, que não são mensuráveis com uma única variável e que são criadas com base em covariâncias entre variáveis observadas. Os MEE minimizam o efeito de confusão residual nas associações, especialmente em estudos observacionais ^20,21.

Este estudo propõe analisar as relações das variáveis com o desfecho nascimento pré-termo via MEE, utilizando um modelo testado com base na regressão logística com seleção de variáveis de forma hierarquizada. Para tanto, utilizou-se a pesquisa realizada em Londrina, Paraná, Brasil, em 2009, conduzida por Silva et al. ². Baseando-se no quadro teórico conceitual do trabalho de origem foi redefinido um modelo (Figura 1) que incorporou variáveis latentes. Além disso, o MEE possibilita a inclusão de variáveis que têm efeito de mediação entre as variáveis de exposição e o desfecho.

Figura 1
Quadro teórico-conceitual de fatores associados a nascimentos pré-termo.

Métodos

Os dados referem-se a estudo do tipo observacional de caso-controle de base populacional de nascimentos hospitalares realizado em Londrina (junho de 2006 a março de 2007), em que os casos eram nascimentos pré-termo (< 37 semanas de gestação) e os controles eram nascimentos não pré-termo (37 semanas e mais). O estudo original de Silva et al. ² utilizou todas as variáveis na forma categórica. Neste trabalho, optou-se por utilizar as variáveis na forma contínua visando a minimizar o fato de que variáveis categóricas não apresentam distribuição aproximada normal e podem apresentar problemas de assimetria que afetam os erros-padrão, as estimativas da variância residual e a estatística qui-quadrado usada nas medidas de ajuste do MEE. Os registros que tinham as variáveis contínuas completas totalizaram 296 casos e 329 controles, representando 90% e 89%, respectivamente, da base original.

A fim de detectar possíveis vieses ocasionados devido à diferença no número de registros do banco de dados original foi empregado o teste de independência, baseado no qui-quadrado de Pearson ²², que mostrou que as distribuições das variáveis presentes nas duas bases são estatisticamente semelhantes, indicando não haver vieses de seleção de eventos.

O quadro teórico de Silva et al. ² estava organizado em blocos e, no presente trabalho, foram transcritos para representar possíveis constructos (Figura 1). O bloco 1 do modelo original está representado como características socioeconômicas maternas e familiares. O bloco 2 foi separado em características pré-concepcionais e história reprodutiva materna. O bloco 3 foi separado em condições psicossociais, hábitos maternos contemplando atividade física e de trabalho e assistência na gestação (pré-natal). O bloco 4, intercorrências maternas durante a gestação e, por fim, o bloco 5: características fetais foram mantidas.

Baseando-se nesse modelo, os passos adotados para fazer a modelagem foram:

(1) Modelo de mensuração: verificação da composição dos fatores via análise fatorial confirmatória (AFC). Criação de variáveis latentes: com base nos constructos elaborados sugeridos pelos blocos de variáveis e revisão da literatura que definiram o modelo teórico, por meio da AFC, considerando-se o número de fatores cujo autovalor é superior a 1 ²³.

(2) Modelo estrutural: relaciona variáveis latentes, observadas e o desfecho nascimento pré-termo, valendo-se dos modelos em que este desfecho é binário. Foram obtidos efeitos diretos e indiretos sobre o desfecho. Foram utilizados índices de modificação para identificar relações não consideradas anteriormente e melhorar o ajuste do modelo ²⁴.

Foram usados estimadores de mínimos quadrados ajustados pela média e variância (WLSMV) com rotação geomin²⁵. Após a obtenção do modelo final e testados os ajustes do modelo, este foi reestimado por máxima verossimilhança robusta (MLR) a fim de obter estimativas de razão de chances (OR). As medidas de ajuste do modelo utilizadas foram: Tucker-Lewis (TLI; referência de bom ajuste: TLI > 0,95), índice de ajustamento comparativo (CFI; referência: > 0,95), raiz do erro quadrático médio de aproximação (RMSEA; referência: < 0,05) e raiz do resíduo quadrático médio ponderado (WRMR; referência: < 0,90). Foi considerado o nível de significância de 5%. Para as análises foi usado o software MPlus (https://www.statmodel.com/) e, para os gráficos, o CmapTools (https://cmap.ihmc.us/).

Características socioeconômicas maternas e familiares

Essas características contaram com as seguintes variáveis: idade da mãe em anos (< 19; 20 a 34; 35 e mais), renda per capita familiar em salários mínimos, número de pessoas por cômodo, anos de estudos da mãe, anos de estudos do responsável pelo sustento familiar, ocupação do responsável pelo sustento familiar (qualificada; manual ou semiqualificada; do lar ou estudante; desempregado), tipo de ocupação do responsável pelo sustento familiar (formal; informal; não inserido), localidade de moradia (não favela; favela), tipo de família (nuclear - composta por marido, mulher e filhos; apenas mulher e filhos; não nuclear; outros arranjos familiares), mora com companheiro há menos de 2 anos (sim; não), responsável pelo sustento familiar (mãe; pai; outros), presença de idosos maiores de 60 anos (sim; não), raça/cor da mãe (branca/amarela/indígena; preta/parda), tipo de acabamento da moradia (completo; incompleto), migração materna (sim; não), número de crianças menores de 10 anos (sim; não) e idade do responsável pelo sustento (em anos).

Aspectos psicossociais

As condições da gestação abrangendo questões psicossociais no início da gravidez e durante a gestação foram: tentativa de aborto (sim; não), gravidez planejada (sim; não), reações da mãe, pai da criança e família com relação à gravidez e preocupações durante a gestação (negativa; positiva), tentativa de aborto (sim; não), brigas com o companheiro (sim; não), separação (sim; não) e preocupações durante a gestação (sim; não).

Intercorrências maternas na gestação

Como intercorrências foram consideradas: sangramento vaginal, hipertensão, eclampsia, alteração do volume amniótico, internação durante a gravidez, anemia, diabetes, infecção genital, corrimento, placenta prévia e infecção urinária (sim; não).

Condição reprodutiva materna

As variáveis utilizadas para medir as características reprodutivas maternas foram: número de gestações, cesáreas anteriores, filho anterior pré-termo, filho anterior baixo peso, idade materna e tratamento para engravidar (sim; não).

Trabalho e atividade física

As variáveis usadas para medir as atividades foram: se a mãe trabalhou durante a gestação, se o trabalho exigia esforços físicos, e se realizava atividade em casa com esforço.

Assistência pré-natal

A inadequação do pré-natal envolveu as seguintes variáveis observadas: início (1º trimestre), número de consultas (mínimo 3), realização de exames (urina, sangue e ultrassonografia), procedimentos e orientações básicas. As categorias foram: pré-natal adequado (todos os itens positivos), inadequado I (uma ou duas respostas negativas), inadequado II (três ou mais respostas negativas) ou não fez pré-natal. Para a modelagem, foram consideradas duas variáveis indicadoras para representar a inadequação ao pré-natal: (i) inadequação ao pré-natal I - se o pré-natal foi Inadequado (recebeu 1) e se foi Adequado (recebeu 0) e (ii) inadequação ao pré-natal II - se o pré-natal foi inadequado do tipo II, que inclui a mãe que não fez pré-natal (recebeu 1) e se foi inadequado do tipo I ou adequado (recebeu 0).

Outras variáveis

A fim de contemplar todas as variáveis selecionadas pelo estudo original ², foram incluídas como variáveis observadas: índice de massa corporal (IMC), prática de atividade física, uso de bebida alcoólica e gravidez múltipla.

Aspectos éticos

O estudo seguiu os preceitos éticos estabelecidos na Resolução nº 196/1996 do Conselho Nacional de Saúde, que regulamenta pesquisas envolvendo seres humanos, sendo aprovado pelo Comitê de Ética da Faculdade de Saúde Pública da Universidade de São Paulo (protocolo nº 404211/2013).

Resultados

Durante a modelagem, alterações foram feitas com base no modelo inicial. Baseando-se nos blocos Características Socioeconômicas Maternas e Familiares e Aspectos Psicossociais foram geradas as variáveis latentes descritas a seguir.

Características socioeconômicas maternas e familiares

Foram selecionados dois fatores cuja porcentagem de variância explicada acumularam 71,5%, com os ajustes do modelo adequados (RMSEA = 0,053; CFI = 0,988; TLI = 0,977) (Tabela 1). Foram excluídas as variáveis cuja correlação com os fatores foi abaixo de 0,30: raça/cor da mãe, tipo de acabamento da moradia, migração materna, número de crianças menores de 10 anos, idade da mãe, idade do responsável pelo sustento e tipo de ocupação do responsável pelo sustento familiar. Esses dois fatores representam, respectivamente, a vulnerabilidade socioeconômica (VSE) e a vulnerabilidade familiar (VF). Com base na forma como foram categorizadas as variáveis, escores mais elevados da variável latente VSE correspondem a menores níveis de escolaridade da mãe e do responsável pela família e maior frequência de moradia em favela, já a renda familiar per capita tem sinal negativo indicando que domicílio com menor valor de renda apresenta maior vulnerabilidade. Já os escores mais elevados da vulnerabilidade familiar correspondem às relações recentes com o companheiro, família não nuclear e presença de idoso.

Tabela 1
Resultado da análise fatorial com características socioeconômicas. Londrina, Paraná, Brasil, 2009.

Variáveis	Fator 1	Fator 2
Baixa escolaridade do responsável pelo sustento familiar (até 8 anos de estudos)	0,987	-0,019
Baixa escolaridade materna (até 8 anos de estudos)	0,795	-0,022
Renda per capita (salários mínimos)	-0,790	-0,088
Número de pessoas por cômodos	0,735	0,056
Tipo de moradia (favela)	0,648	0,011
Família não nuclear	0,158	0,930
Mãe responsável pelo sustento familiar	-0,069	0,812
Presença de idosos	0,005	0,792
Mãe mora com companheiro há menos de 2 anos	-0,103	0,692

Aspectos psicossociais

As condições psicológicas maternas afetadas por momentos de preocupações e estresse durante a gestação são dimensões que esta variável latente tem o objetivo de medir (Tabela 2). As variáveis tentativa de aborto, brigas com companheiro, separação e preocupações durante a gestação foram excluídas devido à baixa correlação (< 0,30) com as demais variáveis. As variáveis selecionadas foram: se houve um planejamento da gravidez e as reações à gravidez da mãe, pai da criança e da família. Foi considerado um único fator cuja porcentagem da variância explicada acumulou 75,3% com valores de ajustes adequados (RMSEA = 0,056; CFI = 0,991; TLI = 0,972).

Com base na categorização das variáveis, escores mais elevados da variável latente não aceitação da gravidez (NAG) correspondem às reações negativas da gravidez, tanto da própria mãe, do pai do bebê e da família. Esse fator também indica a vulnerabilidade emocional quanto ao recebimento de apoio ao engravidar.

Tabela 2
Cargas fatoriais rotacionadas (método geomin) - não aceitação da gravidez. Londrina, Paraná, Brasil, 2009.

Variáveis	Fator 1
Gravidez planejada	0,911
Reação negativa ou de preocupação da mãe à gravidez	0,969
Reação negativa ou de preocupação do pai da criança à gravidez	0,774
Reação negativa ou de preocupação da família	0,628

Condição reprodutiva materna

Este bloco foi representado pela combinação do número de gestações, filho de baixo peso, filho anterior pré-termo e idade materna associados com o desfecho de forma gradual com suas razões de chance e representa a interação destas variáveis (Tabela 3).

Tabela 3
Distribuição de casos e controles segundo variáveis selecionadas. Londrina, Paraná, Brasil, 2009.

Variáveis	Caso		Controle		Total
	n	%	n	%	n	%
Condição reprodutiva materna
a) Uma ou + gestações, não prematuro anterior e idade entre 20 e 34 anos	86	29,1	129	39,2	215	34,4
b) Nulípara, idade entre 13 e 19 anos	35	11,8	44	13,4	79	12,6
c) Uma ou + gestações, não prematuro anterior e 35 anos ou mais	22	7,4	38	11,6	60	9,6
d) Nulípara, idade entre 20 e 34 anos	76	25,7	83	25,2	159	25,4
e) Uma ou + gestações, não prematuro anterior, não baixo peso anterior e idade entre 13 e 19 anos	11	3,7	10	3,0	21	3,4
f) Nulípara, idade 35 anos ou mais	9	3,0	2	0,6	11	1,8
g) Uma ou + gestações, prematuro anterior, baixo peso anterior	57	19,3	23	7,0	80	12,8
Atividades comportamentais da mãe: esforço físico
a) Mãe não trabalhou fora, sem atividades de esforço em casa	77	26,0	77	23,4	154	24,6
b) Mãe não trabalhou fora, com atividades de esforço em casa	55	18,6	57	17,3	112	17,9
c) Mãe trabalhou fora, sem esforços físicos, sem atividades de esforço em casa	94	31,8	92	28,0	186	29,8
d) Mãe trabalhou fora, com esforços físicos, sem atividades de esforço em casa	15	5,1	22	6,7	37	5,9
e) Mãe trabalhou fora, sem esforços físicos, com atividades de esforço em casa	28	9,5	47	14,3	75	12,0
f) Mãe trabalhou fora, com esforços físicos, com atividades de esforço em casa	27	9,1	34	10,3	61	9,8

Intercorrências maternas na gestação

No modelo de regressão logística do estudo original, o bloco relacionado às intercorrências maternas na gestação apresentou um importante ajuste na análise de regressão múltipla ², sugerindo que mais de um destes agravos poderiam estar presentes, bem como alguns deles podem ter levado à internação da mãe durante a gestação. Assim, a variável para representar o número de intercorrências (0; 1 ou mais) (Tabela 3). Ainda, foi possível verificar a presença de internação mais frequente na medida em que aumentava a quantidade de intercorrências - das mães com pelo menos 1 intercorrência, 40% foram internadas.

Trabalho e atividade física

O esforço físico foi composto com as variáveis: se a mãe trabalhou durante a gestação, se o trabalho exigia esforços físicos, e se realizava atividade em casa com esforço (Tabela 3).

Modelo final

O modelo completo foi desenvolvido baseando-se na construção das variáveis latentes e as variáveis observadas (Tabela 4).

Na Tabela 4, são apresentadas as medidas dos coeficientes padronizados e das razões de chance para as relações do modelo estrutural. O ajuste do modelo foi satisfatório (CFI = 0,962; TLI = 0,952; WRMR = 1,005; RMSEA (IC95%) = 0,032 (0,026-0,039)).

Tabela 4
Estimativas dos parâmetros dos modelos para nascimento pré-termo. Resultados utilizando-se os estimadores mínimos quadrados ajustados pela média e variância (WLSMV) e máxima verossimilhança robusta (MLR). Londrina, Paraná, Brasil, 2009.

Variáveis e indicadores dos fatores	WLSMV				MLR
	Estimativa padronizada	Estimativa	EP	Valor de p	Estimativa	EP	Valor de p	OR
Modelo de mensuração
VSE
Pessoas por cômodos	0,652	1,000	0,000	-	1,000	0,000	-	-
Renda per capita	-0,877	-2,483	0,134	0,000	-2,412	0,170	0,000	-
Escolaridade da mãe	0,841	4,798	0,476	0,000	7,274	0,843	0,000	-
Escolaridade do pai	0,741	3,403	0,286	0,000	5,872	0,744	0,000	-
Tipo de moradia	0,661	2,716	0,427	0,000	5,875	0,924	0,000	-
VF
Companheiro há menos de 2 anos	0,669	1,000	0,000	-	1,000		0,000	-		-
Família não nuclear	0,994	9,722	30,591	0,751	2,877		1,204	0,017		-
Presença de idoso	0,768	1,332	0,256	0,000	1,812		0,414	0,000		-
Relação mãe/responsável	0,776	1,367	0,218	0,000	1,503		0,245	0,000		-
NAG
Reação da mãe quando descobriu a gravidez	0,785	1,000	0,000	-	1,000		0,000	-		-
Reação do pai quando soube da gravidez	0,825	1,150	0,302	0,000	1,239		0,278	0,000		-
Reação da família quando descobriu	0,825	1,150	0,272	0,000	1,003		0,253	0,000		-
Modelo estrutural
NAG =
VF	0,365	0,515	0,129	0,000	0,561		0,133	0,000		-
VSE	0,267	1,043	0,241	0,000	2,225		0,524	0,000		-
Condição reprodutiva materna	0,025	0,016	0,032	0,616	0,055		0,058	0,343		-
Frequência de bebida alcoólica na gestação	0,166	0,983	0,368	0,007	1,265		0,590	0,032		-
Intercorrências =
Condição reprodutiva materna	0,039	0,021	0,028	0,454	0,052		0,042	0,215		1,053
Esforço físico	0,067	0,044	0,034	0,203	0,046		0,053	0,380		1,047
Caminhadas	-0,020	-0,047	0,133	0,726	-0,267		0,199	0,179		0,765
IMC	0,099	0,155	0,076	0,041	0,295		0,127	0,021		1,343
Frequência de bebida alcoólica na gestação	0,041	0,198	0,222	0,372	0,439		0,411	0,285		1,552
Gestação múltipla	0,224	0,972	0,228	0,000	1,621		0,392	0,000		5,057
Preocupações	0,107	0,253	0,117	0,031	0,426		0,191	0,026		1,531
Inadequação do pré-natal I =
VSE	0,091	0,323	0,219	0,140	0,713		0,384	0,063		2,041
VF	0,088	0,112	0,099	0,256	0,163		0,116	0,160		1,177
NAG	0,347	0,314	0,094	0,001	0,199		0,093	0,032		1,221
Condição reprodutiva materna	0,149	0,086	0,030	0,005	0,100		0,046	0,028		1,106
Esforço físico	-0,092	-0,066	0,038	0,082	-0,050		0,060	0,399		0,951
Inadequação do pré-natal II =
VSE	0,096	0,367	0,395	0,353	1,258		0,878	0,152		3,519
VF	-0,006	-0,008	0,175	0,963	0,105		0,200	0,601		1,110
NAG	0,467	0,457	0,169	0,007	0,326		0,174	0,061		1,385
Condição reprodutiva materna	0,188	0,118	0,046	0,001	0,127		0,083	0,127		1,158
Esforço físico	-0,159	-0,124	0,059	0,034	-0,154		0,130	0,237		0,857
Idade gestacional =
VSE	-0,083	-0,551	0,611	0,367	0,007		0,383	0,986		1,007
VF	0,073	0,175	0,233	0,452	0,022		0,102	0,827		1,023
NAG	-0,420	-0,714	0,462	0,122	-0,045		0,076	0,556		0,956
Condição reprodutiva materna	0,091	0,098	0,072	0,172	0,213		0,047	0,000		1,237
Inadequação pré-natal I	0,226	0,426	0,227	0,051	0,457		0,224	0,041		1,579
Inadequação pré-natal II	0,555	0,965	0,539	0,044	1,457		0,500	0,004		4,293
Intercorrências	0,584	1,209	0,489	0,013	1,825		0,204	0,000		6,201
Gestação múltipla	0,214	1,924	0,739	0,009	2,876		0,773	0,000		17,747
Esforço físico	-0,034	-0,046	0,090	0,612	-0,146		0,064	0,023		0,864
Preocupações	0,040	0,198	0,227	0,384	0,424		0,218	0,052		1,527
Caminhadas	-0,133	-0,641	0,300	0,032	-0,606		0,235	0,010		0,545
IMC	0,146	0,473	0,188	0,012	0,449		0,157	0,004		1,567
Frequência de bebida alcoólica na gestação	0,114	1,147	0,555	0,039	1,131		0,393	0,004		3,099
Correlações
VF com condição reprodutiva materna	0,065	0,116	0,096	0,225	0,226		0,148	0,127		-
VSE com condição reprodutiva materna	0,170	0,109	0,022	0,000	0,100		0,030	0,001		-
VSE com esforço físico	-0,181	-0,093	0,022	0,000	-0,096		0,021	0,000		-
VSE com VF	0,175	0,051	0,016	0,001	0,079		0,027	0,003		-
VSE com IMC	0,147	0,032	0,009	0,001	0,025		0,010	0,014		-
VSE com caminhadas	-0,228	-0,033	0,015	0,033	-0,033		0,007	0,000		-
IMC com caminhadas	-0,043	-0,013	0,013	0,345	-0,011		0,011	0,024		-
Interceptos
Pessoas por cômodos	1,456	0,724	0,154	0,000	0,852		0,020	0,000		-
Renda per capita	0,093	0,085	0,217	0,693	0,037		0,037	0,931		-

• EP: erro-padrão; IMC: índice de massa corproal; NAG: não aceitação da gravidez; OR: razão de chances; VF: vulnerabilidade familiar; VSE: vulnerabilidade socioeconômica.

Na representação gráfica do modelo final (Figura 2) são apresentadas as estimativas padronizadas. Para algumas relações são apresentadas as estimativas padronizadas e a OR. Os dois valores para algumas variáveis indicam os efeitos para inadequação pré-natal I e II.

Figura 2
Modelo resultante para o nascimento pré-termo.

Os maiores efeitos diretos sobre o nascimento pré-termo foram: intercorrências (0,584), seguido da atenção de pré-natal inadequado II (0,555) e atenção de pré-natal inadequado I (0,226), gestação múltipla (0,214), condição reprodutiva (0,091), IMC (0,146), consumo de bebida alcoólica (0,114), caminhadas (-0,133) e esforço (-0,034). Além de apresentar um efeito direto (0,214), a gravidez múltipla registrou também um efeito indireto sobre o desfecho via intercorrências (0,224). Ao se considerar o efeito total de gravidez múltipla com o efeito direto e via intercorrências seu efeito (0,382) é mais influente.

O IMC durante a gestação também apresentou efeito indireto via intercorrências (0,099), assim, seu efeito total sobre a prematuridade é de 0,205. A frequência de bebida alcoólica, além do efeito direto, apresentou um efeito indireto via NAG e inadequação do pré-natal, o efeito total foi de 0,126. Ter preocupações na gravidez registrou somente um efeito indireto (0,107) sobre a prematuridade via intercorrência.

As variáveis latentes NAG e VSE apresentaram efeitos diretos sobre a assistência pré-natal inadequada. As variáveis latentes registraram sobre o nascimento pré-termo os efeitos indiretos: NAG (0,228), VSE (0,139) e VF (0,123).

Discussão

A atenção pré-natal, no modelo atual, apresenta efeito direito com sinal positivo sobre o nascimento pré-termo, como apontaram outros estudos ^5,6. É no pré-natal que tanto a gestante quanto o bebê têm um acompanhamento preventivo, de orientação, de esclarecimento e de diagnóstico de qualquer alteração da saúde da gestante e/ou do bebê ²⁶. O efeito da VSE, da VF, da condição reprodutiva e da NAG sobre a prematuridade ocorre de modo indireto por meio da inadequação da assistência pré-natal, indicando assim o papel de mediação desta variável sobre a prematuridade. A assistência pré-natal tem um papel mediador importante sobre desfechos negativos da gestação e reflete condições sociais, econômicas e psicológicas das mães. Em uma análise ampla de indicadores de qualidade da assistência pré-natal no Brasil, 15% das mães recebem orientação adequada de pré-natal e 60% recebem todas as orientações e exames complementares ¹⁵. Suas condições de vida e as vulnerabilidades emocional e afetiva podem influenciar o conhecimento sobre saúde durante a gravidez ^8,13,27,28. A não adequação do pré-natal foi observada em mulheres jovens, multíparas, sem parceiro, sem emprego pago, de pouca escolaridade, de mais baixas classes socioeconômicas ²⁹. Com base em nossos resultados, a relação entre esses fatores sugere a existência de fatores de vulnerabilidade socioeconômica, familiar e psicológica, e de condições maternas sobre a adequação do pré-natal.

A presença de intercorrências, que representam a presença de sangramento, pressão alta, eclampsia, alteração do volume do líquido amniótico, infecção do trato genital e diabetes, tem um efeito direto sobre o nascimento pré-termo ¹. Além de efeito direto, a gestação múltipla apresentou também um efeito indireto via intercorrências, fato este que pode estar relacionado à ocorrência de rotura prematura de membranas ou outros agravos materno-fetais ^1,2,5. A variável condições reprodutivas maternas apresentou efeito direto sobre o desfecho. Neste trabalho, essa variável foi composta pela idade materna e sua respectiva história reprodutiva com relação ao número de gestações e características dos filhos anteriormente gerados ou não (baixo peso e pré-termo), estes descritos ^1,5,6,7.

A separação dos fatores socioeconômicos e psicossociais em três fatores auxiliou a compreensão dos vários papéis desta complexa dimensão. Em complemento às características psicossociais, a variável “preocupações” traz a natureza de estresse vivido durante a gravidez que foi estatisticamente significante no estudo original. No presente estudo não foi encontrada associação direta sobre idade gestacional, mas um efeito indireto, via intercorrências, que aponta outra natureza de preocupações, ou seja, a saúde ^1,2. A VSE, sendo uma variável latente distal, apresenta correlações com variáveis tanto distais como proximais, como esforço físico, IMC, caminhadas, condição reprodutiva materna e VF. A VSE está relacionada ao desfecho de forma indireta no presente estudo por dois caminhos. O primeiro, via inadequação do pré-natal, e o outro via NAG. A dimensão socioeconômica é bastante complexa, não se resume somente às privações de necessidades materiais de bem-estar, mas também à negação de oportunidades nas relações sociais, como de acesso ao trabalho e à assistência à saúde ^30,31. Tanto a inadequação pré-natal quanto a não aceitação da gestação capturam tais vulnerabilidades.

Com relação ao esforço fisico, o resultado aqui apresentado apontou uma associação negativa com pré-termo, diferentemente do que se observa em outros estudos, incluindo-se uma revisão sistemática na qual não se observou efeito na maioria dos estudos ou efeito moderado com as condições de trabalho com esforço físico. A variável esforço físico, da forma como foi construída, pode traduzir o estágio de vida da mãe com relação à sua plena atividade produtiva/trabalho, quer seja pela sua idade, quer seja pela renda per capita domiciliar, apontando não VSE (correlação com a variável VSE verificada no modelo). Por outro lado, Behrman & Butler ²⁸ sinalizam que o trabalho com esforço pode ser entendido também como um indicador de circunstâncias socioeconômicas favoráveis, isto é, a capacidade de ter e manter um emprego, os benefícios que recebem como também sinaliza a satisfação psicológica que alguns tipos de trabalho apresentam. Assim, o resultado que foi obtido no presente trabalho não é claramente explicado e outras abordagens devem ser realizadas para se entender melhor o papel do esforço. O mesmo comportamento observou-se para a variável prática de caminhadas ou outra atividade física, isto é, efeito direto e negativo com pré-termo.

O IMC é uma variável que apresentou efeito direto sobre a idade gestacional e efeito indireto via intercorrência, como já apontado em Padilha et al. ³². No presente trabalho também foi encontrada correlação entre IMC e VSE, indicando valores de IMC fora dos níveis referendados - mães com IMC < 19kg/m² (magras) e IMC ≥ 30kg/m² (obesas) em mães com maior VSE.

Esse mesmo comportamento foi observado para a frequência de ingestão de bebida alcoólica durante a gravidez, apresentando efeito direto e efeito indireto sobre a idade gestacional. O efeito indireto foi via NAG e inadequação do pré-natal. Alguns estudos também encontraram associação com o alto consumo de bebida alcoólica e a idade gestacional, porém os mecanismos não são claros ^1,28,33, e outros não encontram associação ³⁴.

O modelo atual também indicou um caminho sobre idade gestacional via NAG e inadequação do pré-natal, e não foram encontradas referências que discutem estas possíveis relações. Outros trabalhos poderiam explorar melhor esses caminhos. Vale ressaltar novamente o papel do pré-natal como canal de cuidado da mãe e do bebê.

O estudo de caso-controle original ² empregou um quadro teórico delineado em blocos hierarquizados com o uso do modelo de regressão logística ³⁵, no qual foram definidos blocos de variáveis em que estas são agrupadas por apresentarem alguma característica comum. Este estudo teve como base esse modelo ao se considerar a estrutura criada incluindo-se variáveis não observadas (latentes). A priori, criar um modelo de equações estruturais valendo-se do modelo hierarquizado é possível, pois a estrutura em blocos guarda um elemento importante das relações dessas variáveis com o desfecho: a temporalidade, caracterizada com os conceitos “distal”, “intermediário” e “proximal”.

No entanto, dificuldades para reescrever o modelo ocorreram: os blocos agrupam variáveis que nem sempre têm correlações entre elas, de forma que cada relação deve ser reavaliada para a inclusão no modelo.

Além das baixas correlações entre as variáveis que compõem os fatores, existe uma série de outros problemas que podem ser enfrentados ao se utilizar o MEE. Problemas de estimação - não convergência, falta de identificação, matriz não positiva definida e estimação de variâncias negativas - não são incomuns ²⁶.

Para algumas variáveis, problemas na estimação ocorreram por conta de caselas vazias (zero) nas combinações bivariadas como, por exemplo, entre “presença de gestação anterior” e “mães nulíparas”, em que estruturalmente há um zero neste cruzamento, o que impediu a estimação da covariância e fez-se necessário usar a solução proposta por Muthén & Muthén ³⁶ que é combinar uma nova variável, isto é, criar uma variável de interação. Levou-se em consideração sua relação com o desfecho (via OR). Uma alternativa é a utilização de uma perturbação na casela, ou seja, troca-se o zero por um valor pequeno ³⁷.

A estimação de OR possibilita a interpretação conhecida na área de epidemiologia. Obter o modelo de OR depois de ter gerado o modelo com WLSMV foi a alternativa encontrada para se fazer a estimação mais adequada aos dados ³⁸. Muthén et al. ³⁹ relataram que WLSMV foi desenvolvido baseando-se em amostras pequenas e moderadas. Beauducel & Herzberg ⁴⁰ também reportaram resultados mais satisfatórios com WLSMV em variáveis com poucas categorias (2 ou 3, comparada a 4 ou mais).

Dentre as limitações há o fato da pesquisa ser do tipo observacional, que não nos permite falar diretamente em causalidade das relações; o número de indivíduos é pequeno frente ao número de relações estudadas, embora isto não seja problema para se obter os ajustes dos modelos; o fato de que o estudo não foi delineado para observar as relações na forma proposta pela modelagem.

Do ponto de vista metodológico, a construção do MEE realizada com base em um modelo já construído (original), valendo-se do modelo de regressão hierarquizada, apresentou vantagens na medida em que pode ser utilizado como um modelo inicial e permitiu incluir caminhos que representam as relações entre as variáveis do modelo (entre elas e o desfecho), além de identificar mediações propostas. As alterações propostas com base em tal modelo inicial visaram a apreender as variáveis envolvidas no trabalho original. O processo de construção do MEE ressaltou: (1) a necessidade de se rediscutir o papel das variáveis quando se escolhe utilizar o modelo de equações estruturais; (2) que modelos complexos requerem estudos que já considerem a criação de variáveis latentes em sua formulação; (3) que a combinação de variáveis para formar uma nova variável é um recurso quando não se tem variáveis latentes bem definidas a priori e que resolvem problemas de caselas vazias nos cruzamentos de variáveis; e (4) as mesmas variáveis foram significantes utilizando-se os diferentes estimadores (WLSMV e MLR) na parte estrutural do modelo.

Conclusão

Foi possível utilizar os resultados do trabalho original para fazer a modelagem com equações estruturais com base na revisão do modelo teórico conhecendo-se as relações com o desfecho e o modelo de regressão logística. A aplicação da metodologia proposta identificou a presença de constructos (VSE, NAG, VF), verificar que a inadequação do pré-natal como mediadora, além de identificar efeitos diretos e indiretos de variáveis sobre o desfecho nascimento pré-termo.

Agradecimentos

À Fundação de Amparo á Pesquisa do Estado de São Paulo (Fapesp; processo nº 2010/20784-8) pelo financiamento.

Referências

• ¹ Goldenberg RL, Culhane JF, Iams JD, Romero R. Epidemiology and causes of preterm birth. Lancet 2008; 371:75-84.
• ² Silva AMR, Almeida MF, Matsuo T, Soares DA. Fatores de risco para nascimentos pré-termo em Londrina, Paraná, Brasil. Cad Saúde Pública 2009; 25:2125-38.
• ³ Hair Jr. JF, Anderson RE, Tatham RL, Black WC. Análise multivariada de dados. 6th Ed. Porto Alegre: Bookman; 2009.
• ⁴ Assunção PL, Novaes HMD, Alencar GP, Melo ASO, Almeida MF. Fatores associados ao nascimento pré-termo em Campina Grande, Paraíba, Brasil: um estudo caso-controle. Cad Saúde Pública 2012; 28:1078-90.
• ⁵ Silveira MF, Victora CG, Barros AJD, Santos IS, Matijasevich A, Barros FC. Determinants of preterm birth: Pelotas, Rio Grande do Sul State, Brazil, 2004 birth cohort. Cad Saúde Pública 2010; 26:185-94.
• ⁶ Aragão VMF, Silva AAM, Aragão LF, Barbieri MA, Bettiol H, Coimbra LC, et al. Risk factors for preterm births in São Luís, Maranhão, Brasil. Cad Saúde Pública 2004; 20:57-63.
• ⁷ Bezerra LC, Oliveira SMJV, Latorre MRDO. Prevalência e fatores associados à prematuridade entre gestantes submetidas à inibição de trabalho de parto prematuro. Rev Bras Saúde Matern Infant 2006; 6:223-9.
• ⁸ Goudard MJF, Simões VMF, Batista RFL, Queiroz RCS, Alves MTSSB, Coimbra LC, et al. Inadequação do conteúdo da assistência pré-natal e fatores associados em uma coorte no nordeste brasileiro. Ciênc Saúde Coletiva 2016; 21:1227-38.
• ⁹ Feijen-de Jong EI, Jansen DEMC, Baarveld F, van der Schans CP, Schellevis FG, Reijneveld SA. Determinants of late and/or inadequate use of prenatal healthcare in high-income countries: a systematic review. Eur J Public Health 2011; 22:904-13.
• ¹⁰ Carvalho VCP, Araújo TVB. Adequação da assistência pré-natal em gestantes atendidas em dois hospitais de referência para gravidez de alto risco do Sistema Único de Saúde, na cidade de Recife, Estado de Pernambuco. Rev Bras Saúde Matern Infant 2007; 7:309-17.
• ¹¹ Ribeiro ERO, Guimarães AMDN, Bettiol H, Lima DDF, Almeida MLD, Souza L, et al. Risk factors for inadequate prenatal care use in the metropolitan area of Aracaju, Northeast Brazil. BMC Pregnancy Childbirth 2009, 9:31.
• ¹² Trevisan MR, Lorenzi DS, Araújo NM, Ésber K. Perfil da assistência pré-natal entre usuárias do Sistema Único de Saúde em Caxias do Sul. Rev Bras Ginecol Obstet 2002; 24:293-9.
• ¹³ Kramer MS, Demissie K, Yang H, Platt RW, Sauvé R, Liston R. The contribution of mild and moderate preterm birth to infant mortality. Fetal and Infant Health Study Group of the Canadian Perinatal Surveillance System. JAMA 2000; 284:843-9.
• ¹⁴ Koga D. Medida de cidades: entre territórios de vida e territórios vividos. São Paulo: Cortez Editora; 2003.
• ¹⁵ Tomasi E, Fernandes PAA, Fischer T, Siqueira FCV, Silveira DS, Thumé E, et al. Qualidade da atenção pré-natal na rede básica de saúde do Brasil: indicadores e desigualdades sociais. Cad Saúde Pública 2017; 33:e00195815.
• ¹⁶ Borges ALV, Cavalhieri FB, Hoga LAK, Fujimori E, Barbosa LR. Planejamento da gravidez: prevalência e aspectos associados. Rev Esc Enferm USP 2011; 45:1679-84.
• ¹⁷ Rothman KJ, Greenland S. Modern epidemiology. 3rd Ed. Philadelphia: Lippincott-Raven; 2012.
• ¹⁸ Susser M. Causal thinking in the health sciences concepts and strategies of epidemiology. New York: Oxford University Press; 1973.
• ¹⁹ Silva AAM, Vasconcelos AGG, Bettiol H, Barbieri MA. Socioeconomic status, birth weight, maternal smoking during pregnancy and adiposity in early adult life: an analysis using structural equation modeling. Cad Saúde Pública 2010; 26:15-29.
• ²⁰ Amorim LDAF, Fiaccone RL, Santos CAST, Santos TN, Morais LTLP, Oliveira NS, et al. Modelagem de equações estruturais em epidemiologia. Cad Saúde Pública 2010; 26:2251-62.
• ²¹ Hoyle RH. Structural equation modeling: concepts, issues, and applications. Thousand Oaks: SAGE; 1995.
• ²² Bussab WO, Morettin PA. Estatística básica. 8ª Ed. São Paulo: Saraiva; 2013.
• ²³ Johnson RA, Wichern DW. Applied multivariate statistical analysis. 6th Ed. New Jersey: Prentice Hall; 2007.
• ²⁴ Schumacker RE, Lomax RG. A beginner's guide to structural equation modeling. 3rd Ed. New York: Routledge; 2010.
• ²⁵ Byrne BM. Structural equation modeling with Mplus: basic concepts, applications, and programming. New York: Routledge; 2011.
• ²⁶ Rios CTF, Vieira NFC. Ações educativas no pré-natal: reflexão sobre a consulta de enfermagem como um espaço para educação em saúde. Ciênc Saúde Coletiva 2007; 12:477-86.
• ²⁷ Líbera BD, Saunders C, Santos MMAS, Rimes KA, Brito FRSS, Balão MR. Avaliação da assistência pré-natal na perspectiva de puérperas e profissionais de saúde. Ciênc Saúde Coletiva 2011; 16:4855-64.
• ²⁸ Behrman RE, Butler AS. Preterm birth: causes, consequences and prevention. Washington DC: National Academies Press; 2007.
• ²⁹ Domingues RMSM, Viellas EF, Dias MAB, Torres JA, Theme-Filha MM, Gama SGN, et al. Adequação da assistência pré-natal segundo as características maternas no Brasil. Rev Panam Salud Pública 2015; 37:140-7.
• ³⁰ Lima AMC. Modelagem de equações estruturais: uma contribuição metodológica para o estudo da pobreza [Dissertação de Mestrado]. Salvador: Universidade Federal da Bahia; 2005.
• ³¹ Instituto de Pesquisa Econômica Aplicada. Pobreza multidimensional no Brasil. Rio de Janeiro: Instituto de Pesquisa Econômica Aplicada; 2006.
• ³² Padilha PC, Saunders C, Machado RCM, Silva CL, Bull A, Sally EOF, et al. Associação entre o estado nutricional pré-gestacional e a predição do risco de intercorrências gestacionais. Rev Bras Ginecol Obstet 2007; 29:511-8.
• ³³ Albertsen K, Andersen AMN, Olsen J, Gronbaek M. Alcohol consumption during pregnancy and the risk of preterm delivery. Am J Epidemiol 2004; 159:155-61.
• ³⁴ Peacock JL, Bland JM, Anderson HR. Preterm delivery: effects of socioeconomic factors, psychological stress, smoking, alcohol, and caffeine. BMJ 1995; 311:531-5.
• ³⁵ Victora CG, Huntly SR, Fuchs SC, Olinto MTA. The role of conceptual frameworks in epidemiological analysis: a hierarchical approach. Int J Epidemiol 1997; 26:224-7.
• ³⁶ Muthén LK, Muthén BO. Mplus user's guide. 6th Ed. Los Angeles: Muthén & Muthén; 2010.
• ³⁷ Oberski DL. lavaan.survey: an R package for complex survey analysis of structural equation models. J Stat Softw 2014; 57:1-27.
• ³⁸ Brown TA. Confirmatory factor analysis for applied research. New York: The Guilford Press; 2015.
• ³⁹ Muthén BO, du Toit SHC, Spisic D. Robust inference using weighted least squares and quadratic estimating equations in latent variable modelling with categorical and continuous outcomes. Los Angeles: University of California, Los Angeles; 1997. (Technical Report).
• ⁴⁰ Beauducel A, Herzberg PY. On the performance of maximum likelihood versus means and variance adjusted weighted least squares estimation in CFA. Struct Equ Modeling 2006; 13:186-203.

Datas de Publicação

• Publicação nesta coleção
  21 Jan 2019
• Data do Fascículo
  2019

Histórico

• Recebido
  05 Dez 2017
• Revisado
  03 Jul 2018
• Aceito
  06 Jul 2018