Izabela: https://rpubs.com/alex_istrate/608122
Ioana: https://rpubs.com/alex_istrate/612312
Pentru variabilele calitative, am folosit pie-uri sau bar-chart-uri și am calculat frecvențele absolute și relative ale categoriilor formate. Pentru a descrie relațiile dintre variabilele calitative am folosit testele Fisher și Odds-Ratio (OR) cu IC95% sau Chi² și Cramer V. Am prezentat grafic rezultatele sub forma unor bar-chart-uri.
Pentru a studia relațiile dintre variabilele cantitative și cele binare, am folosit testul Mann-Whitney (MW). Am prezentat valorile p generate de aceste teste precum și mediile ±DS ale grupurilor și diferența mediilor cu IC95% asociat. Pentru a descrie relațiile dintre variabilele calitative am folosit testul Chi² sau Fisher și indicatorii Cramer V și Odds-Ratio (OR) cu IC95%. Am prezentat grafic rezultatele sub forma unor bar-chart-uri.
Pentru a studia relațiile dintre variabilele cantitative am folosit coeficientul de corelație Spearman (R), cu valoarea p asociata și am prezentat grafic relațiile sub forma unor scatter-plot-uri pe care am adăugat linia de regresie.
Am folosit WPS Office 2019 pentru managementul bazei de date. Pentru toate analizele statistice și graficele ulterioare am folosit R version 4.0.0 (2020-04-24) [1]. Am considerat p < 0.05 ca fiind semnificativ statistic și p < 0.1 ca prezentând doar o tendință spre semnificativitate statistică.
Din totalul de 40 pacienti, 16 (40.0%, IC95% = 24.9% - 56.7%) au fost diagnosticați cu schizofrenie paranoidă și 24 (60.0%) schizofrenie nediferențiată. Diferența dintre dintre aceste proporții nu a fost semnificativă statistic (test binomial: p = 0.268). Scorul CDSS a avut valori între 2 și 8, cu o medie 5.25 și o deviație standard de 2.33. Scorul CDSS nu a fost semnificativ statisitc diferit la pacienții cu schizofrenie paranoidă față de cei cu schizofrenie nediferențiată (5.812 ±2.40 vs 4.875 ±2.25, p=0.295).
Figura 1: Distribuția pacienților și a scorului CDSS în funcție de diagnostic.
Figura 2: Distribuția pacienților în funcție de scorul CDSS și diagnostic.
Majoritatea cazurilor (n=29, 72.5%) au fost femei. Scorul CDSS a fost statistic similar la cele 2 sexe chiar dacă a fost mai mare la bărbați (F: 4.897 ±2.13, mediana=4, M: 6.182 ±2.68, mediana=7, p=0.141).
Vârsta, cu valori de la 20 la 63 ani și cu o medie de 39.6 ±10.8 ani (mediana = 40 ani), nu s-a corelat semnificativ statistic cu scorul CDSS (R=-0.107, p=0.510).
Durata educației a avut valori de la 8 la 18 ani de studiu, cu o medie de 13.50 ±2.34 ani (mediana=13 ani), nu s-a corelat semnificativ statistic cu scorul CDSS (R=0.156, p=0.337).
Majoritatea cazurilor (n = 23, 57.5%) nu aveau un partener stabil, acești pacienți având un scor CDSS semnificativ mai mare decât cei cu un partener stabil (singur: 6.043 ±2.29, mediana = 7; partener stabil: 4.176 ±1.98, mediana = 4; p=0.011).
Majoritatea cazurilor (n=26, 65.0%) nu erau angajați, acești pacienți având un scor CDSS statistic egal cu al celor angajați (angajat: 5.571 ±2.53, mediana = 5; neangajat: 5.077 ±2.24, mediana = 5; p=0.586).
Figura 3: Distribuția pacienților pe grupe de vârstă și sexe. Diagnostic: SP=schizofrenie paranoidă, SN=schizofrenie nediferențiată.
Tabel 1: Descrierea demografică a eșantionului.
Variabila | Detalii | Total |
N= | 40 | |
Diagnostic | schizofrenie paranoidă | 16 (40.0%) |
schizofrenie nediferențiată | 24 (60.0%) | |
CDSS | M (min:max) | 5 (2:9) |
μ ±DS | 5.25 ±2.33 | |
Sex | F | 29 (72.5%) |
M | 11 (27.5%) | |
Vârsta (ani) | M (min:max) | 40 (20:63) |
μ ±DS | 39.60 ±10.8 | |
Durata educației (ani) | M (min:max) | 13 (8:18) |
μ ±DS | 13.50 ±2.34 | |
Status de angajat | angajat | 14 (35.0%) |
neangajat | 26 (65.0%) | |
Status marital | singur | 23 (57.5%) |
partener stabil | 17 (42.5%) | |
μ ±DS = Media (deviația standard); M (min:max) = Mediana (min:max); OR = odds-ratio [cu IC 95%] și p calculat prin testul Fisher); |
Tabel 2: Descrierea demografică a eșantionului în funcție de CDSS.
Subset | N | Mean ±SD | Med (min:max) |
CDSS (Shapiro-Wilk normality test: p=0.002) | |||
(total) | 40 (100.0%) | 5.250 ±2.33 | 5.00 (2.00:9.0) |
Diagnostic (Wilcoxon rank sum test with continuity correction: p=0.295) | |||
schizofrenie paranoidă | 16 (40.0%) | 5.812 ±2.40 | 7.00 (2.00:9.0) |
schizofrenie nediferențiată | 24 (60.0%) | 4.875 ±2.25 | 4.50 (2.00:9.0) |
Sex (Wilcoxon rank sum test with continuity correction: p=0.141) | |||
F | 29 (72.5%) | 4.897 ±2.13 | 4.00 (2.00:9.0) |
M | 11 (27.5%) | 6.182 ±2.68 | 7.00 (2.00:9.0) |
Vârsta (ani) (Spearman's rank correlation rho: R=-0.107, p=0.510) | |||
(total) | 40 (100.0%) | 39.600 ±10.77 | 40.00 (20.00:63.0) |
Durata educației (ani) (Spearman's rank correlation rho: R=0.156, p=0.337) | |||
(total) | 40 (100.0%) | 13.500 ±2.34 | 13.00 (8.00:18.0) |
Status de angajat (Wilcoxon rank sum test with continuity correction: p=0.586) | |||
angajat | 14 (35.0%) | 5.571 ±2.53 | 5.00 (2.00:9.0) |
neangajat | 26 (65.0%) | 5.077 ±2.24 | 5.00 (2.00:9.0) |
Status marital (Wilcoxon rank sum test with continuity correction: p=0.011) | |||
singur | 23 (57.5%) | 6.043 ±2.29 | 7.00 (2.00:9.0) |
partener stabil | 17 (42.5%) | 4.176 ±1.98 | 4.00 (2.00:8.0) |
Figura 4: Distribuția pacienților în funcție de scorul CDSS și sex.
Figura 5: istribuția pacienților în funcție de scorul CDSS și sex.
Figura 6: Distribuția pacienților în funcție de Varsta si CDSS.
Figura 7: Distribuția pacienților în funcție de diagnostic.
Figura 8: Distribuția pacienților în funcție de Statusul marital
Figura 9: Distribuția pacienților și a CDSS în funcție de Statusul de angajat
Figura 10: Distribuția pacienților în funcție de nivelul de educație și CDSS
Figura 11: Distribuția pacienților în funcție de diagnostic.
Fumatul a fost recunoscut de 15 pacienți (37.5%) și nu s-a asociat semnficativ statistic cu scorul CDSS (fumători: 5.267 ±2.43 vs nefumători: 5.240 ±2.31; p=0.877).
Vârsta la debulul bolii a avut valori de la 16 la 42 ani, cu o medie de 28.05 ±7.23 ani, și nu s-a corelat semnificativ statistic cu scorul CDSS (R=0.004, p=0.981).
Durata bolii a avut valori de la 1 la 33 ani, cu o medie de 11.28 ±8.4 ani, și nu s-a corelat semnificativ statistic cu scorul CDSS (R=-0.061, p=0.710).
Numărul de spitalizări a avut valori de la 1 la 15, cu o medie de 5.05 ±3.76 (mediana = 4), și nu s-a corelat semnificativ statistic cu scorul CDSS (R=-0.071, p=0.662).
Scorul GAF a avut valori de la 45 la 75, cu o medie de 61.62 ±8.12 (mediana = 60), și nu s-a corelat semnificativ statistic cu scorul CDSS (R=-0.423, p=0.007).
Tabel 3: Distribuția scorului CDSS, în funcție de Patologie psihiatrică
Subset | N | Mean ±SD | Med (min:max) |
CDSS (Shapiro-Wilk normality test: p=0.002) | |||
(total) | 40 (100.0%) | 5.250 ±2.33 | 5.00 (2.00:9.0) |
Diagnostic (Wilcoxon rank sum test with continuity correction: p=0.295) | |||
schizofrenie paranoidă | 16 (40.0%) | 5.812 ±2.40 | 7.00 (2.00:9.0) |
schizofrenie nediferențiată | 24 (60.0%) | 4.875 ±2.25 | 4.50 (2.00:9.0) |
Vârsta la debut (ani) (Spearman's rank correlation rho: R=0.004, p=0.981) | |||
(total) | 40 (100.0%) | 28.050 ±7.23 | 27.00 (16.00:42.0) |
Durata bolii (ani) (Spearman's rank correlation rho: R=-0.061, p=0.710) | |||
(total) | 40 (100.0%) | 11.275 ±8.40 | 10.00 (1.00:33.0) |
Nr. spitalizări (Spearman's rank correlation rho: R=-0.071, p=0.662) | |||
(total) | 40 (100.0%) | 5.050 ±3.76 | 4.00 (1.00:15.0) |
Fumat (Wilcoxon rank sum test with continuity correction: p=0.877) | |||
da | 15 (37.5%) | 5.267 ±2.43 | 5.00 (2.00:9.0) |
nu | 25 (62.5%) | 5.240 ±2.31 | 4.00 (2.00:9.0) |
GAF (Spearman's rank correlation rho: R=-0.423, p=0.007) | |||
(total) | 40 (100.0%) | 61.625 ±8.12 | 60.00 (45.00:75.0) |
Figura 12: Distribuția pacienților în funcție de Fumat
Figura 13: Distribuția pacienților în funcție de Varsta la debut si CDSS.
Figura 14: Distribuția pacienților în funcție de Varsta la debut
Figura 10: Distribuția pacienților în funcție de nivelul de educație și CDSS
Figura 15: Distribuția pacienților în funcție de diagnostic.
Figura 16: Distribuția pacienților în funcție de Nr. spitalizari și CDSS
Figura 17: Distribuția pacienților în funcție de diagnostic.
Figura 18: Distribuția pacienților în funcție de GAF și CDSS
Figura 19: Distribuția pacienților în funcție de diagnostic.
Dintre comorbidități, cea mai frecventă a fost obezitatea (73.1%), urmată de HTA și suprapondere (câte 5 pacienți, 19.2%). Niciuna dintre comorbiditățile studiate nu s-a asociat semnificativ statistic cu scorul CDSS.
Figura 20: Distribuția comorbidităților în funcție de diagnostic.
Tabel 3: Distribuția scorului CDSS, în funcție de Patologie psihiatrică
Subset | N | Mean ±SD | Med (min:max) |
CDSS (Shapiro-Wilk normality test: p=0.002) | |||
(total) | 40 (100.0%) | 5.250 ±2.33 | 5.00 (2.00:9.0) |
Suprapondere (Wilcoxon rank sum test with continuity correction: p=0.869) | |||
da | 5 (19.2%) | 4.600 ±2.70 | 4.00 (2.00:9.0) |
nu | 21 (80.8%) | 4.905 ±2.39 | 4.00 (2.00:9.0) |
Obezitate (Wilcoxon rank sum test with continuity correction: p=0.838) | |||
da | 19 (73.1%) | 5.000 ±2.49 | 4.00 (2.00:9.0) |
nu | 7 (26.9%) | 4.429 ±2.23 | 4.00 (2.00:9.0) |
Dislipidemie (Wilcoxon rank sum test with continuity correction: p=0.418) | |||
da | 1 (3.8%) | 7.000 ±NA | 7.00 (7.00:7.0) |
nu | 25 (96.2%) | 4.760 ±2.40 | 4.00 (2.00:9.0) |
DZ (Wilcoxon rank sum test with continuity correction: p>0.999) | |||
da | 1 (3.8%) | 4.000 ±NA | 4.00 (4.00:4.0) |
nu | 25 (96.2%) | 4.880 ±2.44 | 4.00 (2.00:9.0) |
HTA (Wilcoxon rank sum test with continuity correction: p>0.999) | |||
da | 5 (19.2%) | 5.000 ±2.92 | 4.00 (2.00:9.0) |
nu | 21 (80.8%) | 4.810 ±2.34 | 4.00 (2.00:9.0) |
Tiroidită (Wilcoxon rank sum test with continuity correction: p>0.999) | |||
da | 1 (3.8%) | 4.000 ±NA | 4.00 (4.00:4.0) |
nu | 25 (96.2%) | 4.880 ±2.44 | 4.00 (2.00:9.0) |
Figura 21: Asocierile dintre scorul CDSS și Comorbidități.
Dintre tratamente, cel mai frecvent a fost RC (35.0%, 14 pacienți), urmat de OLANZ (30.0%) și RISP (12.5%). Niciunul dintre tratamentele studiate nu s-a asociat semnificativ statistic cu scorul CDSS.
Figura 22: Distribuția Tratamentelor în funcție de diagnostic.
Tabel 4: Distribuția scorului CDSS, în funcție de Tratament.
Subset | N | Mean ±SD | Med (min:max) |
CDSS (Shapiro-Wilk normality test: p=0.002) | |||
(total) | 40 (100.0%) | 5.250 ±2.33 | 5.00 (2.00:9.0) |
ABIL (Wilcoxon rank sum test with continuity correction: p=0.327) | |||
da | 4 (10.0%) | 6.250 ±2.63 | 6.00 (4.00:9.0) |
nu | 36 (90.0%) | 5.139 ±2.31 | 5.00 (2.00:9.0) |
DEPA (Wilcoxon rank sum test with continuity correction: p=0.540) | |||
da | 1 (2.5%) | 7.000 ±NA | 7.00 (7.00:7.0) |
nu | 39 (97.5%) | 5.205 ±2.34 | 5.00 (2.00:9.0) |
INVEGA (Wilcoxon rank sum test with continuity correction: p=0.125) | |||
da | 1 (2.5%) | 2.000 ±NA | 2.00 (2.00:2.0) |
nu | 39 (97.5%) | 5.333 ±2.30 | 5.00 (2.00:9.0) |
LEPO (Wilcoxon rank sum test with continuity correction: p=0.726) | |||
da | 1 (2.5%) | 4.000 ±NA | 4.00 (4.00:4.0) |
nu | 39 (97.5%) | 5.282 ±2.35 | 5.00 (2.00:9.0) |
OLANZ (Wilcoxon rank sum test with continuity correction: p=0.492) | |||
da | 12 (30.0%) | 5.750 ±2.80 | 6.50 (2.00:9.0) |
nu | 28 (70.0%) | 5.036 ±2.12 | 4.50 (2.00:9.0) |
PALI (Wilcoxon rank sum test with continuity correction: p=0.540) | |||
da | 1 (2.5%) | 7.000 ±NA | 7.00 (7.00:7.0) |
nu | 39 (97.5%) | 5.205 ±2.34 | 5.00 (2.00:9.0) |
RC (Wilcoxon rank sum test with continuity correction: p=0.473) | |||
da | 14 (35.0%) | 4.786 ±1.93 | 4.50 (2.00:9.0) |
nu | 26 (65.0%) | 5.500 ±2.52 | 5.00 (2.00:9.0) |
RISP (Wilcoxon rank sum test with continuity correction: p>0.999) | |||
da | 5 (12.5%) | 5.200 ±1.79 | 5.00 (3.00:7.0) |
nu | 35 (87.5%) | 5.257 ±2.42 | 5.00 (2.00:9.0) |
ROMP (Wilcoxon rank sum test with continuity correction: p=0.930) | |||
da | 1 (2.5%) | 5.000 ±NA | 5.00 (5.00:5.0) |
nu | 39 (97.5%) | 5.256 ±2.36 | 5.00 (2.00:9.0) |
SERO (Wilcoxon rank sum test with continuity correction: p=0.917) | |||
da | 3 (7.5%) | 5.333 ±2.52 | 5.00 (3.00:8.0) |
nu | 37 (92.5%) | 5.243 ±2.35 | 5.00 (2.00:9.0) |
SOLI (Wilcoxon rank sum test with continuity correction: p=0.778) | |||
da | 2 (5.0%) | 5.000 ±4.24 | 5.00 (2.00:8.0) |
nu | 38 (95.0%) | 5.263 ±2.29 | 5.00 (2.00:9.0) |
XEPLION (Wilcoxon rank sum test with continuity correction: p=0.274) | |||
da | 1 (2.5%) | 8.000 ±NA | 8.00 (8.00:8.0) |
nu | 39 (97.5%) | 5.179 ±2.32 | 5.00 (2.00:9.0) |
Figura 23: Asocierile dintre scorul CDSS și Tratament (partea 1).
Figura 24: Asocierile dintre scorul CDSS și Tratament (partea 2).
Cu toate că am identificat puține asocieri semnificative statistic între scorul CDSS și restul variabilelor considerate separat, am identificat un număr mare de asocieri în analiza multivariată. Forest-plot-ul și tabelul de mai jos prezintă aceste asocieri, rezultate din regresia lineară multivariată a scorului CDSS cu restul variabilelor. Dintre acestea, am selectat prim metoda Stepwise doar acei predictori care influiențează semnfificativ statistic performanța regresiei. În grafic, am ales să reprezint coeficienții beta standardizați, deci regresia nu a fost calculată cu valorile originale ale variabilelor continue cu scorurile Z ale acestora (Z = media / deviația standard).
În mod semnificativ statistic, ajustând pentru restul predictorilor, creșterea CDSS s-a corelat cu: sexul masculin, statusul de angajat, statusul marital (singur), diagnosticul de schizofrenie paranoidă, durata mai mare a bolii și tratamentul cu PALI. Scăderea CDSS s-a corelat cu: sexul feminin, numărul mai mare de spitalizări, și scorul GAF mai mare.
Fără semnificație statistică (dar relevante pentru performața modelului statistic), scorul CDSS mai mre s-a corelat cu: Vârsta la debut, fumatul și tratamentul cu SOLI sau XEPLION, în timp ce scorul CDSS mai mic s-a corelat cu tratamentul cu INVEGA.
Figura 25: Asocierile multivariate dintre scorul CDSS restul variabilelor (selecție Stepwise). Coeficienți standardizați cu IC95%. * p<0.05, ** p<0.01, *** p<0.001.
Tabel 5: Asocierile multivariate dintre scorul CDSS restul variabilelor (selecție Stepwise).
Predictor | Beta (SE) | IC 95% | Beta std (SE) | IC 95% std | p |
(Intercept) | 9.91 (3.03) | 3.69, 16.14 | -0.95 (0.36) | -1.69, -0.22 | 0.003 |
Sex: F | -1.59 (0.59) | -2.81, -0.37 | -0.68 (0.26) | -1.21, -0.16 | 0.013 |
Status de angajat: angajat | 1.40 (0.65) | 0.06, 2.74 | 0.60 (0.28) | 0.02, 1.18 | 0.042 |
Status marital: singur | 3.05 (0.64) | 1.74, 4.36 | 1.31 (0.27) | 0.75, 1.87 | <0.001 |
Diagnostic: schizofrenie paranoidă | 1.36 (0.58) | 0.16, 2.56 | 0.59 (0.25) | 0.07, 1.10 | 0.028 |
Vârsta la debut (ani) | 0.72 (0.38) | -0.05, 1.49 | 0.22 (0.12) | -0.02, 0.46 | 0.067 |
Durata bolii (ani) | 1.80 (0.48) | 0.81, 2.79 | 0.65 (0.17) | 0.29, 1.01 | <0.001 |
Nr. spitalizări | -0.32 (0.10) | -0.53, -0.12 | -0.52 (0.16) | -0.85, -0.19 | 0.003 |
Fumat | 1.08 (0.64) | -0.23, 2.39 | 0.46 (0.27) | -0.10, 1.03 | 0.102 |
INVEGA | -2.79 (1.77) | -6.42, 0.84 | -1.20 (0.76) | -2.76, 0.36 | 0.126 |
PALI | 4.19 (1.64) | 0.81, 7.57 | 1.80 (0.71) | 0.35, 3.25 | 0.017 |
SOLI | 1.87 (1.20) | -0.60, 4.34 | 0.80 (0.52) | -0.26, 1.86 | 0.131 |
XEPLION | 2.02 (1.71) | -1.50, 5.54 | 0.87 (0.74) | -0.64, 2.38 | 0.249 |
GAF | -0.15 (0.04) | -0.23, -0.07 | -0.53 (0.14) | -0.81, -0.24 | <0.001 |
.