Descriptive statistics
Table 1 shows that among the deceased elderly, 51.7% were women and 48.3% were men. A significant proportion of the sample passed away at an advanced age, with 26.6% aged between 80-89 years and 41.0% aged between 90-99 years. In terms of educational attainment, the majority had no schooling (60.4%), while 31.7% had primary education. The data shows that the majority reside in towns or villages (87.8%), with only a small percentage living in cities (12.2%). Additionally, 23.6% were married at the time of death, while 76.4% had other marital statuses. Roughly equal proportions participated in old-age insurance (48.2%) and did not (51.8%). Most (73.6%) lived with others, while 26.4% lived alone. The majority (88.3%) passed away at home, with the remaining 11.7% passing away in other locations. Among those who passed away, 64.5% suffered serious illness before death, 18.4% were completely self-sufficient, and 24.5% were unable to take care of themselves.
Table 2 shows a significant change in the severity of loneliness among older individuals from 2008 to 2014. In 2008, 37.8% reported no feelings of loneliness, a percentage that decreased to 29.2% in 2014. The majority of those reporting loneliness felt 'rarely lonely' or 'sometimes lonely', accounting for 54.1% in 2008 and increasing to 62% in 2014. The minimum proportion (8.2%) of older people reported frequent or constant feelings of loneliness and showed little change (8.8%) by 2014. Table 3 shows a decline in the proportion of people rating their health as 'very good' and 'good' from 52.5% in 2008 to 37.9% in 2014. Conversely, the proportion of those rating their self-assessed health as 'poor' and 'very poor' increased from 14.6% in 2008 to 22.2% in 2014.
Table 1 Descriptive statistics of the sample
|
subcategory
|
N
|
Percentage/mean
|
N
|
/
|
1423
|
/
|
Live alone or not
|
Yes
|
376
|
26.4
|
|
No
|
1047
|
73.6
|
place of death
|
Home
|
1257
|
88.3
|
|
Others
|
166
|
11.7
|
Seriously ill before death
|
Yes
|
505
|
35.5
|
|
No
|
918
|
64.5
|
Self-care ability
before death
|
/
|
1423
|
3.6
|
Sex
|
Male
|
687
|
48.3
|
|
Female
|
736
|
51.7
|
Age of Death
|
65-69
|
2
|
0.2
|
|
70-79
|
152
|
10.7
|
|
80-89
|
376
|
26.6
|
|
90-99
|
580
|
41.0
|
|
100 above
|
304
|
21.5
|
Number of years of schooling
|
No Schooling
|
859
|
60.4
|
|
Primary School
|
451
|
31.7
|
|
Junior Middle School
|
76
|
5.5
|
|
High School
|
25
|
1.8
|
|
University and above
|
10
|
0.7
|
Residence
|
City
|
173
|
12.2
|
|
Others
|
1250
|
87.8
|
Marital Status
|
Married
|
336
|
23.6
|
|
Others
|
1087
|
76.4
|
Pension insurance
before death
|
Participated
|
686
|
48.2
|
|
Non-participation
|
737
|
51.8
|
Amount of pension
|
Less than¥1,000
|
708
|
72.7
|
|
¥1000-¥1999
|
96
|
10.5
|
|
¥2000-¥2999
|
89
|
8.5
|
|
¥3000and above
|
81
|
8.3
|
Table 2 Change in loneliness among older persons, 2008-2014 (%)
Degree of loneliness
|
2008
|
2011
|
2014
|
Always lonely
|
2.0
|
1.6
|
1.6
|
Often lonely
|
6.2
|
5.6
|
7.2
|
Sometimes lonely
|
21.4
|
23.4
|
28.2
|
Rarely lonely
|
32.7
|
36.5
|
33.8
|
Never lonely
|
37.8
|
32.8
|
29.2
|
Table 3 Change in self-rated health of older persons, 2008-2014 (%)
Self-rated health
|
2008
|
2011
|
2014
|
Very good
|
12.5
|
10.4
|
5.8
|
Good
|
40.0
|
33.2
|
32.1
|
Average
|
32.9
|
36.3
|
39.9
|
Bad
|
13.3
|
19.0
|
19.2
|
Very bad
|
1.3
|
1.2
|
3.0
|
Group-based growth model
To better reflect the internal heterogeneity of changes in loneliness, this section introduces the group-based growth model to investigate the trend of loneliness among older adults. Table 4 shows that the three-category fit outperformed the two-category one, as indicated by the significant LMR and BLMR, along with increasing Entropy, while the AIC, BIC, and aBIC continued to decrease. The AIC, BIC, and aBIC values increased in category four compared to category three, indicating that category four was less well-fitted. In category five, although the AIC, BIC, and aBIC values are on the decline, the LMR of the model in category five was not significant. Therefore, the three-category model in LCGM was selected as the analytical tool for subsequent analyses in this paper.
Table 4 LCGM model fitting information
Categories
|
LL
|
AIC
|
BIC
|
aBIC
|
LMR
|
BLRT
|
Entropy
|
2C
|
-5839.62
|
11695.25
|
11737.33
|
11711.92
|
0.0013
|
0.0000
|
0.61
|
3C
|
-5799.18
|
11620.37
|
11678.23
|
11643.29
|
0.0404
|
0.0000
|
0.68
|
4C
|
-5764.35
|
11626.37
|
11700.01
|
11655.54
|
0.0384
|
0.0000
|
0.75
|
5C
|
-5743.79
|
11521.58
|
11611.00
|
11557.00
|
0.0564
|
0.0000
|
0.75
|
Fig. 1 illustrates the classification of changes in loneliness among older people into three groups using the LCGM model. The largest group (61.9%) is characterized as the 'not lonely group,' which has maintained low levels of loneliness from the first measurement in 2008 to the final measurement in 2014. The second group, the 'rising loneliness group', accounts for the second-largest share (30.8%). Loneliness was low in the first survey (2008), but it gradually increased in both the second (2011) and third (2014) surveys, i.e. older people were becoming 'lonelier'. The third and smallest group (7.3%) is the decreasing loneliness group. Members of this group had a high level of loneliness in the first survey (2008) but saw a gradual decrease in loneliness gradually over the second survey (2011) and the third survey (2014). It is evident that the majority do not feel lonely. However, there is an increase over time in the number of older people experiencing loneliness, with minimal reduction among those whose loneliness decreases.
Regression analysis
Table 5 shows that increased loneliness among older adults correlates with a higher likelihood of being chronically bedridden at the end of life compared to those in the non-lonely group. Surprisingly, falling loneliness did not decrease but rather increased the probability of chronic bedriddenness in older adults. Conversely, those not lonely exhibited a significantly lower probability of chronic bedriddenness at the end of life compared to the other groups. Among the control variables, the probability of being bedridden at the end of life was significantly higher for those seriously ill compared to those not. Model 2's -2log likelihood decreased upon the addition of the variable 'serious illness before dying', indicating an enhanced explanatory power compared to Model 1. The significance levels for the prediction of long-term bedriddenness decreased in both the groups that experienced an increase and a decrease in loneliness. The significance of loneliness in both groups could be partially explained by the component of 'serious illness at the end of life'. The -2loglikelihood of Model 3 continued to decrease with the inclusion of 'self-rated health', indicating a further enhancement in the explanatory power of the model.
In Model 4, older adults categorized in the non-lonely group underwent a significant reduction in their end-of-life self-care inability when 'serious illness before dying' and 'self-rated health' were not factored in. In Model 5, the inclusion of 'serious illness before dying' gave rise to a significant decrease in the -2 log likelihood, indicating an enhancement in the model's explanatory power. Similarly, in Model 6, the inclusion of 'self-rated health' led to a continued decrease in -2log likelihood and a decrease in the significance level of the non-loneliness group's prediction of long-term bed rest. The significance level for the non-lonely group's prediction of long-term bedriddenness also declined. Among the control variables, older adults who lived alone exhibited a lower probability of chronic bedriddenness at the end of life. In conclusion, the impact of changes in loneliness on chronic bedriddenness at the end of life became less significant when the influence of 'self-rated health' or 'serious illness at the end of life' was factored in. This suggests that the physical health of older adults may serve as a potential pathway through which changes in loneliness contribute to chronic bedriddenness at the end of life.
Table 5 Effect of changes in loneliness on prolonged bedriddenness at end of life
|
Model 1
|
Model 2
|
Model 3
|
Model 4
|
Model 5
|
Model 6
|
Sex
(Female)
|
-0.205
|
-0.223
|
-0.219
|
-0.236
|
-0.236
|
-0.249
|
Years of Schooling
|
0.001
|
0.001
|
0.001
|
-0.001
|
-0.001
|
0.001
|
Urban residence
(town and rural area)
|
0.173
|
0.172
|
0.189
|
0.155
|
0.155
|
0.172
|
Age of Death
|
0.005
|
0.008
|
0.009
|
0.004
|
0.004
|
0.009
|
Married
(Other
marital status)
|
0.106
|
0.072
|
0.073
|
0.11
|
0.078
|
0.079
|
Participation in
pension insurance
(non-participation)
|
-0.035
|
-0.061
|
-0.066
|
-0.04
|
-0.065
|
-0.07
|
Amount of pension
|
0.000
|
0.000
|
0.000
|
0.000
|
0.000
|
0.000
|
Live alone
(Non-living alone)
|
-0.275+
|
-0.257+
|
-0.271+
|
-0.268+
|
-0.251+
|
-0.266+
|
Died at home
(Other places)
|
0.197
|
0.235
|
0.247
|
0.178
|
0.214
|
0.228
|
Illness before dying
(No illnesses)
|
|
0.392*
|
0.372**
|
|
0.391**
|
0.371**
|
Self-rated health
|
|
|
0.124.
|
|
|
0.133
|
Loneliness Rising Group
(Not Alone Group)
|
0.229**
|
0.229*
|
0.204*
|
|
|
|
Loneliness Decline Group
(Not Alone Group)
|
1.214***
|
1.216**
|
1.177**
|
|
|
|
Not Alone Group
(Other Groups)
|
|
|
|
-0.364**
|
-0.365**
|
-0.335*
|
Constant
|
-0.91
|
-1.334
|
-1.797
|
-0.622
|
-0.884
|
-1.417
|
-2loglikelihood
|
1303.98
|
1295.96
|
1294.23
|
1313.62
|
1305.53
|
1303.48
|
Note:+ P<0.1;*P<0.05;**P<0.01;***P<0.001
As shown in Table 6, increased loneliness and decreased loneliness significantly weakened the end-of-life self-care ability of the elderly, and the absence of loneliness significantly increased this ability. Compared with Model1 and Model2, the R of the model increased significantly upon the inclusion of "serious illness before dying," and the significance level of the effect of increasing loneliness and decreasing loneliness on the elderly's self-care ability began to decrease. Model 3 showed results similar to Model 2 without significant changes. In terms of the other control variables, men had a better end-of-life self-care ability compared to women, and older people who lived alone displayed a superior end-of-life self-care ability.
In Model 4, the probability of end-of-life self-care inability is significantly lower for older adults in the non-loneliness group when "serious illness before dying" and "self-rated health" are not included. The addition of "serious illness before dying" in Model 5 significantly increases R, suggesting an increased explanatory power of the model. Model 6 showed that the significance level of the prediction of loneliness on the end-of-life self-care ability of older adults decreased significantly after the addition of "self-rated health". In conclusion, after adding "self-rated health" or "serious illness before dying", the significance of changes in loneliness on elderly people's end-of-life self-care ability is weakening, i.e., elderly people's physical health may serve as a potential pathway through which changes in loneliness affect their end-of-life self-care ability.
Table 6 Effect of changes in loneliness on self-care at the end of life
|
Model 1
|
Model 2
|
Model 3
|
Model 4
|
Model 5
|
Model 6
|
Sex
(Female)
|
-0.378*
|
-0.408*
|
-0.402*
|
-0.396*
|
-0.425**
|
-0.418**
|
Years of Schooling
|
0.023
|
0.024
|
0.025
|
0.022
|
0.023
|
0.024
|
Urban residence
(town and rural area)
|
0.085
|
0.082
|
0.102
|
0.076
|
0.072
|
0.094
|
Age of Death
|
0.002
|
0.007
|
0.010
|
0.002
|
0.007
|
0.010
|
Married
(Other marital status)
|
-0.119
|
-0.176
|
-0.175
|
-0.116
|
-0.174
|
-0.173
|
Participation in
pension insurance
(non-participation)
|
0.286*
|
0.242
|
0.236
|
0.284+
|
0.239
|
0.233
|
Amount of pension
|
-0.000
|
-0.000
|
-0.000
|
-0.000
|
-0.000
|
-0.000
|
Live alone
(Non-living alone)
|
-0.507**
|
-0.474**
|
-0.491**
|
-0.506**
|
-0.473**
|
-0.490**
|
Died at home
(Other places)
|
0.035
|
0.098
|
0.114
|
0.024
|
0.088
|
0.105
|
Illness before dying
(No illnesses)
|
|
0.691***
|
0.663***
|
|
0.693***
|
0.665***
|
Self-rated health
|
|
|
0.167+
|
|
|
0.173+
|
Loneliness Rising Group
(Not Alone Group)
|
0.224**
|
0.222*
|
0.187*
|
|
|
|
Loneliness Decline Group
(Not Alone Group)
|
0.736*
|
0.718*
|
0.661*
|
|
|
|
Not Alone Group
(Other Groups)
|
|
|
|
-0.299*
|
-0.294*
|
-0.255+
|
Constant
|
3.305
|
2.561
|
1.940
|
3.652
|
2.900
|
2.216
|
R
|
0.182
|
0.231
|
0.237
|
0.176
|
0.226
|
0.233
|
Note:+ P<0.1;*P<0.05;**P<0.01;***P<0.001
Path analysis
The above analyses suggest that changes in loneliness may have an impact on older people's end-of-life self-care ability, with this influence potentially mediated by health status. Therefore, in this section, the variable of "self-rated health" was introduced to further explore the underlying mechanism through path analysis. First, as shown in Fig. 2, the highest proportion of older adults (61.9%) belongs to the non-lonely group, characterized by better self-rated health. This significantly reduces the probability of end-of-life bedriddenness and significantly improves the end-of-life self-care ability of older adults. Secondly, the previous analysis revealed that a large proportion (30.8%) of older persons were "becoming lonely", which, as shown in Fig. 3, led to a decline in self-rated health, which further increased the probability of bedriddenness and end-of-life self-care inability.
Finally, older people with declining loneliness are the least represented of the three categories (7.3%), and this group of older people has a unique pattern of changes in loneliness compared to the other two categories. Fig. 4 shows a positive correlation between self-rated health and declining loneliness. The analysis of the previous two categories of older adults suggests that declining loneliness in older adults is accompanied by declining self-rated health, a conclusion inconsistent with established facts. One possible explanation for the decline in loneliness among older people is a passive decline due to much attention from family members as a result of their poor self-rated health. This form of care, although passive, reduces older people's feelings of loneliness. Regarding end-of-life autonomy, the reduction in loneliness does not positively impact older individuals due to disadvantages resulting from a prolonged history of poor self-assessed health. As a result, older individuals may still experience low end-of-life autonomy.