Noua comandă mediate extinde puternica suită Stata de inferență cauzală pentru a susține analiza de mediere cauzală. Analiza cauzală identifică și cuantifică efectele cauzale. Analiza de mediere cauzală le descâlcește. Sunt aceste efecte mediate prin intermediul unei alte variabile, un mediator?

Alegeți una dintre cele 23 de combinații de modele de rezultat și mediator, inclusiv modele liniare, logit și Poisson, pentru a estima efectul total și a-l descompune în efecte directe și indirecte (prin mediator). Calculați efectele directe controlate și proporția mediată. Reconstituiți efectele în rapoarte de cote, de risc și de rată de incidență. Reprezentați grafic efectele estimate. Obțineți predicții. Și multe altele.

. mediate (wellbeing, logit) (bonotonin, logit) (exercise)

Iteration 0:  EE criterion =  2.047e-17
Iteration 1:  EE criterion =  1.570e-32

Causal mediation analysis                                Number of obs = 2,000

Outcome model:     Logit
Mediator model:    Logit
Mediator variable: bonotonin
Treatment type:    Binary

Robust

Coefficient std. err. z P>|z| [95% conf. interval]

.1130778 .0287055 3.94 0.000 .056816 .1693397

.1457939 .0357194 4.08 0.000 .0757851 .2158027

.2588717 .0213767 12.11 0.000 .2169742 .3007692

Note: Outcome equation includes treatment–mediator interaction.

Efectul total (TE) estimat este de 0,26. Deoarece variabila noastră de rezultat este binară, acest efect este măsurat pe o scară de probabilitate. Interpretăm acest efect la fel ca un efect mediu al tratamentului: dacă fiecare individ din populație ar face exerciții fizice, probabilitatea unei bunăstări mai mari ar crește cu 0,26 puncte pe scara probabilităților, comparativ cu cazul în care nimeni nu ar face exerciții fizice.

Celelalte rezultate prezentate în tabel sunt estimările efectului natural indirect (NIE) și ale efectului natural direct (NDE). NIE ne spune în ce măsură exercițiul fizic afectează bunăstarea prin producerea de bonotonină. NDE surprinde efectul exercițiilor fizice asupra bunăstării prin alte mecanisme decât bonotonina. În acest caz, NIE de 0,11 înseamnă că o creștere de 0,11 a probabilității unei stări de bine mai mari se datorează efectului exercițiilor fizice prin producerea de bonotonină. Și pe baza NDE, o creștere de 0.15 se datorează altor mecanisme.

Mai sus, interpretăm efectele estimate pe scara de probabilitate ca diferențe de risc. Dacă am dori, am putea să le interpretăm și în termeni de rapoarte de risc sau rapoarte de cote. Pentru a obține rapoarte de probabilitate, de exemplu, folosim comanda de postestimare estat sau:

. estat or
estat or requires potential-outcome means; refitting model ...

Transformed treatment effects                            Number of obs = 2,000

Robust

Odds ratio std. err. z P>|z| [95% conf. interval]

1.575108 .1827157 3.92 0.000 1.254785 1.977204

1.87189 .2785201 4.21 0.000 1.398393 2.505713

2.948429 .2768771 11.51 0.000 2.452772 3.544249

Efectul total corespunde unui odds ratio de 2,95, care, în acest caz, este produsul dintre odds ratio-ul efectului direct și cel al efectului indirect. Pentru a calcula ratele de risc, am putea utiliza estat rr, iar dacă am fi adaptat un model Poisson pentru rezultat, am fi putut utiliza estat irr pentru a calcula ratele de incidență.

De asemenea, putem estima efectele directe controlate (CDE) folosind estat cde. De exemplu, am putea fi interesați de efectul direct în ipoteza contrafactuală că fie fiecare individ din populație experimentează o creștere a nivelului de bonotonină, fie nimeni nu o face. Pentru a face acest lucru, specificăm estat cde cu opțiunea mvalue(0 1) pentru a estima efectul direct controlat mediu al tratamentului cu mediatorul (binar) bonotonină setat la 0 sau 1:

. estat cde, mvalue(0 1)

Controlled direct effect                                 Number of obs = 2,000

Mediator variable: bonotonin
Mediator values:
  1._at: bonotonin = 0
  2._at: bonotonin = 1

Delta-method

CDE std. err. z P>|z| [95% conf. interval]

.1391299 .039573 3.52 0.000 .0615682 .2166916

.200756 .0505278 3.97 0.000 .1017234 .2997887

Efectul direct al exercițiilor fizice asupra creșterii stării de bine este de 0,14 pe scara probabilității dacă nimeni nu experimentează o creștere a nivelului de bonotonină și putem observa că efectul este de 0,2 dacă toți membrii populației ar experimenta o creștere a nivelului de bonotonină.

Din nou, dacă am dori să exprimăm aceste efecte în termeni de rapoarte de risc sau rapoarte de șanse, putem folosi opțiunile rr sau, respectiv, or. Aici folosim opțiunea or pentru a estima efectele directe controlate pe scara odds-ratio:

. estat cde, mvalue(0 1) or

Controlled direct effect                                 Number of obs = 2,000

Mediator variable: bonotonin
Mediator values:
  1._at: bonotonin = 0
  2._at: bonotonin = 1

Delta-method

Odds ratio std. err. z P>|z| [95% conf. interval]

1.835771 .3049524 3.66 0.000 1.325621 2.542244

2.257759 .4785054 3.84 0.000 1.490306 3.420422

În cele din urmă, am putea fi interesați să atribuim un număr la valoarea medierii. Putem folosi estat proportion pentru a calcula proporția de mediere:

. estat proportion

Proportion mediated                                      Number of obs = 2,000

Robust

Proportion std. err. z P>|z| [95% conf. interval]

.4368103 .1164639 3.75 0.000 .2085453 .6650752

Proporția estimată se bazează pe estimările efectului pe scara de probabilitate. Efectul indirect reprezintă 44% din efectul total.