The study was conducted in accordance with the Declaration of Helsinki and was approved by the institutional ethics committee of Jinhua Municipal Central Hospital (No. 2022129). The requirement for written informed consent was waived because only retrospective, de-identified data were analyzed. The research tools used in this protocol are listed in the Table of Materials.
1. Study design and source population
A retrospective observational cohort study was conducted at the comprehensive stroke center of Jinhua Municipal Central Hospital, Zhejiang, China. Consecutive admissions from June 2024 through September 2025 were screened. Patients were included if they had been admitted within 24 h of symptom onset, had a baseline NIHSS score of 5 or less, and had received a treating-team diagnosis of minor ischemic stroke compatible with a small subcortical or perforator-territory process. Patients were excluded for other documented etiologic subtypes, a baseline NIHSS score greater than 5, an onset-to-admission time greater than 24 h, severe pre-existing neurological dysfunction or severe systemic disease, or incomplete core clinical or imaging data. Every exclusion was recorded in the flow diagram.
2. Etiologic classification and brain imaging
The etiologic classification assigned during routine care by the treating stroke team was abstracted. Because this was a retrospective electronic-record study, the assignment was not considered prospective, centralized, blinded, or independently re-adjudicated. For MRI-confirmed cases, an acute subcortical lesion compatible with the neurological syndrome on diffusion-weighted imaging, with a corresponding apparent diffusion coefficient abnormality and no cortical infarction, was required. For CT-only cases, exclusion of hemorrhage and a cortical territorial infarct, together with a compatible lacunar clinical syndrome, was required. CT alone was acknowledged as potentially unable to directly confirm a small acute subcortical infarct.
Imaging modality, onset-to-imaging time, maximum lesion diameter, slice thickness, number of involved axial slices, ventral pontine surface involvement, parent-artery stenosis, and nonculprit stenosis were recorded. Suspected BAD was operationalized as any of the following: a maximum lesion diameter of at least 15 mm, involvement of at least three axial slices, or ventral pontine surface involvement. Detailed lesion side, anatomic location, responsible artery segment, and a traceable, independent, blinded reread were unavailable and were not imputed.
3. Angiographic assessment of stenosis
CTA or MRA stenosis measurements were abstracted from the clinical imaging reports. The percentage of intracranial stenosis was calculated using the WASID formula:
[1 - (diameter at the narrowest lumen
/diameter of the adjacent normal reference artery)] x 10015
Parent-artery stenosis was defined as narrowing of at least 50% in the artery supplying the infarct territory. Nonculprit stenosis was defined as narrowing of at least 50% outside that territory. The primary model variable, any vascular stenosis of at least 50%, was considered positive if either component was positive. Because the responsible artery segment was not retained, segment-specific measurements could not be independently verified.
4. Baseline clinical, laboratory, and treatment variables
Age, sex, vascular history, smoking, alcohol use, pre-stroke antiplatelet use, pre-stroke modified Rankin Scale score, onset-to-admission time, admission SBP and diastolic BP (DBP), baseline NIHSS, platelet count, blood glucose, neutrophil count, lymphocyte count, and blood-draw time were extracted. NLR was calculated as the absolute neutrophil count divided by the absolute lymphocyte count.
Intravenous thrombolysis, single versus dual antiplatelet strategy, treatment-start time after admission, statin intensity, in-hospital BP-lowering treatment, glucose management, anticoagulation, and the recorded timing of antithrombotic exposure relative to blood collection were reported. No participant received in-hospital anticoagulation. Intravenous fluid volume was not included because it was not reliably recoverable. These interventions were treated as potential confounders or consequences of the clinical state rather than randomized exposures.
5. Neurological assessments and END
NIHSS scores documented at admission, approximately 24 and 72 h were used, and additional assessments triggered by clinical worsening were included. The scores were recorded by the routine stroke team. Predictor blinding and centralized duplicate scoring were not documented. Primary END was defined as an increase of at least 2 total NIHSS points within 72 h after admission relative to baseline, consistent with contemporary lacunar-stroke literature2,3. The time and documented cause of deterioration were recorded. All-cause qualifying deterioration was included in the primary outcome, and the model was repeated using ischemic-cause END only.
6. Missing-data handling
Missingness was audited before analysis. For three records with a missing raw NLR value but available neutrophil and lymphocyte counts, NLR was calculated deterministically. For four records with a missing raw 24-h SBP SD but available mean SBP and coefficient of variation, SD was calculated as coefficient of variation × mean SBP. After these deterministic derivations, the four primary predictors and the outcome were complete, and all 240 patients were included in the primary model without multiple imputation. The single record with a missing pre-blood-draw antithrombotic status was excluded only from the treatment-adjusted exploratory model.
7. Statistical analysis and model development
Approximately symmetric continuous variables were summarized as mean ± SD, and skewed or ordinal variables were summarized as median (interquartile range). Welch t tests or Mann–Whitney U tests were used as appropriate. Categorical variables were summarized as n (%) and compared using chi-square or Fisher's exact tests.
The primary logistic regression model included continuous admission SBP scaled per 10 mmHg, baseline NIHSS per point, NLR per unit, and any vascular stenosis of at least 50%. The intercept, coefficients, standard errors, ORs, 95% CIs, and two-sided P values were reported. SBP nonlinearity was examined using a restricted cubic spline with three knots at the 10th, 50th, and 90th percentiles (130.9, 157.0, and 194.0 mmHg), adjusted for NIHSS, NLR, and stenosis. A likelihood-ratio test was used for the nonlinear term. The linear SBP term was used when nonlinearity was not supported.
An SBP threshold of 169 mmHg was retained only as a secondary simplified analysis. The threshold was selected using the maximum Youden index in the development data, and cutoff uncertainty was quantified using 1,000 bootstrap resamples. The unvalidated 0–12 score and low-, intermediate-, and high-risk categories were not reinstated.
Discrimination was estimated using the AUC and 95% bootstrap CI. The combined model was compared with each component using paired bootstrap differences, and superiority was not claimed when the difference CI included zero. Optimism was estimated using 1,000 bootstrap resamples. Out-of-bag bootstrap predictions were generated for calibration and calculation of the Brier score. Decision-curve analysis was conducted only across threshold probabilities of 5%–30%, for which the modeled action was enhanced neurological observation and clinician reassessment rather than automatic treatment.
Sensitivity analyses were performed in the MRI-only cohort, after exclusion of parent-artery stenosis of at least 50%, after exclusion of suspected BAD, and using ischemic-cause END only. Exploratory ridge-penalized expanded models containing imaging, onset-time, and treatment variables were fitted to reduce instability. In a 24-hlandmark analysis, patients who deteriorated before 24 hwere removed, and a landmark model was compared with models that added SBP SD, coefficient of variation (CV), or average real variability (ARV). These models were interpreted as dynamic post-admission predictions. Transparent prediction-model reporting principles were followed16. The reproducible revision analysis was performed using Python 3.13.5, pandas 2.3.3, NumPy 2.3.3, and SciPy 1.16.2.