É necessária uma assinatura do JoVE para visualizar este conteúdo. Faça login ou inicie seu teste gratuito.
Artigo de investigação
Uma pontuação baseada em aprendizado profundo extraída de radiografias laterais toracolombares, combinada com fatores de risco clínicos, permitiu a predição precisa de fraturas vertebrais incidentes em até dois anos. O modelo validado internamente apresentou melhor discriminação, calibração, reclassificação e benefício decisório do que o modelo clínico, apoiando estratégias de estratificação de risco individualizada e de manejo preventivo precoce.
A identificação precoce de pacientes com risco de fratura vertebral incidente permanece desafiadora porque a avaliação clínica rotineira de risco não capta completamente a fragilidade espinhal local. Este estudo de coorte retrospectivo unicêntrico avaliou se características de aprendizado profundo (DL) extraídas de radiografias laterais toracolombares basais melhoram a predição de fratura vertebral incidente em até 2 anos quando combinadas com fatores de risco clínicos. Um total de 2.173 pacientes foi incluído e dividido cronologicamente em uma coorte de derivação (n = 1.449) e uma coorte de validação interna (n = 724). As características de DL foram derivadas das radiografias basais, e a regressão LASSO-Cox foi utilizada para selecionar preditores e construir um modelo clínico, um modelo de DL e um modelo combinado. O desempenho foi avaliado por correção de otimismo bootstrap, validação interna temporal, calibração, análise de curva de decisão, melhoria na reclassificação líquida dependente do tempo (NRI), melhoria na discriminação integrada (IDI) e análises de sensibilidade. De 2.048 características candidatas de DL, 5 foram mantidas para gerar um escore de DL, que permaneceu um preditor independente no modelo combinado (HR 1,64, IC 95% 1,34–2,01; P < 0,001). Na validação interna, o modelo combinado alcançou um índice C de 0,759, uma AUC de 2 anos de 0,774 e um escore de Brier de 2 anos de 0,077, todos superiores ao modelo clínico, com boa calibração (intercepto 0,012; inclinação 0,972). Em comparação com o modelo clínico, o modelo combinado também melhorou a reclassificação (NRI de 2 anos de 0,316 na derivação e 0,241 na validação) e a discriminação (IDI de 2 anos de 0,047 e 0,033, respectivamente; todos P < 0,01) e proporcionou maior benefício líquido na análise de curva de decisão. As análises de sensibilidade foram consistentes com os resultados principais. A combinação de características de DL de radiografias laterais toracolombares com fatores de risco clínicos pode permitir uma predição individualizada mais precisa de fratura vertebral incidente em até 2 anos.
A fratura vertebral é um dos tipos mais comuns de fraturas por fragilidade osteoporótica e é particularmente frequente na região toracolombar. Pode levar a dor crônica, perda de altura, deformidade cifótica, mobilidade limitada e aumentar o risco de refratura e prognóstico adverso1. Na prática clínica, uma proporção considerável de pacientes não apresenta sintomas típicos antes da ocorrência da fratura, e muitos casos são identificados apenas em exames de imagem de acompanhamento, o que sugere que confiar exclusivamente em sintomas ou diagnóstico retrospectivo dificulta a realização de triagem oportuna em populações de alto risco2,3. A avaliação de risco existente baseia-se principalmente em informações como idade, sexo, índice de massa corporal, fratura por fragilidade prévia, diabetes, exposição a glicocorticoides e densidade mineral óssea, que podem refletir o quadro de fragilidade óssea sistêmica, mas é difícil caracterizar plenamente a fragilidade estrutural local e as anormalidades mecânicas da coluna toracolombar, o que também representa uma dificuldade central persistente na predição do risco de nova fratura vertebral4. A radiografia lateral toracolombar é um dos exames de imagem da coluna mais utilizados e acessíveis na prática clínica. Ela não apenas pode mostrar a morfologia vertebral, mas também pode conter fenótipos ocultos relacionados a fraturas futuras, como alterações da placa terminal, textura óssea rarefeita, cunhificação leve e desalinhamento5. Estudos anteriores concentraram-se principalmente na detecção de fraturas vertebrais já existentes, no diagnóstico de osteoporose ou na avaliação de risco por meio de indicadores medidos manualmente6,7. Evidências recentes mostraram ainda que a fratura vertebral prevalente e a osteoporose identificadas por aprendizado profundo em imagens laterais da coluna, juntamente com fatores clínicos de risco, podem melhorar a predição de fratura incidente5; no entanto, as evidências ainda são limitadas quanto à predição de fratura vertebral incidente especificamente em pacientes sem fratura vertebral-alvo na linha de base, utilizando radiografias toracolombares laterais de rotina e características locais de aprendizado profundo (DL). Métodos de inteligência artificial têm sido usados na análise de imagens da coluna, mas estudos diretamente voltados a esse cenário clínico específico ainda são escassos, e a avaliação sistemática da calibração, do benefício líquido a partir da análise de decisão e da validação por divisão temporal ainda é insuficiente nesse contexto8.
Portanto, é difícil responder a uma pergunta clinicamente mais relevante: os recursos extraídos por aprendizado profundo de radiografias laterais toracolombares de rotina podem fornecer informações incrementais independentes e significativas com base na avaliação clínica de risco9? Com base no exposto, este estudo adotou um delineamento de coorte retrospectivo unicêntrico, extraiu características de aprendizado profundo de radiografias laterais toracolombares e as combinou com fatores de risco clínicos para construir um modelo de predição de risco de fratura vertebral incidente em até 2 anos, avaliando a discriminação, calibração, robustez e valor clínico do modelo por meio de validação interna temporal, correção de otimismo por bootstrap e análise de sensibilidade. Este estudo focou-se no alerta individualizado de risco com base em radiografia de rotina, integrando fenótipos ocultos de fragilidade local por imagem e informações sistêmicas de suscetibilidade clínica em uma ferramenta preditiva interpretável, com o objetivo de fornecer base para a identificação de alto risco, acompanhamento intensificado e intervenção preventiva.
Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.
Este estudo foi revisado e aprovado pelo Comitê de Ética Médica do Oitavo Hospital Popular de Xangai, Xangai, China (número de aprovação 2026-102-03-02). Como este estudo foi retrospectivo e todos os dados foram desidentificados antes da análise, o comitê de ética isentou o consentimento informado dos pacientes.
Desenho do estudo:
Tipo de estudo
Este estudo foi um estudo de coorte retrospectivo unicêntrico, e o banco de dados do estudo foi criado utilizando dados do sistema de arquivamento e comunicação de imagens (PACS), do sistema de informações radiológicas (RIS) e do sistema de prontuário eletrônico do hospital. A população do estudo compreendeu pacientes consecutivos que realizaram exame de radiografia digital lateral toracolombar no hospital. O período de inclusão estendeu-se de 1º de janeiro de 2018 a 31 de dezembro de 2023, e o prazo final de acompanhamento foi 31 de dezembro de 2025. O relatório do estudo seguiu as recomendações TRIPOD+AI e STROBE para garantir a padronização da apresentação de estudos de modelos preditivos envolvendo inteligência artificial e estudos observacionais.
Local do estudo e fonte dos casos
Os casos foram obtidos a partir do processo rotineiro de diagnóstico clínico e tratamento de pacientes ambulatoriais, de emergência e internados no hospital. Os dados de imagem foram todos obtidos a partir dos arquivos DICOM originais no PACS, e os dados clínicos foram obtidos a partir de prontuários eletrônicos estruturados, sistemas laboratoriais e registros de prescrições. A data do primeiro exame radiográfico lateral do tórax e da coluna lombar que atendeu aos critérios de inclusão durante o período do estudo foi definida como a data basal; quando o mesmo paciente tinha múltiplos exames que atendiam aos critérios, apenas o mais antigo foi mantido como exame basal para evitar inclusão repetida. Todos os dados foram anonimizados antes da análise, e as informações de imagem e clínicas foram pareadas utilizando um número de identificação único do estudo.
População do estudo:
Critérios de inclusão
Os critérios de inclusão foram os seguintes: idade igual ou superior a 50 anos; realização de exame radiográfico digital lateral toracolombar em posição ortostática padrão no hospital durante o período do estudo; imagens basais em formato DICOM rastreável; visualização completa das vértebras T10 a L4 nas imagens basais; ausência de fratura vertebral pré-existente de T10 a L4 na análise das imagens basais; variáveis clínicas basais previamente especificadas extraídas dos prontuários eletrônicos; pelo menos um exame de acompanhamento com radiografia toracolombar, tomografia computadorizada (TC) ou ressonância magnética (RM) dentro de 24 meses após o basal, ou ocorrência de nova fratura vertebral confirmada por imagem dentro de 24 meses após o basal.
Critérios de exclusão
Os critérios de exclusão foram os seguintes: fratura vertebral de T10 a L4 na linha de base; histórico definitivo de lesão traumática de alta energia na linha de base ou durante o acompanhamento; tumor primário ou metastático da coluna vertebral, infecção da coluna vertebral ou doença óssea destrutiva; cirurgia prévia de fixação interna toracolombar, vertebroplastia ou cifoplastia; escoliose com ângulo de Cobb maior que 30° ou deformidade cifótica evidente (incluindo deformidade cifótica do tipo Scheuermann, se presente), resultando na incapacidade de identificar com precisão as placas terminais de T10 a L4; artefato de movimento evidente, exposição anormal, oclusão metálica ou alcance insuficiente de exibição na imagem; incapacidade de confirmar variáveis-chave na linha de base ou informações sobre os desfechos a partir dos prontuários eletrônicos.
Processo de construção de coorte retrospectiva
A triagem da população do estudo foi realizada independentemente por dois pesquisadores de acordo com critérios pré-especificados, e as discordâncias foram resolvidas por meio de discussão para alcançar consenso. Após a conclusão da triagem dos casos, foi realizada a agrupamento em séries temporais com base na data basal: os pacientes incluídos de 1º de janeiro de 2018 a 31 de dezembro de 2021 constituíram a coorte de derivação para seleção de características e construção do modelo; os pacientes incluídos de 1º de janeiro de 2022 a 31 de dezembro de 2023 constituíram a coorte de validação interna para avaliação do desempenho do modelo. A divisão temporal, em vez da divisão aleatória, pode reduzir o risco de vazamento de informações e está mais próxima do cenário real de aplicação do modelo em pacientes subsequentes. O processo de triagem da população do estudo é apresentado na forma de um fluxograma.
Desfecho primário e sua determinação:
Definição do desfecho primário
O desfecho primário deste estudo foi a primeira fratura vertebral por fragilidade de T10 a L4 dentro de 24 meses após a linha de base. A janela temporal de predição do estudo foi previamente definida como 2 anos, e a saída do modelo foi a probabilidade individual de risco de fratura vertebral incidente nos próximos 2 anos.
Critérios para a determinação de uma fratura vertebral incidente
A fratura vertebral incidente foi definida da seguinte forma: em comparação com as imagens basais, exames de acompanhamento mostraram uma redução de 20% ou mais na altura anterior, média ou posterior de qualquer corpo vertebral de T10 a L4, com uma redução absoluta de altura de pelo menos 4 mm, ou o aparecimento de novo colapso da placa terminal ou interrupção cortical10. A determinação do desfecho foi realizada de forma abrangente com base em radiografias de tórax e coluna lombar, tomografia computadorizada (TC) e ressonância magnética (RM) de acompanhamento. A análise das imagens foi realizada independentemente por dois radiologistas musculoesqueléticos, com 8 e 12 anos de experiência diagnóstica relevante, respectivamente, e nenhum deles teve acesso aos dados clínicos ou aos resultados da saída do modelo durante a análise das imagens; em caso de discordância, a decisão foi tomada por um radiologista musculoesquelético sênior com 18 anos de experiência. Fraturas vertebrais causadas por tumor, infecção ou trauma de alta energia não foram contabilizadas como eventos de desfecho.
Ponto inicial, ponto final e janela de observação do acompanhamento
O ponto inicial do acompanhamento foi a data do exame radiográfico lateral de tórax e lombar basal. O ponto final do acompanhamento foi definido como o mais precoce entre os seguintes momentos: a data da primeira fratura vertebral incidente, 24 meses após o exame basal, a data do último exame de imagem da coluna que confirmou ausência de fratura vertebral ou a data do óbito. Fraturas que apareceram pela primeira vez após 24 meses não foram incluídas no desfecho primário. Pacientes sem eventos de desfecho foram considerados censurados.
Coleta de dados clínicos e definição de variáveis clínicas candidatas:
Dados demográficos e clínicos gerais
Os dados clínicos de referência foram extraídos do sistema de prontuário eletrônico por dois pesquisadores de acordo com um formulário padronizado de relato de caso, sem revisar os resultados da determinação do desfecho durante a extração. Os dados demográficos e clínicos gerais coletados incluíram idade, sexo, altura, peso e índice de massa corporal. A idade foi definida como a idade real na data de referência; peso e altura foram obtidos do registro mais próximo à data de referência, dentro de 30 dias antes ou depois dessa data; o índice de massa corporal foi calculado como o peso dividido pelo quadrado da altura, em quilogramas por metro quadrado.
História médica, uso de medicamentos e dados relacionados ao metabolismo ósseo
Com base na disponibilidade clínica e na generalização do modelo, os seguintes fatores de risco clínicos candidatos foram pré-especificados para inclusão: histórico prévio de fratura por fragilidade, diabetes melito tipo 2, artrite reumatoide, uso crônico de glicocorticoides orais e tratamento antiosteoporose na linha de base. As medições padronizadas da densidade mineral óssea na linha de base e o escore FRAX não foram pré-especificados como preditores candidatos, pois não estavam uniformemente disponíveis como variáveis padronizadas na linha de base em toda a coorte; em vez disso, diversos fatores clínicos relacionados ao FRAX foram considerados separadamente como variáveis candidatas individuais. O histórico prévio de fratura por fragilidade, o diagnóstico de doenças subjacentes e as informações sobre medicação foram todos derivados de prontuários eletrônicos, registros de alta hospitalar e sistemas de prescrição antes da linha de base, sendo exigido que todas as variáveis já existissem antes da linha de base para garantir que os preditores precedessem temporalmente o evento desfecho.
Critérios de definição para variáveis clínicas
A história prévia de fratura de fragilidade foi definida como uma fratura ocorrida após os 40 anos de idade, causada por lesão de baixa energia e claramente registrada no prontuário médico; fraturas do crânio, ossos faciais, ossos dos dedos das mãos e ossos dos dedos dos pés não foram incluídas nesta definição. O diabetes melito tipo 2 foi definido como um diagnóstico claro registrado antes do momento basal ou pelo uso prolongado de medicamentos hipoglicemiantes. A artrite reumatoide foi definida como um diagnóstico claro feito por um especialista em reumatologia no prontuário médico. O uso crônico de glicocorticoides por via oral foi definido como uma dose equivalente de prednisona de não menos que 5 mg/d por um período de não menos que 3 meses dentro do ano anterior ao momento basal. O tratamento antiosteoporose no momento basal foi definido como o uso contínuo de qualquer um dos seguintes medicamentos — bifosfonatos, denosumabe, teriparatida, raloxifeno, calcitonina, alfacalcidol ou calcitriol — nos 3 meses anteriores ao momento basal, por um período de duração de não menos que 8 semanas. A idade e o índice de massa corporal foram tratados como variáveis contínuas na modelagem e não foram categorizados artificialmente.
Aquisição de dados de imagem e pré-processamento de imagens
Protocolo de aquisição de radiografia lateral toracolombar
Todas as imagens basais foram radiografias laterais toracolombares padrão obtidas pelo sistema de radiografia digital do hospital. Durante o exame, os pacientes assumiram uma posição ortostática natural, com ambos os membros superiores flexionados para frente para reduzir a sobreposição dos ombros, e a área de imagem cobriu a região de T10 a L4. Utilizou-se controle automático de exposição para o exame, com uma faixa de tensão do tubo de 80–95 kV e uma distância fonte-imagem de 110 cm. Para o mesmo paciente, quando estavam disponíveis múltiplas radiografias laterais elegíveis na data basal, aquela com alcance de exibição completo e melhor qualidade de imagem foi selecionada como objeto de análise.
Critérios de inclusão de imagens e controle de qualidade
Exigia-se que as imagens basais atendessem aos seguintes critérios de qualidade: visualização completa das vértebras T10 a L4 e de suas placas superior e inferior; margens anterior e posterior vertebrais, placas terminais e limites corticais bem definidos; ausência de artefato evidente de movimento; ausência de superexposição ou subexposição grave; ausência de oclusão por metal em grande área; e ausência de distorção morfológica evidente causada por rotação da posição corporal. Também foram excluídas imagens com alterações degenerativas graves ou osteófitos que impediam a identificação confiável das margens vertebrais ou das placas terminais. Dois radiologistas musculoesqueléticos realizaram a avaliação de qualidade de todas as imagens basais, e qualquer imagem que não atendesse a algum dos critérios-chave de qualidade foi excluída.
Pré-processamento e padronização de imagens
Todas as imagens DICOM foram anonimizadas antes da análise. Os passos de pré-processamento incluíram a padronização da orientação das imagens, redimensionamento para uma resolução espacial de 0,30 mm × 0,30 mm, truncamento dos valores de escala de cinza entre o percentil 0,5o e o percentil 99,5o e a padronização dos valores dos pixels para o intervalo 0–1 utilizando o método de normalização min-max. O fluxo de trabalho de pré-processamento descrito foi mantido consistente entre a coorte de derivação e a coorte de validação, sendo totalmente concluído automaticamente por scripts previamente especificados, a fim de reduzir o viés causado por operações manuais.
Extração de características de imagem por aprendizado profundo:
Determinação da região de interesse
A região de interesse foi a região de projeção lateral da coluna entre a placa superior de T10 e a placa inferior de L4. Um radiologista musculoesquelético com 8 anos de experiência realizou a anotação em caixa retangular de todas as imagens basais no software ITK-SNAP, com o limite anterior definido a 5 mm anterior à margem anterior do corpo vertebral e o limite posterior definido a 5 mm posterior à margem posterior do corpo vertebral11; outro radiologista musculoesquelético com 12 anos de experiência revisou as imagens caso a caso. A ROI era uma caixa retangular em nível regional, e não uma segmentação rigorosa do contorno vertebral; portanto, osteófitos marginais comuns não foram removidos separadamente e podiam ser parcialmente incluídos se estivessem dentro do limite previamente especificado, enquanto os casos com alterações degenerativas graves o suficiente para obscurecer as margens vertebrais ou as placas vertebrais já haviam sido excluídos durante a avaliação da qualidade da imagem. Para avaliar a reprodutibilidade da anotação da região, 50 imagens foram selecionadas aleatoriamente e reanotadas pelo mesmo radiologista após 4 semanas, e independentemente reanotadas pelo segundo radiologista, para análise subsequente da estabilidade das características. Após o recorte da ROI, todas as imagens foram redimensionadas uniformemente para 224 × 224 pixels.
Arquitetura do modelo de aprendizado profundo e processo de extração de características
Este estudo utilizou a rede neural convolucional ResNet50 como extrator de características de aprendizado profundo. Os parâmetros da rede foram inicializados com pesos pré-treinados do ImageNet, e foi realizada uma adaptação de domínio auto-supervisionada em todas as imagens de ROI de referência da coorte de derivação, sem utilizar rótulos de desfecho durante o processo de adaptação. Especificamente, foi utilizada uma tarefa auto-supervisionada contrastiva, na qual duas visualizações independentemente aumentadas geradas a partir da mesma imagem de ROI foram tratadas como um par positivo, enquanto visualizações de pacientes diferentes dentro do mesmo mini-lote foram tratadas como pares negativos, de modo que o codificador pudesse adaptar-se à distribuição das imagens do estudo. O treinamento do modelo utilizou o otimizador AdamW, com uma taxa de aprendizado inicial definida em 1 × 10^-4, tamanho de lote de 64 e 200 épocas de treinamento; durante o treinamento, foi realizada ampliação dos dados com rotação de ±5°, redimensionamento de 0,9 a 1,1 vezes, translação de no máximo 10 pixels e perturbação de contraste de ±10%12. Essas ampliações foram utilizadas para gerar visualizações em pares para a tarefa auto-supervisionada, sendo utilizadas apenas imagens não rotuladas da coorte de derivação nesta etapa. Após a adaptação de domínio, nenhuma sintonização fina supervisionada por desfecho foi realizada, e o codificador principal adaptado foi fixado para a extração de características. Após a conclusão da adaptação de domínio, o vetor de 2.048 dimensões proveniente da camada de agrupamento médio global foi extraído como as características candidatas de aprendizado profundo para cada paciente.
Triagem de características de imagem e redução da dimensionalidade
Primeiro, o coeficiente de correlação intraclasse das características foi calculado com base nas 50 imagens com anotação repetida, e as características com ICC intraobservador e interobservador não inferiores a 0,80 foram mantidas para garantir a estabilidade das características frente a pequenas variações na ROI. Em seguida, as características mantidas foram padronizadas pelo escore Z na coorte de derivação, características com variância nula foram removidas e, para características com coeficiente de correlação par a par absoluto superior a 0,90, apenas uma delas foi mantida. Por fim, a regressão LASSO-Cox foi utilizada para seleção de características, e o parâmetro de penalização foi determinado por validação cruzada de 10 dobras de acordo com o critério de 1-SE. As características com coeficientes de regressão não nulos foram ponderadas e somadas de acordo com seus coeficientes para construir o escore de aprendizado profundo (escore DL)13. Após a fórmula de pontuação ser determinada na coorte de derivação, ela foi fixada inalterada e aplicada diretamente à coorte de validação interna.
Pré-processamento e integração de preditores candidatos:
Tratamento de dados ausentes e padronização de dados
Todas as variáveis clínicas candidatas foram obtidas a partir de campos estruturados do prontuário médico. As variáveis com taxa de omissão superior a 20% foram excluídas do processo de modelagem. Os valores ausentes restantes foram tratados por meio de imputação múltipla por equações encadeadas, gerando 10 conjuntos de dados imputados; o modelo de imputação incorporou todos os preditores candidatos, a variável indicadora do desfecho e a estimativa de risco cumulativo de Nelson-Aalen para preservar ao máximo a informação sobre o tempo até o evento. As variáveis clínicas contínuas e o escore DL foram padronizados utilizando a média e o desvio padrão da coorte de derivação, e os mesmos parâmetros de transformação foram aplicados à coorte de validação; as variáveis binárias foram codificadas uniformemente como 0 ou 1.
Seleção de fatores de risco clínicos
A pré-especificação dos fatores de risco clínicos candidatos baseou-se na interpretabilidade clínica, evidências prévias e disponibilidade de dados, não tendo sido utilizada triagem com base no valor P univariável. As variáveis clínicas candidatas incluídas na seleção por meio da regressão LASSO-Cox foram idade, sexo, índice de massa corporal, histórico prévio de fratura por fragilidade, diabetes melito tipo 2, artrite reumatoide, uso crônico de glicocorticoides orais e tratamento antiosteoporose na linha de base; altura e peso foram coletados de forma descritiva e utilizados para calcular o índice de massa corporal, mas não foram inseridos separadamente no modelo. A regressão LASSO-Cox foi realizada separadamente nos 10 conjuntos de dados imputados da coorte de derivação, sendo o parâmetro de penalização selecionado por validação cruzada com 10 dobras; as variáveis com coeficientes não nulos em pelo menos 7 conjuntos de dados imputados foram incluídas no modelo clínico final. Idade e índice de massa corporal foram testadas quanto a relações não lineares utilizando splines cúbicos restritos; caso o termo não linear não fosse estatisticamente significativo, mantinha-se a forma linear. A multicolinearidade foi avaliada por meio do fator de inflação da variância, e as variáveis com fator de inflação da variância superior a 5 não foram mantidas simultaneamente.
Construção do conjunto preditor combinado
Para evitar o sobreajuste causado pela inserção direta de características de imagem de alta dimensão no modelo, as informações de aprendizado profundo foram primeiro comprimidas em uma única variável contínua, o escore DL, e então inseridas conjuntamente em uma modelagem combinada, juntamente com os fatores de risco clínicos selecionados. Nenhum termo de interação foi pré-especificado no modelo combinado, a fim de manter a parcimônia e a interpretabilidade do modelo. O conjunto final de preditores combinados consistiu no escore DL e nas variáveis clínicas mantidas.
Construção do modelo de predição de risco:
Estratégia de modelagem
Na coorte de derivação, o modelo clínico, o modelo de aprendizado profundo e o modelo combinado foram estabelecidos separadamente. Os modelos utilizaram regressão de riscos proporcionais de Cox, com a primeira fratura vertebral por fragilidade incidente nos 24 meses após a linha de base como desfecho do estudo, e as regras de censura são descritas na definição de acompanhamento acima. Para controlar o sobreajuste, a complexidade do modelo combinado foi restringida antes da modelagem, e manteve-se uma razão relativamente alta de eventos por parâmetro tanto quanto possível. Os coeficientes de regressão finais e os erros padrão de cada modelo foram estimados separadamente nos 10 conjuntos de dados imputados e depois combinados utilizando as regras de Rubin. A função de risco basal foi estimada de acordo com o método de Breslow, e a probabilidade individual de risco em 2 anos foi calculada.
Construção do modelo clínico
O modelo clínico incluiu os fatores de risco clínicos mantidos após a seleção LASSO-Cox. Todas as variáveis contínuas foram mantidas na forma contínua e não foram dicotomizadas. Após o ajuste do modelo, a suposição de riscos proporcionais foi testada utilizando os resíduos de Schoenfeld; para as variáveis que não satisfizeram a suposição de riscos proporcionais, um termo de interação com ln(tempo) foi adicionado para correção. O modelo clínico foi utilizado para caracterizar a capacidade preditiva das informações clínicas tradicionais para fratura vertebral incidente.
Construção do modelo de imagem de aprendizado profundo
O modelo de aprendizado profundo foi estabelecido como um modelo de riscos proporcionais de Cox utilizando a pontuação de DL como único preditor, com o objetivo de quantificar a capacidade preditiva das características de aprendizado profundo obtidas a partir de radiografias laterais toracolombares basais para o risco de fratura vertebral incidente dentro de 2 anos. Este modelo não introduziu nenhuma informação clínica e, portanto, serviu como um modelo unimodal de imagem para comparação com os demais modelos.
Construção do modelo combinado
O modelo combinado acrescentou ainda o escore de DL com base no modelo clínico, construindo um modelo abrangente de predição com base em características de aprendizado profundo provenientes de radiografias laterais toracolombares combinadas com fatores de risco clínicos. Após o estabelecimento do modelo combinado, foi elaborado um nomograma de risco de 2 anos com base em seus coeficientes de regressão, para estimativa individualizada do risco e exibição da aplicação clínica.
Validação interna e avaliação do desempenho do modelo:
Método de validação interna
A validação interna adotou uma estratégia de validação interna única e separada no tempo. Todos os modelos estabelecidos na coorte de derivação foram diretamente aplicados à coorte de validação recrutada entre 1º de janeiro de 2022 e 31 de dezembro de 2023, após os parâmetros serem fixados, sem novo ajuste. Além disso, foram realizadas 1.000 reamostragens por bootstrap dentro da coorte de derivação para obter estimativas de desempenho corrigidas para otimismo, com o objetivo de avaliar a estabilidade do modelo.
Avaliação da discriminação
A discriminação do modelo foi avaliada pelo índice de concordância de Harrell e pela AUC dependente do tempo de 2 anos, calculada com base no método de ponderação pela probabilidade inversa de censura, ambos com intervalos de confiança de 95%. Uma discriminação mais alta indica que o modelo é mais capaz de distinguir indivíduos que desenvolverão ou não fraturas vertebrais incidentes no futuro. As diferenças na discriminação entre modelos foram calculadas utilizando o método de reamostragem bootstrap com intervalos de confiança de 95%.
Avaliação da calibração
A calibração do modelo foi avaliada utilizando a curva de calibração do risco em 2 anos, o intercepto de calibração, a inclinação da calibração e a pontuação de Brier em 2 anos. A curva de calibração foi traçada com base nos decís do risco previsto e foi corrigida por reamostragem bootstrap. Um intercepto de calibração próximo a 0, uma inclinação de calibração próxima a 1 e uma pontuação de Brier mais baixa indicam boa concordância entre o risco previsto e o risco realmente observado.
Avaliação do valor de aplicação clínica
O valor de aplicação clínica do modelo foi avaliado por meio de uma análise de curva de decisão de 2 anos, comparando o benefício líquido sob diferentes probabilidades de limiar. A faixa de probabilidade de limiar foi previamente definida como 0,05–0,30 para abranger o intervalo de risco que pode ser utilizado clinicamente para acompanhamento intensificado, avaliação óssea adicional ou manejo de intervenção14. Considerou-se que um modelo com maior benefício líquido possui melhor valor de suporte à decisão clínica.
Comparação de modelos e determinação do melhor modelo
O modelo clínico, o modelo de aprendizado profundo e o modelo combinado foram comparados de forma abrangente por meio de discriminação, calibração, pontuação de Brier e curva de decisão. O ganho do modelo combinado em relação ao modelo clínico foi quantificado adicionalmente utilizando a melhoria na reclassificação líquida dependente do tempo de 2 anos e a melhoria na discriminação integrada. O melhor modelo foi pré-especificado como aquele que simultaneamente apresentou maior discriminação, boa calibração, menor erro de predição e maior benefício líquido.
Análise Estatística:
As variáveis contínuas foram inicialmente avaliadas quanto ao padrão de distribuição utilizando o teste de Shapiro-Wilk; aquelas que apresentaram distribuição normal foram expressas como média ± desvio padrão, enquanto as com distribuição assimétrica foram relatadas como mediana e intervalo interquartílico; as variáveis categóricas foram apresentadas como número de casos e porcentagem. As comparações das características basais entre a coorte de derivação e a coorte de validação foram realizadas utilizando o teste t para amostras independentes, teste de Mann-Whitney U, teste χ2 ou teste exato de Fisher, conforme apropriado. As comparações basais foram usadas apenas para descrever as características das coortes e não serviram como base para seleção de variáveis. Todos os testes estatísticos foram bicaudais, e um valor de P < 0,05 foi considerado estatisticamente significativo. As análises estatísticas foram realizadas no software R, utilizando principalmente os pacotes survival, glmnet, mice, rms, timeROC e rmda; o pré-processamento das imagens e a análise de aprendizado profundo foram realizados no ambiente Python e PyTorch. Para avaliar a robustez dos resultados, foi realizada adicionalmente uma análise de casos completos como análise de sensibilidade.
Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.
Processo de construção do coorte retrospectivo e características basais dos coortes
Durante o período do estudo, foram recuperados registros radiográficos laterais toracolombares, e 6.114 pacientes foram incluídos para triagem após a remoção de duplicatas. Após a exclusão progressiva de pacientes com idade < 50 anos, daqueles com fraturas pré-existentes na linha de base e daqueles com acompanhamento insuficiente, um total de 2.173 pacientes foram finalmente incluídos, sendo...
Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.
O modelo combinado manteve o desempenho ideal mesmo após a correção para otimismo e validação interna temporal, sugerindo que as características de aprendizado profundo obtidas a partir da radiografia lateral toracolombar não representam uma simples repetição das informações clínicas, mas podem fornecer informações incrementais independentes e verificáveis para a avaliação de risco de fratura vertebral incidente nos próximos 2 anos. Sua importância reside na incorporação do contexto de fragilidade sistêmica e da fragilid...
Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.
Os autores declaram que não possuem conflitos de interesse.
Os autores agradecem ao pessoal do hospital do estudo pelo apoio na recuperação de imagens, extração e gerenciamento de dados. Os autores também agradecem a todos os clínicos e tecnólogos radiológicos envolvidos na assistência ao paciente e na aquisição de imagens. Este estudo foi financeiramente apoiado pelo Projeto de Pesquisa Médica no Distrito de Xuhui em 2024 (SHXH202405).
Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.
| Nome | Empresa | Número de catálogo | Comentários |
|---|---|---|---|
| pacote glmnet | CRAN | N/A | Utilizado para análise de regressão LASSO-Cox. |
| ITK-SNAP | Universidade da Pensilvânia / Projeto ITK-SNAP | N/A | Utilizado para anotação de ROI em imagens basais. |
| pacote mice | CRAN | N/A | Utilizado para imputação múltipla. |
| Python | Fundação Python de Software | versão 3.10 | Utilizado para pré-processamento de imagens e análise de aprendizado profundo. |
| PyTorch | Fundação PyTorch / Fundação Linux | versão 2.1 | Utilizado para desenvolvimento de modelos de aprendizado profundo e extração de características. |
| versão do R | Fundação R para Computação Estatística | versão 4.3.2 | Utilizado para análise estatística. |
| pacote rmda | CRAN | N/A | Utilizado para análise de curva de decisão. |
| pacote rms | CRAN | N/A | Utilizado para desenvolvimento do modelo e análise de calibração. |
| pacote survival | CRAN | N/A | Utilizado para análise de regressão de riscos proporcionais de Cox. |
| pacote timeROC | CRAN | N/A | Utilizado para análise de AUC dependente do tempo. |
Acesso restrito. Inicie sessão ou comece um teste para visualizar este conteúdo.