The data consist of 400 stroke patients who were treated under follow-up in FHRH, Bahir Dar, Ethiopia. Among four hundred patients who had stroke taken from the repository of FHRH, 188(47%) were females. The patients mean age and corresponding 95%CI was (58; 95%CI = 55–58). Around 210(53%) of patients were from rural. Regarding the previous history of other diseases(comorbidities), around 166(42%) of the patients were hypertensive, 127(32%) were diabetic, 104 (26%) had heart diseases, and 122 (31%) of the patients had an ischemic stroke respectively (See Table 2).
Table 2: Frequency distribution of predictors variables
Characteristic
|
Label
|
n (%)
|
Gender
|
Gender
|
|
|
female (0)
|
188 (47%)
|
|
male (1)
|
212 (53%)
|
Residence
|
Residence of patients
|
|
|
rural (0)
|
210 (53%)
|
|
urban (1)
|
190 (48%)
|
Hypertension
|
History of hypertension
|
|
|
no (0)
|
234 (59%)
|
|
yes (1)
|
166 (42%)
|
Diabetes
|
History of diabetes
|
|
|
no (0)
|
273 (68%)
|
|
yes (1)
|
127 (32%)
|
AF
|
History of atrial fibrillation
|
|
|
no (0)
|
285 (71%)
|
|
yes (1)
|
115 (29%)
|
Hdiseases
|
History of heart disease
|
|
|
no (0)
|
296 (74%)
|
|
yes (1)
|
104 (26%)
|
Tstroke
|
type of stroke
|
|
|
hemorrhage (0)
|
278 (70%)
|
|
ischemic (1)
|
122 (31%)
|
|
|
|
STcomp
|
stroke complications
|
|
|
no (0)
|
290 (73%)
|
|
yes (1)
|
110 (28%)
|
TIA
|
History of transit ischemic attack
|
|
|
no (0)
|
362 (91%)
|
|
yes (1)
|
38 (9.5%)
|
HIV
|
History of HIV/AIDs
|
|
|
no (0)
|
391 (98%)
|
|
yes (1)
|
9 (2.3%)
|
Glasgow
|
Glasgow coma scale score
|
|
|
13-15 (0)
|
254 (64%)
|
|
3-8 (1)
|
100 (25%)
|
|
9-12 (2)
|
46 (12%)
|
Characteristic
|
Label
|
Mean(95%CI)
|
Age
|
age of the patients
|
56(55-58)
|
Hemoglobin
|
Hemoglobin level
|
14.30(14.09-14.51)
|
SBP
|
systolic blood pressure
|
133(131-135)
|
DBP
|
diastolic blood pressure
|
84(82-85)
|
BUN
|
blood urea nitrogen
|
30.56(29.10-32.03)
|
WBC
|
white blood cell count
|
9.24(8.79-9.68)
|
Cholesterol
|
Total cholesterol level
|
155.48(149.20-161.76)
|
Exploring Longitudinal Data using area percentage
Figure 1 shows that at baseline 14% were sample patients who had stroke had no symptoms and/or no disability despite symptoms and the percentage increased at every visit time, whereas the percent of slight and severe disability were 24% and 25% respectively and decrease over time. But the percentage of patients who had stroke having a moderate disability at baseline were 37% and changes through time with some degree of similarity.
Exploring number and percent of observations for functional status
Figure 2 shows the total number of observations for each ordinal outcome in terms of the stroke patient’s functional status. The x axis of the bar graph displays a cumulative number of treatment observations within each functional status status group. The bar graph showed that patients with a moderate disability had the highest cumulative number of treatment observations 1335(37.04%) and slight disability had lowest cumulative number of treatment observation 758(21.05%). And also, the study observed that patient with no symptom and sever/dead had 670(18.6%) and 841(23.35%) number of treatments.
Frequency Profile Plot for Functional status
Figure 3 shows that the total number of patients with a moderate disability was higher at the baseline, and then it strictly decreased after 3 months. Because they have more probability of experiencing dropout and death overtime and the same is true for slight disability However, the number of patients with no symptoms or slight disability despite symptoms increased until 12 months, and then it decreased monotonically.
Partial Proportional Odds Model (PPOM)
Univariate analysis for PPOM
In this sub-section a univariate analysis was done and then Age, diabetes, Stroke complication, blood urea nitrogen, white blood cell count, and Glasgow coma scale were selected.
Proportional Odds Assumption
Based on POM assumption checking Age, Stroke complication, blood urea nitrogen, and Glasgow coma scale. Had significant interaction with threshold value This indicates that the proportional odds assumption is violated for those variables.
Multivariable Analysis for PPOM
Multivariable analysis was done by using all the significant predictors together by assuming a partial proportional odds model (See Table 3).
Table 3 showed that age, blood urea nitrogen, Glasgow coma scale, and stroke complication had statistically significant interaction effects with thresholds. indicated that the effects of the predictors can’t be assumed identical across the three cumulative logits. By holding this assumption, the result can be interpreted as follows.
age has a significant effect on the odds of patients who had stroke being in the “no symptoms or no disability despite symptoms” functional status (aOR = 0.97; 95%CI: 0.96–0.97), meaning for a 1-year increase in patient's age the chance of being in the “no symptoms or no disability despite symptoms” functional status reduced by 3%. Similarly, for a 1-year increase in patient's age the chance of being in the “no symptoms or no disability despite symptoms or slight disability” functional status increased by 2%(aOR = 1.02;95%CI = 1.01–1.03), controlling the effect of all other variables in the model. Furthermore, age has different effects and directions across thresholds, indicating that as age increases the level of disability also increases or decreases the level of functional status.
The odds of being in “no symptoms or no disability despite symptoms” of functional status for 1 unit increase in blood urea nitrogen is 0.98(aOR = 0.98;95%CI = 0.97–0.99) and 1.03(aOR = 1.03;95%CI = 1.02–1.04) for being in the “no symptoms or no disability despite symptoms or slight disability”, keeping all other variables in the model as constant.
For the Glasgow coma scale, the odds of being no symptoms or no disability despite symptoms versus slight disability, moderate, or severe/dead with GCS 3–8 is 1.62(aOR = 1.62;95%CI = 1.25–2.1) times higher than that of patients whose Glasgow coma scale 13–15. The odds of being no symptoms or no disability despite symptoms or slight disability versus moderate or severe/dead with GCS 3–8 is 1.36(aOR = 1.36;95%CI = 1.01–1.83) times higher than that of patients whose
Table 3: Estimate and 95% CI for the PPOM.
Effect
|
Proportional odds
|
Non-proportional odds
|
|
Exp (β) (95%CI)
|
Exp (β) (95%CI)
No vs slight, moderate, severe/death
|
Exp (β) (95%CI)
No, slight vs moderate, severe/death
|
Exp (β) (95%CI)
No, slight, moderate vs severe/death
|
Thresholds
|
|
|
|
|
No Symptom
|
0.09(0.01,0.99)*
|
0.09(0.01,0.99)*
|
-
|
-
|
Slight disability
|
8.04(5.08, 12.73)***
|
-
|
8.04(5.08,12.73)***
|
-
|
Moderate disability
|
194269.2(74500.42,506581.25)***
|
-
|
-
|
194269.2(74500.42,506581.25)***
|
Predictors
|
|
|
|
|
Age
|
-
|
1.11(1.06, 1.13)***
|
1(0.99, 1.01)
|
1.06(1.04, 1.07)***
|
Bun
|
-
|
1.04(1.01, 1.07)*
|
0.97(0.96, 0.98)***
|
1(0.99, 1.01)
|
glasgow:3-8 vs 13-15
|
-
|
4.15(1.08, 15.99)*
|
0.64(0.45, 0.9)*
|
3.86(2.44, 6.17)***
|
glasgow:9-12 vs 13-15
|
-
|
10.86(1.84, 63.94)**
|
0.61(0.44, 0.85)**
|
0.8(0.45, 1.4)
|
time slopes
|
|
|
|
|
diabetes: yes vs no
|
1.78(1.69,1.88)***
|
-
|
-
|
-
|
WBC
|
0.99(0.98,0.995)***
|
-
|
-
|
-
|
stroke complications: yes vs no
|
-
|
0.88(0.84,0.92)***
|
1.17(1.08,1.15)***
|
1.24(1.16,1.30)***
|
*p < 0.05, **p < 0.01, ***p < 0.001.
BUN, Blood Urea Nitrogen; PPOM, Partial Proportional Odds Model; WCC, White Cell Count
Glasgow coma scale 13–15. The odds of being no symptoms or no disability despite symptoms, slight disability, moderate versus severe/dead with GCS 3–8 is 0.28(aOR = 0.28;95%CI = 0.19–0.39) times lower than that of patients whose Glasgow coma scale 13–15. Similarly, 0.49(aOR = 0.49;95%CI = 0.34–0.72) times lower to be being no symptoms or no disability despite symptoms versus slight disability, moderate, or severe/dead, 2.44(aOR = 2.44;95%CI = 1.84–3.22) times higher to be being no symptoms or no disability despite symptoms or slight disability versus moderate or severe/dead, and 3.78(aOR = 3.78;95%CI = 2.36–6.11) times higher to be being no symptoms or no disability despite symptoms, slight disability, or moderate versus severe/dead comparing with patients whose Glasgow coma scale 13–15.
Considering the interaction effects of the covariates with time, since we didn’t reject the proportional odds assumption for WBC and diabetes, therefore, we assumed an identical effect across three cumulative logit models for those covariates. for every white blood cell increase with follow-up time the odds of being in the first, second, and third category are the same and it is1.01(1.01;95%CI = 1.007–1.02). the odds of diabetic patients who had stroke with follow-up time being in “no symptoms or no disability despite symptoms versus slight disability, moderate, or severe/dead”, in “no symptoms or no disability despite symptoms, slight disability versus moderate, severe/dead”, and in “no symptoms or no disability despite symptoms, slight disability, moderate versus severe/dead” are the same(aOR = 1.96;95%CI = 1.89–2.1).
Finally, by considering the time slope with nonproportionality, the odds of patients with stroke complications being in the “no symptoms or no disability despite symptoms versus slight disability, moderate, or severe/dead” is 1.1(aOR = 1.1;95%CI = 1.06–1.13), being in the “no symptoms or no disability despite symptoms slight disability versus moderate, or severe/dead” is 0.93(aOR = 0.93;95%CI = 0.9–0.96), and being in the “no symptoms or no disability despite symptoms, slight disability, or moderate versus severe/dead” is 0.89(aOR = 0.89;95%CI = 0.85–0.93) with follow-up time.
model Diagnosis
In fact, a crucial step in the modeling process is the application of diagnostic procedures for model checking. The adequacy of the fitted model was checked for the possible presence of outliers and influential values were presented in Fig. 4 for the appropriateness of the PPOM using Q-Q plot and Empirical Bayes (EB) estimates of the random effects in multilevel models to represent how individuals deviate from the population averages and are often extracted to detect outliers or used as predictors in the follow-up analysis. To sum up, figures given below declared the PPOM was adequately fitted the data (Fig. 4).
The global assumption of mixed proportionality of odds can be checked by LRT, where the deviance of the model assuming proportional odds \(\:({D}_{ev}PO\)) is compared to the deviance of the partial proportional odds model \(\:({D}_{ev}PPO\))
Under the null hypothesis, the asymptotic distribution of statistics in the above equation is chi-square with degrees of freedom given by l = (d.f PPO − d.f PO) [38].
Table 4 shows that: The LRT = 116.8091 (p < 0.001), indicating that the proportional odds assumption for the full model was not satisfied. This suggested that the effect of one or more of the explanatory variables was likely to differ across separate binary models fit for the cumulative thresholds.
Table 4: Checking the fit of the model for longitudinal model
Model
|
AIC
|
BIC
|
Deviance
|
chisq
|
P(>chisq)
|
Proportional odds
|
4892.416
|
4936.323
|
4870.416
|
|
|
Partial-proportional odds
|
4795.607
|
4879.428
|
4753.607
|
116.8091
|
<0.001
|