logo

1 To-do list

  • Comparar pacientes com e sem osteoporose ( coluna T: com osteoporose 1; sem osteoporose 2). Variáveis para comparação: B, C, F,H, I, J, K, L, M, N,O, P, U,V, W, X, Y, Z, AA, AB, AC, AD, AH,AI, AM,AO.

  • Encontrar fatores de risco independentes (análise multi)

2 Descritiva Geral

Characteristic N = 381
idade 57 (49, 67)
sexo
    Feminino 23 (60.53%)
    Masculino 15 (39.47%)
tempo_de_dialise 71 (24, 107)
tipo_de_tx
    1 27 (71.05%)
    2 11 (28.95%)
tempo_de_tx 181 (78, 288)
doenca_de_base
    DRPA 2 (5.71%)
    FABRY 1 (2.86%)
    Hipertensão 6 (17.14%)
    Indeterminada 26 (74.29%)
ptx_previa
    Não 33 (86.84%)
    Sim 5 (13.16%)
fraqueza_muscular
    Não 35 (92.11%)
    Sim 3 (7.89%)
fratura
    Não 36 (94.74%)
    Sim 2 (5.26%)
dor_ossea
    Não 18 (47.37%)
    Sim 20 (52.63%)
prurido
    Não 38 (100.00%)
calcificacao_vascular
    Não 22 (68.75%)
    Sim 10 (31.25%)
calcio 10.40 (9.53, 11.38)
fosforo 2.75 (2.10, 3.13)
fa 106 (65, 155)
vit_d 33 (29, 36)
alb 4.25 (4.00, 4.43)
pth 206 (113, 284)
cr 1.30 (1.00, 1.60)
tfg 54 (41, 68)
colecalciferol
    Não 31 (81.58%)
    Sim 7 (18.42%)
dose_colecalciferol 20,000 (17,000, 32,500)
dose_pred 5.00 (5.00, 5.00)
aza
    Não 24 (63.16%)
    Sim 14 (36.84%)
ciclosporina
    Não 29 (76.32%)
    Sim 9 (23.68%)
tacrolimo
    Não 18 (47.37%)
    Sim 20 (52.63%)
osteoporose_total
    Não 18 (47.37%)
    Sim 20 (52.63%)
mmf
    Não 19 (50.00%)
    Sim 19 (50.00%)
dose_mmf 1,040 (720, 1,440)
t_scdore_colo -1.90 (-2.50, -0.55)
t_score_femur_total -0.80 (-1.60, 0.00)
t_score_coluna -2.50 (-3.50, -0.50)
1 Median (IQR); n (%)

3 Comparações para osteoporose total

Characteristic Não, N = 181 Sim, N = 201 p-value2
idade 50 (48, 57) 62 (55, 71) 0.002
sexo 0.2
    Feminino 9 (50.00%) 14 (70.00%)
    Masculino 9 (50.00%) 6 (30.00%)
tempo_de_dialise 96 (45, 116) 48 (18, 101) 0.2
tipo_de_tx 0.11
    1 15 (83.33%) 12 (60.00%)
    2 3 (16.67%) 8 (40.00%)
tempo_de_tx 95 (53, 189) 252 (138, 312) 0.009
doenca_de_base 0.6
    DRPA 1 (5.88%) 1 (5.56%)
    FABRY 1 (5.88%) 0 (0.00%)
    Hipertensão 4 (23.53%) 2 (11.11%)
    Indeterminada 11 (64.71%) 15 (83.33%)
ptx_previa 0.2
    Não 14 (77.78%) 19 (95.00%)
    Sim 4 (22.22%) 1 (5.00%)
fraqueza_muscular >0.9
    Não 17 (94.44%) 18 (90.00%)
    Sim 1 (5.56%) 2 (10.00%)
fratura 0.5
    Não 18 (100.00%) 18 (90.00%)
    Sim 0 (0.00%) 2 (10.00%)
dor_ossea 0.11
    Não 11 (61.11%) 7 (35.00%)
    Sim 7 (38.89%) 13 (65.00%)
prurido
    Não 18 (100.00%) 20 (100.00%)
calcificacao_vascular 0.7
    Não 11 (73.33%) 11 (64.71%)
    Sim 4 (26.67%) 6 (35.29%)
calcio 10.90 (9.63, 11.50) 9.75 (9.48, 10.78) 0.081
fosforo 2.90 (2.40, 3.40) 2.50 (2.10, 3.00) 0.2
fa 112 (74, 136) 88 (48, 155) 0.3
vit_d 29 (26, 37) 34 (32, 36) 0.6
alb 4.40 (4.20, 4.60) 4.10 (3.90, 4.35) 0.12
pth 234 (151, 357) 131 (93, 228) 0.034
cr 1.50 (1.00, 1.70) 1.15 (1.00, 1.40) 0.11
tfg 50 (39, 70) 59 (43, 64) 0.6
colecalciferol >0.9
    Não 15 (83.33%) 16 (80.00%)
    Sim 3 (16.67%) 4 (20.00%)
dose_colecalciferol 20,000 (17,000, 25,000) 27,500 (17,500, 41,250) 0.7
dose_pred 5.00 (5.00, 5.00) 5.00 (5.00, 5.00) 0.3
aza 0.076
    Não 14 (77.78%) 10 (50.00%)
    Sim 4 (22.22%) 10 (50.00%)
ciclosporina 0.7
    Não 13 (72.22%) 16 (80.00%)
    Sim 5 (27.78%) 4 (20.00%)
tacrolimo 0.022
    Não 5 (27.78%) 13 (65.00%)
    Sim 13 (72.22%) 7 (35.00%)
mmf 0.2
    Não 7 (38.89%) 12 (60.00%)
    Sim 11 (61.11%) 8 (40.00%)
dose_mmf 1,080 (720, 1,440) 1,000 (720, 1,440) 0.7
t_scdore_colo -0.30 (-1.25, 0.40) -2.50 (-2.75, -1.90) <0.001
t_score_femur_total 0.00 (-0.80, 1.00) -1.60 (-2.63, -0.85) <0.001
t_score_coluna -0.50 (-1.70, 0.20) -3.35 (-4.23, -3.00) <0.001
1 Median (IQR); n (%)
2 Wilcoxon rank sum test; Pearson’s Chi-squared test; Fisher’s exact test

4 Buscando fatores de risco para a osteoporose

4.1 Análise Univariada

Characteristic N OR1 95% CI1 p-value
idade 38 1.13 1.05, 1.25 0.004
sexo 38
    Feminino — —
    Masculino 0.43 0.11, 1.59 0.2
tempo_de_dialise 38 0.99 0.98, 1.00 0.11
tipo_de_tx 38 3.33 0.77, 17.9 0.12
tempo_de_tx 38 1.01 1.00, 1.02 0.010
doenca_de_base 35
    DRPA — —
    FABRY 0.00 >0.9
    Hipertensão 0.50 0.01, 17.5 0.7
    Indeterminada 1.36 0.05, 37.1 0.8
ptx_previa 38
    Não — —
    Sim 0.18 0.01, 1.41 0.15
dor_ossea 38
    Não — —
    Sim 2.92 0.80, 11.5 0.11
calcio 38 0.50 0.23, 0.98 0.054
fosforo 36 0.44 0.14, 1.11 0.11
alb 24 0.08 0.00, 1.25 0.11
pth 38 1.00 0.99, 1.00 0.2
cr 37 0.24 0.04, 1.04 0.078
dose_pred 37 1.52 0.89, 4.35 0.2
aza 38
    Não — —
    Sim 3.50 0.89, 15.9 0.083
tacrolimo 38
    Não — —
    Sim 0.21 0.05, 0.79 0.025
mmf 38
    Não — —
    Sim 0.42 0.11, 1.53 0.2
dose_mmf 18 1.00 1.00, 1.00 0.6
t_scdore_colo 35 0.22 0.06, 0.50 0.003
t_score_femur_total 33 0.25 0.08, 0.55 0.005
t_score_coluna 37 0.09 0.01, 0.32 0.006
1 OR = Odds Ratio, CI = Confidence Interval

Na análise univariada somente as variáveis idade (p = 0.004), tempo de tx (p = 0.01), e tacrolimo (p = 0.025) foram estatisticamente significativas para osteoporose. Para cada ano adicional na idade, há uma chance 1.13 (IC 95%: 1.05-1.25) maior de desenvolver osteoporose. Para cada mês adicional de tx há uma chance de 1.01 (C 95%: 1.00-1.02) maior de desenvolver osteoporose. Os pacientes que usam tacrolimo têm um risco reduzido de desenvolver osteoporose 0.21 (IC 95%: 0.05-0.79), ou, de forma similiar, podemos dizer que os pacientes que não fazem uso de tacrolimo, possuem 4.76 vezes mais risco de desenvolver osteoporose do que as que usam, com IC 95%: 1.26-20.

4.2 Análise Multivariável

Na análise multivariável foram testados alguns modelos com todas as variáveis anteriores, porém, algumas delas forneceram resultados incoerentes (valores de OR muito grandes) e, portanto, foram excluídas do modelo. A manutenção ou exclusão de variáveis também pode ser feita por importância clínica das mesmas. Se alguma das variáveis listadas abaixo não forem coerentes, podemos excluí-la(s) e rodar uma nova análise multivariável.

Characteristic OR1 95% CI1 p-value
tacrolimo
    Não — —
    Sim 0.22 0.00, 8.50 0.4
idade 1.29 1.07, 1.81 0.040
tempo_de_tx 1.00 0.98, 1.03 0.8
mmf
    Não — —
    Sim 0.43 0.00, 315 0.8
aza
    Não — —
    Sim 0.02 0.00, 49.7 0.3
sexo
    Feminino — —
    Masculino 0.88 0.02, 29.5 >0.9
tempo_de_dialise 1.0 0.95, 1.04 0.8
dose_pred 1.31 0.21, 399 0.8
cr 0.07 0.00, 1.45 0.14
pth 0.99 0.98, 1.00 0.15
fosforo 0.22 0.00, 6.35 0.4
ptx_previa
    Não — —
    Sim 0.06 0.00, 48.3 0.4
tipo_de_tx 2.88 0.05, 439 0.6
1 OR = Odds Ratio, CI = Confidence Interval

4.3 Modelo somente com as variáveis que foram significativas na univariada

Characteristic OR1 95% CI1 p-value
idade 1.12 1.02, 1.26 0.025
tempo_de_tx 1.00 0.99, 1.01 0.5
tacrolimo
    Não — —
    Sim 0.39 0.05, 2.80 0.3
1 OR = Odds Ratio, CI = Confidence Interval

Na análise multivariável somente a idade foi um fator estatisticamente significativo para a osteoporose, quando ajustada parao tempo de tx e uso de tacrolimo.

4.4 Modelo incluindo as variáveis: T-score colo, T-score-coluna e T-score femur

Estas variáveis foram significativas, assim, vamos incluí-las no modelo multivariável com as outras variáveis também significativas.

Characteristic OR1 95% CI1 p-value
idade 0.61 0.00, Inf >0.9
tempo_de_tx 0.48 0.00, Inf >0.9
tacrolimo
    Não — —
    Sim 0.00 0.00, Inf >0.9
t_scdore_colo 0.00 0.00, Inf >0.9
t_score_femur_total 6,871,869,371,826 0.00, Inf >0.9
t_score_coluna 0.00 0.00, Inf >0.9
1 OR = Odds Ratio, CI = Confidence Interval

Como visto na tabela acima, essas variáveis em conjunto geram um modelo inconsistente, com OR zero e valores infinitos no intervalo de confiança. Isto pode ocorrer por vários motivos como:

  • Tamanho da amostra pequeno

  • Multicolinearidade

  • Relações não-lineares entre preditores

  • entre outros.

4.4.1 Testando um outro modelo somente com os t-score do colo e femur total

Characteristic OR1 95% CI1 p-value
idade 1.04 0.92, 1.18 0.5
tempo_de_tx 1.00 0.99, 1.02 0.7
tacrolimo
    Não — —
    Sim 0.43 0.02, 6.32 0.5
t_scdore_colo 0.07 0.00, 0.76 0.077
t_score_femur_total 2.86 0.28, 71.3 0.4
1 OR = Odds Ratio, CI = Confidence Interval

5 Correlações entre Pth e os t-scores

## [1] "Coeficiente de correlação entre pth e t-score colo: 0.35"
## [1] "Coeficiente de correlação entre pth e t-score femur total: 0.24"
## [1] "Coeficiente de correlação entre pth e t-score coluna: 0.21"
  • Os coeficientes indicam uma correlação fraca entre as variáveis, os testes estatísticos abaixo, informam se há uma correlação estatisticamente significativa.
## [1] "PTH e t-score colo"
## 
##  Kendall's rank correlation tau
## 
## data:  df_mdl_osteoporose$pth and df_mdl_osteoporose$t_scdore_colo
## z = 2.9044, p-value = 0.00368
## alternative hypothesis: true tau is not equal to 0
## sample estimates:
##      tau 
## 0.349045
## [1] "PTH e t-score femur total"
## 
##  Kendall's rank correlation tau
## 
## data:  df_mdl_osteoporose$pth and df_mdl_osteoporose$t_score_femur_total
## z = 1.7706, p-value = 0.07663
## alternative hypothesis: true tau is not equal to 0
## sample estimates:
##       tau 
## 0.2192409
## [1] "PTH e t-score coluna"
## 
##  Kendall's rank correlation tau
## 
## data:  df_mdl_osteoporose$pth and df_mdl_osteoporose$t_score_coluna
## z = 2.1991, p-value = 0.02787
## alternative hypothesis: true tau is not equal to 0
## sample estimates:
##       tau 
## 0.2541615
  • Para t-score da coluna (p = 0.028) e do colo (p = 0.0037) houve uma correlação estatisticamente significativa com o pth, porém, como mencionado acima, esta correlação é fraca. Assim, a predição de tais escores dado um valor do pth é limitada.

5.1 Visualização