Larger social networks may increase stigma against vocal illness: An integrated empirical and computational study of deciphering help-seeking behaviors and vocal stigma

Theoretical Framework of Health Stigma

Stigma is considered a serious threat to population health inequalities as it disrupts resource availability, social relationships and coping behaviors within society [19]. Individuals with a stigmatized medical condition may become reluctant to seek medical help [20]. For example, a meta-analysis on mental health stigma found a negative correlation between stigma and help-seeking [21]. Similarly, a meta-analysis of stigma and health outcomes in HIV/AIDS found that those who experience symptoms of stigma were 21% less likely to access health and social services [22]. Individual studies also found associations between stigma and avoidance or delaying of help-seeking in other health issues, including alcoholism [23] and cancers [24]. Similarly, stigma against vocal illness, i.e., vocal stigma, is a known barrier to vocal performers seeking help for work-related illness.

In 2019, Stangl et al. published a seminal paper “The Health Stigma and Discrimination Framework” in BMC Medicine [25]. The paper proposed a hierarchy of social levels on which stigma operates, and a series of processes leading from the creation of a stigma to its measurable outcomes. In brief, the proposed framework consists of a multi-level social hierarchy, namely, individual, interpersonal, organizational, community, and public policy. Within this hierarchy, the stigmatization process starts with constituent factors that drive or facilitate health stigma. Then, individuals with a specific health condition are “marked” or “labeled” with a stigma, followed by manifestations, i.e., attitudes, behaviors, and experiences that arise from the stigma in those who are marked by the stigma and those around them. Finally, the stigma manifestations affect a range of measurable outcomes on affected populations (e.g., access to healthcare) and the society (e.g., public policies against discrimination).

By adopting Stangl et al.’s framework to vocal stigma, vocal performers with vocal impairments can entail catastrophic fears about losing the employment, social status, and creative fulfilment that their art provides [26]. In other words, they fear the consequences of not being able to meet the vocal demands of performance. This fear is a strong candidate for a driving factor in vocal stigma among vocal performers. However, individuals face different vocal demands and consequences for being unable to perform. For example, actors and singers face different vocal demands, and full-time performers face greater consequences from voice disorders than those who have other sources of income. (Figure 1)

• Download figure
• Open in new tab

Figure 1. Proposed Health Stigma and Discrimination Framework for Vocal Stigma.

Adapted from Stangl et al. (2019)

Predicting Help-Seeking Behavior

Predicting help-seeking and other behaviors that influence health is a major topic in the field of health psychology [27]. One prominent psychological model for making these predictions is the Information Motivation Behavioral Skills (IMB) model [28]. Information represents a person’s knowledge and beliefs about the behavior, and about the issue that the behavior is intended to address. Motivation represents a person’s attitudes toward the predicted behavior, and their perception of social norms around it, i.e., what attitudes does the person believe other people hold? An example can be perceived costs and benefits of the behavior. Behavioral Skills represents a person’s ability to perform the behavior. This component can include objective skills, psychological skills and more importantly self-efficacy, i.e., a person’s intention to seek help [29].

By applying this IMB model to vocal stigma, information includes a person’s knowledge and belief about voice disorders and related issues. Motivation includes a person’s attitudes and perception of social norms about seeking professional medical help for a voice disorder. Finally, behavioral skills represent the intention to seek professional medical help if they acquire a voice disorder (Figure 2).

• Download figure
• Open in new tab

Figure 2. Proposed Information Motivation Behavioral Skills Model of Predicting Help-Seeking Behavior for Voice Disorders.

Adopted from Fisher et al. (2002).

Computer Simulation of Social Networks

In response to the well-recognized health challenges presented by stigma, various strategies have been proposed to reduce the stigma around conditions such as HIV and mental illness. At the level of social-stigma, one common approach is to provide information about a health condition to the general public, emphasizing that affected individuals are not to blame for their condition. Other interventions includ promoting empathy and understanding toward affected individuals, for example via testimonials or social contact between affected and unaffected individuals. For self-stigma, a provision of counselling or coping strategies can also be implemented to affected individuals [30–32].

Of these, social contact has emerged as the most successful strategy for reducing stigma in adults [30, 33]. Computer simulation has emerged as a powerful, quantitative way to study social networks and their effects on human behavior and public policy. In fact, social networks have been investigated as a structural determinant of health on their role in the spread of social behaviors and attitudes affecting health outcomes [34]. For instance, testable models have been developed to simulate the mechanisms of misinformation, bias, and polarization in social networks [35, 36], characterize the effects of misinformation during the COVID-19 pandemic [37], and help understand stigma and its impact on obesity, depression, mental illness, heart disease, and HIV-prevention [38–41].

In particular, agent-based models (ABM) have become a pivotal simulation platform in social and behavioral sciences in the fast-growing field of social computing [42–44]. ABM have been applied for building and testing theories in stock markets [45–47], consumer behaviors [48–53], and psychological cooperation [54–57]. For driver behaviors in particular [58–66], ABMs are used to simulate how the behavior of a driver (an individual agent) will affect other drivers (other agents) and the evolving traffic patterns (emergent population behavior). Rather than individual isolated decisions, social behaviors often resulted from interactions among people with diverse backgrounds over time. Such characteristics match very well with the framework of ABMs.

In essence, an ABM is built up by a collection of entities called agents, whose state and behavior are governed by a set of computational rules. Agents “live” in the virtual world, which is discretized into grids, or patches. Each patch is a certain situation or a physical environment that an agent lives in. In the application of social behaviors, each agent can represent an individual person with distinctive features with respect to demography, psychology and social history. Second, an individual agent makes its own decision and adapts the behavior to the environment over time that effectively represents the dynamics of social and physical influences on human behavior in the real world [43]. Lastly, ABM is technically flexible for users to explore a wide parameter space and a long temporal scale, which otherwise would be very costly or at times prohibitive with empirical experiments. For example, ABM can simulate human interactions in social networks with the ability to control for network sizes, type and frequency of interaction, and observe individual and group outcomes from weeks to years. Given that the process of stigmatization is multiscale and time-evolving, ABM presents a new, exciting opportunity in the research of health stigma.

Study Objectives

The objective of this study was to investigate the experience of vocal stigma using empirical survey questionnaires and computer network models. The primary focus was to predict how peer interactions could affect an individual’s help-seeking behaviors and vocal stigma experience especially in vocal performers. With a better understanding of vocal stigma, clinicians, social epidemiologists and public health policy makers will be more informed of devising interventions aimed at reducing vocal stigma and improving access to voice healthcare.

2.1. Empirical Study of Vocal Stigma Profiling

2.2. Computational Study of Simulating Social Interactions and Vocal Stigma

2.2.1. Vocal Stigma Agent-Based Models: Conceptual Framework

In this vocal stigma agent-based model (VS-ABM), each agent represents an avatar of a vocal performer (VP) or a non-vocal performer (non-VP), who functions as a social actor and interacts with their peers to share information and their motivations for help-seeking related to voice health. Depending on an agent’s geolocation, they have access to a small or a large peer network in the VS-ABM. When two agents interact, they will exchange information and motivations related to vocal health.

An agent’s vocal health (history and frequency of disordered voice), stigma experience (self- and social-stigma) and IMB are connected via two feedback loops (social and personal), which in turn affects an agent’s tendency to seek help (Figure 3). The social loop simulates the interrelationship between social-stigma and IMB. Also, an agent’s motivation is linked to its own frequency of voice disorder. For instance, if an individual has frequent relapses in voice disorders, they will likely have less motivation to seek help from social partners and professionals. In contrast, the personal loop simulates the interrelationship between self-stigma and IMB. Also, an agent’s information related to vocal health is contingent on its history of voice disorder. For instance, if an individual has prior experience with a voice disorder, they will likely have more information about the nature of it and may have reduced their self-stigma as a result.

• Download figure
• Open in new tab

Figure 3. The Framework of Vocal Stigma Agent-Based Models (VS-ABM).

An agent’s vocal health, vocal stigma, information, motivation, and behavioral skills (IMB) are presumably interrelated and modulated via the social and personal feedback loops, which collectively affects an agent’s tendency to seek help.

3.1 Overall Relation between Total Stigma and Social Network Size

Our original hypothesis for Research Question #1 was that vocal stigma would decrease with larger social networks. Unexpectedly, this hypothesis was not supported by our ABM simulation results in looking at changes in total stigma relative to the number of interactions with network sizes (50-400 agents) and initial social-stigma profiles (low vs average vs high) for each VP, non-VP and general populations (Figure 4A).

• Download figure
• Open in new tab

Figure 4. Total Stigma Experience as Functions of Network Size and Number of Interactions in Agents with Low, Average and High Initial Social-Stigma.

(A) Aggregated trends of mean total stigma scores. Each trend aggregated from 100 simulations initialized with social-stigma and network size. (B) Varying outcomes of mean total stigma score from individual simulations for VP agents with high risk of voice disorders, high initial social-stigma, and no intervention within the 400-agent network size. The aggregated trend, same as VP group in Fig 4A, for all 100 simulations is shown in bold. (C) A representative snapshot of simulated social networks with 400 VP agents near the maximal total stigma after 100,000 interactions in the 400-VP network. An agent’s social-stigma is color-coded from low (Green) to high (Red).

First on the VP networks, agents in the two smallest network sizes (50 and 100) had their total stigma stabilized as early as 25,000 social interactions. Their final total stigma values were also lower than initial values, regardless of initial social-stigma. In contrast, when the networks became sizable (> 150 agents), the trajectories of total stigma started to follow a log-normal distribution pattern. That is, agents’ total stigma quickly increased to a peak and then slowly leveled off but did not always return to initial values. Interestingly, VP agents with high initial social-stigma took 25% fewer social interactions (50,000 to 100,000) to reach the peak of total stigma, compared to those with low initial social-stigma (the general and non-VP groups). This initial observation suggests that total stigma likely worsens over time with social interactions in larger network groups, particularly for individuals prone to higher social-stigma such as VP at high risk for voice disorder and a tendency to not seek intervention.

To further examine this paradox, we proceeded to examine individual-specific variations in the evolution of total stigma in the largest network (400 agents) from the VP group with initial high stigma up to 500,000 interactions (Figure 4B). Three distinct patterns emerged from this case simulation. (1) First, one set of simulations quickly leveled off towards a total stigma value lower than the initial value, which is denoted as the “decreasing cohort”. (2) Another set of simulations, which is denoted as the “increasing cohort”, increased rapidly towards a maximum between 50,000 and 200,000 interactions and then slowly leveled off to a total stigma value larger than the initial value. (3) Lastly, the “chaotic cohort” is a set of agent simulations that appeared to fluctuate randomly and did not appear to stabilize even after 500,000 social interactions. An extended set of simulations were further run for the “chaotic cohort” and confirmed that the total stigma value continued fluctuating randomly up to 10,000,000 interactions (data not shown).

Except for the “chaotic cohort”, most agents’ total stigma appeared to approach maximum around 100,000 social interactions with large social networks. To understand one such maxima, we further visualized an established social network from the “increasing cohort” in 400 VP agents after 100,000 social interactions (Figure 4C). Multiple clusters or social groups with similar social-stigma values were formed within the global network, suggesting a homogenization with these highly interconnected social groups. Also, these clusters were not isolated but rather fairly connected to neighboring groups. In other words, peers were connecting with the whole set of the neighboring group, rather than simply bottlenecking through a few nearby peers.

In particular, a few social groups appeared to have a heterogeneous mix of agents with initial low and initial high social-stigma profiles (Figure 4C: networks with mixed green and red color), confirming that agents could interact with opposite social-stigma values. Especially in this “increasing cohort”, the emergence of mixed social groups in the global network coincides with the stabilization of total stigma around the critical point of 100,000 social interactions. We speculated that, at the beginning of interactions, agents would establish social groups with peers around them, becoming more like their neighbors. Then, when social groups became homogenous, the total stigma would approach to its peak. Eventually, when social groups with polarized views interacted, a mix of agents with low and high social-stigma emerged and, paradoxically, stabilized the population’s social-stigma.

To answer Research Question #1, our simulation results indicated that, on average, vocal stigma increases with larger social groups. However, there are non-trivial trends where total stigma diminishes as social groups become less polarized. In addition, our results suggested that interactions between groups are not bottlenecked with a small set of agents but distributed among the entire social group. The clusters can exchange information/motivation with multiple members of neighboring groups, which may ultimately result in less polarized neighboring groups.

3.2 Social Network Size in relation to Each IMB and Vocal Stigma Components

For Research Question #2, the original hypothesis was that social-stigma would decrease when agents were more willing to socially interact, i.e., lower social-stigma and higher IMB. This hypothesis was partially supported by our ABM simulation results in comparing the probability distributions of final IMB and final total stigma scores between the smallest and largest social networks (50 vs 400 agents) (Figure 5A-B).

• Download figure
• Open in new tab

Figure 5. Probability Distribution of Total IMB and Total Stigma Experience.

A total of 500,000 social interactions were simulated for VP, Non-VP and general populations with low, average or high initial social-stigma. (A) Total IMB and its components (information, motivation, and behavioral skills) at 50 and 400 network sizes (B) Total stigma and its components (self-stigma and social-stigma) at 50 and 400 network sizes (C) Average change in total IMB (diamonds) and total stigma (circles) from their initial values. Each row of subplots represents a specific magnitude of change, i.e., from less than 1 SD (σ) to more than 2 SD increase (inc.) or decrease (dec.), in total IMB and total stigma. % = percentage of simulations with corresponding changes in total stigma across all network sizes. (D) Ratio of simulation with total stigma changes more than 1 SD (Panel C row 1 + 2) relative to those with little or no change (Panel C row 3).

Following a total of 500,000 social interactions, a clear bimodal distribution of total IMB score was noted in large but not in small networks across all three agent populations (Figure 5A). Of the total IMB score, the first mode (local maximum) appeared around 60 and 75, whereas the second mode appeared around 110 in the probability density function. Each IMB element (information, motivation and behavioral skills) also appeared to be highly coupled. That is, an agent population, regardless of its initial social-stigma, can be sub-grouped into either with low IMB and high IMB. Such sub-group distinction is much more pronounced in large social networks than small ones.

For the total stigma, a bimodal distribution also appeared in large but not small social networks across all three agent populations (Figure 5B). Of the total stigma score, the first mode was centered around 10, whereas the second mode was around 20 within the data. Each total stigma element (social- and self-stigma) also displayed two distinctive peaks, but their distance was more separated in social-stigma score. Additionally, comparing high-risk VP groups, those with social intervention had overall lower social-stigma, whereas the maximal mode for the no intervention group showed higher social-stigma. Our results indicated that social-stigma experience has more weight than self-stigma in the total stigma experience.

Further analysis was performed to evaluate the magnitude of changes in total IMB and total stigma by aggregating VP, non-VP, and general populations before and after the 500,000 interactions (Figure 5C-D). For total IMB, about 35% and 63% of agents showed increases or no changes of their scores in smaller networks (50 and 100 agents) at the end of simulations. With larger network sizes (>150 agents), about 36% of agents showed significant increases in their total IMB scores with more than two standard deviations (SD, σ), i.e., between 30 and 40 points increase from initial values. At the same time, about 39% and 21% of agents showed reduction in their IMB scores more than 1 and 2 SD respectively in these larger networks (Figure 5C). For total stigma, all agents in smaller networks had negligible changes of their scores, whilst 25% of them had total stigma increased by more than 6 points when they were in larger networks. For those with high risk and high social-stigma, VP agents were 1.5 times more likely to experience more total stigma if they had no intervention, compared to those with social intervention, when they are in the largest social network (Figure 5D).

All in all, the aforesaid simulation results confirmed that higher total IMB was associated with lower total stigma for small networks. Interestingly, IMB outcomes vary across network sizes. Despite minimal change to total stigma, IMB tended to decrease more for larger social networks, while IMB for smaller networks tended to increase. The bimodal distribution noted in large social networks (i.e., 400) may be attributed to the increased likelihood that extreme values of social-stigma arise in the population. Also, within a large social network, their social groups tend to be more polarized with group size as noted in Figure 4C, leading to a wider separation between the two modes of the bimodal distributions.

To address Research Question #2, our results supported the hypothesis in cases of small networks, i.e., lower social-stigma leads to an overall decrease in total stigma. Larger networks are far more likely to have significantly increased total stigma. For high-risk VPs, social interventions might help mitigate high social-stigma and eventually a reduction in total stigma over time.

3.3 Network Structure in relation to Total IMB and Total Stigma

To understand the connections among agents in the social network, degree of node and clustering coefficient were computed in relation to the changes in total IMB and total stigma when the social network was established (i.e., after 100,000 social interactions) in all 4000 simulations.

First, the degree of node represents the number of social links that an agent has in a social network. (Figure 6A) The average number of links ranged from 2.5 in the smallest social network (50 agents) to 13.5 in the largest social network (400 agents). This positive linear relationship is expected because each agent has more peers to interact within their vicinity in larger social networks. For instance, the relative patch coverage increases from 5% in the 50-agent network to nearly 36% in the 400-agent network.

• Download figure
• Open in new tab

Figure 6. Network Structure Analysis of Network Size Effects on IMB and Vocal Stigma.

(A) Distribution of average number of peer links (degree of nodes) by network size. (B) Distribution of average local clustering coefficient by network size. (C) Total IMB values decrease as a function of average number of peer links, which is highly correlated with network size. A clear bifurcation appears around 4.5 peer links. (D) Total stigma values increase as a function of average number of peer links. No bifurcation appears. However, distributions show a larger range of changing stigma values for more peer links.

Second, the clustering coefficient represents how close-knit an agent’s social group (i.e., cluster) i interconnected in a social network. The clustering coefficient is measured on a scale of 0 to 1. A lower cluster coefficient is when peers are only connected through the agent, whereas a higher cluster coefficient is when agents’ peers connected independently and create denser and tighter connections.

Interestingly, the average clustering coefficient for the network does not show a linear trend with network size (Figure 6B). A wide spread of clustering coefficients was observed in small social networks wherea the probability distribution approximated to a normal distribution towards a larger social network. The cluster coefficient is also higher in larger networks (400 agents; mean (SD) = 0.61 (0.009)) than that of smaller networks (50 agents; mean (SD) = 0.31 (0.086)). In small networks, agents were sparsely randomized in the simulation space, which geographically favored agents being connected with their neighboring peers only but not reaching out to other social groups that are further away.

Finally, we analyzed how the number of peer links (i.e., social partners) affects an agent’s total IMB and total stigma scores when the social network is stabilized, i.e., all links are maximally connected. For total IMB, more peer links might relate to a gradual decrease in IMB. (Figure 6C) Upon the stabilization of social networks, the IMB increases by 7.6 points in the smallest 50-agent network (1.5 social links); whereas it decreases by 10.6 points in the largest 400-agent network (12.5 social links). For total stigma, more peer links might contribute to a slight increase in total stigma (Figure 6D). Upon network stabilization, the total stigma decreases by 1.2 points in the 50-agent network yet increases by 3.5 points in the 400-agent network.

4.1. Diverging population outcomes and heterogeneity within social groups

Three population behaviors, namely decreasing, increasing, and chaotic cohorts, were identified across all network sizes. The “decreasing cohort” had characteristics that resembled the small network trends, such that total stigma decreases and quickly stabilizes at a value less than the original population stigma. The “increasing cohort” had characteristics that resembled the large network trends, where total stigma increases to a maximum and then levels off towards a value larger than the initial stigma value. Such stigma maxima are seen when the social network of the population is fully established, and social groups are largely homogenous. When polarized homogenous groups start to interact and exert more influence on each other, greater heterogeneity of the groups emerges at a population level. Lastly, the “chaotic set” is referred to an oscillatory behavior when total stigma does not stabilize. Such complex patterns of vocal stigma are not programmed (i.e., hard-coded) in VF-ABM but rather emerge de novo (i.e., emergency) from dynamical systems. Our simulation outcomes may also be in line with a new notion of “stigma mutation”. [69] For instance, stigma associated with COVID-19 has a social-historical context and changes over time in response to the local and global environment. In our case, agents’ social-personal backgrounds (e.g., occupation, medical history, initial social-stigma) modulate the evolution of their total stigma in a network environment.

4.2 Social-stigma drives change in vocal stigma

We further investigated which specific IMB and stigma components drove changes in population outcomes in relation to network size. Overall, information and behavioral skills tend to be coupled more strongly with motivation being inversely related, i.e., higher information and behavioral skills relative to a lower motivation. Social-stigma has more impact than self-stigma on total stigma. The variance of IMB changes is considerably large across all network sizes, ranging from −6 to +5 SD changes on average. In contrast, the variance of total stigma changes remains quite consistent – around 75%-85% of all simulations show less than 1 SD change on average. Moreover, in larger network sizes, social interactions led to three distinctive patterns of social-stigma, namely, (1) small decreases in low/ average risk groups, (2) small increases in high-risk groups; or (3) large increases regardless of risk group. In particular, larger networks with higher risk groups of voice disorders have considerably higher odds of increasing the total stigma than those of lower risk groups. Among these higher risk individuals, those without medical interventions have higher odds (4:4) to increase total stigma after social interaction than those who received treatment (3:4). That is, for large networks (i.e., 400 agents), about 51% of VP with no intervention are likely to increase their stigma after interacting with their peers; this likelihood drops to 41% if they receive intervention.

Another key observation is that IMB and stigma do not always have a directional relationship. For instance, if we looked at simulated populations whose total stigma changed more than 1 SD with social interactions, their IMB scores could increase by 20 points at most. For those with changes more than 2 SD in total stigma, IMB decreased significantly by 30 points – in this case, vocal stigma might get worse despite informational and social intervention. Similar notions were also confirmed in reducing public stigma of mental illness, in which stigma is difficult to eliminate especially if it is deep rooted in social relationships [70].

4.3 Larger social groups may promote vocal stigma

Two standard network measures were analyzed to evaluate how connected VPs were (i.e., degree of node) and how close their social groups were linked (i.e., clustering coefficient). As the network size increased, a VP agent had more peers to interact and connect with in their vicinity. That is, for larger networks, peers had consistently larger social groups. However, as we simulated hundreds of thousands of social interactions, all agents were able to eventually interact and connect with all their potential peers and all social groups were maximally interconnected or clustered. As social group networks stabilize and become maximally connected, a gradual increase in total stigma and decrease in IMB was evident in larger but not smaller social groups. These analytics also provided part of the evidence to answer Research Question #1, i.e., larger networks stabilized towards more polarized social groups. This observation aligns with similar social network findings of polarization of political groups and how negativity spreads through larger social networks faster than positivity [71]-[72].

4.4. Limitations and Future Directions

This study investigates vocal stigma, for the very first time, in the context of social computing with expected technical challenges. To start, our vocal-stigma specific questionnaire was deployed to Canadians only for this initial round of analysis. Additional data from other countries can be collected to evaluate the sensitivity and specificity of VS-ABM or if additional organization or institutional factors (e.g., reimbursement system) will be needed to be included in future simulations. Second, our initial model results suggest the presence of a “chaotic set” in the population behavior. A more rigorous analysis of systems dynamics can be performed to evaluate when the system is in equilibrium or cycling in parameter space, which is out of the scope of this paper. Third, causal models can be considered to further examine the direction and causal relationships between components of IMB and stigma, since we primarily focused on correlations between these measures. Lastly, for the network analysis, we took a measured approach to establish a general trend for the population dynamics. However, more complex network measures such as information transmission and node similarity can be employed to characterize the details of stigma spread through the social network, instead of focusing on population averages in IMB and stigma.