Este estudo é uma análise retrospectiva aprovada pelo Comitê de Ética Institucional do Nantong First People's Hospital. O estudo foi registrado em ClinicalTrials.gov (Número de registro: NCT06833099). Como as informações de saúde de todos os participantes foram anonimizadas, o consentimento informado não foi necessário. Os consumíveis e equipamentos utilizados estão listados na Tabela de Materiais.
1. População do estudo
O estudo incluiu dados clínicos e registros de imagem pré-operatórios de 436 pacientes submetidos à discectomia interlaminar endoscópica percutânea (PEID) para dor lombar e dor nas pernas devido à hérnia de disco L5-S1 no Nantong First People's Hospital entre 1º de janeiro de 2020 e 30 de junho de 2024. Com base nos critérios de exclusão, 309 pacientes foram incluídos, entre os quais a hérnia de disco lombar recorrente (rLDH) foi confirmada por meio de escores pós-operatórios da Escala Visual Analógica (EVA) e imagens de acompanhamento. Dos 309 pacientes, 33 apresentaram rLDH pós-cirurgia, enquanto os 276 restantes tiveram alívio significativo da dor lombar e nas pernas, com escores EVA reduzidos em mais de 60%. As variáveis de imagem para todos os participantes foram derivadas de exames pré-operatórios de raios-X, tomografia computadorizada e ressonância magnética, combinados com informações clínicas detalhadas, incluindo sexo, idade, altura, peso, IMC, escores EVA e outras variáveis relevantes (Figura 1).
2. Critérios de inclusão e exclusão
Critérios de inclusão para rLDH: (1) Pacientes com hérnia de disco lombar L5-S1 submetidos a PEID de segmento único. (2) Exames de imagem abrangentes concluídos dentro de um mês antes da cirurgia. (3) Redução do escore EVA pós-operatória ≥60%, seguida de aumento do escore e confirmação por imagem. (4) Nenhuma outra anormalidade detectada na imagem. (5) Período mínimo de acompanhamento de 6 meses.
Critérios de inclusão para não rLDH: (1) Pacientes com hérnia de disco lombar L5-S1 submetidos a PEID de segmento único. (2) Exames de imagem abrangentes concluídos dentro de um mês antes da cirurgia. (3) Redução do escore EVA pós-operatória ≥60%, sem recidiva. (4) Nenhuma outra anormalidade detectada na imagem. (5) Período mínimo de acompanhamento de 6 meses.
Critério de exclusão: (1) Presença de outras condições patológicas que causam lombalgia, como infecção discal, tumores espinhais, doenças ósseas metabólicas ou osteoporose. (2) História de disco lombar ou outras cirurgias da coluna vertebral. (3) Baixa qualidade de imagem ou dados de exame incompletos. (4) Perda de acompanhamento.
3. Variáveis categóricas e contínuas
As características clínicas pré-operatórias e os dados de imagem dos pacientes foram analisados estatisticamente e medidos (Tabela 1 e Tabela 2). As variáveis categóricas foram usadas para distinguir características básicas da doença, fatores de estilo de vida e outras variáveis, enquanto as variáveis contínuas representaram medidas específicas que descrevem o estado fisiológico dos pacientes e as alterações de imagem pré-operatórias. Para reduzir o viés, foram implementadas medidas rigorosas de controle de qualidade. Dois radiologistas e cirurgiões de coluna com mais de 10 anos de experiência clínica foram responsáveis pela análise estatística e medição dos dados de imagem. Para casos complexos, os dois médicos resolveram problemas por meio de consulta conjunta para garantir precisão e consistência no processamento de dados.
As variáveis categóricas incluíram: sexo, diabetes, hipertensão, doenças cardiovasculares e cerebrovasculares (CCD), escoliose, estenose espinhal, fatores desencadeantes (por exemplo, atividade extenuante, entorse, exposição ao frio, impacto ou nenhum gatilho claro), duração da doença (mais de 6 meses ou menos de 6 meses), dormência ou fraqueza, tipo de protrusão (contida ou não contida), degeneração discal (Graus I, II, III e IV, V), degeneração discal adjacente (Graus I, II, III e IV, V), alterações modic e calcificação discal.
As variáveis contínuas incluíram: idade, duração da cirurgia, altura, peso, índice de massa corporal (IMC), diâmetro máximo da hérnia de disco (MDHD), altura posterior do disco (PDH), distância entre os dois centros das vértebras (DBTVC), ângulo de orientação da articulação facetária (FJOA), ângulo do disco (AD), ângulo de inclinação sacral (SSA), ângulo de lordose lombar (LLA) e índice de altura posterior do disco (PDHI, PDHI = PDH / DBTVC).
4. Limpeza de dados e seleção de variáveis
Primeiro, o arquivo de dados "zjkj3" foi lido em R (versão 4.3.1, https://www.r-project.org/, Plataforma: x86_64-w64-mingw32/x64 (64 bits), Copyright (C) 2023 The R Foundation for Statistical Computing) usando a função read_excel e armazenado em uma variável de dados. A codificação de caracteres selecionada foi UTF-8 (padrão do sistema), que é um método de codificação amplamente utilizado capaz de suportar os idiomas escritos da maioria dos países ao redor do mundo. A função select foi então usada para separar a variável de destino das variáveis de recurso. Para garantir a reprodutibilidade da divisão dos dados, uma semente aleatória fixa de 3 foi definida e o conjunto de dados foi dividido em um conjunto de treinamento (80%) e um conjunto de teste (20%). Finalmente, uma função de regressão LASSO foi definida e uma parte dos dados de treinamento foi amostrada aleatoriamente usando a função glmnet para ajustar a regressão LASSO. O parâmetro de regularização ideal λ foi determinado usando regularização L1 e validação cruzada (padrão de 10 vezes), garantindo que nenhum dado duplicado do paciente fosse misturado entre os conjuntos de treinamento e teste. O erro de validação cruzada para cada valor λ foi calculado para determinar o valor λ ideal.
5. Desenvolvimento e avaliação de modelos
SVM (Support Vector Machine): Este modelo mapeia dados de treinamento para um hiperplano, maximizando a margem entre duas classes para prever o alvo. As configurações específicas de hiperparâmetros para o modelo SVM neste estudo são: kernel = "linear", custo = 1. A Figura 1 suplementar mostra as matrizes de confusão para SVM.
DT (Árvore de Decisão): Este modelo divide recursivamente os dados para criar uma estrutura de árvore para prever a variável de destino. A configuração de hiperparâmetro específica para o modelo DT neste estudo é: max_depth = 3. A Figura 2 suplementar mostra as matrizes de confusão para DT.
ADA (AdaBoost): Este modelo combina vários alunos fracos (normalmente árvores de decisão) para criar um aluno forte, melhorando o desempenho da classificação. As configurações específicas de hiperparâmetros para o modelo ADA neste estudo são: n_estimators = 150, learning_rate = 0,1, semente = 80. A Figura 3 suplementar mostra as matrizes de confusão para ADA.
LGBM (Light Gradient Boosting Machine): Este modelo usa um algoritmo de histograma para dividir recursos contínuos, acelerando o processo de treinamento e reduzindo a carga de memória, tornando-o uma estrutura de aumento de gradiente eficiente e rápida, particularmente adequada para conjuntos de dados em grande escala. As configurações específicas de hiperparâmetros para o modelo LGBM neste estudo são: num_leaves = 5, learning_rate = 0,05, n_estimators = 50. A Figura 4 suplementar mostra as matrizes de confusão para LGBM.
RF (Random Forest): Este modelo constrói várias árvores de decisão e combina suas previsões para melhorar a precisão e controlar o sobreajuste. As configurações específicas de hiperparâmetros para o modelo de RF neste estudo são: ntree = 310, mtry = 1, maxnodes = 10, max_depth = 1, seed = 80. A Figura 5 suplementar mostra as matrizes de confusão para RF.
XGB (Extreme Gradient Boosting): Este modelo é um algoritmo aprimorado baseado em árvores de decisão, utilizando um novo algoritmo de árvore de decisão de aumento de gradiente generalizado para acelerar a construção do modelo, com forte aplicabilidade para tarefas de classificação e regressão. As configurações específicas de hiperparâmetros para o modelo XGB neste estudo são: nrounds = 100, max_depth = 2, eta = 0,39, gamma = 0, colsample_bytree = 0,88, seed = 80. A Figura 6 suplementar mostra as matrizes de confusão para XGB.
O desempenho dos seis modelos foi avaliado usando as seguintes métricas: conjunto de teste ROC AUC, acurácia, sensibilidade, especificidade, valor preditivo positivo (VPP), valor preditivo negativo (VPN) e escore F1. A fórmula de cálculo é a seguinte22:

TP: Verdadeiro Positivo; TN: Verdadeiro Negativo; FP: Falso positivo; FN: Falso negativo.
6. Pesquisa em grade e ajuste de hiperparâmetros
Para otimizar o desempenho do modelo, a pesquisa em grade foi utilizada para explorar várias combinações de hiperparâmetros. Por meio de ajuste e avaliação iterativos, a combinação ideal de hiperparâmetros foi identificada para maximizar os valores de AUC ROC para os conjuntos de treinamento e teste, melhorando assim o desempenho geral do modelo. Os intervalos de hiperparâmetros ajustados são apresentados na Tabela 3. A pesquisa em grade é fornecida na Figura Suplementar 7, Figura Suplementar 8 e Tabela Suplementar 1, Tabela Suplementar 2, Tabela Suplementar 3 e Tabela Suplementar 4.
7. Classificação de importância variável
Os modelos de melhor desempenho - Random Forest (RF) e Extreme Gradient Boosting (XGB) - foram selecionados para classificação de importância variável. Essas classificações ajudam a identificar os preditores mais críticos de L5-S1 rLDH após PEID.
8. Análise estatística
Todas as análises estatísticas foram realizadas usando o software R. As variáveis categóricas foram relatadas como porcentagens, enquanto as variáveis contínuas foram expressas como média ± desvio padrão. As características clínicas e os parâmetros de imagem foram comparados entre o grupo não recorrente (não rLDH, n=276) e o grupo recorrente (rLDH, n=33).