Modelling Social Care Provision in An Agent-Based Framework with Kinship Networks

TL;DR

This paper introduces an agent-based model simulating UK social care, incorporating kinship networks and socioeconomic factors, to better understand informal and formal care dynamics amid demographic shifts.

Contribution

The paper presents a novel agent-based framework that models social care decisions using kinship networks and socioeconomic data, capturing complex care provision patterns.

Findings

Model replicates real-world care patterns

Kinship networks influence care distribution

Socioeconomic factors impact care decisions

Abstract

Current demographic trends in the UK include a fast-growing elderly population and dropping birth rates, and demand for social care amongst the aged is rising. The UK depends on informal social care -- family members or friends providing care -- for some 50\% of care provision. However, lower birth rates and a graying population mean that care availability is becoming a significant problem, causing concern amongst policy-makers that substantial public investment in formal care will be required in decades to come. In this paper we present an agent-based simulation of care provision in the UK, in which individual agents can decide to provide informal care, or pay for private care, for their loved ones. Agents base these decisions on factors including their own health, employment status, financial resources, relationship to the individual in need, and geographical location. Results demonstrate that the model can produce similar patterns of care need and availability as is observed in the real world, despite the model containing minimal empirical data. We propose that our model better captures the complexities of social care provision than other methods, due to the socioeconomic details present and the use of kinship networks to distribute care amongst family members.

Tables

Table 1. Table 1: The different care need categories, with the number of hours of care required per week

Care need category	Weekly hours of care required
None	0
Low	8
Moderate	16
Substantial	32
Critical	80

Equations

$$R=K\sum _{i=1}^n{y}_i^p$$

$$w=F{e}^{c{e}^{-rt}}$$

$$c=ln\frac{I}{F}$$

$$e=\frac{C_p-C_b}{U_b-U_p}$$

Full text

K****eywords Social Care, Kinship Networks, Agent-Based Modelling

1 Introduction

As human lifespans continue to lengthen and birth-rates drop throughout much of the developed world, many nations are experiencing an increase in demand for social care – the provision of personal and medical care for people in need of assistance due to age, disability or other factors. In the UK, the elderly consume the largest share of social care, and lower birth-rates mean that the supply of available carers is decreasing over time even as demand is growing rapidly [1]. As a result, social care is a frequent topic in UK policy debates, with widespread concern that the country will be unable to afford the significant public investment needed to fulfill this growing care need.

Age UK notes that nearly 50% of those aged 75 and over are living with a long-term limiting health condition – and the fastest-growing age group in the UK is the over-85s, so the problem will only get worse in the coming decades [2]. As a result, social care need is rising at a pace that outstrips the growth of public and private care supply. Unmet social care need is therefore an increasingly important and widespread social problem in the UK. The 2017 Ipsos MORI report Unmet Need for Care showed that “more than half of those with care needs had unmet need for at least some of their needs” [3]. This means that less than half of those elderly individuals in need of assistance with activities of daily living (ADLs) were able to receive sufficient care. Large numbers of UK citizens are thus living without sufficient care; an estimated 1.2 million people did not receive sufficient care for their needs in 2017, an increase of 48% since 2010 [2].

The increased pressure on the social care system also has an impact on the health care system, with delayed discharges from hospital being a particularly expensive consequence of unmet care need. Patients in need of social care frequently stay in hospital longer than necessary due to a lack of care availability, leading to a shortage of beds for other patients and increased costs to hospitals. According to Age UK, between 2010 and 2016 the number of additional bed-days attributable a lack of available home-care packages increased by 181.7%, while wait times for residential care placements increased by 40% [2]. In 2016 the National Audit Office estimated the total delays attributable to care shortages amounted to 2.7 million bed-days per year, resulting in an annual cost to the National Health Service (NHS) of approximately £820 million [4].

The UK is largely dependent on informal social care, or care provided free-of-charge by family members and loved ones, in order to meet the needs of the population. Informal care is enormously widespread in the UK and is much larger than the formal care infrastructure. The Family Resources Survey 2013/14 showed that there were 5.3 million informal carers in the UK [5], while projections indicate that the number of people receiving informal care will increase by 60% in the period 2015-2035 [6]. The increase in demand is such that Carers UK proposes that carer numbers would need to increase 40% over the next two decades to meet demand [7]. At the same time, the number of people using privately-funded social care is expected to rise by almost 50% by 2035, and private expenditure for social care is projected to increase from £6.8 billion to almost £20 billion in the same period, almost a three-fold increase [6].

Understanding how the need for social care evolves and the process through which informal and formal care is provided is a vital component in developing sustainable social care policies. In this regard, the importance of support and care-giving networks has long been recognized [8]. In particular, research has shown that informal care is provided mostly through care networks with an average of three to five members [9]. These networks are predominantly composed of an individual’s close relatives. Aldridge and Huges [5] report that 72% of carers provide care to a member of their immediate family, whether a parent (40%), partner (18%), son or daughter (14%), while Petrie and Kirkup [10] show that 51% of carers provide care to someone in their household. Wettstein and Zulkarnain[11] show that, in 2011, 31% of informal care in the US was provided by spouses, 47% by children and 18% by other relatives (sons-in-law, daughters-in-law and grandchildren), with only 4% being provided by non-relatives.

Moreover, empirical research has shown that the kind of social care provided is affected by socioeconomic status. Petrie and Kirkup [10] report that people working in routine occupations and those with lower qualifications are more likely to provide informal care. Moreover, Laing and Buisson [12] find that, while across the UK an estimated 24% of care home residents are funded through private ‘top-ups’, in North-East England only 18% are funded privately, as opposed to 54% in the wealthier South-East. Overall, these findings suggest that informal care becomes less common in higher socioeconomic status groups than in lower groups, while formal care becomes more common.

Finally, the social care literature also reveals a significant gender gap in social care provision. Wettstein et al. [11] reported that in the US daughters are almost twice as likely (31%) to provide care as compared to sons (16%). These findings are confirmed by Petrie and Kirkup [10], who find that 59% of family carers in the UK are women. At the population level, 16% of women provide informal care as compared to 12% of men.

In this paper, we propose an agent-based model of the UK informal and privately-funded formal social care system. This model reflects the complexity of a system where demographic, social and economic processes interact to determine the dynamics of social care demand and supply. Our aim is to provide a theoretical framework that allows us to to improve our understanding of the mechanisms driving unmet social care need. Moreover, using ABMs enables us to model scenarios of economic and social policy change in virtual populations, providing a means to test complex policies targeted at multiple levels of society. Such models can allow policy-makers to experiment with potential policy interventions and reveal any possible unintended side-effects of those policies prior to implementing them in the real world.

Previous work has attempted to address the social care problem using agent-based simulation approaches [13] [14]. Here we present a simulation that significantly extends these previous efforts. The model provides a more comprehensive simulation of social care provision behaviour, via the inclusion of a detailed socioeconomic model and the representation of care provision as not just a simple one-to-one exchange of resources, but a complex negotiation taking place between members of the care receivers’ kinship networks.

2 The Model

2.1 Motivations

Here we present an agent-based model as a formal representation of the complex demographic and social processes affecting informal social care demand and supply, and the dynamics of the social care outcomes resulting from their interaction over time. Our primary concern at this stage, was not the precise replication of empirical, real-world data but, rather, the development of a theoretical framework potentially capable of representing the full complexity of the social care system. In our results we aim for qualitative similarity to real-world social care trends, rather than precise numerical replication, in order to determine if our modelling of the underlying processes is producing appropriate outcomes.

With these motivations in mind, at this initial stage the behavioural assumptions made in this paper should not be considered necessarily accurate, but rather should be seen as a first approximation, i.e. a necessary point of departure allowing us to provide proof-of-concept results and demonstrate the model’s potential as a policy-making tool. Further work and more specialized expertise will be needed to refine and revise the model’s behavioural assumptions in order for the model to be used in this way.

Ultimately, our aim is for this simulation to facilitate the development and evaluation of alternative social care policies, including complex interventions aimed at behavioural change. Our addition of health care costs and a detailed economic model to the simulation can also clarify the impact of interventions on other key areas of policy, enable more robust policy evaluation, and reduce the incidence of unintended consequences derived from otherwise well-meaning policy prescriptions. Future versions of this framework can also be applied to nations other than the UK, simply by altering the simulated geography and population data.

2.2 Basics of the model

The model itself is complex and contains many detailed economic and social processes and sub-processes; this section provides a high-level overview of the model’s functionality. For details of the previous models which formed the initial basis for this simulation, please refer to Noble et al. [13] and Silverman et al. [14]. 111For those who wish to examine the model more closely, or run it for yourself, you can find the annotated Python code available at https://github.com/UmbertoGostoli/ABM-for-social-care/releases/tag/v.09

The simulated agents occupy a space roughly based on UK geography. Agents live in houses which form towns containing clusters of up to 1225 houses. These clusters vary in size in rough proportion to real-world population density, which varies across the 8$\times$12 grid composing the model’s geography. The agent population is scaled down from real UK levels at a factor of roughly 1:10,000.

The model updates in one-year time steps. Initial populations are generated and randomly distributed in the year 1860, and the model then runs until 2040 when final social care costs are calculated and recorded. We start the model in 1860 to ensure the population dynamics have time to settle before empirical population data (UK Census data) is integrated into the model in 1951.

2.3 Agent Life-Course

Newborn agents are classed as dependent children until age 16, at which point they reach adulthood (the minimum working age in the UK). The agents then decide whether to continue to study or to start looking for a job. This choice is repeated every two years until the age of 24. This is a probabilistic choice that depends on the household’s income level and on the parent’s education level. When the agent decides to look for a job (or when it reaches age 24 at the latest), the agent then enters the workforce. During their working life, agents can be hired, fired and can change jobs. When an agent is unemployed, they find a job with a certain probability which depends on the unemployment rate (which is an input of the model), then start earning an income and paying tax. At a certain age, agents retire from the workforce (at age 65 when default parameter settings are used), at which point they stop paying tax and begin receiving a pension.

This version of the simulation uses a Gompertz-Makeham mortality model to approximate mortality rates until 1951, as in Noble et al. [13]. From 1951 the simulation uses mortality rates from the Human Mortality Database [15]. After 2009 a Lee-Carter model is used to generate future mortality rates, as in Silverman et al. [14].

2.3.1 Partnership Formation

Upon reaching adulthood, agents can form partnerships with other agents 222Hereafter we use ‘partnership’ as shorthand for relationships capable of producing children.. Employed male agents and adult female agents are randomly paired with probabilities based on the agents’ socioeconomic difference, age difference and geographical distance (with the relative weight of these three factors depending on the model’s parameters). Age-specific annual divorce probabilities determine whether a couple dissolves their partnership. Fertility rates follow the procedures outlined in Silverman et al. [14], in which data from [16] and the Office for National Statistics [17] are used from 1950–2009, at which point Lee-Carter projections are used.

2.3.2 Migration

As in the real world, agents can migrate for a variety of reasons. An agent becomes independent when they leave the parental home. This can happen when they form a partnership, when they find a job in a different town from the family home, or with a certain probability when they find a job in the same town of their parental home. When a partnership dissolves, the male agent will move elsewhere on the map, while any dependent children resulting from that partnership will stay with the mother. A family will relocate to a new house if one of the two parents finds a job in another town or when the family needs a larger house due to the increased size of the family. Retired agents with social care needs may elect to move into the household of one of their children, with a probability directly proportional to the amount of care provided by that household. Rarely, dependent children will be orphaned before reaching adulthood; in that case the agent will be adopted by a household of their kinship network or by a randomly-selected couple if no such household is available.

2.4 Health status and care need

Agents begin in a healthy state and have no need of additional care. They may enter different categories of care need depending on age- and gender-specific probabilities. Table 1 shows the five possible categories of care need and the amount of care required per week at each level of need. Agents who enter a state of care need do not recover and return to normal health, but instead continue to progress to higher levels of need over time. The probability of progressing to higher care need levels depends positively on the agent’s age and on the discounted sum of the agent’s unmet care need in past periods 333The assumption is that a prolonged period of unmet care will increase the agents’ frailty and, therefore, the probability that their health will further deteriorate..

Social care provision is linked to informal care availability in an agent’s kinship network. An agent’s kinship network’s nodes consist of the households of agents with a familial relationship with the agent; the degree of kinship is defined as the network distance between the household and the agent. If they have time or available income, agents will provide informal or formal care to anyone in their kinship network with care need. The amount of care agents are available to provide depends on their status (through their income), their kinship relationship with the receiver and their geographical distance from the receiver (see the Kinship Network sub-section below for details).

2.5 Model Enhancements

This updated version of the Linked Lives model presented in Silverman et al. [14] has been rewritten from the ground up and substantially extended for greater detail and realism. The following features have been introduced:

• •

The population is composed of five socioeconomic status groups.

• •

The care supply is provided by the agent’s kinship network, (a network of households which have a kinship relationship to the agent).

• •

Households allocate part of their income to care provision, which can be in the form of both informal and formal care.

• •

A salary function implying an inverse relationship between time taken off work to provide informal care and the agent’s hourly wage.

• •

Unmet social care needs affect the agents’ hospitalization probability (and the associated health care costs).

2.5.1 Socioeconomic Status Groups

The population is composed of 5 distinct socioeconomic status (SES) groups. These categories follow the Approximated Social Grade, a socioeconomic classification produced by the Office for National Statistics, which is composed of 6 categories (A, B, C1, C2, D and E). For convenience, we redistributed category E (state pensioners, casual and lowest grade workers, unemployed with state benefits only) into the categories D (semi-skilled and unskilled manual workers) and C2 (skilled manual workers), to maintain a unimodal distribution. We initialised the groups’ distribution to roughly reflect the 2016 UK distribution. The SES groups are characterized by different education levels, different career paths (represented by income growth curves) and different unemployment rates.

The introduction of SES groups has a number of effects on the various stages of agent life-courses. A higher socioeconomic position is associated with lower mortality and fertility rates, and to a lower probability of developing care need. In the marriage market, the probability that two opposite-sex individuals will form a couple depends on their SES distance (which determine the marriage probability together with the partners’ geographical distance and age difference)444The relationship between SES distance and probability of forming a couple is asymmetric: the probability of getting married decreases less rapidly with the SES distance if the higher-status individual is male rather than female.. In the job market, a higher SES is associated with higher starting and maximum salaries (but a lower salary growth rate); a higher probability to find a job and a lower probability to be fired (probabilities which are reflected in a lower unemployment rate); and a higher probability to change jobs if the job offer comes from a different town from their current hometown.

Given that in this model care supply and relocation decisions depend on income level, the socioeconomic position of an agent affects its behaviour (and that of the household it belongs to) through the agent’s income, as a higher socioeconomic position is associated with a higher income level. For example, the share of income the household allocates to care supply depends positively on the household’s per capita income (see the Formal Care section below).

We included an inter-generational mobility process which allows agents to move to a different SES group from their parents’. Each SES group is associated with an education level. From the age of 16, an agent can decide whether to continue its studies (a choice that will allow the agent to reach a higher education level and therefore a higher SES group) or start searching for a job (in which case the agent is assigned the SES group associated with the education level reached). This choice is made by the agents every two years, until the age of 24 (i.e. at ages 16, 18, 20 and 22) 555We assume each education step lasts 2 years. The educational stages correspond to A-level, Higher National Diploma, Degree and Higher Degree.. The probability that an agent keeps studying depends on three factors: the per capita available income (i.e, net of social care costs) of the agent’s household; the difference between the maximum education level reached by the agent’s parents and the agent’s current education level; and the amount of time the agent allocates to informal social care supply. In this highly stylized inter-generational mobility process, an agent’s SES group is determined by the education level the agent reaches666Meaning the SES group of agents ending their studies when they are 16, 18, 20, 22 and 24 year-old, will be, respectively, D, C2, C1, B and A..

2.5.2 Kinship Networks

Each agent is associated with a kinship network – a network of households containing at least one agent with a kinship relation to them. This network includes:

• •

the agent’s household (distance 0);

• •

the households of the agents’ parents and of children who are not part of the agent’s household (distance I);

• •

the households of the agents’ grandparents, grandchildren, brothers and sisters who are not part of the agent’s household nor of the households at distance I (distance II);

• •

the households of the agents’ uncles, aunts, nephews and nieces who are not part of the agent’s household nor of the households at distance I and II (distance III).

The kinship network has two functions. First, the network defines the total care supply available for a particular agent in need. The care receiver’s total care supply is the sum of the available care supply of all the members of all the households that are part of the kinship network which met certain conditions. These conditions represent kinship- and space-specific limitations in the supply of care. Kinship relationship and distance determine the quantity of informal care that each member of the receiver’s kinship network is available to supply, as shown in Table 2.5.2:

[htb] Amount of care agents can provide depending on their status and distance from the receiver

Agent status Household (D-0) D-I D-II D-III

Teenager 16 0 0 0

Student 24 16 8 4

Employed 28* 20* 12* 8*

Unemployed 32 24 16 8

Retired 48 36 20 10

• •

*These are the minimum amounts the employed can offer, representing the informal care provided outside of the working hours. If needed, additional informal care can be supplied by these agents taking time off their working hours. Moreover they can use their income to pay for formal care. See the Formal Care section for details.

In particular, we assume first that the amount of care which an agent is available to provide to another agent depends positively on their kinship’s closeness. Agents living with the care receiver are an exception, as these agents’ available care supply is assumed to be equal to that of the receiver’s next of kin (i.e. spouse, children and parents), independently from the kinship degree (although most of the time the members of the receiver’s household are his/her next of kin). The second factor affecting the informal care supply availability is the physical distance from the receiver’s household. With regards to the physical distance we can distinguish three classes of care suppliers:

• •

agents living in the receiver’s household;

• •

relatives living in another household in the same town as the care receiver;

• •

relatives living in a different town.

We assume that only households living in the same town of the care receiver can provide informal care. Moreover, in the case of formal care we assume that only the care receiver’s household and the households at distance I (i.e., the care receiver’s parents or children) are available to provide formal care.

The care allocation process proceeds in a series of steps in which a 4-hour ‘quantum’ of care supply is transferred to an agent needing care from one of the households with available supply in its kinship network. The care allocation function first samples a care receiver from the pool of people with unmet social care need who have a kinship network with available care supply. The probability of each care receiver being sampled is directly proportional to the care receivers’ quantity of unmet care need. Then the allocation function samples a household from the selected care receiver’s kinship network with a probability that is directly proportional to a distance-weighted measure of the household’s available supply 777For an household at kinship network’s distance x from the care receiver, the weight is the reciprocal of the exponential function of the product of x and a parameter.. Once the supplying household has been selected, a ’quantum’ of care is transferred from one of the household’s members with available supply to the care receiver.

The member who is to provide care within the selected household is determined in two steps. First, one of the household’s six possible care sources is selected with a probability proportional to the residual available care of each source. The six care sources consist of the five groups that can provide informal care as shown in Table 2: teenager, student, retired, unemployed and employed; plus a sixth source which represents the amount of care which the household is available to supply by allocating part of the household income (a category we call out-of-income care). The household’s out-of-income care supply is the share of income that the household has available to allocate to care, either directly in the form of formal care, or indirectly in the form of informal care provided by employed household members (who provide care by taking unpaid time off work). If one of the first five sources is selected, the household member with the greater residual available care within that category will provide care. If the out-of-income category is selected, the household will provide formal care if the lowest hourly wage among the hourly wages of the employed household members is higher than the hourly cost of social care. Otherwise, the employed household member with the lowest hourly wage will provide the ‘quantum’ of care888Given that the quantity of out-of-income care supply depends on the household’s per-capita income, this care allocation mechanism implies that the probability that members of the households will spend time providing informal care is inversely related to the household’s per-capita income. In other words, the wealthier the supplying household, the more likely that it will provide formal care.

At the end of this step, both the care receiver’s care need and the supplying household’s availability are reduced by the amount of care transferred (i.e. the 4-hour ‘quantum’ of care). The allocation process is repeated until the set of care receivers with unmet care need and available care supply is empty (i.e. there are no more care receivers with outstanding care need and available care supply). The ultimate result of this care allocation function is that units of social care are distributed from potential providers to receivers across the receiving agents’ kinship network, and decisions about who provides care and choosing informal or formal care provision are made according to the composition, socioeconomic position and employment status of potential care providers. We suggest this detailed decision-making process better represents the complexities of care decisions that families need to navigate. The specifics of this complex process can be adjusted further by incorporating insights received from qualitative and quantitative data on informal care-givers and their decision-making.

The kinship network also allows households to compute the informal care attraction associated with each town, which represents a rough measure of the informal care that a household expects to get or supply in a given town. The households use this town-specific attraction to make relocation decisions 999In this model, agents face the choice of relocating to other towns either if they receive a job offer or partner with someone from another town.. More precisely, for a particular household $H_{i}$, the informal care attraction of a town $T$, is a function of the sum of the members of the households with a kin relationship with H, living in the town $T$ 101010The household’s kinship network is obtained by joining the kinship networks of the household’s members.. The contribution of a household $H_{j}$ to the social care attraction is weighted according to the degree of kinship between $H_{i}$ and $H_{j}$. This weighted sum is then multiplied by the complement to one of the share of care provided by the government. Two assumptions characterise the towns’ informal care attraction: first, the larger the local kinship network in a town (in terms of number of people who are part of it), the higher the amount of social care the agent can expect to receive (or to supply) and the higher the social care ‘value’ associated with a town; second, the higher the share of care supplied by the government, the lower the importance of the potential informal care in the relocation decision.

Apart from this ‘network size’ factor, a household’s probability of relocating depends on the relocation cost, which increases with the household’s size and the number of years the household’s members lived in the current town 111111The relocation cost R can be thought of as a measure of the social capital developed by the household in their current town. This social capital would be lost by relocating to another town, so it acts as a barrier to relocation. Formally, relocation cost is computed as:

$R=K\sum\limits_{i=1}^{n}y_{i}^{p}$

where n is the set of the household’s members, yi is the number of years the household member i spent in the current town and p is a parameter with value smaller than 1 (as additional years will increase the member’s social capital by increasingly smaller amounts). K is a scaling parameter., and the town’s homophily attraction, which depends on the town’s share of people belonging to other SES groups. Towns with a larger share of unoccupied houses are more likely to be chosen for relocation.

2.5.3 Formal Care

Both informal and formal care are allocated through the care recipient’s kinship network. Each household allocates a share of its income to care, which increases with the household’s per capita income 121212More precisely, with x being the household’s per-capita income, the share allocated to care supply is the complement to one of the reciprocal of the exponential function of the product of x and a parameter.. We assume that formal care can be provided only by the care receiver’s own household or from the households in the care receiver’s kinship network with distance equal to 1131313In other words, only the care receiver’s parents and children are available to pay for formal care if they cannot provide informal care.. As explained above, the choice between the informal and formal care is stochastic, with probabilities equal to the relative availability of the two kinds of care. In order for the household to supply formal care, the hourly wage of the working member with the lowest wage (and available time left) must be higher than the price of formal social care. If the hourly wage of this member is lower than the price of formal social care, they will take time off work to provide informal care (as in this case it is cheaper for the household to give up the agent’s salary than pay for formal care). However, the supplying household can provide informal care only if it is in the same town of the care receiver, otherwise it will provide only formal care.

2.5.4 Salary function.

We model the hourly salary that an agent receives (on average) when it finds a job using the Gompertz function. This function is a double exponential that takes the following three SES-specific arguments:

• •

the initial salary level;

• •

the final salary level;

• •

the salary growth rate.

The salary growth rate is multiplied by the agent’s cumulative work experience, which is the discounted sum of all the fractions of the working week allocated to work (if an agent works full time, this fraction is equal to 1) 141414Formally, the salary function is:

$w=Fe^{ce^{-rt}}$

where:

$c=ln{\frac{I}{F}}$

and F is the maximum hourly wage, I is the initial hourly wage, r is the wage growth rate and t is the discounted cumulative work experience..

This equation implies that if an agent takes time off work to provide informal care, this will result in less work experience and, therefore a lower hourly salary. Given the properties of the care allocation mechanism, the agents employed with the lowest hourly salary will be more likely to provide informal care in the future (if there are not retired/students/unemployed agents or if the household’s retired/students/unemployed agents do not have available informal care supply). In this model new mothers devote their time to caring for their newborn, meaning that we see a gender pay gap emerge due to the interaction between:

• •

the allocation of part of the household’s income to care supply;

• •

the choice of the care giver within the supplying household;

• •

the salary function.

Retired agents receive a pension which is proportional to their final income level. We assume that when an agent needs care they retire due to sickness and thereafter receive a pension. If their care need is low or moderate, their ill-health pension is their normal pension reduced in accordance with the ratio between lost working years and the maximum number of working years (in other words, the working years of a person in good health). If the agent’s care need is substantial or critical, the ratio is computed by reducing the number of lost working years by 50%.

2.5.5 Hospitalisation

We assume that agents with care need will spend additional time in the hospital, time whose duration is a function of the agents’ care need level and the average discounted share of unmet care need. The higher the agent’s care need level and the higher her average share of unmet care need, the greater the number of days the agent is expected to spend in the hospital each year. By multiplying the sum of the hospitalisation’s duration by the hospitalisation cost per day, we can determine the cost of unmet care need for the public health service.

2.5.6 Simulation steps

More specifically, the simulation unfolds through the following eleven steps:

deaths: with a given probability which depends on age, SES and care need level, some agents are removed from the population; 2. 2.

adoptions: children without parents are adopted; 3. 3.

births: with a probability which depends on age and SES, married woman give birth to new agents; 4. 4.

divorces: some couples dissolve (and the male relocates); 5. 5.

marriages: with a certain probability which depends on age, SES and geographical location, some singles get married and go to live together. 6. 6.

social care allocation: social care is transferred from agents with available hours (or income) to agents with social care needs; 7. 7.

age transitions: the age of the agents is incremented, and their age-related status is updated. 8. 8.

social transitions: students decide whether to start working or keep studying, depending on their family income per capita, their parents’ SES and their care responsibilities. In case they start looking for a job, they are assiged the SES associated to the education level they have reached. 9. 9.

job market: employed and unemployed agents receive new job offers (and eventually accept them) and employed agents are fired, with probabilities which depend on the age and SES-specific unemployment rate. 10. 10.

relocations: agents who have accepted job offers from towns different from their current town, relocate. 11. 11.

care transitions: with a probability which depends on age, SES, current care need level and unmet care need, agents pass to higher care need levels or are hospitalised.

3 Social Policy Experiments

ABM allows us to investigate the effects of virtual social or economic policies. This model is characterised by various policy-related parameters, the value of which can be directly or indirectly related to specific measures of social care policy. Therefore, by varying these parameters, we can examine ‘what if’ scenarios and simulate social care outcomes and costs under alternative social care policies.

For illustrative purposes, we investigated the effect on social care of two policies: the introduction of tax-deductible social care expenses, and direct governmental funding of care for people above a certain level of care need. For the first policy, we assume that the carers are allowed to deduct 100% of social care expenses from their tax. In the second policy experiment, we assume that the government directly pays for the social care expenses of those agents with the highest social care need (i.e., care need level at ‘critical’, see Table 1).

We assume that the two policies are implemented from simulation year 2020 and compare the outputs of these two policy scenarios with the benchmark no-policy scenario in the period 2020-2040. We present two results: the hours of unmet care need in each scenarios and the incremental cost-effectiveness ratio (ICER) of the two policies considered 151515In this paper we formally define the ICER as:

$e=\frac{C_{p}-C_{b}}{U_{b}-U_{p}}$

where C is the total additional cost of policy implementation, U is the total unmet care need and the subscripts p and b indicate two scenarios (i.e. the policy scenario and the benchmark scenario, respectively). In this case, being b the no-policy scenario, Cb is equal to [math]. .

4 Results

As described above, this simulation is highly complex with numerous processes at play, and a number of parameters governing system behaviour. For these early-stage results, we present figures from a representative single run at default parameter values plus a comparison of this run, taken as the benchmark, with two simulations where we introduce two alternative policy interventions in the year 2020. The single-run charts displayed here were chosen to highlight the key features of this new simulation framework, namely the modelling of formal care, additional economic and labour market details such as socioeconomic status groups, interaction between social care need and health care demand, and agent kinship networks.

Figure 2 shows the evolution of hours of informal and formal care delivered and of unmet social care need per week for the whole population over the period 1990-2040. In our simulations, the total population reaches a peak of around 10,000 agents in the period 2025, and then decreases slowly to around 8000 by 2040. As in Bijak et al. [18], the simulation’s population projections roughly match those of the Office for National Statistics in the UK if we assume no international migration.

In Figure 2 we can see that from 2010 informal care cannot keep pace with the growth in care need, with a consequent rapid increase of unmet care need, which reaches a peak in 2025. In Figure 2 we show the weekly hours of informal and formal care delivered and the unmet social care need per recipient. Here we can see that the negative trend of the unmet care need after 2025 in Figure 2 is entirely due to the decrease of population: in fact, the unmet care need per recipient keeps growing quite steadily up to 2033.

Figure 4 shows the mean informal and formal care received and the unmet care need by care need level over the period 2020-2040 (note that the heights of the bars correspond to the weekly hours of care required as shown in Table 1). We can see that most of the unmet care need is due to the agents with the highest care need. Figure 4 shows the mean informal and formal care received and the unmet care need by SES group (with I being the poorest SES group). We can see that the relative weight of informal care with respect to formal care decreases from the poorest to the wealthiest SES group, a result which reflects empirical findings. As expected, by comparing the green bars in Figure 4, we can see that the poorest SES group has a somewhat higher level of unmet care need compared to the wealthiest group but also that unmet care need is a significant share of total care need across SESs.

Figure 6 shows the share of total informal care provided by women. At the population level this figure starts from just below 75% in 1990, then decreases steadily until 2020, when it starts fluctuating just above 60%, a value which is in line with empirical findings. Moreover, there is a marked difference between SES groups, with the poorest SES group having the highest (and growing) gender inequality. As shown in Figure 6, this inequality in care provision affects income inequality between genders. At the population level, female agents’ income is about 5% lower than the male’s income. Therefore our model can explain at least part of the gender pay gap.

Figure 8 shows the total informal and formal care provided in the period 2020-2040 according to the kinship distance between caregiver and receiver: household (distance 0); parents and children (distance 1); grandparents, grandchildren and brothers (distance 2); aunts, uncles and nephews (distance 3). We can see that most social care is supplied within the household, however a significant part comes from outside the household, especially formal care. Figure 8 shows annual per capita health care cost due to the hospitalization of people with care needs. We can see that although the total unmet care need levels off after 2020 (as shown in Figure 2), the per capita health care cost keeps growing until the late 2030s. This is due to the decrease in population after 2015, resulting in the per capita health care cost reaching its peak 15 years after the unmet social care need.

In the last two charts we investigate and compare the efficacy of two policy interventions. These two policy experiments are illustrative examples that show how this model may be used as a tool for the design and evaluation of social care policies. In these examples we consider a tax deduction policy, in which the households are allowed to deduct all the social care expenses from their tax base, and a direct funding policy, in which the social care needs of people with the highest need level (the ‘critical’ level in Table 1) are taken care of directly by the state. The policies are implemented starting in simulation year 2020, and the two scenarios are compared to the benchmark no-policy scenario.

In Figure 10 we can see that the two policies have a positive impact on total unmet social care need, as expected. However, while the tax deduction policy has a marginal effect, decreasing the total unmet care need by about 11%, the direct funding policy has a drastic effect, reducing the total unmet care to 23% of the original level. Although the tax deduction policy has a smaller effect on the unmet social care need, Figure 10 shows that the tax deduction policy is more cost-effective than the direct funding policy.

This is because, while the effect of direct funding on total unmet social care is seven times the effect of the tax deduction, the cost of the former policy compared to the latter is much higher, as shown in Figure 12. We note however that notwithstanding its cost-effectiveness, the effect of the tax deduction policy on unmet care demand is limited and thus it is not a solution to the growing social care demand on its own. Finally, Figure 12 shows the discounted value of the hospitalization costs in the period 2020-2040 in the three policy scenarios. We can see that, apart from the direct costs associated with each policy, our model allows us to estimate the policies’ spillover effects (in this case, the effect on hospitalization cost), effects which we need to take into account to assess the relative advantages and disadvantages of alternative policies.

5 Discussion

While the results presented here are still early, the outcomes of these simulation runs suggest that this model can produce broadly realistic portraits of the coming trends in UK social care. The overall population dynamics largely mirror those in evidence in the real world, with the notable exception of international migration which is not modelled here. The simulation also produces inequalities in care provision by both SES and gender; while the data presently available is not sufficient to determine if the socioeconomic inequalities are accurate, the gender pay gap is reflective of current UK norms.

The simulation results on unmet care need lend credence to the worries expressed by UK policy-makers. The simulation shows that unmet care need will continue to grow over time, and given that 1.2 million older people in England alone are not receiving sufficient care [2], further growth in these figures could lead to severe consequences for public health. The results also show marked inequalities in care provision, with the wealthy capable of supporting their aged relatives in need through privately-funded care while still staying in work and receiving high wages. Among the lower-income groups women are providing an overwhelmingly larger share of informal care as compared to men, meaning that women at the lower end of the socioeconomic scale are more likely to be pushed out of work and toward unpaid care provision.

Our use of kinship networks as the mechanism for distributing care illustrates that kinship distance impacts care provision behaviour, and thus is an important aspect of care decision-making that should be taken into account in future research. We observed a significant difference in caring behaviours between within-household and kinship distance I agents, with the latter providing much more formal care than informal. This suggests that future models in this area should consider kinship networks and their impact on caring behaviours. Understanding the negotiation process within families regarding care provision will be an important aspect for policy-makers to examine as well, as policies put into place to support and encourage informal care may need to take account of these complex social aspects of care.

Finally, the speculative policy scenarios we provide here demonstrate the efficacy of this platform as an aid to policy-makers who wish to examine the impact and possible unintended side-effects (spillover effects) of their planned policy interventions. The scenarios showed that while some policies such as tax deductions may seem like an easy ‘win’ for the policy-maker, the actual impact on social care need is minimal. For significant reductions in care need and related health care costs, more expansive – and expensive – interventions are needed. The direct funding scenario here is a relatively simplistic example used to demonstrate this point; the model is capable of simulating the results of much more nuanced policy interventions, due to the detailed modelling of key social care mechanisms and related economic and social processes.

6 Future Work

While the current simulation does broadly reflect the expected trends in social care provision, validating the results is difficult. Ongoing projects in Scotland are linking administrative data sources to develop a clearer picture of social care provision and receipt across the country. In future revisions of the model we intend to make use of these data by focusing the simulation on Scotland rather than the whole UK. As the framework matures and incorporates more real-world data, we will use Gaussian process emulators to perform detailed sensitivity analyses, following the example of Silverman et al. [14].

Apart from the validation process, this model can be expanded in various directions. First, according to the Office for National Statistics the proportion of informal childcare to GDP increased from 13.8% to 17.6% in the period 2005-2014 [19]. Child care represents a significant part of the total household’s informal care need and, therefore, it affects the care supply which can be allocated to social care needs. Our next planned update to the model will include a child care mechanism which will allow us to model the interaction with social care supply.

Second, according to the Office for National Statistics, in 2016, 28.2% of births in England and Wales were to women who were not born in the UK [20]. Moreover, projections show that post-2016 immigration accounts for 77% of total population growth until 2041 [21]. The addition of international migration to the model will allow us to understand UK demographic dynamics in future decades.

Following on from the simplistic policy comparison presented here, in future work we will consult with social care policy-makers in Scotland to model proposed social care policies in more detail. Social care policy is very complex, with the programmes available varying not only by region (England, Wales, Scotland and Northern Ireland), but by individual local councils. The model is capable of representing these details by varying relevant policy parameters across the simulated UK geography. Once we are able to replicate the current state of social care policy in the UK, we will then be able to develop detailed evaluations of the possible impacts of new policy interventions. Our model’s ability to capture detailed individual-level care decision-making will enable policy-makers to examine policies aimed at behavioural change as well as broader economic and social policies.

Finally, while the current simulation is very UK-centric, the core simulation engine can be altered easily to examine the situation in other countries. Ultimately we hope to produce a simulation framework that is capable of modelling informal and formal social care across a variety of cultural and economic contexts.

7 Data Availability

The code is available on GitHub:

https://github.com/UmbertoGostoli/ABM-for-social-care/tree/v.09

The output data used for the figures is available in this Dryad repository:

https://datadryad.org/review?doi=doi:10.5061/dryad.6tm6183

8 Authors’ Contributions

Umberto Gostoli developed the model (based on a previous version, developed by Eric Silverman), run the simulations, gathered the data and produced the charts. Eric Silverman conceived the research and helped develop the model. Both authors contributed to draft the paper. All authors gave final approval for publication.

9 Funding

The authors are part of the Complexity in Health Improvement Programme supported by the Medical Research Council (MC_UU_12017/14) and the Chief Scientist Office (SPHSU14).

10 Acknowledgements

We would like to thank Chris Patterson and Lauren White from the MRC/CSO Social and Public Health Sciences Unit for their helpful comments.