Affectionate touch and diurnal oxytocin levels: An ecological momentary assessment study

  1. Ekaterina Schneider
  2. Dora Hopf
  3. Corina Aguilar-Raab
  4. Dirk Scheele
  5. Andreas B Neubauer
  6. Uta Sailer
  7. René Hurlemann
  8. Monika Eckstein  Is a corresponding author
  9. Beate Ditzen  Is a corresponding author
  1. Institute of Medical Psychology, Center for Psychosocial Medicine, Heidelberg University Hospital, Germany
  2. Heidelberg University, Germany
  3. Department of Social Neuroscience, Faculty of Psychology, Ruhr University Bochum, Germany
  4. Department for Education and Human Development, DIPF|Leibniz Institute for Research and Information in Education, Germany
  5. Center for Research on Individual Development and Adaptive Education of Children at Risk, Germany
  6. Department of Behavioural Medicine, Faculty of Medicine, Institute of Basic Medical Sciences, University of Oslo, Norway
  7. Department of Psychiatry, University of Oldenburg, Germany
4 figures, 3 tables and 1 additional file

Figures

Flowchart of the recruitment process.

Figure 1 depicts the recruitment stages of both the online and the ecologically momentary assessments (EMA) study. Participants were recruited between April 1 and July 30, 2020, via online media and local newspapers. Inclusion criteria: fluency in German, minimum age of 18 y, and willingness to participate voluntarily. In total, 1483 individuals agreed to participate, of whom 1050 participants filled out the online questionnaires of interest. Out of the 472 participants who were interested in the EMA study, 247 finished the assessments.

Figure 2 with 1 supplement
Diurnal oxytocin and cortisol trajectories.

Panels (A) and (B) illustrate the daily oxytocin (pg/ml) and cortisol (ng/ml) trajectories across 2 d and all participants. Gray area indicates cortisol and oxytocin awakening response. Error bars indicate 95% confidence intervals.

Figure 2—figure supplement 1
Diurnal oxytocin levels depending on relationship status (A) and living arrangements (B).
Associations between occurrence of affectionate touch and psychological and hormonal state.

Panels (A) to (G) illustrate violin plots with density distributions of subjective ratings of general and Covid-19-related burden, stress, anxiety, happiness, cortisol, and oxytocin, depending on whether touch occurred or not. Each dot represents one assessment. Central dots (black) represent each mean. Black lines represent the standard deviations. * indicates statistically significant results (p<0.05). +indicates a statistical trend (p<0.1). ♠ indicates statistically significant between-person effect. ♣ indicates statistically significant within-person effect.

Associations between momentary affectionate touch intensity and psychological and hormonal state.

Panels (A) to (G) illustrate the results of random intercept and fixed slopes models depicting associations of momentary intensity of touch with self-reported general burden and Covid-19 related burden, stress, anxiety, happiness, cortisol, and oxytocin. Gray lines indicate the overall predicted slope, whereas the blue lines indicate the individual’s predicted slopes with their minimum and maximum predicted values as endpoints. Gray areas depict the 95% confidence band. * indicates statistically significant results (p<0.05).+indicates a statistical trend (p<0.1). ♠ indicates statistically significant between-person effect. ♣ indicates statistically significant within-person effect.

Tables

Table 1
Sample characteristics of online survey and ecological momentary assessment.
Sample characteristics of online survey participants
Men (n = 227)Women (n = 815)Non-binary (n = 4)Missing (n = 4)
M (SD)M (SD)M (SD)M (SD)
Age (years)34.67 (15.18)36.74 (14.58)45.50 (24.73)40.25 (15.39)
General psychological distress*10.22 (6.75)13.18 (7.49)21.00 (5.89)19.50 (14.66)
Anxiety5.15 (3.68)6.95 (4.12)10.50 (3.11)9.25 (7.68)
Depression 5.07 (3.55)6.22 (4.10)10.50 (3.32)10.25 (8.18)
Loneliness §37.18 (10.15)39.33 (10.95)53.00 (13.24)47.50 (19.50)
Attitude toward social touch 33.58 (10.18)34.76 (12.27)38.50 (23.39)40.75 (16.92)
Sample characteristics of ecological momentary assessment participants
Men (n = 74)Women (n = 173)
MSDMSD
Age (years)30.9913.6233.0512.41
Cortisol (ng/ml) **8.402.028.682.31
Oxytocin (pg/ml)**176.12106.18164.5496.74
Covid-19-related burden**36.9824.6141.7823.49
General burden39.9623.8947.2021.78
Stress levels††29.4915.8835.6217.00
Anxiety levels††18.3915.9824.1420.08
Happiness levels††71.1317.0967.8718.42
Intensity of affectionate touch††65.2120.0056.5723.13
  1. Table depicts means (M) and standard deviations (SD). Number of participants indicated as (n).

  2. *

    Hospital Anxiety and Depression Scale (HADS total score).

  3. HADS Anxiety subscale.

  4. HADS Depression subscale.

  5. §

    University of California, Los Angeles Loneliness Scale (UCLA Loneliness).

  6. Social Touch Questionnaire (STQ).

  7. **

    Out of 2964 possible data points, n = 2724 remained for analysis after excluding outliers, samples that were not stored as instructed or below detection limit, sampling problems.

  8. ††

    Momentary self-reported state.

Table 2
Results of the associations between affectionate touch and self-reported psychological affective states.
‍(A) Random intercept and random slopes models
EffectsGeneral burdenCovid-19 burdenStressAnxietyHappiness
Fixed effectsWithin-person
 Intercept47.734 (4.668); p<0.00143.145 (5.009); p<0.00134.431 (3.596); p<0.00121.096 (3.958); p<0.00171.080 (3.778); p<0.001
 ‍Touch*0.462 (1.561); p=0.767–2.660 (1.424); p=0.062–4.187 (1.994); p=0.036–0.217 (1.510); p=0.8861.557 (1.599); p=0.331
‍Between-person
 Touch*–5.560 (3.791); p=0.144–7.478 (4.068); p=0.067–7.534 (2.907); p=0.010–1.483 (3.210); p=0.64512.420 (3.068); p<0.001
Covariates
 ‍Age–0.186 (0.122); p=0.128–0.158 (0.131); p=0.228–0.136 (0.094); p=0.150–0.090 (0.103); p=0.3880.051 (0.099); p=0.605
 ‍Sex4.709 (3.325); p=0.1583.761 (3.565); p=0.2934.986 (2.524); p=0.0506.568 (2.814); p=0.021–3.318 (2.681); p=0.217
 Day–2.196 (0.864); p=0.011–2.602 (0.981); p=0.008–4.189 (1.228); p<0.001–3.145 (0.875); p<0.0011.637 (0.901); p=0.070
‍Random effects (SD)
 ‍Intercept21.48722.79413.69417.64716.597
 ‍Touch*9.1410.6428.0947.9068.945
 Residual12.82114.86418.74113.07613.478
(B) Random intercept and fixed slopes models
EffectsGeneral burdenCovid-19 burdenStressAnxietyHappiness
Fixed effects
Within-person
 Intercept44.439 (6.130); p<0.00139.748 (6.114); p<0.00132.966 (4.609); p<0.00124.466 (5.277); p<0.00171.437 (4.626); p<0.001
 ‍Touch intensity–0.077 (0.028); p=0.008–0.068 (0.036); p=0.058–0.148 (0.044); p<0.001–0.065 (0.029); p=0.0260.085 (0.030); p=0.005
‍Between-person
 ‍‍Touch intensity–0.121 (0.086); p=0.163–0.138 (0.087); p=0.115–0.223 (0.067); p=0.001–0.102 (0.074); p=0.1710.314 (0.066); p<0.001
Covariates
 Age–0.031 (0.161); p=0.847–0.082 (0.161); p=0.610–0.090 (0.125); p=0.475–0.158 (0.143); p=0.2700.007 (0.122); p=0.952
 ‍Sex2.647 (4.253); p=0.5354.143 (4.214); p=0.3272.650 (3.069); p=3895.932 (3.645); p=0.106–0.030 (3.171); p=0.993
 ‍‍Day–3.788 (1.056); p<0.001–4.930 (1.302); p<0.001–4.791 (1.583); p=0.003–3.718 (1.078); p<0.0013.695 (1.112); p=0.001
Random effects (SD)
 Intercept22.46721.60113.20618.62215.842
 ‍Residual11.81214.62618.30012.12112.591
  1. Table depicts coefficients (standard errors in parentheses) and p-values of associations between (A) the presence and (B) intensity of affectionate touch and psychological variables. Number of observations = 593–1023, Number of participants 162–227.

  2. *

    0 = no, 1 = yes.

  3. 0 = male, 1 = female.

Table 3
Results of the associations between affectionate touch and hormonal levels.
(A) Random intercept and random slopes models(B) Random intercept and fixed slopes models
EffectsCortisolOxytocinEffectsCortisolOxytocin
Fixed effectsWithin-personFixed effectsWithin-person
Intercept2.941 (0.165); p<0.0014.973 (0.405); p<0.001Intercept2.744 (0.370); p<0.0014.657 (0.798); p<0.001
‍Touch*–0.019 (0.060); p=0.756–0.030 (0.076); p=0.688Touch intensity–0.001 (0.001); p=0.3670.006 (0.002); p=0.003
‍Between-personBetween-person
Touch*–0.121 (0.057); p=0.036–0.145 (0.147); p=0.329Touch intensity–0.001 (0.002); p=0.5040.002 (0.003); p=0.489
‍CovariatesCovariates
‍Age–0.001 (0.002); p=0.633–0.013 (0.005); p=0.011Age–0.001 (0.003); p=0.776–0.019 (0.007); p=0.007
‍Sex–0.020 (0.044); p=0.647–0.129 (0.123); p=0.293Sex–0.001 (0.063); p=0.989–0.167 (0.151); p=0.272
‍Day–0.050 (0.035); p=0.155–0.010 (0.059); p=0.863Day–0.032 (0.051); p=0.5280.077 (0.086); p=0.370
‍Time-fall –0.476 (0.023); p<0.0010.051 (0.035); p=0.145Time-fall –0.471 (0.032); p<0.0010.025 (0.049); p=0.611
‍Body mass index–0.014 (0.005); p=0.0090.012 (0.015); p=0.422Body mass index0.005 (0.009); p=0.5850.014 (0.022); p=0.528
‍Eating*–0.077 (0.075); p=0.3020.032 (0.117); p=0.788Eating*–0.125 (0.099); p=0.2100.099 (0.149); p=0.506
‍Drinking*–0.007 (0.081); p=0.9330.045 (0.126); p=0.722‍Drinking*–0.011 (0.108); p=0.919–0.039 (0.161); p=0.810
‍Caffeine*0.127 (0.043); p=0.004–0.099 (0.070); p=0.158‍Caffeine*0.131 (0.064); p=0.0430.157 (0.102); p=0.125
‍Alcohol*–0.030 (0.063); p=0.628–0.174 (0.096); p=0.070‍Alcohol*–0.012 (0.072); p=0.865–0.154 (0.112); p=0.171
‍Cigarettes*0.104 (0.072); p=0.1530.053 (0.139); p=0.705Cigarettes*0.089 (0.098); p=0.368–0.054 (0.179); p=0.766
‍Physical activity*0.039 (0.040); p=0.3320.108 (0.065); p=0.098‍Physical activity*–0.072 (0.056); p=0.1990.014 (0.086); p=0.873
‍Sleep duration §–0.003 (0.004); p=0.371–0.009 (0.008); p=0.272‍Sleep duration §–0.029 (0.029); p=0.3210.046 (0.060); p=0.438
‍Sleep quality –0.001 (0.001); p=0.289–0.000 (0.002); p=0.912‍Sleep quality –0.001 (0.001); p=0.477–0.000 (0.002); p=0.960
Problem falling asleep*0.031 (0.053); p=0.552–0.245 (0.101); p=0.015Problem falling asleep*0.030 (0.081); p=0.709–0.133 (0.148); p=0.370
Sleeping pills*0.032 (0.109); p=0.7720.031 (0.254); p=0.904Sleeping pills*–0.203 (0.290); p=0.4870.771 (0.453); p=0.093
Forced awake*–0.008 (0.041); p=0.8380.012 (0.088); p=0.894Forced awake*0.016 (0.058); p=0.7880.075 (0.117); p=0.520
Brushing teeth*0.045 (0.036); p=0.2170.096 (0.056); p=0.090Brushing teeth*0.029 (0.051); p=0.5680.016 (0.077); p=0.839
Random effects (SD)Random effects (SD)
‍Intercept0.1490.542Intercept0.1790.546
‍Touch*0.2730.166Touch intensity
‍Residual0.3340.477Residual0.3390.461
  1. Table depicts unstandardized coefficients (standard errors in parentheses) and p-values of hormonal associations with (A) the presence and (B) intensity of affectionate touch. Number of observations = 251–545. Number of participants = 88–152.

  2. *

    0 = no, 1 = yes.

  3. 0 = male, 1 = female.

  4. 0 = time point 1–3, 1 = time point 4, 2 = time point 5, 3 = time point 6.

  5. §

    In hours.

  6. 1 = very bad, 101 = very good.

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ekaterina Schneider
  2. Dora Hopf
  3. Corina Aguilar-Raab
  4. Dirk Scheele
  5. Andreas B Neubauer
  6. Uta Sailer
  7. René Hurlemann
  8. Monika Eckstein
  9. Beate Ditzen
(2023)
Affectionate touch and diurnal oxytocin levels: An ecological momentary assessment study
eLife 12:e81241.
https://doi.org/10.7554/eLife.81241