1 FUNCTIONS

2 DATA PREPARATION

2.1 Load all sheets of an Excel

X: Leader delegation=LDLV Level-1 Moderator [Mo1]: Leader emotional exhaustion =LEDV Level-2 Moderator [Mo2]: LMX=LMX

Mediator [M]: Follower work engagement =MWEV Y: Follower task performance(t+1)=LPERFV2

Controls: [C1] Interaction Freuencey=MIFV, [C2] Follower work engagement (t-1)=MWEV0, [C4] Leader delegation(t-1)=LDLV0, [C5] Leader emotional exhaustion (t-1)=LEDV0, [C3] Follower task performance(t-1)=LPERFV0

2.2 Reproduce HRs (f - mod lmx.out): All group mean centered

    X=LDLV0; !leader delegation
    C0=MIFV0;
    C1=MIFV;
    C2=MWEV0;
    C3=LPERFV0;
    C4=LDLV;
    C5=LEDV;
    !C6=LPERFV;
    !Mo1=LEEHV0; !leader emotionale exhaustion
    Mo2=MLMXP; !LSABP;
    Me=MWEV;
    Y=LPERFV2;
    CENTER X(GROUPMEAN);
    !CENTER X Mo1(GROUPMEAN);
    !Inter1 = X*Mo1;
    CENTER Mo2 (GRANDMEAN);

2.3 Same day: All group mean centered

3 BASED ON C1

## 
## Hierarchical Linear Model (HLM)
## (also known as) Linear Mixed Model (LMM)
## (also known as) Multilevel Linear Model (MLM)
## 
## Model Information:
## Formula: MWEV ~ LDLV.GroC * LEEHV.GroC + C1 + (LDLV.GroC + LEEHV.GroC + Int1 | MNum)
## Level-1 Observations: N = 859
## Level-2 Groups/Clusters: MNum, 106
## 
## Model Fit:
## AIC = 2561.051
## BIC = 2637.143
## R_(m)² = 0.01014  (Marginal R²: fixed effects)
## R_(c)² = 0.62124  (Conditional R²: fixed + random effects)
## Omega² = NA  (= 1 - proportion of unexplained variance)
## 
## ANOVA Table:
## ─────────────────────────────────────────────────────────────────
##                       Sum Sq Mean Sq NumDF  DenDF     F     p    
## ─────────────────────────────────────────────────────────────────
## LDLV.GroC               8.23    8.23  1.00  30.34 10.34  .003 ** 
## LEEHV.GroC              0.37    0.37  1.00  53.43  0.47  .496    
## C1                      2.12    2.12  1.00 737.25  2.67  .103    
## LDLV.GroC:LEEHV.GroC    4.43    4.43  1.00  18.92  5.56  .029 *  
## ─────────────────────────────────────────────────────────────────
## 
## Fixed Effects:
## Unstandardized Coefficients (b or γ):
## Outcome Variable: MWEV
## ───────────────────────────────────────────────────────────────────────────
##                          b/γ    S.E.     t    df     p      [95% CI of b/γ]
## ───────────────────────────────────────────────────────────────────────────
## (Intercept)            5.120 (0.114) 44.87 104.3 <.001 *** [ 4.894,  5.347]
## LDLV.GroC              0.143 (0.045)  3.22  30.3  .003 **  [ 0.052,  0.234]
## LEEHV.GroC             0.022 (0.032)  0.69  53.4  .496     [-0.042,  0.085]
## C1                     0.043 (0.027)  1.63 737.3  .103     [-0.009,  0.096]
## LDLV.GroC:LEEHV.GroC  -0.075 (0.032) -2.36  18.9  .029 *   [-0.142, -0.008]
## ───────────────────────────────────────────────────────────────────────────
## 'df' is estimated by Satterthwaite approximation.
## 
## Standardized Coefficients (β):
## Outcome Variable: MWEV
## ───────────────────────────────────────────────────────────────────────────
##                            β    S.E.     t    df     p        [95% CI of β]
## ───────────────────────────────────────────────────────────────────────────
## LDLV.GroC              0.081 (0.025)  3.22  30.3  .003 **  [ 0.030,  0.132]
## LEEHV.GroC             0.015 (0.022)  0.69  53.4  .496     [-0.029,  0.060]
## C1                     0.035 (0.021)  1.63 737.3  .103     [-0.007,  0.077]
## LDLV.GroC:LEEHV.GroC  -0.062 (0.026) -2.36  18.9  .029 *   [-0.118, -0.007]
## ───────────────────────────────────────────────────────────────────────────
## 
## Random Effects:
## ──────────────────────────────────────────
##  Cluster  K   Parameter   Variance     ICC
## ──────────────────────────────────────────
##  MNum     106 (Intercept)  1.26064 0.61302
##               LDLV.GroC    0.02066        
##               LEEHV.GroC   0.00147        
##               Int1         0.00421        
##  Residual                  0.79580        
## ──────────────────────────────────────────