Home Weight Loss and Sexual Activity in Adult Obese Individuals: Establishing a Causal Link
Article Open Access

Weight Loss and Sexual Activity in Adult Obese Individuals: Establishing a Causal Link

  • Harald Tauchmann ORCID logo EMAIL logo and Ansgar Wübker
Published/Copyright: August 15, 2023

Abstract

Obesity may not only be linked to undesirable health outcomes but also to limitations in sexual life. The present paper aims to assess whether there is a causal relationship between weight loss and sexual activity in adult obese individuals. To address the endogeneity of weight loss that is likely to result in biased estimation results, the analysis is based on data from a randomized field experiment. In this experiment, financial weight-loss rewards were offered to a random subgroup of participants and can be used as an exogenous source of weight variation in an instrumental variables approach. Estimation results indicate that for obese males losing weight, both the probability of being involved in a sexual relationship and the frequency of sexual intercourse increases. The same pattern of results is found when pooling both sexes in the sample. Due to the small share of women in the sample, the analysis yields no reliable results specifically for females.

JEL Classification: I10; I18; C93

1 Introduction

A close link between a fulfilling sex life and general life satisfaction is well documented in the literature (e.g. Woloski-Wruble et al. 2010; Schafer et al. 2013; Soysal and Smith 2022). A growing body of literature also demonstrates that a higher frequency of sexual activity is associated with better health. For example, Kloner et al. (2016) and Soysal and Smith (2022) document reductions in cardiovascular events later in life, lower risk of fatal coronary events, prostate and breast cancer, and better quality of life. Impairments in sexual life thus are likely to have a strong negative impact on individual happiness and life satisfaction and likewise deteriorate the general health condition.

One possible reason for limitations in sexual life is obesity; see Larsen et al. (2007) and Sarwer et al. (2018) for reviews of the medical literature focussing on the linkage of obesity and sexual dysfunction such as erectile dysfunction. Though the majority of analyses reviewed by Larsen et al. (2007) and Sarwer et al. (2018) find a positive connection, the evidence is still mixed.[1] Taking a less physiologically oriented perspective and considering deficits in sexual life beyond sexual dysfunctions, Kolotkin et al. (2006) list ‘lack of enjoyment’, ‘lack of desire’, ‘difficulties with performance’, and ‘avoidance of sexual encounters’ as examples of obesity related impairments of sexual life. Rowland et al. (2017) state that the relationship between obesity and sexual function is complex and multivariate, with at least three different likely pathways: direct effects of adipose tissue, effects of comorbidities, and psychologically driven effects.[2] In general, besides undesirable general health outcomes, obesity may have a more direct detrimental effect on happiness and life satisfaction that operates through an unsatisfactory sex life.

The non-medical, health economics oriented literature has taken a very different view on the link of body weight to sexual life. This relatively small literature (Ali et al. 2014; Cawley 2001; Cawley et al. 2006; Halpern et al. 2005; Neymotin and Downing-Matibag 2014; Sabia and Rees 2011) focusses on individuals of very young age, adolescents in particular. Due to the focus on this special age group, sexual activity is associated with premature initiation to sex, potentially exerting detrimental effects on adolescents’ later lives. For female adolescents, Sabia and Rees (2008) find a causal and detrimental effect of early sexual intercourse on psychological well-being. Ali et al. (2014) allude to disease and pregnancy risks teenagers take when having early sex, suggesting that a postponed initiation to sexual activities might be a beneficial sideeffect of teenage overweight.

The contribution of the present paper is to (i) conduct an analysis focusing on adult individuals, for whom the above line of argument does not apply and whose sexual behavior, most likely, differs from the behavior of the age group studied by the literature cited above. Unlike the majority of the medical literature, our focus is (ii) on (self-reported) sexual activity rather than specific sexual dysfunctions such as erectile dysfunction. This is important as obesity related limitations in sex life may often be unrelated to any physiological deficit but can possibly be attributed to social and psychological body-weight related factors, such as being less attractive to a potential partner, a lack of self-confidence when dating, or less enjoyment in and desire for sex. Finally, the present analysis contributes to the existing evidence by (iii) establishing a causal effect of weight loss on sexual activity as opposed to a mere correlation in the data.

Causality in the link between body weight and sex life has not always attracted much attention. Halpern et al. (2005), for instance, estimate that the probability for adolescent girls of being involved in a romantic relationship (with or without sexual intercourse) increases by 6–7 percent with one BMI unit less. Using Add Health (National Longitudinal Survey of Adolescent Health; Harris 2013) data,[3] a US longitudinal school-based survey, Cawley et al. (2006) find that the odds of initiation to sexual intercourse for obese girls are just 32 percent of those of normal-weight girls. However, both analyses do not account for unobserved heterogeneity. Hence, the striking results may well be attributed to unobserved confounding factors, such as personal character traits being relevant for both obesity related behavior (eating habits, physical activity) and mating behavior. Numerous medical studies (see Larsen et al. 2007; Sarwer et al. 2018) likewise focus on the association of overweight and impairments in sex life rather than on a causal relationship.

More recent papers address possible endogeneity by using instrumental variables (IV) estimation. Analyzing data from Add Health and using siblings’ BMI and mother’s obesity status as instruments for the respondents’ body weight, Sabia and Rees (2011) find that for girls aged 14 to 17 one BMI unit more reduces the probability of sexual activity by 3.5 percentage points. In contrast, no significant effect is found for boys. Rather, the point estimate for boys bares the opposite sign in several specifications. Ali et al. (2014) take a similar approach by instrumenting body weight with maternal obesity status. For white girls – but not for blacks – they find that higher body weight or obesity significantly lowers the probability of having had sex or having been involved in a romantic relationship.

Though both papers carefully discuss the validity of the instruments used, one may still doubt whether mothers’ – and other relatives’ – body weight is a valid instrument: body weight as well as sexual behavior – or at least attitudes towards it – may be intergenerationally transmitted from mother to child (e.g. Taris 2000). If so, maternal body weight could influence children’s sex life through mothers’ own sexual behavior, given that body weight matters for sexual behavior. Another concern with respect to instrument validity is the social environment children share with their mothers and siblings. It may strongly influence both, body weight and attitudes towards sex for both mothers and children and, most likely, cannot entirely be controlled for by including covariates in the analysis.

While following Ali et al. (2014) and Sabia and Rees (2011) in relying on instrumental variables estimation, we do not use maternal body weight as instrument. Rather, in our analysis identification is based on exogenous variation in BMI having been artificially induced in a randomized field experiment. More precisely, we use data (Augurzky et al. 2017) originating from a randomized trial that was conducted to evaluate the effectiveness of financial incentives for weight loss (Augurzky et al. 2012, 2018; Eilers and Pilny 2018). In this experiment, individuals in the intervention groups were exposed to financial incentives for losing body weight, while control group members were not. Following Hafner et al. (2021) and Reichert (2015), treatment status can serve as an instrument for the change in body weight. Besides using different instruments and focussing on adult individuals, our analysis materially deviates from Ali et al. (2014) and Sabia and Rees (2011) by only considering obese[4] individuals. That is, we are interested in whether – even moderate – weight loss in obese individuals makes a difference to their sexual lives.

The identification strategy of the present analysis is related to the approach of Esposito et al. (2004). This medical study also relies on data from a randomized field experiment with obese participants. In this study, the intervention was not exposition to financial incentives but guidance on how to lose weight combined with weight monitoring through monthly group sessions. Esposito et al. (2004) estimate a reduced form effect by comparing post-intervention IIEF (International Index of Erectile Function) scores between the intervention and the control group. The statistically significant advantage in IIEF found for the intervention group nevertheless provides striking evidence for body weight mattering for sexual dysfunctions in obese persons. However, by focusing on erectile dysfunction alone, Esposito et al. (2004) may miss out on other obesity related impairments of sexual life. By considering two different measures of sexual activity, our analysis takes a broader perspective on the link between obesity and sex life.

The remainder of this paper is organized as follows. Section 2 introduces the data and describes the underlying experiment. Section 3 discusses the econometric model and Section 4 presents the estimation results for the basic model. In Section 5, alternative model specifications and estimation results are discussed. Section 6 summarizes and discusses the main findings.

2 Data

2.1 The Experiment

The data (Augurzky et al. 2017) used in this analysis originate from a large scale, four-phase, randomized field experiment that is described in more detail elsewhere (Augurzky et al. 2012, 2018; Eilers and Pilny 2018; Hafner et al. 2021). The experiment was conducted between March 2010 and July 2013. Roughly 700 obese[5] individuals were recruited in four rehabilitation hospitals in Germany (For details about rehab in Germany, see e.g. Ziebarth 2014; Krauth and Bartling 2017). After a typically three weeks stay, but shortly before discharge, the physician in charge set all participants an individual weight-loss target that they were prompted to realize within four months. Subsequently, the participants were randomly assigned to two incentive groups (incentive 1501–4 and incentive 3001–4)[6] and one control group (incentive 01–4). Covariate balancedness across the groups is documented in Augurzky et al. (2018, Table 2). Members of the former two groups were promised € 150 and € 300, respectively, for realizing (or exceeding) their individual weight-loss target within four months, while members of the control group were not exposed to any financial incentives for reducing body weight.

The weigh-in at the end of the weight-reduction phase – and all further weigh-ins – were carried out at assigned pharmacies. As the pharmacy was selected by the experimenter and not by the test person, the location of weigh-in was exogenous to the participants. Shortly after the weigh-in, all participants – irrespective of group membership in the weight-reduction phase and irrespective of success – were prompted to comply with their target weight by the end of a subsequent six months weight-maintenance phase.[7] Yet, conditional on success in the weight reduction phase, a second randomization took place at the same time. Participants who had lost at least 50 percent of their contractual weight-loss target were randomly assigned to two intervention groups (incentive 2505–10, incentive 5005–10) and one control group (incentive 05–10). The former two could gain up to € 250 and € 500, respectively, if they achieved their target weight or less ten months after rehab discharge. Members of the control group were not promised any rewards. Individuals who were not successful in the weight-reduction phase were hence effectively assigned to the control group. Six months later, i.e. ten months after rehab discharge, another weigh-in took place under identical conditions. Yet, at this point in time, the incentive scheme was terminated and no further randomization took place. However, all participants were still requested to comply with their target weights and another weigh-in was announced to take place one year later. This final year serves as the follow-up phase. See Figure 1 for the time line of the entire experiment.

Figure 1: 
Time line of the four-phase experiment. N indicates sample sizes at the four points in time at which data was collected. The differences in N (177, 106, and 96) are the number of drop-outs in the respective phase of the experiment. Source: Augurzky et al. (2018), subject to minor modifications.
Figure 1:

Time line of the four-phase experiment. N indicates sample sizes at the four points in time at which data was collected. The differences in N (177, 106, and 96) are the number of drop-outs in the respective phase of the experiment. Source: Augurzky et al. (2018), subject to minor modifications.

By the end of each of the four phases, the participants had to answer a detailed questionnaire. The final one, being sent 22 months after rehab discharge, is the only one to include questions about the participants’ sex life. Yet, only for participants with even identification numbers – i.e. for a fifty percent random sample[8] – the questionnaire included these touchy questions. The experimenters wanted to reduce the risk of too many individuals dropping out by asking them – possibly displeasing – questions about their sex lives. However, this concern proved immaterial as the drop-out rate virtually did not differ between individuals with even and with odd identification numbers.

The experiment population was subject to significant sample attrition. 695 individuals started the weight-reduction phase after rehab discharge. 177 participants dropped out during the weight-reduction phase, i.e. they did not show up at the weigh-in after four months, leading to only 518 individuals entering the weight-maintenance phase. In this phase, another 106 individuals dropped out. The follow-up phase was also subject to some sample attrition reducing the number to 316 individuals for whom weight information is available by the end of the experiment. Among these, 174 were asked questions regarding sex life. Only 17 denied such information. This corresponds to a rate of item-non-response of less than 10 percent, which is remarkably low compared to many surveys addressing sexual behavior (cf. Fenton et al. 2001). Attrition usually occurred in terms of just not showing up at the next scheduled weigh-in.

2.2 Outcome Variables

In the final questionnaire 22 months after rehab discharge, the participants were asked two questions regarding their sex lives, (i) if they were involved in a sexual relationship and (ii) how frequently they had sexual intercourse. Both questions referred to the previous twelve months, in other words to the follow-up phase of the experiment. While the first question was a yes-no question, the latter allowed for five different answers: ‘never’, ‘occasionally’, ‘at least once a month’, ‘at least once a week’, ‘(almost) daily’. Hence, we can use two different measures of sexual activity as dependent variables, a binary indicator (sexpartner11–22) and an ordered categorial one (sexfrequency11–22). 156 and 152 participants, respectively, answered these questions. For estimation, we reduced the number of categories for the latter variable to three (‘never’, ‘occasionally’, ‘regularly’) with the third category combining having sex monthly, weekly or (almost) daily. Table 1 displays the joint distribution of sexpartner11–22 and sexfrequency11–22 (full set of categories). None of the respondents who report they have not had sex in the last year report they have had a sexual partnership in the same period. The data also indicates that the respondents distinguish between living together with a partner and having a sexual relationship.

Table 1:

Joint and marginal distribution of dependent variables.

frequency of sexual intercourse Marginal distribution+
Never Occasionally Monthly Weekly Daily All Males Females
Sexual relationship: no 32 9 1 0 0 45 20 25
 Yes 0 25 32 48 4 111 83 28
Marginal distribution+: all 32 35 33 48 4 157 103 53
 Males 11 26 22 37 3 99
 Females 21 9 11 11 1 53
  1. Notes: +Due to item non-response, values do not exactly sum up. 157 responses to (at least) one of the two questions regarding sex life; 156 responses to question about sexual relationship; 152 responses to question about frequency of intercourse.

2.3 Explanatory and Instrumental Variables

In our preferred model, we use the change in BMI between the start of month 5 and the end of month 10 (weight-maintenance phase) as the key explanatory variable (denoted ΔBMI5–10).[9] It is, therefore, measured earlier in time than the outcome variables.

The key concern to the present analysis is the possible endogeneity of both body weight and its change over time. First, unobserved individual characteristics may have an effect on either variable. One may hypothesize that individuals with high self-esteem are more successful in finding a partner and are also less vulnerable to overeating, obesity and failure to stick to intentions. Self-confidence and self-esteem are possibly decisive for being successful when seriously trying to become less obese. Preferences for physical pleasures represent a further channel through which individual heterogeneity may generate a spurious correlation. If such preferences are strong, one is likely to have more desire for both sex and calorie intake. Second, reverse causality might also be a source of endogeneity bias. One such channel is overeating induced by sexual frustration. Yet, one may also argue that being satisfied with one’s sex life reduces the pressure to comply with ideals of beauty, rendering obesity less costly in psychological terms.

Given the lack of information about how the respondents’ sex life develops over time, instrumental variables estimation is the first choice for identifying a causal relation. As the data used in this analysis originates from a randomized field experiment, it includes promising variables for being used as instruments. More precisely, these variables are incentive 1501–4, incentive 3001–4, incentive 2505–10, and incentive 5005–10, i.e. the indicators for randomly assigned group membership in the experiment’s two intervention phases. In principle, exogeneity is guaranteed by the experimental design. Moreover, substantial explanatory power of these indicators for weight change has been established elsewhere in the literature (Augurzky et al. 2012, 2018; Hafner et al. 2021; Reichert 2015). Yet, since we are interested in weight change in temporal proximity to the period for which information on sex life is available, the second intervention (incentive 2505–10, incentive 5005–10) appears to be the more promising exogenous source of variation in weight loss;[10] see also footnote 9.

With respect to the second incentive, group membership – though randomly assigned – is not purely exogenous since only successful weight losers had a chance of being assigned to one of the incentive groups. In our preferred specification, we thus condition on having qualified for the second randomization and use only incentive 2505–10 and incentive 5005–10 as instruments for Δ BMI5–10. This guarantees – at the cost of a reduced sample size – randomness of group membership in the estimation sample.

Yet, another threat to identifying the effect of treatment remains, which is selective attrition from the experiment. We address this issue by tabbing an exogenous source of variation in attrition. That is, we use an exogenous component of the (travel) costs of attending the weigh-in as an instrument for attrition. More specifically, a dummy (pharmacy nearby1–22), indicating that the place of the weigh-in is located within the same zip-code area as the participant’s place of residence, enters the model as instrument for attrition from the experiment. Since the pharmacy to be visited for the weigh-ins was assigned by the experimenter, this variable is almost as good as randomly assigned. Moreover, it proves to have substantial explanatory power for showing-up at the weigh-ins.[11]

2.4 Control Variables

The small sample size precludes specifying a rich regression model with numerous controls. Hence, only age and sex enter the model as basic socioeconomic controls. In addition, BMI at the end of the weight-loss phase (BMI4) and an indicator for living together with a partner at that point in time (together4) enter the regression model. We include the former to capture the effect of pre-intervention body weight which is most likely endogenous and for which no instruments are available. The latter we include as the closest substitute available for pre-intervention values of the dependent variable, which is not observed in the data. Table 2 displays key descriptives for the explanatory variables and the instruments in the sample used for estimating our preferred specification. Besides descriptives for the pooled sample, we also provide them for the sub-sample of men.[12]

Table 2:

Descriptive statistics for estimation sample.

Mean S.D. Median Min. Max.
Males
Dependent variables:

 sexpartner11–22 0.815 0.391 1.000 0.000 1.000
 sexfrequency11–22 1.476 0.715 2.000 0.000 2.000
  t o g e t h e r 22 + 0.682 0.469 1.000 0.000 1.000

Explanatory variables:

 Δ BMI1–4+ −2.734 1.651 −2.196 −7.748 −0.906
 Δ BMI1–10+ −2.502 2.475 −2.179 −15.347 1.614
 Δ BMI5–10 0.232 1.711 0.311 −8.424 3.720
 age 50.500 7.917 51.000 21.000 63.000
 together0+ 0.742 0.441 1.000 0.000 1.000
 BMI0+ 35.177 5.058 34.273 28.441 50.039
together 4 0.697 0.463 1.000 0.000 1.000
 BMI4 32.443 4.438 31.997 22.857 43.226

Instrumental variables:

 incentive 1501–4+ 0.273 0.449 0.000 0.000 1.000
 incentive 3001–4+ 0.439 0.500 0.000 0.000 1.000
 incentive 2505–10 0.348 0.480 0.000 0.000 1.000
 incentive 5005–10 0.394 0.492 0.000 0.000 1.000
 pharmacy nearby122 0.621 0.489 1.000 0.000 1.000
Males & females
Dependent variables:

 sexpartner11–22 0.763 0.428 1.000 0.000 1.000
 sexfrequency11–22 1.400 0.791 2.000 0.000 2.000
  t o g e t h e r 22 + 0.653 0.478 1.000 0.000 1.000

Explanatory variables:

 Δ BMI1–4+ −2.660 1.488 −2.108 −7.748 −0.906
 Δ BMI1–10+ −2.425 2.434 −2.236 −15.347 2.384
 Δ BMI510 0.235 1.746 0.238 −8.424 4.147
 age 50.061 8.705 51.000 21.000 63.000
 female 0.327 0.471 0.000 0.000 1.000
 together0+ 0.704 0.459 1.000 0.000 1.000
 BMI0+ 35.670 4.897 34.732 28.441 50.039
 together4 0.673 0.471 1.000 0.000 1.000
 BMI4 33.010 4.524 32.387 22.857 48.111

Instrumental variables:

 incentive 1501–4+ 0.327 0.471 0.000 0.000 1.000
 incentive 3001–4+ 0.378 0.487 0.000 0.000 1.000
 incentive 2505–10 0.296 0.459 0.000 0.000 1.000
 incentive 5005–10 0.439 0.499 0.000 0.000 1.000
 pharmacy nearby1–22 0.643 0.482 1.000 0.000 1.000
  1. Notes: Statistics for 98 individuals who provide information on sexpartner11–22 or sexfrequency11–22 and qualified for the second randomization (i.e. observations entering the outcome equation of the preferred specification). +Variable used in alternative model specifications. Descriptive statistics for the women sub-sample reported in the Appendix Table 5.

3 Estimation Procedure

The econometric model rests on formulating a likelihood function that jointly considers the three model components (outcome equation, ‘first stage equation’ for weight loss, selection equation) and hence allows for simultaneous estimation by maximum likelihood (ML). The estimation approach incorporates the assumption of joint normality, which makes us use a probit model with an endogenous regressor (Wooldridge 2002, p. 476) augmented by a Heckman-type selection equation. Though the model components are estimated simultaneously, the number of observations effectively contributing to the identification may still be heterogeneous across the components. For more details on simultaneous estimation of multi-equation models with jointly normally distributed errors, see Roodman (2011). This paper also introduces an accompanying software component.[13] Simple two-stage least squares combined with the selection correction proposed by Olsen (1980) is a fully linear, multi-step – arguably simpler, yet rather inefficient – alternative to joint ML estimation. However, the ordered categorial nature of the outcome sexfrequency11–22 renders a fully linear model questionable for the present application, making us stick to ML estimation.[14]

Weak identification has become a major concern in instrumental variables estimation. The usually reported F-statistic of the test on the instruments – and the corresponding rule-of-thump threshold value of 10 – do not, however, have a statistical basis in the context of non-linear IV-models (Denzer and Weiser 2021). Therefore, we follow the procedure as suggested in Denzer and Weiser (2021) and report the Shapiro–Wilk W-statistic as a measure of instrument strength not limited to linear models. The value of this statistic is compared to the critical values (0.994 for 5 percent relative bias, 0.952 for 10 percent relative bias) provided by Denzer and Weiser (2021, see p. 28) in the context of a probit model with endogenous regressors.[15]

4 Estimation Results for the Preferred Specification

In this section, we discuss estimation results for our preferred model specifications described above. Table 3 displays results for the model variant using sexpartner11–22 as dependent variable, while Table 4 displays results for the model explaining sexfrequency11–22. Since the determinants of sex life are likely to differ substantially between genders, we also estimated separate models for males and females. Yet, when using the rather small all-women sample, the ML estimator suffered from severely weak instruments (W-statistics of about 0.1) and convergence issues, which made us abstain from reporting results for that specification. On the contrary, with the outcome sexpartner11–22 the ShapiroWilk W-statistic exceeds the threshold value of 0.952, both for the pooled and the all-males sample, which argues against weak identification. This does not in the same way apply to the estimations considering the frequency of sexual intercourse as the outcome variable. Yet, the W-statistics still take rather large values of about 0.9; see Table 4.

Table 3:

Sexual relationship, preferred model (estimated coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexpartner11–22)

Δ BMI5–10 −0.359** 0.164 −0.322* 0.178
age 0.011 0.025 −0.000 0.012
female −0.059 0.225
together 4 0.260 0.299 0.482** 0.235
BMI 4 −0.062* 0.037 −0.058* 0.034
constant 1.258 1.802 1.436 1.201

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.399** 0.500 −1.090** 0.389
incentive 500 5–10 −1.068** 0.479 −0.750** 0.365
age 0.007 0.028 −0.013 0.018
female 0.112 0.334
together 4 −0.111 0.418 0.289 0.323
BMI 4 −0.032 0.030 0.012 0.021
constant 2.038 2.013 1.004 1.311

Selection equation (dependent variable: sexinfo22)

pharmacy nearby 1 22 0.880** 0.205 0.858** 0.166
incentive 250 5–10 −0.251 0.249 −0.010 0.205
incentive 500 5–10 0.141 0.252 0.162 0.196
age 0.020 0.014 0.014 0.011
female 0.110 0.177
together 4 −0.384+ 0.262 −0.305+ 0.200
BMI 4 −0.061** 0.027 −0.050** 0.019
constant 0.486 1.261 0.279 0.844

Ancillary parameters

σ (instrumental equ. error S.D.) 2.133** 0.099 2.167** 0.072
ρ bmi sex 0.689** 0.347 0.555 0.414
ρ info sex 0.691* 0.389 0.809** 0.302
ρ info bmi −0.047 0.144 −0.039 0.106
# of observations (over all) 212 307
# of observations (outcome equation) 65 97
log likelihood −515.8 −763.5
Shapiro–Wilk W-statistic (weak identification) 0.957 0.958
Sample split (p-value, LR-test) 0.390
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. All model equations estimated using only individuals who successfully lost body weight during month one to four. Simultaneous ML estimation assuming joint normality (binary probit, linear, binary probit). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. Test on sample split refers to the null model that pools males and females. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

Table 4:

Frequency of Intercourse, selection correction (estimated coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Main equation (dependent variable: sexfrequency11–22)

Δ BMI5–10 −0.316** 0.131 −0.331** 0.122
age −0.014 0.025 −0.015 0.012
female 0.013 0.210
together 4 0.356 0.318 0.563** 0.272
BMI 4 −0.054 0.039 −0.054+ 0.033

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.183** 0.495 −0.922** 0.379
incentive 500 5–10 −1.210** 0.473 −0.878** 0.355
age 0.003 0.028 −0.014 0.018
female 0.105 0.335
together 4 −0.110 0.432 0.260 0.321
BMI 4 −0.034 0.029 0.011 0.020
constant 2.290 1.948 1.089 1.279

Selection equation (dependent variable: sexinfo22)

pharmacy nearby 1 22 0.815** 0.220 0.843** 0.172
incentive 250 5–10 −0.085 0.280 0.126 0.209
incentive 500 5–10 0.213 0.266 0.125 0.202
age 0.022+ 0.014 0.014 0.010
female 0.129 0.181
together 4 −0.296 0.278 −0.229 0.207
BMI 4 −0.063** 0.028 −0.053** 0.019
constant 0.312 1.334 0.247 0.855

Ancillary parameters

ordered probit threshold 1 −2.776 1.952 −2.371* 1.346
ordered probit threshold 2 −2.174 1.815 −1.912+ 1.228
σ (instrumental equ. error S.D.) 2.124** 0.095 2.154** 0.069
ρ bmi sex 0.805** 0.186 0.747** 0.246
ρ info sex 0.243 0.371 0.434+ 0.296
ρ info bmi −0.061 0.149 −0.045 0.109
# of observations (over all) 212 307
# of observations (outcome equation) 63 95
log likelihood −547.2 −808.6
Shapiro–Wilk W-statistic (weak identification) 0.920 0.866
Sample split (p-value, LR-test) 0.047
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. All model equations estimated using only individuals who successfully lost body weight during month one to four. Simultaneous ML estimation assuming joint normality (ordered probit, linear, binary probit). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. Test on sample split refers to the null model that pools males and females. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

The estimates for the cross-equation error correlation ρ bmi sex are positive throughout and statistically significant (except for sexpartner11–22, pooled sample) suggesting that endogeneity of weight loss needs to be taken into account. From a non-technical perspective, this finding conflicts with our earlier reasoning that self-confidence is a major source of unobserved heterogeneity. It rather argues in favor of the desire for physical pleasures driving the correlation. For the estimates of ρ info sex the pattern is less clear. Though also positive and of substantial magnitude, they are only statistically significant (at conventional levels) for the models with binary outcome. But still, the rather strong estimated cross-equation error correlations suggest that selection cannot be ignored. For the outcome sexpartner11–22 pooling the sample over the sexes is clearly not rejected by a likelihood-ratio test. In line with this, the estimated coefficients exhibit little difference between the all-male and the pooled sample. The corresponding test for the outcome sexfrequency11–22 is at the margin of rejecting the pooled model. Yet, the estimated coefficients still do not deviate strikingly.

Turning to the question of what makes a sexual relationship more likely, Table 3 indicates that being involved in a sexual relationship is not significantly associated with age. Not surprisingly, those who lived with a partner by the end of the weight-reduction phase are more likely to have a sexual partner 6–18 months later (statistically significant only for the pooled sample). Model estimation also yields a significantly negative association with initial BMI. Yet, this coefficient does not isolate a causal relationship but is likely to capture the influence of unobserved individual heterogeneity. Most importantly, the coefficient of Δ BMI1–10 is statistically significant (in the pooled sample only at the 10 percent level) and negative for both the pooled and the all-males model. This indicates that – conditional on initial weight and initial relationship status – losing body weight increases the likelihood of being involved in a sexual relationship. The coefficient estimates correspond to marginal effects of substantial size. At its respective maximum value ( β ̂ Δ BMI ϕ ( 0 ) ) the corresponding marginal effects are as high as −0.143 (all male) and −0.128 (pooled), meaning one BMI unit reduction in body weight increases the probability of having a sexual partner by more than ten percentage points. The mean marginal effects in the estimation sample are of similar magnitude (−0.112 and −0.099, respectively).

To dig deeper into the interdependence of body weight and sexual activity, we re-estimated the model with an indicator for living together with a partner (together22) as an alternative dependent variable; see Table 6 and Figure 3 in the Appendix. The estimated coefficients were in absolute terms smaller and by all standards statistically insignificant (and of the opposite sign). In other words, it is not the relationship status per se, but its sexual nature that is significantly affected by losing body weight.

In line with this result, we find statistically significant and negative effects of weight gain on the outcome frequency of intercourse, too. Table 4 shows that the point estimates are very similar to their counterparts from using the indicator for a ‘sexual relationship’ as the outcome variable. This implies that – according to these estimates – losing one BMI unit has quite a big effect on the probability of having sex regularly. More specifically, for the most affected individual the effect amounts to 12.6 percentage points (all-males) and 13.2 percentage points (pooled), respectively.

To sum up, our preferred model shows that even moderate weight loss in obese individuals has remarkably strong effects on their sex lives, both at the extensive and the intensive margin. Nevertheless, this result is driven by male individuals. Albeit including females in the estimation sample has little impact on the estimation results, the rather small size of the sub-sample of women does not allow for a reliable analysis of the link between weight change and sex life separately for women.

5 Robustness Checks

We consider conditioning on success in the weight-reduction phase and simultaneously accounting for sample attrition as the most appropriate approach to dealing with the limitations of the data. This approach, however, results in estimating a rather complex model – which involves a trivariate normal error distribution – using a rather small sample. Therefore, we also estimate several alternative model specifications, which are technically simpler and/or allow for using a larger estimation sample.

More specifically: (i) We estimate the reference model just without a selection equation. (ii) We condition on success in the weight-loss phase not by restricting the sample to successful participants, but by including an indicator for success as control, which is instrumented by the weight-loss phase incentive indicators incentive 1501–4 and incentive 3001–4. This specification, which does not involve a selection equation either, allows for using a substantially larger sample than the reference. (iii) We do not condition on success but use the total weight change over both, the weight-loss and the weight-maintenance phase (Δ BMI1–10) as explanatory variable, which is instrumented by all four incentive indicators incentive 1501–4, incentive 3001–4, incentive 2505–10, and incentive 5005–10. This specification is rather simple as it does not involve a selection equation and allows using a relatively large sample compared to the reference model. (iv) We estimate a fully linear, one-to-one counterpart to the reference model combining two-stage-least squares with the selection correction suggested by Olsen (1980); cf. Section 3.

Figure 2 graphically displays the results from these alternative models for the coefficient/marginal effect of Δ BMI and contrasts them with the results of the preferred model. Full regression output is provided in the Appendix; see Tables 714. Starting with specification (i), Figure 2 (upper left panel) indicates that it makes little difference to the estimated effects of weight loss when attrition is not taken into account. Without a selection equation, we estimate slightly stronger effects on the probability of living in a sexual relationship and slightly smaller effects on the frequency of intercourse. Yet, the confidence intervals overlap heavily for the two specifications. One exemption is the pooled model specification. Here the effect on sexfrequency11–22 turns statistically insignificant at the 10 percent level for the simpler specification without selection equation. The pattern of weight-change coefficients found for specification (ii) exhibits almost exactly the same pattern found for specification (i); see Figure 2, upper right panel. That is, including a control indicating success in the weight-loss phase makes little difference compared to excluding unsuccessful participants. As before, without modeling attrition from the experiment in the sample that pools men and women, the effect on the frequency of intercourse turns insignificant. Yet, from a statistical perspective the estimated coefficients still hardly deviate from their counterparts in the reference specification. Things are different for specification (iii) (Figure 2, lower left panel) which does not condition on success in the weight-loss phase, considers the weight change over both intervention phases, and uses all incentive-group identifiers as instruments for the latter. Although in qualitative terms the results do not deviate from the reference model, the estimated coefficients become much smaller. For the outcome sexfrequency11–22, they are clearly statistically insignificant, though the confidence intervals are tighter than for the reference model (presumably because of the substantially larger sample used for estimation). Thus, specification (iii) questions the effect of weight change at least on the outcome frequency of intercourse. However, this specification uses instruments whose validity is questionable by the design of the experiment (cf. Section 2.3). The deviation from the reference may, hence, not indicate a lack of robustness but might just be due to using invalid instruments. In order to make the results of the fully linear multi-step model (iv) comparable to the reference, the lower right panel of Figure 2, depicts (average) marginal effects rather than coefficient estimates. In qualitative terms, the results are similar to their counterpart from the reference model. Yet, the estimated effects are bigger in absolute terms. The most striking difference to reference is, however, the much wider confidence intervals. This lack of precision can presumably be attributed to the rather inefficient multi-step estimation procedure. Based on the linear model, one can hardly make any reliable statement regarding the link between weight loss and sex life. In turn, the results from the linear model neither confirm nor challenge their counterparts from the non-linear specification of reference.

Figure 2: 
Grey: reference specification (Tables 3 and 4); black: (i) reference model without selection equation, (ii) model that includes an indicator for success in weight-loss phase as control, (iii) model that does not condition an success in weight-loss phase, (iv) 2SLS combined with linear Olsen (1980) selection correction. Spikes indicate 90 percent confidence intervals; confidence intervals bootstrapped for the linear model. Markers: ,  outcome sexpartner11–22; ,  outcome sexfrequency11–22. Source: own calculations.
Figure 2:

Grey: reference specification (Tables 3 and 4); black: (i) reference model without selection equation, (ii) model that includes an indicator for success in weight-loss phase as control, (iii) model that does not condition an success in weight-loss phase, (iv) 2SLS combined with linear Olsen (1980) selection correction. Spikes indicate 90 percent confidence intervals; confidence intervals bootstrapped for the linear model. Markers: , outcome sexpartner11–22; , outcome sexfrequency11–22. Source: own calculations.

6 Discussion and Conclusions

This paper analyzed whether weight reduction in adult obese individuals has an impact on their sex lives, namely whether or not they live in a sexual partnership and how frequently they have sex. Based on data generated by a randomized field experiment, where randomization provides an exogenous source of variation in body weight exploited for identification, we find that obese males are substantially more likely to be involved in a sexual relationship when losing weight. This effect is of substantial magnitude, amounting to about 14 percentage points for the most affected individual in our preferred model. In the preferred specification, the results for the frequency of sexual intercourse are very similar. Yet, in some alternative specifications, the effect on the frequency of intercourse is not statistically significant. This applies in particular to estimation samples not restricted to males. Moreover, weak identification seems to be more of an issue if the intensive margin of sex life is analyzed.

In contrast to those findings for men, the analysis does not yield any significant and reliable effect for women. This can be explained by the small number of women in our data that seriously hampers the empirical analysis. Another explanation might be gender differences in sex life. According to Baumeister et al. (2001) there is ample empirical evidence for men – as compared to women – having more frequent and more intense sexual desires. The strength of sexual motivation, also referred to as sex drive, is for instance reflected in the desired number of partners and the desired frequency of intercourse, which are closely related to the outcome variables of the present study. The stronger male sex drive may facilitate detecting a causal effect of weight loss for men. Based on an extensive literature survey, Petersen and Hyde (2010) state only small gender differences in many sexual behaviors. But casual sex and attitudes toward casual sex are among the exceptions to this general finding. Men report a more active sexual behavior or express more permissive attitudes towards them than women do. This diverging predisposition could be another explanation for the significant effects found for men but not for women, concerning the probability of being involved in a sexual relationship and the frequency of sexual intercourse.

For many individuals, a fulfilling sex life is an important facet of general life satisfaction and happiness as well as better general health. Therefore, our key results may provide additional motivation for obese males to reduce overweight, similar to tobacco education where preventing impairments in sexual life is an important argument for smoking cessation (Linnebur 2006). From a health policy perspective, improved sexuality might be seen as one of the benefits of weight loss. We cannot say much about the mechanism behind fulfilling sex and better mental and physical health, but scientists in public health and medicine (see e.g. Soysal and Smith 2022) argue that sexual activity can be considered as a form of physical exercise (Smith et al. 2019), so those who engage in regular sexual activity are likely to enjoy the mental and physical benefits of a physically active lifestyle. Moreover, during sexual activity or at the peak of intercourse, endorphins are released, endogenous opioid peptides that act as neurotransmitters and produce a feeling of pleasure or bliss (see Ponzi and Dandy 2019). As argued by Soysal and Smith (2022), circulating endorphins and their association with higher natural killer cell activity have also proved to prevent lung infections and to play an important role in ameliorating lung cancer and many other diseases (Ponzi and Dandy 2019; Zhang et al. 2021).

Moreover, our key result provides further evidence for substantial social costs of obesity that are not necessarily linked to pathological health problems. This raises the question about the channels through which the established effects operate, which cannot be answered on basis of the data used in this study. One possibility is recovery from sexual dysfunction on which the medical literature has focussed. In their review Sarwer et al. (2018), for instance, stress that weight loss, whether from bariatric surgery or less intensive interventions, is associated with significant and clinically meaningful changes in sexual functioning and relevant reproductive hormones. Yet, given the strong effects of a moderate weight reduction found in the data, other channels likely play a role. Greater self-confidence when initiating sexual contact, improved attractiveness to potential sexual partners, feeling more comfortable with their own body, and more desire for sex may act as such channels. One limitation of the present analysis is that it does not allow for disentangling the effect of a reduced body weight from the effect of the process of weight reduction. The strong effects found may therefore rather be attributed to being successful in meeting a major individual challenge than to weight loss itself. One possible route to address this question is to evaluate the effects of bariatric surgeries on sexual activity.

While the existing non-medical literature has a strong focus on the effect of overweight on young individuals’ sex lives, our result explicitly refers to adult obese individuals. This can be interpreted in terms of another desirable consequence of weight reduction. In fact, the argument of overweight having the desirable side effect of counteracting premature initiation to sex, put forward elsewhere in the literature, does not apply given the age group considered in the present analysis. Our results rather provide further evidence for obesity imposing restrictions on private lives, but also for even moderate weight reduction mitigating these limitations.


Corresponding author: Harald Tauchmann, Professur für Gesundheitsökonomie, Friedrich-Alexander-Universität Erlangen-Nürnberg, Findelgasse 7/9, 90402 Nürnberg, Germany; RWI – Leibniz-Institut für Wirtschaftsforschung, Essen, Germany; and CINCH (Health Economics Research Centre), Essen, Germany, E-mail:
The authors are grateful for valuable comments to Arndt Reichert, Joachim Winter, Astrid Schürmann, Nils Heinrich, Christiane Wuckel, Simon Reif, and to the participants of the Economics Research Seminar at the University of Mainz, the Tilburg University TiSEM Seminar: Econometrics and Statistics, and the 2015 dggö Health Econometrics Workshop. Excellent research assistance from Elena Yurkevich and Zekun Zhang is gratefully acknowledged.

A Appendix

A.1 Descriptives for Women Sub-Sample

Table 5:

Descriptive statistics for women sup-sample.

Mean S.D. Median Min. Max.
Dependent variables:

 sexpartner11–22 0.656 0.483 1.000 0.000 1.000
 sexfrequency11–22 1.250 0.916 2.000 0.000 2.000
  t o g e t h e r 22 + 0.594 0.499 1.000 0.000 1.000

Explanatory variables:

 Δ BMI1–4+ −2.507 1.084 −2.022 −5.688 −1.327
 Δ BMI1–10+ −2.265 2.377 −2.462 −7.249 2.384
 Δ BMI510 0.242 1.844 0.000 −3.366 4.147
 age 49.156 10.217 51.000 21.000 63.000
 together0+ 0.625 0.492 1.000 0.000 1.000
 BMI0+ 36.689 4.452 35.339 29.714 49.671
 together4 0.625 0.492 1.000 0.000 1.000
 BMI4 34.182 4.544 33.369 27.690 48.111

Instrumental variables:

 incentive 1501–4+ 0.438 0.504 0.000 0.000 1.000
 incentive 3001–4+ 0.250 0.440 0.000 0.000 1.000
 incentive 2505–10 0.188 0.397 0.000 0.000 1.000
 incentive 5005–10 0.531 0.507 1.000 0.000 1.000
 pharmacy nearby122 0.688 0.471 1.000 0.000 1.000
  1. Notes: Statistics for 32 women who provided information on sexpartner11–22 or sexfrequency11–22 and qualified for the second randomization.

A.2 Results for Alternative Outcome together22

Figure 3: 
Estimated coefficients of weight change for reference specification (grey) and for specification using alternative outcome together22 (black). Spikes indicate 90 percent confidence intervals. Markers:  sexpartner11–22;  sexfrequency11–22;  together22. Source: own calculations.
Figure 3:

Estimated coefficients of weight change for reference specification (grey) and for specification using alternative outcome together22 (black). Spikes indicate 90 percent confidence intervals. Markers: sexpartner11–22; sexfrequency11–22; together22. Source: own calculations.

Table 6:

Living together, preferred model (estimated coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: together22)

Δ BMI5–10 0.164 0.603 0.332 0.301
age 0.017 0.063 0.033 0.054
female −0.433 0.549
together 4 3.199 2.597 2.653 3.679
BMI 4 0.039 0.039 0.026 0.048
constant −3.520 3.787 −3.657 5.409

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.420** 0.490 −1.014** 0.380
incentive 500 5–10 −1.030** 0.496 −0.794** 0.365
age 0.004 0.028 −0.014 0.018
female 0.117 0.360
together 4 −0.112 0.745 0.267 0.418
BMI 4 −0.032 0.033 0.011 0.021
constant 2.152 1.979 1.073 1.234

Selection equation (dependent variable: sexinfo22)

pharmacy nearby 1 22 0.818** 0.206 0.856* 0.167
incentive 250 5–10 −0.284 0.248 −0.016 0.201
incentive 500 5–10 0.174 0.244 0.129 0.198
age 0.022* 0.012 0.015+ 0.010
female 0.098 0.178
together 4 −0.409+ 0.271 −0.305+ 0.207
BMI 4 −0.060** 0.027 −0.050** 0.019
constant 0.386 1.233 0.280 0.838

Ancillary parameters

σ (instrumental equ. Error S.D.) 0.753** 0.043 0.768** 0.031
ρ bmi sex −0.365 2.590 −0.808 1.959
ρ info sex −0.327 0.697 −0.206 0.449
ρ info bmi −0.060 0.154 −0.041 0.107
# of observations (over all) 212 307
# of observations (outcome equation)
log likelihood −523.04 −763.5
Shapiro–Wilk W-statistic (weak identification) 0.967 0.867
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. All model equations estimated using only individuals who successfully lost body weight during month one to four. Simultaneous ML estimation assuming joint normality (binary probit, linear, binary probit). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

A.3 Full Results for Alternative Model Specifications

Table 7:

Sexual relationship, ref. spec. without selection equation (est. coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexpartner11–22)

Δ BMI5–10 −0.403** 0.087 −0.368** 0.154
age −0.003 0.016 −0.015 0.014
female −0.325 0.310
together 4 0.738** 0.384 1.133** 0.474
BMI 4 −0.034 0.029 −0.045 0.033
constant 1.549 1.383 2.369* 1.372

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.409** 0.396 −1.040** 0.325
incentive 500 5–10 −1.015** 0.369 −0.719** 0.326
age 0.004 0.020 −0.014 0.016
female 0.117 0.291
together 4 −0.103 0.378 0.266 0.310
BMI 4 −0.031 0.029 0.012 0.023
constant 2.095 1.501 −1.003 1.177

Ancillary parameters

σ (instrumental equation error S.D.) 2.116** 0.112 2.150** 0.094
ρ (error correlation) 0.851** 0.142 0.687** 0.347
# of observations (over all) 179 260
# of observations (outcome equation) 65 97
log likelihood −410.1 −605.7
Shapiro–Wilk W-statistic (weak identification) 0.896 0.890
Sample split (p-value, LR-test) 0.170
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. All model equations estimated using only individuals who successfully lost body weight during month one to four. Simultaneous ML estimation assuming joint normality (binary probit, linear). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. Test on sample split refers to the null model that pools males and females. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

Table 8:

Frequency of intercourse, ref. spec. without selection equation (est. coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexfrequency11–22)

Δ BMI5–10 −0.275* 0.146 −0.268+ 0.186
age −0.023 0.019 −0.026* 0.014
female −0.115 0.249
together 4 0.519* 0.313 0.945a 0.348
BMI 4 −0.045+ 0.030 −0.049+ 0.031

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.280** 0.440 −0.980** 0.346
incentive 500 5–10 −1.170** 0.358 −0.833** 0.310
age 0.003 0.020 −0.014 0.016
female 0.108 0.291
together 4 −0.106 0.379 0.262 0.311
BMI 4 −0.033 0.029 0.011 0.023
constant 2.236+ 1.502 1.056 1.176

Ancillary parameters

ordered probit threshold 1 −3.245* 1.730 −3.200** 1.436
ordered probit threshold 2 −2.528+ 1.546 −2.558** 1.273
σ (instrumental equation error S.D.) 2.120** 0.112 2.151** 0.094
ρ (error correlation) 0.750** 0.233 0.627* 0.368
# of observations (over all) 179 260
# of observations (outcome equation) 63 95
log likelihood −437.9 −645.6
Shapiro–Wilk W-statistic (weak identification) 0.878 0.875
Sample split (p-value, LR-test) 0.015
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. All model equations estimated using only individuals who successfully lost body weight during month one to four. Simultaneous ML estimation assuming joint normality (ordered probit, linear). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. Test on sample split refers to the null model that pools males and females. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

Table 9:

Sexual relationship, control for success in weight-loss phase (est. coef.).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexpartner11–22)

Δ BMI5–10 −0.425** 0.166 −0.386** 0.150
age −0.007 0.016 −0.020* 0.012
female −0.690** 0.289
together 4 0.579** 0.289 1.233** 0.484
BMI 4 −0.039 0.029 −0.027 0.023
success 4 0.549 0.951 −0.612 0.627
constant 1.732 1.240 2.391** 0.957

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.382** 0.396 −1.018** 0.295
incentive 500 5–10 −0.962** 0.372 −0.731** 0.300
age −0.012 0.020 −0.028** 0.013
female –0.288 0.222
together 4 −0.077 0.311 0.276 0.237
BMI 4 −0.002 0.023 0.009 0.017
success 4 −0.924 1.022 −1.086 0.816
constant 2.281+ 1.444 2.384** 1.039

Instrumental equation (dependent variable: success4)

incentive 150 1–4 0.599** 0.165 0.559** 0.145
incentive 300 1–4 0.562** 0.166 0.650** 0.146
age 0.000 0.009 −0.006 0.007
female −0.179+ 0.122
together 0 0.373** 0.165 0.260** 0.130
BMI 0 0.013 0.013 0.006 0.010
constant −0.840 0.725 −0.245 0.532

Ancillary parameters

σ (instrumental equation error S.D.) 1.986** 0.141 2.096** 0.159
ρ bmi sex 0.641* 0.379 0.707** 0.311
ρ success sex −0.413 0.540 0.498 0.365
ρ success bmi 0.435* 0.258 0.486** 0.205
# of observations (over all) 347 519
# of observations (outcome equation) 103 156
log likelihood −808.3 −1246.2
Shapiro–Wilk W-statistic (weak identification) 0.940 0.902
Sample split (p-value, LR-test) 0.017
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. Simultaneous ML estimation assuming joint normality (binary probit, linear, binary probit). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. Test on sample split refers to the null model that pools males and females. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

Table 10:

Frequency of intercourse, control for success in weight-loss phase (est. coef.).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexfrequency11–22)

Δ BMI5–10 −0.296* 0.175 −0.268+ 0.183
age −0.031** 0.015 −0.031** 0.012
female −0.533** 0.211
together 4 0.416+ 0.262 0.979** 0.308
BMI 4 −0.040+ 0.025 −0.031+ 0.021
success 4 0.084 0.794 −0.829+ 0.531

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.322** 0.376 −0.973** 0.308
incentive 500 5–10 −1.075** 0.328 −0.822** 0.288
age −0.013 0.016 −0.028** 0.013
female −0.294 0.222
together 4 −0.075 0.286 0.277 0.236
BMI 4 −0.004 0.022 0.009 0.017
success 4 −0.772 0.883 −1.097 0.788
constant 2.294* 1.233 2.408** 1.034

Instrumental equation (dependent variable: success4)

incentive 150 1-4 0.661** 0.167 0.553** 0.143
incentive 300 1-4 0.555** 0.166 0.647** 0.144
age 0.000 0.009 −0.006 0.007
female −0.177+ 0.122
together 0 0.324** 0.162 0.264** 0.129
BMI 0 0.013 0.013 0.007 0.010
constant −0.823 0.703 −0.255 0.531

Ancillary parameters

ordered probit threshold 1 −3.717** 1.378 −3.475** 0.976
ordered probit threshold 2 −2.956** 1.271 −2.829** 0.880
σ (instrumental equation error S.D.) 1.977** 0.146 2.101** 0.156
ρ bmi sex 0.558+ 0.356 0.607* 0.357
ρ success sex −0.286 0.479 0.519+ 0.334
ρ success bmi 0.384+ 0.246 0.493** 0.196
# of observations (over all) 347 519
# of observations (outcome equation) 99 152
log likelihood −851.6 −1310.6
Shapiro–Wilk W-statistic (weak identification) 0.925 0.920
sample split (p-value, LR-test) 0.000
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. Simultaneous ML estimation assuming joint normality (ordered probit, linear, binary probit). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. Test on sample split refers to the null model that pools males and females. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

Table 11:

Sexual relationship, unconditional on success (estimated coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Main equation (dependent variable: sexpartner1122)

Δ BMI110 −0.194** 0.086 −0.166* 0.096
age 0.009 0.020 −0.007 0.015
female −0.672** 0.254
together 0 1.313** 0.381 1.508** 0.280
BMI 0 −0.053* 0.030 −0.045* 0.025
constant 1.212 1.533 1.675 1.194

Instrumental equation (dependent variable: ΔBMI1–10)

incentive 150 1–4 −0.701* 0.410 −0.762** 0.328
incentive 300 1–4 −0.960** 0.408 −1.077** 0.322
incentive 250 5–10 −2.753** 0.412 −2.289** 0.322
incentive 500 5–10 −2.049** 0.402 −1.877** 0.326
age 0.004 0.021 0.001 0.016
female 0.076 0.271
together 0 −0.485 0.381 −0.167 0.290
BMI 0 −0.082** 0.028 −0.046** 0.021
constant 3.510** 1.640 2.052* 1.194

Ancillary parameters

σ (control function error S.D.) 2.654** 0.114 2.556** 0.090
ρ (error correlation) 0.536** 0.248 0.268 0.276
# of observations (over all) 271 412
# of observations (main equation) 103 156
log likelihood −684.0 −1034.1
Shapiro–Wilk W-statistic (weak identification) 0.963 0.988
Sample split (p-value, LR-test) 0.069
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. Simultaneous ML estimation assuming joint normality (binary probit, linear). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. Test on sample split refers to the null model that pools males and females. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

Table 12:

Frequency of Intercourse, unconditional on success (estimated coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Main equation (dependent variable: sexfrequency11–22)

Δ BMI1–10 −0.033 0.090 −0.057 0.083
age −0.036** 0.018 −0.025* 0.013
female −0.503** 0.212
together 0 0.677** 0.276 1.065** 0.221
BMI 0 −0.050* 0.026 −0.045** 0.021

Instrumental equation (dependent variable: ΔBMI1–10)

incentive 150 1–4 −0.817** 0.408 −0.781** 0.326
incentive 300 1–4 −0.987** 0.412 −1.083** 0.322
incentive 250 5–10 −2.680** 0.421 −2.262** 0.322
incentive 500 5–10 −2.077** 0.405 −1.903** 0.326
age 0.004 0.021 0.001 0.016
female 0.075 0.271
together 0 −0.491 0.381 −0.169 0.290
BMI 0 −0.083** 0.028 −0.047** 0.021
constant 3.607** 1.639 2.070* 1.194

Ancillary parameters

threshold 1 −4.364** 1.420 −3.212** 1.022
threshold 2 −3.362** 1.388 −2.394** 1.012
σ (control function error S.D.) 2.654** 0.114 2.556** 0.089
ρ (error correlation) 0.266 0.254 0.173 0.229
# of observations (over all) 271 412
# of observations (main equation) 99 152
log likelihood −727.5 −1098.7
Shapiro–Wilk W-statistic (weak identification) 0.991 0.997
Sample split (p-value, LR-test) 0.005
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. Simultaneous ML estimation assuming joint normality (ordered probit, linear). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. Test on sample split refers to the null model that pools males and females. No results are reported for the females’ sub-sample because of weak instruments and convergence issues.

Table 13:

Sexual relationship, 2SLS with Olsen correction (estimated coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexpartner11–22)

Δ BMI5–10 −0.173+ 0.117 −0.155 0.137
age 0.001 0.008 −0.002 0.007
female −0.060 0.116
together 4 0.224 0.164 0.302** 0.135
BMI 4 −0.030* 0.016 −0.031** 0.012
Olsen correction −0.706** 0.356 −0.881** 0.289
constant 1.285* 0.662 1.241** 0.482

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.417** 0.411 −1.034** 0.325
incentive 500 5–10 −0.992** 0.392 −0.738** 0.321
age 0.004 0.018 −0.014 0.018
female 0.116 0.296
together 4 −0.102 0.414 0.265 0.336
BMI 4 −0.030 0.041 0.012 0.033
constant 2.073 1.568 1.012 1.389

Selection equation (dependent variable: sexinfo22)

pharmacy nearby 1 22 0.267** 0.064 0.290** 0.054
incentive 250 5–10 −0.076 0.069 −0.004 0.061
incentive 500 5–10 0.074 0.079 0.051 0.061
age 0.007* 0.004 0.004+ 0.003
female 0.029 0.057
together 4 −0.120* 0.069 −0.098* 0.058
BMI 4 −0.016** 0.004 −0.014** 0.004
constant 0.528** 0.234 0.525** 0.209
# of observations (over all) 212 307
# of observations (outcome equation) 65 97
F-statistic, inst. equation (weak identification) 6.442 5.127
F-statistic, sel. equation (weak identification) 8.094 11.108
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. Bootstrapped standard errors. Model estimated by two-stage least squares. Selection model specified as linear probability model (Olsen 1980). All model equations estimated using only individuals who successfully lost body weight during month one to four. No results are reported for the females’ sub-sample because of weak instruments and convergence issues. F-statistic refer to (joint) significance of instruments in respective first-stage regression.

Table 14:

Frequency of intercourse, 2SLS with Olsen correction (estimated coefficients).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexfrequency11–22)

Δ BMI5–10 −0.238 0.233 −0.256 0.245
age −0.016 0.013 −0.014 0.013
female −0.047 0.208
together 4 0.376 0.293 0.614** 0.236
BMI 4 −0.049* 0.029 −0.053** 0.022
Olsen correction −0.432 0.693 −0.965* 0.571
constant 3.467** 1.000 2.983** 0.853

Instrumental equation (dependent variable: Δ BMI5–10)

incentive 250 5–10 −1.417** 0.388 −1.034** 0.319
incentive 500 5–10 −0.992** 0.395 −0.738** 0.336
age 0.004 0.019 −0.014 0.017
female 0.116 0.296
together 4 −0.102 0.423 0.265 0.337
BMI 4 −0.030 0.041 0.012 0.032
constant 2.073 1.615 1.012 1.347

Selection equation (dependent variable: sexinfo22)

pharmacy nearby 1 22 0.267** 0.064 0.290** 0.053
incentive 250 5–10 −0.076 0.072 −0.004 0.060
incentive 500 5–10 0.074 0.075 0.051 0.060
age 0.007* 0.004 0.004+ 0.003
female 0.029 0.053
together 4 −0.120* 0.069 −0.098* 0.056
BMI 4 −0.016** 0.004 −0.014** 0.004
constant 0.528** 0.235 0.525** 0.206
# of observations (over all) 212 307
# of observations (outcome equation) 63 95
F-statistic, inst. equation (weak identification) 6.442 5.127
F-statistic, sel. equation (weak identification) 8.094 11.108
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. Bootstrapped standard errors. Model estimated by two-stage least squares. Selection model specified as linear probability model (Olsen 1980). All model equations estimated using only individuals who successfully lost body weight during month one to four. No results are reported for the females’ sub-sample because of weak instruments and convergence issues. F-statistic refer to (joint) significance of instruments in respective first-stage regression.

A.4 Results for Spec. that uses Weight Change in Weight-Loss Phase

Figure 4: 
Estimated coefficients of weight change for reference specification (grey) and for specification considering weight change in weight-loss phase (Δ BMI1–4, black). Spikes indicate 90 percent confidence intervals; confidence intervals bootstrapped for the linear model. Markers: ,  outcome sexpartner11–22; ,  outcome sexfrequency11–22. See also footnote 9. Source: own calculations.
Figure 4:

Estimated coefficients of weight change for reference specification (grey) and for specification considering weight change in weight-loss phase (Δ BMI1–4, black). Spikes indicate 90 percent confidence intervals; confidence intervals bootstrapped for the linear model. Markers: , outcome sexpartner11–22; , outcome sexfrequency11–22. See also footnote 9. Source: own calculations.

Table 15:

Sexual relationship, only weight change in weight-loss phase considered (est. coefs.).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexpartner11-22)

Δ BMI1–4 −0.189 0.281 −0.051 0.214
age 0.010 0.021 0.003 0.013
female −0.463* 0.262
together 0 0.952* 0.542 0.960** 0.268
BMI 0 −0.073** 0.029 −0.051** 0.023
constant 1.034 1.369 0.766 1.042

Instrumental equation (dependent variable: Δ BMI1–4)

incentive 150 1–4 −0.916** 0.286 −0.951** 0.233
incentive 300 1–4 −0.916** 0.274 −1.053** 0.210
age 0.012 0.012 0.018* 0.010
female 0.359* 0.202
together 0 −0.547** 0.276 −0.502** 0.206
BMI 0 −0.067** 0.016 −0.055** 0.013
constant 1.238 0.986 0.561 0.759

Selection equation (dependent variable: sexinfo 22 )

pharmacy nearby 1 22 1.282** 0.148 1.248** 0.119
incentive 150 1–4 −0.107 0.178 −0.013 0.140
incentive 300 1–4 −0.046 0.185 −0.008 0.150
age 0.011 0.009 0.012+ 0.007
female 0.010 0.127
together 0 −0.094 0.175 −0.175 0.138
BMI 0 −0.046** 0.017 −0.037** 0.012
constant 0.009 0.913 −0.357 0.640

Ancillary parameters

σ (instrumental equ. Error S.D.) 1.954** 0.053 1.899** 0.042
ρ bmi sex 0.270 0.632 −0.023 0.411
ρ info sex 0.792+ 0.278 0.884** 0.132
ρ info bmi −0.144* 0.077 −0.124* 0.064
# of observations (over all) 469 694
# of observations (outcome equation) 103 156
log likelihood −949.4 −1416.5
Shapiro–Wilk W-statistic (weak identification) 0.963 0.988
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. Simultaneous ML estimation assuming joint normality (binary probit, linear, binary probit). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. See also footnote 9.

Table 16:

Freq. of intercourse, only weight change in weight-loss phase considered (est. coefs.).

Males Males & females
Est. coef. S.E. Est. coef. S.E.
Outcome equation (dependent variable: sexfrequency11–22)

Δ BMI1-4 −0.063 0.276 0.113 0.197
age −0.033+ 0.023 −0.022* 0.013
female −0.484** 0.220
together 0 0.578+ 0.355 0.958** 0.229
BMI 0 −0.067** 0.031 −0.047* 0.025

Instrumental equation (dependent variable: Δ BMI1–4)

incentive 150 1–4 −0.940** 0.288 −0.958** 0.232
incentive 300 1–4 −0.891** 0.273 −1.048** 0.210
age 0.013 0.012 0.018* 0.010
female 0.358* 0.202
together 0 −0.547** 0.272 −0.502** 0.204
BMI 0 −0.068** 0.016 −0.056** 0.013
constant 1.243 0.985 0.559 0.750

Selection equation (dependent variable: sexinfo22)

pharmacy nearby 1 22 1.244** 0.151 1.234** 0.122
incentive 150 1–4 −0.148 0.177 −0.050 0.145
incentive 300 1–4 −0.016 0.182 0.027 0.149
age 0.010 0.009 0.011+ 0.007
female 0.044 0.130
together 0 −0.056 0.176 −0.141 0.140
BMI 0 −0.051** 0.017 −0.039** 0.012
constant 0.212 0.893 −0.281 0.629

Ancillary parameters

ordered probit threshold 1 −4.344** 1.679 −2.942** 1.053
ordered probit threshold 2 −3.379** 1.633 −2.158** 1.043
σ (instrumental equ. Error S.D.) 1.954** 0.053 1.899** 0.042
ρ bmi sex 0.210 0.567 −0.193 0.380
ρ info sex 0.381 0.297 0.420* 0.212
ρ info bmi −0.141* 0.078 −0.124* 0.064
# of observations (over all) 469 694
# of observations (outcome equation) 99 152
log likelihood −991.9 −1482.4
Shapiro–Wilk W-statistic (weak identification) 0.975 0.980
  1. Notes: **Significant at 5 %; *significant at 10 %; +significant at 15 %. Simultaneous ML estimation assuming joint normality (ordered probit, linear, binary probit). ρ denotes cross-equation error correlations. Thresholds for Shapiro–Wilk W-statistic (Denzer and Weiser 2021): 0.994 (5 percent relative bias), 0.952 (10 percent relative bias); 499 bootstrap replications (suggested by Denzer and Weiser 2021); ‘W < threshold’ indicates weak identification. See also footnote 9.

References

Adolfsson, B., Elofsson, S., Rössner, S., and Undén, A.-L. (2004). Are sexual dissatisfaction and sexual abuse associated with obesity? A population-based study Obes. Res. 12: 1702–1709, https://doi.org/10.1038/oby.2004.211.Search in Google Scholar

Ali, M.M., Rizzo, J.A., Amialchuk, A., and Heiland, F. (2014). Racial differences in the influence of female adolescents’ body size on dating and sex. Econ. Hum. Biol. 12: 140–152, https://doi.org/10.1016/j.ehb.2013.11.001.Search in Google Scholar

Augurzky, B., Bauer, T.K., Reichert, A.R., Schmidt, C.M., and Tauchmann, H. (2012). Does money burn fat? Evidence from a randomized experiment. Ruhr Econ. Pap.: 368, https://doi.org/10.4419/86788423.Search in Google Scholar

Augurzky, B., Bauer, T.K., Reichert, A.R., Schmidt, C.M., and Tauchmann, H. (2018). Habit Formation, obesity, and cash rewards. Ruhr Econ. Pap. 750, https://doi.org/10.4419/86788871.Search in Google Scholar

Augurzky, B., Bauer, T., Reichert, A., and Tauchmann, H., RWI (2017). RWI-Obesity. Version: 1. RWI – Leibniz institute for economic research. Data, https://doi.org/10.7807/Obesity:2015:v1.Search in Google Scholar

Baumeister, R.F., Catanese, K.R., and Vohs, K.D. (2001). Is there a gender difference in strength of sex drive? Theoretical views, conceptual distinctions, and a review of relevant evidence. Pers. Soc. Psychol. Rev. 5: 242–273, https://doi.org/10.1207/s15327957pspr0503_5.Search in Google Scholar

Cawley, J. (2001). Body weight and the dating and sexual behaviors of young adolescents. In: Michael, R.T. (Ed.), Social awakening: adolescent behavior as adulthood approaches. Russell Sage, pp. 174–198.Search in Google Scholar

Cawley, J., Joyner, K., and Sobal, J. (2006). Size matters: the influence of adolescents’ weight and height on dating and sex. Ration. Soc. 18: 67–94, https://doi.org/10.1177/1043463106060153.Search in Google Scholar

Denzer, M. and Weiser, C. (2021). Beyond F-statistic – a general approach for assessing weak identification. Gutenberg School of Management and Economics Discussion Paper Series, Vol. 2107.Search in Google Scholar

Eilers, L. and Pilny, A. (2018). Data from a randomized experiment: financial incentives on weight loss (RWI-Obesity). Jahrb. Natl. Stat. 238: 591–600, https://doi.org/10.1515/jbnst-2017-0138.Search in Google Scholar

Esposito, K., Giugliano, F., Palo, C.D., Giugliano, G., Marfella, R., D’Andrea, F., D’Armiento, M., and Giugliano, D. (2004). Effect of lifestyle changes on erectile dysfunction in obese men: a randomized controlled trial. JAMA, J. Am. Med. Assoc. 291: 2978–2984, https://doi.org/10.1001/jama.291.24.2978.Search in Google Scholar

Fenton, K., Johnson, A., McManus, S., and Erens, B. (2001). Measuring sexual behaviour: methodological challenges in survey research. Sex. Transm. Infect. 77: 84–92, https://doi.org/10.1136/sti.77.2.84.Search in Google Scholar

Hafner, L., Tauchmann, H., and Wübker, A. (2021). Does moderate weight loss affect subjective health perception in obese individuals? Evidence from field experimental data. Empir. Econ. 61: 2293–2333, https://doi.org/10.1007/s00181-020-01971-8.Search in Google Scholar

Halpern, C.T., King, R.B., Oslak, S.G., and Udry, J.R. (2005). Body mass index, dieting, romance, and sexual activity in adolescent girls: relationships over time. J. Res. Adolesc. 15: 535–559, https://doi.org/10.1111/j.1532-7795.2005.00110.x.Search in Google Scholar

Harris, K.M. (2013). The add health study: design and accomplishments. Carolina Population Center, University of North Carolina at Chapel Hill. Available at: <http://www.cpc.unc.edu/projects/addhealth/data/guides>.Search in Google Scholar

Kloner, R.A., Carson, C., Dobs, A., Kopecky, S., and Mohler, E.R. (2016). Testosterone and cardiovascular disease. J. Am. Coll. Cardiol. 67: 545–557, https://doi.org/10.1016/j.jacc.2015.12.005.Search in Google Scholar

Kolotkin, R.L., Binks, M., Crosby, R.D., Østbye, T., Gress, R.E., and Adams, T.D. (2006). Obesity and sexual quality of life. Obesity 14: 472–479, https://doi.org/10.1038/oby.2006.62.Search in Google Scholar

Krauth, C. and Bartling, T. (2017). Is rehabilitation worth it?: review of economic evaluations of rehabilitation in Germany [Lohnt sich Rehabilitation? : Überblick über gesundheitsökonomische Analysen der Rehabilitation in Deutschland]. Bundesgesundheitsblatt – Gesundheitsforsch. – Gesundheitsschutz 60: 394–401, https://doi.org/10.1007/s00103-017-2514-0.Search in Google Scholar

Larsen, S.H., Wagner, G., and Heitmann, B.L. (2007). Sexual function and obesity. Int. J. Obes. 31: 1189–1198, https://doi.org/10.1038/sj.ijo.0803604.Search in Google Scholar

Linnebur, S.A. (2006). Tobacco education: emphasizing impotence as a consequence of smoking. Am. J. Health Syst. Pharm. 63: 2509–2512, https://doi.org/10.2146/ajhp050288.Search in Google Scholar

Neymotin, F. and Downing-Matibag, T. (2014). Adolescent obesity: a barrier to mate selection? J. Child Adolesc. Behav. 2: 1–3, https://doi.org/10.4172/jcalb.1000157.Search in Google Scholar

Olsen, R.J. (1980). A least squares correction for selectivity bias. Econometrica 48: 1815–1820, https://doi.org/10.2307/1911938.Search in Google Scholar

Petersen, J.L. and Hyde, J.S. (2010). A meta-analytic review of research on gender differences in sexuality, 1993–2007. Psychol. Bull. 136: 21–38, https://doi.org/10.1037/a0017504.Search in Google Scholar

Ponzi, D. and Dandy, M. (2019). The influence of endogenous opioids on the relationship between testosterone and romantic bonding. Hum. Nat. 30: 98–116, https://doi.org/10.1007/s12110-018-9332-4.Search in Google Scholar

Reichert, A.R. (2015). Obesity, weight loss, and employment prospects – evidence from a randomized trial. J. Hum. Resour. 50: 759–810, https://doi.org/10.3368/jhr.50.3.759.Search in Google Scholar

Roodman, D. (2011). Fitting fully observed recursive mixed-process models with cmp. Stata J. 11: 159–206, https://doi.org/10.1177/1536867x1101100202.Search in Google Scholar

Rowland, D.L., McNabney, S.M., and Mann, A.R. (2017). Sexual function, obesity, and weight loss in men and women. Sex. Med. Rev. 5: 323–338, https://doi.org/10.1016/j.sxmr.2017.03.006.Search in Google Scholar

Sabia, J.J. and Rees, D.I. (2008). The effect of adolescent virginity status on psychological well-being. J. Health Econ. 27: 1368–1381, https://doi.org/10.1016/j.jhealeco.2008.05.008.Search in Google Scholar

Sabia, J.J. and Rees, D.I. (2011). The effect of body weight on adolescent sexual activity. Health Econ. 20: 1330–1348, https://doi.org/10.1002/hec.1674.Search in Google Scholar

Sarwer, D.B., Hanson, A.J., Voeller, J., and Steffen, K. (2018). Obesity and sexual functioning. Curr. Opin. Rep. 7: 301–307, https://doi.org/10.1007/s13679-018-0319-6.Search in Google Scholar

Schafer, M.H., Mustillo, S.A., and Ferraro, K.F. (2013). Age and the tenses of life satisfaction. J. Gerontol. B Psychol. Sci. Soc. Sci. 68: 571–579, https://doi.org/10.1093/geronb/gbt038.Search in Google Scholar

Smith, L., Grabovac, I., Yang, L., Veronese, N., Koyanagi, A., and Jackson, S.E. (2019). Participation in physical activity is associated with sexual activity in older English adults. Int. J. Environ. Res. Publ. Health 16: 489, https://doi.org/10.3390/ijerph16030489.Search in Google Scholar

Soysal, P. and Smith, L. (2022). Sexual activity and successful aging. Eur. Geriatr. Med. 13: 1235–1238, https://doi.org/10.1007/s41999-022-00694-6.Search in Google Scholar

Taris, T.W. (2000). Quality of mother-child interaction and the intergenerational transmission of sexual values: a panel study. J. Genet. Psychol. 161: 169–181, https://doi.org/10.1080/00221320009596704.Search in Google Scholar

Woloski-Wruble, A.C., Oliel, Y., Leefsma, M., and Hochner-Celnikier, D. (2010). Sexual activities, sexual and life satisfaction, and successful aging in women. J. Sex. Med. 7: 2401–2410, https://doi.org/10.1111/j.1743-6109.2010.01747.x.Search in Google Scholar

Wooldridge, J.M. (2002). Econometric analysis of cross section and panel data. MIT Press, Cambridge Massachusetts.Search in Google Scholar

World Health Organization (2000). Obesity: preventing and managing the global epidemic, Technical Report 894, Geneva.Search in Google Scholar

Zhang, S., Liu, N., Ma, M., Huang, H., Handley, M., Bai, X., and Shan, F. (2021). Methionine enkephalin (MENK) suppresses lung cancer by regulating the Bcl-2/Bax/caspase-3 signaling pathway and enhancing natural killer cell-driven tumor immunity. Int. Immunopharm. 98: 107837, https://doi.org/10.1016/j.intimp.2021.107837.Search in Google Scholar

Ziebarth, N.R. (2014). Assessing the effectiveness of health care cost containment measures: evidence from the market for rehabilitation care. Int. J. Health Care Finance Econ. 14: 41–67, https://doi.org/10.1007/s10754-013-9138-1.Search in Google Scholar

Received: 2022-06-13
Accepted: 2023-06-09
Published Online: 2023-08-15
Published in Print: 2023-12-15

© 2023 the author(s), published by De Gruyter, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 13.10.2025 from https://www.degruyterbrill.com/document/doi/10.1515/jbnst-2022-0039/html
Scroll to top button