Hospital capacity for medical tourists and local patients
Amirkabir University of Technology · paper under review at Health Care Management Science
With patient behaviour modelled, the mean inpatient admission wait is 14.8 days; the otherwise identical discrete-event-only model shows 25.8 days (30 and 36 replications, Welch t-test).
Problem
A multi-specialty international-patient department serves medical tourists alongside local patients, and how those patients behave (abandoning, switching doctors, choosing a channel) affects the queues.
Approach
A hybrid agent-based and discrete-event simulation in AnyLogic that models channel choice, physician switching, treatment abandonment and deterioration. A 256-run fractional factorial design screened 16 controllable factors across six performance measures, followed by a replicated comparison against an otherwise identical discrete-event-only model.
Result
Bed capacity is the most consistent driver of waiting; priority rules, channel design and appointment timing also matter. Leaving patient behaviour out changes the answer: the discrete-event-only model shows a 25.8-day admission wait against 14.8 days in the hybrid model (42% lower, p < 0.001).
Healthcare operationsSimulationAnyLogicDesign of experiments
Constrained boosting for reduced-order models
Georgia Institute of Technology
An unconstrained leaf value falls outside the valid region; cXGBoost projects it into the intersection of the balls during fitting.
Problem
Predicting a simulation's reduced basis at unseen parameters means regressing a point on a manifold, valid only inside a bounded region.
Approach
cXGBoost enforces that bound inside every leaf of the boosting ensemble rather than correcting afterwards, so predictions are valid by construction and no general-purpose optimiser is called at any split.
Result
Lower mean reconstruction error than the projected Gaussian process and Grassmann interpolation on four parametric problems (cylinder wake, Kolmogorov flow, the Kuramoto–Sivashinsky equation and a forced beam). Released as an open-source Python package.
Actions lead to states; a loss seen after a delay is credited back only to the actions that could have produced that state.
Problem
When an action's outcome only arrives after a delay, most algorithms wait or discard the signal.
Approach
State-EXP3 charges one delayed outcome to every action that could have produced it, and pays for that pooling with an effective dimension that reads the overlap between actions rather than counting them, so the cost tracks how much the actions actually share.
Result
Regret reduced by up to 79% versus action-level weighting and by 32–68% versus a rate-optimal action-only delayed-bandit baseline, with a lower bound showing the cost of delayed feedback cannot be eliminated by pooling.
From the observed modalities, either certify the prediction and act, or abstain and acquire the next modality.
Problem
Clinical models often have to act on incomplete inputs: imaging but no labs, notes but no genomics.
Approach
RouteCert certifies a prediction when the evidence already in hand determines the decision, and otherwise abstains and says which modality to acquire next. Acquisition is separated from calibration, so conformal coverage holds conditional on the terminal modality pattern rather than on average.
Result
Distribution-free coverage guarantees, evaluated on four datasets spanning ECG, speech and wearable sensing (PTB-XL, CMU-MOSEI, MHEALTH, IEMOCAP).
Conformal predictionMultimodal learningMissing dataClinical decision support
Surrogates that keep the data's shape
Georgia Institute of Technology · submitted to the Journal of Computational Physics
Simulation output is multiway; instead of flattening it, each mode keeps its own subspace, interpolated separately, with cores aligned before fitting.
Problem
Simulation output is naturally multiway: space by space by time. The usual move is to flatten it into one long vector, which throws that structure away.
Approach
Every mode gets its own low-dimensional subspace, each interpolated on its own manifold, with the cores aligned before a single low-rank surrogate is fitted.
Result
Across five PDE benchmarks in two and three dimensions, including public CFD data, the interpolation size is 13× to 482× smaller and cached queries run 3.3× to 8.5× faster than a single-subspace surrogate built on the flattened field, with mean errors within 0.22 percentage points. In a 3D grid study at 64³ the size ratio reaches 2,260× and the speedup 11.9×.
Tensor decompositionsSurrogate modelingRiemannian geometryScientific ML
Other projects
LLM agents as simulated patients: prompted three dated OpenAI model snapshots, with a fixed
prompt, a strict JSON schema and cached replies, to act as 36 synthetic patient profiles inside a
queueing simulation of eight capacity policies (120 paired seeds). Calibrated to the same baseline
rates, the models agreed at baseline but not on the best policy under a demand surge; a follow-up with
2,160 sampled actions found that simple behavioral rules picked the same policy, so the study stopped
at its go/no-go gate. Python, OpenAI API · 2026
Retrieval-augmented claim verification: a local RAG pipeline over 5,183 SciFact abstracts
comparing BM25, two dense retrievers, rank fusion and reranking, with an open 1.5B LLM returning JSON
verdicts with citations. Hybrid retrieval raised nDCG@10 from 0.678 to 0.724 over BM25 (p = 0.0008);
retrieval lifted accuracy on claims with evidence but almost removed correct abstentions, and one
distractor abstract cut faithful rationales from 0.251 to 0.142.
Python, PyTorch, Hugging Face · 2026
Probing what language models encode: trained probes at three depths on residual-stream
activations of Pythia models from 70M to 1.4B parameters, testing every finding against shuffled-label
nulls and matched controls; day-of-week order appeared from 160M parameters up (rank correlation about
0.91) and never in the controls. Python, PyTorch, Transformers · 2026
Predictive maintenance and anomaly detection: end-to-end pipeline for sensor-based
fault detection: ARIMA and LSTM forecasting, ensemble classifiers over multivariate telemetry, and
hybrid Weibull and state-space degradation modelling with online retraining, giving over 15 hours of
advance fault warning. Deployed on Spark for low-latency inference across sensor networks.
Python, PySpark · 2025
Multimodal phenotype modeling: built representation-learning pipelines over single-cell
transcriptomics, longitudinal EHR and imaging-derived phenotypes, using contrastive learning for
patient and cellular embeddings and transformer language models to extract structured phenotypes
from unstructured clinical notes. Python, PyTorch · 2024–25
Drug safety and clinical risk: detected adverse drug events from longitudinal EHR with
temporally-aware patient embeddings over medication trajectories, labs and comorbidities; compared
regularized regression, tree ensembles and neural models, and audited calibration and fairness
across demographic strata. Python · 2024–25
Sensitivity-guided diffusion models: constrained diffusion generative models to physically
meaningful latent subspaces, deriving sensitivity-informed perturbation schedules that steer the
forward SDE along high-variance directions to improve sample efficiency in data-scarce scientific
regimes. Python, PyTorch · 2024–25
Anomaly detection in sensor streams: built ARMA-based statistical monitoring with
sliding-window control charts to flag faults in industrial sensor streams early, trading off
detection delay against false-alarm rate. Python · 2023
Loan default prediction: trained and compared ensemble classifiers for consumer credit
risk, evaluating ranking performance alongside probability calibration rather than accuracy alone.
Python · 2023
Safe operating range under hybrid inputs: built a six-degree-of-freedom physics-based
simulator in MATLAB/Simulink from system velocity equations and operational constraints, generated
data across operating scenarios, and trained neural network models in Python to predict system
behaviour and expand the feasible operating range.
MATLAB, Simulink, Python · 2018–20
Additional projects
Internet access and female economic empowerment: cross-national nonparametric analysis of
internet usage against female employment and representation in senior management across 92 countries
(World Bank Gender Data Portal, 2022); Mann–Whitney U tests with Bonferroni correction and
Anderson–Darling / Shapiro–Wilk normality diagnostics. Group project
Breast cancer prediction: built and statistically evaluated diagnostic classification
pipelines, comparing generalization across models rather than training accuracy.
Python, scikit-learn · 2021–22
Nurse scheduling: formulated and solved mixed-integer programs for healthcare workforce
scheduling under demand uncertainty. GAMS · 2020–21
Healthcare operations quality: applied Six Sigma to locate and reduce service quality
deficiencies in a healthcare organization. 2021–22
Food recommender expert system: knowledge-based recommendation tailored to individual
health characteristics and dietary needs. Python · 2021–22
Tensor-structured learning pipelines: deep and tensor-structured architectures for
parameterized scientific fields under geometric constraints.
Python, PyTorch · 2024
Supersonic airfoil CFD: meshed a supersonic airfoil in Gambit and simulated the flow
in Ansys Fluent to capture the shock.
Gambit, Ansys Fluent · 2021–22