#Results of an ANOVA examining differences in emotional exhaustion by physical activity group (low, medium, high)
summary(anova1)
##              Df Sum Sq Mean Sq F value Pr(>F)  
## PA_level      3   1866   622.0   3.133 0.0253 *
## Residuals   516 102449   198.5                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
TukeyHSD(anova1)
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = Burnout_EE ~ PA_level, data = b)
## 
## $PA_level
##                diff        lwr         upr     p adj
## Low-High  3.6206897  -1.034918  8.27629704 0.1876141
## Mod-High -0.4509091  -4.336854  3.43503552 0.9906960
## NA-High  22.0000000 -14.459045 58.45904466 0.4053362
## Mod-Low  -4.0715987  -8.092362 -0.05083596 0.0458708
## NA-Low   18.3793103 -18.094350 54.85297078 0.5641062
## NA-Mod   22.4509091 -13.932518 58.83433651 0.3849595
#Results of an ANOVA examining differences in depersonalization by physical activity group (low, medium, high)
summary(anova2)
##              Df Sum Sq Mean Sq F value Pr(>F)
## PA_level      3     95   31.76   0.713  0.545
## Residuals   516  22982   44.54
TukeyHSD(anova2)
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = Burnout_DP ~ PA_level, data = b)
## 
## $PA_level
##                 diff       lwr       upr     p adj
## Low-High  0.02774784 -2.177292  2.232788 0.9999877
## Mod-High -0.46156250 -2.302066  1.378941 0.9168009
## NA-High   7.89843750 -9.369698 25.166573 0.6404638
## Mod-Low  -0.48931034 -2.393668  1.415048 0.9111902
## NA-Low    7.87068966 -9.404368 25.145747 0.6433754
## NA-Mod    8.36000000 -8.872320 25.592320 0.5950168
#Results of an ANOVA examining differences in personal achievement by physical activity group (low, medium, high)summary(anova3)
summary(anova3)
##              Df Sum Sq Mean Sq F value  Pr(>F)   
## PA_level      3   1218   406.1    5.38 0.00118 **
## Residuals   516  38950    75.5                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
TukeyHSD(anova3)
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = Burnout_PA ~ PA_level, data = b)
## 
## $PA_level
##                 diff         lwr         upr     p adj
## Low-High  -2.9684806  -5.8391073 -0.09785386 0.0395279
## Mod-High  -0.5237784  -2.9198345  1.87227770 0.9428355
## NA-High  -25.0546875 -47.5351718 -2.57420323 0.0219729
## Mod-Low    2.4447022  -0.0344822  4.92388659 0.0548997
## NA-Low   -22.0862069 -44.5757032  0.40328939 0.0563899
## NA-Mod   -24.5309091 -46.9647681 -2.09705008 0.0257329

#ANCOVA—-

##Results of ANCOVA using gender as a covariate
##Female (1), Male (2), TGNC (3)

##Emotional Exhaustion
anova(ancova1)
## Analysis of Variance Table
## 
## Response: Burnout_EE
##            Df Sum Sq Mean Sq F value  Pr(>F)  
## PA_level    3   1866  622.04  3.1485 0.02477 *
## Gender      2    900  449.95  2.2775 0.10358  
## Residuals 514 101549  197.57                  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Posthoc test in emotional exhaustion by gender
posthoc_EE_Gender
## $lsmeans
##  Gender lsmean   SE  df lower.CL upper.CL
##  1        28.9 3.64 514     21.8     36.1
##  2        27.7 3.58 514     20.7     34.8
##  3        37.4 5.86 514     25.9     48.9
## 
## Results are averaged over the levels of: PA_level 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast          estimate   SE  df t.ratio p.value
##  Gender1 - Gender2     1.16 1.27 514   0.915  0.6310
##  Gender1 - Gender3    -8.51 4.77 514  -1.782  0.1765
##  Gender2 - Gender3    -9.67 4.77 514  -2.029  0.1063
## 
## Results are averaged over the levels of: PA_level 
## P value adjustment: tukey method for comparing a family of 3 estimates
#Posthoc test in emotional exhaustion by physical activity group
posthoc_EE_PA
## $lsmeans
##  PA_level lsmean    SE  df lower.CL upper.CL
##  High       25.0  1.93 514     21.2     28.8
##  Low        28.4  1.99 514     24.4     32.3
##  Mod        24.4  1.73 514     21.0     27.8
##  NA         47.6 14.16 514     19.8     75.4
## 
## Results are averaged over the levels of: Gender 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast   estimate    SE  df t.ratio p.value
##  High - Low   -3.341  1.84 514  -1.817  0.2666
##  High - Mod    0.589  1.52 514   0.387  0.9804
##  High - NA   -22.600 14.12 514  -1.601  0.3790
##  Low - Mod     3.930  1.56 514   2.517  0.0585
##  Low - NA    -19.259 14.14 514  -1.362  0.5238
##  Mod - NA    -23.189 14.10 514  -1.645  0.3544
## 
## Results are averaged over the levels of: Gender 
## P value adjustment: tukey method for comparing a family of 4 estimates
##Depersonalization
anova(ancova2)
## Analysis of Variance Table
## 
## Response: Burnout_DP
##            Df  Sum Sq Mean Sq F value  Pr(>F)  
## PA_level    3    95.3  31.756  0.7217 0.53933  
## Gender      2   366.0 183.007  4.1592 0.01615 *
## Residuals 514 22616.1  44.000                  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Posthoc test in depersonalization by gender
posthoc_DP_Gender
## $lsmeans
##  Gender lsmean   SE  df lower.CL upper.CL
##  1        8.91 1.72 514     5.54     12.3
##  2       10.53 1.69 514     7.21     13.9
##  3       12.01 2.77 514     6.57     17.4
## 
## Results are averaged over the levels of: PA_level 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast          estimate   SE  df t.ratio p.value
##  Gender1 - Gender2    -1.62 0.60 514  -2.702  0.0194
##  Gender1 - Gender3    -3.09 2.25 514  -1.373  0.3558
##  Gender2 - Gender3    -1.47 2.25 514  -0.654  0.7901
## 
## Results are averaged over the levels of: PA_level 
## P value adjustment: tukey method for comparing a family of 3 estimates
#Posthoc test in depersonalization by PA group
posthoc_DP_PA
## $lsmeans
##  PA_level lsmean    SE  df lower.CL upper.CL
##  High       8.54 0.912 514     6.75     10.3
##  Low        9.05 0.939 514     7.20     10.9
##  Mod        8.40 0.816 514     6.80     10.0
##  NA        15.95 6.681 514     2.82     29.1
## 
## Results are averaged over the levels of: Gender 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast   estimate    SE  df t.ratio p.value
##  High - Low   -0.510 0.868 514  -0.588  0.9358
##  High - Mod    0.136 0.719 514   0.190  0.9976
##  High - NA    -7.413 6.662 514  -1.113  0.6819
##  Low - Mod     0.647 0.737 514   0.877  0.8166
##  Low - NA     -6.903 6.672 514  -1.035  0.7292
##  Mod - NA     -7.550 6.653 514  -1.135  0.6681
## 
## Results are averaged over the levels of: Gender 
## P value adjustment: tukey method for comparing a family of 4 estimates
##Personal Achievement by Gender and PA Group
anova(ancova3)
## Analysis of Variance Table
## 
## Response: Burnout_PA
##            Df Sum Sq Mean Sq F value   Pr(>F)   
## PA_level    3   1218  406.13  5.4773 0.001036 **
## Gender      2    838  419.07  5.6518 0.003733 **
## Residuals 514  38112   74.15                    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Posthoc test in gender by PA group
posthoc_PA_Gender
## $lsmeans
##  Gender lsmean   SE  df lower.CL upper.CL
##  1        28.2 2.23 514     23.8     32.6
##  2        26.3 2.19 514     22.0     30.6
##  3        20.6 3.59 514     13.6     27.7
## 
## Results are averaged over the levels of: PA_level 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast          estimate    SE  df t.ratio p.value
##  Gender1 - Gender2     1.92 0.779 514   2.467  0.0370
##  Gender1 - Gender3     7.62 2.925 514   2.604  0.0256
##  Gender2 - Gender3     5.69 2.921 514   1.949  0.1261
## 
## Results are averaged over the levels of: PA_level 
## P value adjustment: tukey method for comparing a family of 3 estimates
#Posthoc test in personal achievement by PA group
posthoc_PA_PA
## $lsmeans
##  PA_level lsmean   SE  df lower.CL upper.CL
##  High      32.32 1.18 514     30.0     34.6
##  Low       28.75 1.22 514     26.4     31.1
##  Mod       31.37 1.06 514     29.3     33.5
##  NA         7.74 8.67 514     -9.3     24.8
## 
## Results are averaged over the levels of: Gender 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast   estimate    SE  df t.ratio p.value
##  High - Low    3.565 1.127 514   3.164  0.0089
##  High - Mod    0.944 0.933 514   1.012  0.7426
##  High - NA    24.572 8.649 514   2.841  0.0241
##  Low - Mod    -2.620 0.957 514  -2.739  0.0322
##  Low - NA     21.007 8.662 514   2.425  0.0736
##  Mod - NA     23.628 8.636 514   2.736  0.0325
## 
## Results are averaged over the levels of: Gender 
## P value adjustment: tukey method for comparing a family of 4 estimates