We develop a conservative continuous-time stochastic control framework for treatment optimization from irregularly sampled patient trajectories. The unknown patient dynamics are modeled as a controlled stochastic differential equation with treatment as a continuous-time control. Naive model-based optimization can exploit model errors and propose out-of-support controls, so optimizing the estimated dynamics may not optimize the true dynamics. To limit extrapolation, we add a consistent signature-based MMD regularizer on path space that penalizes treatment plans whose induced trajectory distribution deviates from observed trajectories. The resulting objective minimizes a computable upper bound on the true cost. Experiments on benchmark datasets show improved robustness and performance compared to non-conservative baselines.
Figures & tables
Figure 1: Effect of increasing regularization strength λ . (Left) Improvement of the ground truth costs for λ={1,10,100} relative to the ground truth costs with no regularization, i.e., λ=0 on three optimization problems across 15 different initial conditions (Cancer (Explicit), Cancer (Relative), Covid-19 (Target Trajectory). Boxplots show medians and interquartile ranges (IQR). Whiskers extend to the farthest point within 1.5×IQR from boxes. (Right) Example of predicted optimal treatment plans for λ∈{0,1,10,100} together with treatment plans seen during training and ground truth optimal treatment plan.
Table 1: Benchmarking the optimized treatment plans on all three tasks; Cancer (Explicit), Cancer (Relative) and Covid-19 (Target Trajectory).
Figure 2: Example patient trajectories under treatment plans optimized with different conservatism levels ( λ ). Tumor volume over time for two optimized treatment plans. Shaded regions indicate variability across simulated trajectories, the bold curve denotes the mean trajectory, and the dot marks the mean final tumor volume.
Appendix figures & tables2 assets
Supplementary material from the paper’s appendix.
Appendix
Symbol
Meaning
I
I:=[ts,tf] , the considered finite time horizon
N
the number of training samples
D
the dataset D:={(z(i),y(i),u(i))}i∈[N]
dz
dimension of the stochastic covariate process Z:Ω×I→Rdz
dy
dimension of the stochastic outcome process Y:Ω×I→Rdy
du
dimension of the stochastic treatment process U:Ω×I→Rdu
Appendix
Table 2: Summary of notation
Figure 3: Effects of regularization strength parameter λ .