Skill: Pharmacodynamics & Dose-Response Recovery

Context

You are analyzing drug effects (e.g., Anticonvulsants, Psychedelics) on neural activity. Your goal is to reverse-engineer the "Effect Function" $E(C)$ from the dynamical changes observed in the EEG.

Strategy

We assume the brain dynamics obey $\dot{x} = A x + B(u) x$, where $u$ is the Drug Concentration (proxy) or Dose. We use SINDy with control inputs to identify the structure of $B(u)$.

Pipeline steps

1. Data Setup

• Dataset: Mendeley Pharmaco-EEG (Rat) or PsiConnect.
• Input:
  • $x(t)$: Neural features (e.g., Total Power, Complexity/Entropy).
  • $u(t)$: Drug concentration curve (often modeled as a simple decay $u(t) = D_0 e^{-kt}$ if blood samples unavailable).