Training the next generation of pharmacometric modelers: a multisector perspective.

The current demand for pharmacometricians outmatches the supply provided by academic institutions and considerable investments are made to develop the competencies of these scientists on-the-job. Even with the observed increase in academic programs related to pharmacometrics, this need is unlikely to change in the foreseeable future, as the demand and scope of pharmacometrics applications keep expanding. Further, the field of pharmacometrics is changing. The field largely started when Lewis Sheiner and Stuart Beal published their seminal papers on population pharmacokinetics in the late 1970's and early 1980's and has continued to grow in impact and use since its inception. Physiological-based pharmacokinetics and systems pharmacology have grown rapidly in scope and impact in the last decade and machine learning is just on the horizon. While all these methodologies are categorized as pharmacometrics, no one person can be an expert in everything. So how do you train future pharmacometricians? Leading experts in academia, industry, contract research organizations, clinical medicine, and regulatory gave their opinions on how to best train future pharmacometricians. Their opinions were collected and synthesized to create some general recommendations.

Characterization of anti-drug antibody dynamics using a bivariate mixed hidden-markov model by nonlinear-mixed effects approach.

Biological therapies may act as immunogenic triggers leading to the formation of anti-drug antibodies (ADAs). Population pharmacokinetic (PK) models can be used to characterize the relationship between ADA and drug disposition but often rely on the ADA bioassay results, which may not be sufficiently sensitive to inform on this characterization.In this work, a methodology that could help to further elucidate the underlying ADA production and impact on the drug disposition was explored. A mixed hidden-Markov model (MHMM) was developed to characterize the underlying (hidden) formation of ADA against the biologic, using certolizumab pegol (CZP), as a test drug. CZP is a PEGylated Fc free TNF-inhibitor used in the treatment of rheumatoid arthritis and other chronic inflammatory diseases.The bivariate MHMM used information from plasma drug concentrations and ADA measurements, from six clinical studies (n = 845), that were correlated through a bivariate Gaussian function to infer about two hidden states; production and no-production of ADA influencing PK. Estimation of inter-individual variability was not supported in this case. Parameters associated with the observed part of the model were reasonably well estimated while parameters associated with the hidden part were less precise. Individual state sequences obtained using a Viterbi algorithm suggested that the model was able to determine the start of ADA production for each individual, being a more assay-independent methodology than traditional population PK. The model serves as a basis for identification of covariates influencing the ADA formation, and thus has the potential to identify aspects that minimize its impact on PK and/or efficacy.

A disease model predicting placebo response and remission status of patients with ulcerative colitis using modified Mayo score.

The Mayo Clinical Score is used in clinical trials to describe the clinical status of patients with ulcerative colitis (UC). It comprises four subscores: rectal bleeding (RB), stool frequency (SF), physician's global assessment, and endoscopy (ENDO). According to recent US Food and Drug Administration guidelines (Ulcerative colitis: developing drugs for treatment, Guidance Document, https://www.fda.gov/regulatory-information/s. 2022), clinical response and remission should be based on modified Mayo Score (mMS) relying on RB, SF, and ENDO. Typically, ENDO is performed at the beginning and end of each phase, whereas RB and SF are more frequently available. Item response theory (IRT) models allow the shared information to be used for prediction of all subscores at each observation time; therefore, it leverages information from RB and SF to predict ENDO. A UC disease IRT model was developed based on four etrolizumab phase III studies to describe the longitudinal mMS subscores, placebo response, and remission at the end of induction and maintenance. For each subscore, a bounded integer model was developed. The placebo response was characterized by a mono-exponential function acting on all mMS subscores similarly. The final model reliably predicted longitudinal mMS data. In addition, remission was well-predicted by the model, with only 5% overprediction at the end of induction and 3% underprediction at the end of maintenance. External evaluation of the final model using placebo arms from five different studies indicated adequate performance for both longitudinal mMS subscores and remission status. These results suggest utility of the current disease model for informed decision making in UC clinical development, such as assisting future clinical trial designs and evaluations.

Joint longitudinal model-based meta-analysis of FEV1 and exacerbation rate in randomized COPD trials.

Model-based meta-analysis (MBMA) is an approach that integrates relevant summary level data from heterogeneously designed randomized controlled trials (RCTs). This study not only evaluated the predictability of a published MBMA for forced expiratory volume in one second (FEV1) and its link to annual exacerbation rate in patients with chronic obstructive pulmonary disease (COPD) but also included data from new RCTs. A comparative effectiveness analysis across all drugs was also performed. Aggregated level data were collected from RCTs published between July 2013 and November 2020 (n = 132 references comprising 156 studies) and combined with data used in the legacy MBMA (published RCTs up to July 2013 - n = 142). The augmented data (n = 298) were used to evaluate the predictive performance of the published MBMA using goodness-of-fit plots for assessment. Furthermore, the model was extended including drugs that were not available before July 2013, estimating a new set of parameters. The legacy MBMA model predicted the post-2013 FEV1 data well, and new estimated parameters were similar to those of drugs in the same class. However, the exacerbation model overpredicted the post-2013 mean annual exacerbation rate data. Inclusion of year when the study started on the pre-treatment placebo rate improved the model predictive performance perhaps explaining potential improvements in the disease management over time. The addition of new data to the legacy COPD MBMA enabled a more robust model with increased predictability performance for both endpoints FEV1 and mean annual exacerbation rate.

Improved Confidence in a Confirmatory Stage by Application of Item-Based Pharmacometrics Model: Illustration with a Phase III Active Comparator-Controlled Trial in COPD Patients.

PURPOSE: The current study aimed to illustrate how a non-linear mixed effect (NLME) model-based analysis may improve confidence in a Phase III trial through more precise estimates of the drug effect. METHODS: The FULFIL clinical trial was a Phase III study that compared 24 weeks of once daily inhaled triple therapy with twice daily inhaled dual therapy in patients with chronic obstructive pulmonary disease (COPD). Patient reported outcome data, obtained by using The Evaluating Respiratory Symptoms in COPD (E-RS:COPD) questionnaire, from the FULFIL study were analyzed using an NLME item-based response theory model (IRT). The change from baseline (CFB) in E-RS:COPD total score over 4-week intervals for each treatment arm was obtained using the IRT and compared with published results obtained with a mixed model repeated measures (MMRM) analysis. RESULTS: The IRT included a graded response model characterizing item parameters and a Weibull function combined with an offset function to describe the COPD symptoms-time course in patients receiving either triple therapy (n = 907) or dual therapy (n = 894). The IRT improved precision of the estimated drug effect compared to MMRM, resulting in a sample size of at least 3.64 times larger for the MMRM analysis to achieve the IRT precision in the CFB estimate. CONCLUSION: This study shows the advantage of IRT over MMRM with a direct comparison of the same primary endpoint for the two analyses using the same observed clinical trial data, resulting in an increased confidence in Phase III.

Improved Decision-Making Confidence Using Item-Based Pharmacometric Model: Illustration with a Phase II Placebo-Controlled Trial.

This study aimed to illustrate how a new methodology to assess clinical trial outcome measures using a longitudinal item response theory-based model (IRM) could serve as an alternative to mixed model repeated measures (MMRM). Data from the EXACT (Exacerbation of chronic pulmonary disease tool) which is used to capture frequency, severity, and duration of exacerbations in COPD were analyzed using an IRM. The IRM included a graded response model characterizing item parameters and functions describing symptom-time course. Total scores were simulated (month 12) using uncertainty in parameter estimates. The 50th (2.5th, 97.5th) percentiles of the resulting simulated differences in average total score (drug minus placebo) represented the estimated drug effect (95%CI), which was compared with published MMRM results. Furthermore, differences in sample size, sensitivity, specificity, and type I and II errors between approaches were explored. Patients received either oral danirixin 75 mg twice daily (n = 45) or placebo (n = 48) on top of standard of care over 52 weeks. A step function best described the COPD symptoms-time course in both trial arms. The IRM improved precision of the estimated drug effect compared to MMRM, resulting in a sample size of 2.5 times larger for the MMRM analysis to achieve the IRM precision. The IRM showed a higher probability of a positive predictive value (34%) than MMRM (22%). An item model-based analysis data gave more precise estimates of drug effect than MMRM analysis for the same endpoint in this one case study.

Detecting placebo and drug effects on Parkinson's disease symptoms by longitudinal item-score models.

This study tested the hypothesis that analyzing longitudinal item scores of the Unified Parkinson's Disease Rating Scale could allow a smaller trial size and describe a drug's effect on symptom progression. Two historical studies of the dopaminergic drug ropinirole were analyzed: a cross-over formulation comparison trial in 161 patients with early-stage Parkinson's disease, and a 24-week, parallel-group, placebo-controlled efficacy trial in 393 patients with advanced-stage Parkinson's disease. We applied item response theory to estimate the patients' symptom severity and developed a longitudinal model using the symptom severity to describe the time course of the placebo response and the drug effect on the time course. Similarly, we developed a longitudinal model using the total score. We then compared sample size needs for drug effect detection using these two different models. Total score modeling estimated median changes from baseline at 24 weeks (90% confidence interval) of -3.7 (-5.4 to -2.0) and -9.3 (-11 to -7.3) points by placebo and ropinirole. Comparable changes were estimated (with slightly higher precision) by item-score modeling as -2.0 (-4.0 to -1.0) and -9.0 (-11 to -8.0) points. The treatment duration was insufficient to estimate the symptom progression rate; hence the drug effect on the progression could not be assessed. The trial sizes to detect a drug effect with 80% power on total score and on symptom severity were estimated (at the type I error level of 0.05) as 88 and 58, respectively. Longitudinal item response analysis could markedly reduce sample size; it also has the potential for assessing drug effects on disease progression in longer trials.

Item Response Model Adaptation for Analyzing Data from Different Versions of Parkinson's Disease Rating Scales.

PURPOSE: The aim of this work was to allow combination of information from recent and historical trials in Parkinson's Disease (PD) by developing bridging methodology between two versions of the clinical endpoint. METHODS: A previously developed Item Response Model (IRM), that described longitudinal changes in Movement Disorder Society (MDS) sponsored revision of Unified Parkinson's Disease Rating Scale (UPDRS) [MDS-UPDRS] data from the De Novo PD cohort in Parkinson's Progression Markers Initiative, was first adapted to describe baseline UPDRS data from two clinical trials, one in subjects with early PD and another in subjects with advanced PD. Assuming similar IRM structure, items of the UPDRS version were mapped to those in the MDS-UPDRS version. Subsequently, the longitudinal changes in the placebo arm of the advanced PD study were characterized. RESULTS: The parameters reflecting differences in the shared items between endpoints were successfully estimated, and the model diagnostics indicated that mapping was better for early PD subjects (closer to De Novo cohort) than for advanced PD subjects. Disease progression for placebo in advanced PD patients was relatively shallow. CONCLUSION: An IRM able to handle two variants of clinical PD endpoints was developed; it can improve the utilization of data from diverse sources and diverse disease populations.

Sensitivity of Pegfilgrastim Pharmacokinetic and Pharmacodynamic Parameters to Product Differences in Similarity Studies.

In this work, a previously developed pegfilgrastim (PG) population pharmacokinetic-pharmacodynamic (PKPD) model was used to evaluate potential factors of importance in the assessment of PG PK and PD similarity. Absolute neutrophil count (ANC) was the modelled PD variable. A two-way cross-over study was simulated where a reference PG and a potentially biosimilar test product were administered to healthy volunteers. Differences in delivered dose amounts or potency between the products were simulated. A different baseline absolute neutrophil count (ANC) was also considered. Additionally, the power to conclude PK or PD similarity based on areas under the PG concentration-time curve (AUC) and ANC-time curve (AUEC) were calculated. Delivered dose differences between the products led to a greater than dose proportional differences in AUC but not in AUEC, respectively. A 10% dose difference from a 6 mg dose resulted in 51% and 7% differences in AUC and AUEC, respectively. These differences were more pronounced with low baseline ANC. Potency differences up to 50% were not associated with large differences in either AUCs or AUECs. The power to conclude PK similarity was affected by the simulated dose difference; with a 4% dose difference from 6 mg the power was approximately 29% with 250 subjects. The power to conclude PD similarity was high for all delivered dose differences and sample sizes.

A Novel Method for Analysing Frequent Observations from Questionnaires in Order to Model Patient-Reported Outcomes: Application to EXACT® Daily Diary Data from COPD Patients.

Chronic obstructive pulmonary disease (COPD) is a progressive lung disease with approximately 174 million cases worldwide. Electronic questionnaires are increasingly used for collecting patient-reported-outcome (PRO) data about disease symptoms. Our aim was to leverage PRO data, collected to record COPD disease symptoms, in a general modelling framework to enable interpretation of PRO observations in relation to disease progression and potential to predict exacerbations. The data were collected daily over a year, in a prospective, observational study. The e-questionnaire, the EXAcerbations of COPD Tool (EXACT®) included 14 items (i.e. questions) with 4 or 5 ordered categorical response options. An item response theory (IRT) model was used to relate the responses from each item to the underlying latent variable (which we refer to as disease severity), and on each item level, Markov models (MM) with 4 or 5 categories were applied to describe the dependence between consecutive observations. Minimal continuous time MMs were used and parameterised using ordinary differential equations. One hundred twenty-seven COPD patients were included (median age 67 years, 54% male, 39% current smokers), providing approximately 40,000 observations per EXACT® item. The final model suggested that, with time, patients more often reported the same scores as the previous day, i.e. the scores were more stable. The modelled COPD disease severity change over time varied markedly between subjects, but was small in the typical individual. This is the first IRT model with Markovian properties; our analysis proved them necessary for predicting symptom-defined exacerbations.

A Population Pharmacokinetic-Pharmacodynamic Model of Pegfilgrastim.

Neutropenia and febrile neutropenia (FN) are serious side effects of cytotoxic chemotherapy which may be alleviated with the administration of recombinant granulocyte colony-stimulating factor (GCSF) derivatives, such as pegfilgrastim (PG) which increases absolute neutrophil count (ANC). In this work, a population pharmacokinetic-pharmacodynamic (PKPD) model was developed based on data obtained from healthy volunteers receiving multiple administrations of PG. The developed model was a bidirectional PKPD model, where PG stimulated the proliferation, maturation, and margination of neutrophils and where circulating neutrophils in turn increased the elimination of PG. Simulations from the developed model show disproportionate changes in response with changes in dose. A dose increase of 10% from the 6 mg therapeutic dose taken as a reference leads to area under the curve (AUC) increases of ~50 and ~5% for PK and PD, respectively. A full random effects covariate model showed that little of the parameter variability could be explained by sex, age, body size, and race. As a consequence, little of the secondary parameter variability (Cmax and AUC of PG and ANC) could be explained by these covariates.

Dosing and Switching Strategies for Paliperidone Palmitate 3-Month Formulation in Patients with Schizophrenia Based on Population Pharmacokinetic Modeling and Simulation, and Clinical Trial Data.

Paliperidone palmitate 3-month formulation (PP3M), a long-acting injectable atypical antipsychotic, was recently approved in the US and Europe for the treatment of schizophrenia in adult patients who have already been treated with paliperidone palmitate 1-month formulation (PP1M) for ≥4 months. This article reviews the pharmacokinetic rationale for the approved dosing regimens for PP3M, dosing windows, management of missed doses and treatment discontinuation, switching to other formulations, and dosing in special populations. Approved PP3M dosing regimens are based on the comparisons of simulations with predefined dosing regimens using paliperidone palmitate and oral paliperidone extended release (ER) population pharmacokinetic models (one-compartment model with two saturable absorption processes for PP3M; one-compartment model with parallel zero- and first-order absorption for PP1M; two-compartment model with sequential zero- and first-order absorption for ER) versus clinical trial data. Covariates were obtained by resampling subject covariates from the pharmacokinetics database for PP1M and PP3M. Simulation scenarios with varying doses and covariate values were generated. The population median and 90% prediction interval of the simulated concentration-time profiles were plotted for simulation outcomes evaluation. Simulations described in this paper provide (a) simulated plasma exposures for switching from PP1M to PP3M, (b) support for a once-every-3-months injection cycle, (c) information on dosing windows and managing missed doses of PP3M, (d) important guidance on PP3M dosing in special patient populations, and (e) key PP3M pharmacokinetic exposure metrics based on the population pharmacokinetic PP3M model. Population pharmacokinetics provided practical guidance to establish dosing regimens for PP3M.

Modeling a Composite Score in Parkinson's Disease Using Item Response Theory.

In the current work, we present the methodology for development of an Item Response Theory model within a non-linear mixed effects framework to characterize the longitudinal changes of the Movement Disorder Society (sponsored revision) of Unified Parkinson's Disease Rating Scale (MDS-UPDRS) endpoint in Parkinson's disease (PD). The data were obtained from Parkinson's Progression Markers Initiative database and included 163,070 observations up to 48 months from 430 subjects belonging to De Novo PD cohort. The probability of obtaining a score, reported for each of the items in the questionnaire, was modeled as a function of the subject's disability. Initially, a single latent variable model was explored to characterize the disease progression over time. However, based on the understanding of the questionnaire set-up and the results of a residuals-based diagnostic tool, a three latent variable model with a mixture implementation was able to adequately describe longitudinal changes not only at the total score level but also at each individual item level. The linear progression rates obtained for the patient-reported items and the non-sided items were similar, each of which roughly take about 50 months for a typical subject to progress linearly from the baseline by one standard deviation. However for the sided items, it was found that the better side deteriorates quicker than the disabled side. This study presents a framework for analyzing MDS-UPDRS data, which can be adapted to more traditional UPDRS data collected in PD clinical trials and result in more efficient designs and analyses of such studies.

Population Pharmacokinetics of a Novel Once-Every 3 Months Intramuscular Formulation of Paliperidone Palmitate in Patients with Schizophrenia.

OBJECTIVES: Our objective was to characterize the population pharmacokinetics of paliperidone after intramuscular administration of its long-acting 3-month formulation palmitate ester at various doses and at different injection sites (deltoid and gluteal muscles). METHODS: This retrospective analysis included pooled data from 651 subjects from one phase I study (single injection of the 3-month formulation) and one phase III study (multiple injections of both 1- and 3-month formulations). A total of 8990 pharmacokinetic samples with valid concentration time points were available for this analysis. Nonlinear mixed-effects modelling of the pooled data was conducted using NONMEM software. Knowledge from a previously developed 1-month formulation model was used as a starting point to build the 3-month formulation model. RESULTS: The final model describing the plasma concentrations after administration of the 3-month formulation was a one-compartment model with first-order elimination and two saturable absorption processes (rapid and slow). The apparent volume of distribution estimated for the 3-month formulation was not the same as for the previously modelled 1-month formulation. Apparent clearance (CL), apparent volume of distribution (V), and fraction of the absorbed dose (F3) were estimated to be 3.84 l/h, 1960 L, and 20.9 %. For slow absorption, the maximum absorption rate constant (k a1 max), amount of paliperidone at the absorption site when half of the maximum absorption rate was achieved (k amt1 50), and Hill factor (γ) were estimated to be 90.4 µg/h, 120 mg, and 1.44, respectively. For rapid absorption, the maximum absorption rate constant (k a3 max) and amount of paliperidone at the absorption site when half of the maximum absorption rate was achieved (k amt3 50) were estimated to be 164 µg/h and 21.4 mg, respectively. CONCLUSION: The final model with two saturable absorption processes provided a good description of the pharmacokinetic characteristics of paliperidone after intramuscular administration of its long-acting 3-month formulation palmitate ester. In addition to the structural covariates (creatinine clearance on CL, body mass index on V, and injection volume on both absorption rates), injection site and sex were identified as covariates on k a max of the slow absorption process (k a1 max). Clinical trial registration numbers: NCT01559272, NCT01529515, and NCT01515423.

Approaches for modeling within subject variability in pharmacometric count data analysis: dynamic inter-occasion variability and stochastic differential equations.

Parameter variation in pharmacometric analysis studies can be characterized as within subject parameter variability (WSV) in pharmacometric models. WSV has previously been successfully modeled using inter-occasion variability (IOV), but also stochastic differential equations (SDEs). In this study, two approaches, dynamic inter-occasion variability (dIOV) and adapted stochastic differential equations, were proposed to investigate WSV in pharmacometric count data analysis. These approaches were applied to published count models for seizure counts and Likert pain scores. Both approaches improved the model fits significantly. In addition, stochastic simulation and estimation were used to explore further the capability of the two approaches to diagnose and improve models where existing WSV is not recognized. The results of simulations confirmed the gain in introducing WSV as dIOV and SDEs when parameters vary randomly over time. Further, the approaches were also informative as diagnostics of model misspecification, when parameters changed systematically over time but this was not recognized in the structural model. The proposed approaches in this study offer strategies to characterize WSV and are not restricted to count data.

Evaluation of bias, precision, robustness and runtime for estimation methods in NONMEM 7.

NONMEM is the most widely used software for population pharmacokinetic (PK)-pharmacodynamic (PD) analyses. The latest version, NONMEM 7 (NM7), includes several sampling-based estimation methods in addition to the classical methods. In this study, performance of the estimation methods available in NM7 was investigated with respect to bias, precision, robustness and runtime for a diverse set of PD models. Simulations of 500 data sets from each PD model were reanalyzed with the available estimation methods to investigate bias and precision. Simulations of 100 data sets were used to investigate robustness by comparing final estimates obtained after estimations starting from the true parameter values and initial estimates randomly generated using the CHAIN feature in NM7. Average estimation time for each algorithm and each model was calculated from the runtimes reported by NM7. The method giving the lowest bias and highest precision across models was importance sampling, closely followed by FOCE/LAPLACE and stochastic approximation expectation-maximization. The methods relative robustness differed between models and no method showed clear superior performance. FOCE/LAPLACE was the method with the shortest runtime for all models, followed by iterative two-stage. The Bayesian Markov Chain Monte Carlo method, used in this study for point estimation, performed worst in all tested metrics.

Improved utilization of ADAS-cog assessment data through item response theory based pharmacometric modeling.

PURPOSE: This work investigates improved utilization of ADAS-cog data (the primary outcome in Alzheimer's disease (AD) trials of mild and moderate AD) by combining pharmacometric modeling and item response theory (IRT). METHODS: A baseline IRT model characterizing the ADAS-cog was built based on data from 2,744 individuals. Pharmacometric methods were used to extend the baseline IRT model to describe longitudinal ADAS-cog scores from an 18-month clinical study with 322 patients. Sensitivity of the ADAS-cog items in different patient populations as well as the power to detect a drug effect in relation to total score based methods were assessed with the IRT based model. RESULTS: IRT analysis was able to describe both total and item level baseline ADAS-cog data. Longitudinal data were also well described. Differences in the information content of the item level components could be quantitatively characterized and ranked for mild cognitively impairment and mild AD populations. Based on clinical trial simulations with a theoretical drug effect, the IRT method demonstrated a significantly higher power to detect drug effect compared to the traditional method of analysis. CONCLUSION: A combined framework of IRT and pharmacometric modeling permits a more effective and precise analysis than total score based methods and therefore increases the value of ADAS-cog data.

Performance comparison of various maximum likelihood nonlinear mixed-effects estimation methods for dose-response models.

Estimation methods for nonlinear mixed-effects modelling have considerably improved over the last decades. Nowadays, several algorithms implemented in different software are used. The present study aimed at comparing their performance for dose-response models. Eight scenarios were considered using a sigmoid E(max) model, with varying sigmoidicity and residual error models. One hundred simulated datasets for each scenario were generated. One hundred individuals with observations at four doses constituted the rich design and at two doses, the sparse design. Nine parametric approaches for maximum likelihood estimation were studied: first-order conditional estimation (FOCE) in NONMEM and R, LAPLACE in NONMEM and SAS, adaptive Gaussian quadrature (AGQ) in SAS, and stochastic approximation expectation maximization (SAEM) in NONMEM and MONOLIX (both SAEM approaches with default and modified settings). All approaches started first from initial estimates set to the true values and second, using altered values. Results were examined through relative root mean squared error (RRMSE) of the estimates. With true initial conditions, full completion rate was obtained with all approaches except FOCE in R. Runtimes were shortest with FOCE and LAPLACE and longest with AGQ. Under the rich design, all approaches performed well except FOCE in R. When starting from altered initial conditions, AGQ, and then FOCE in NONMEM, LAPLACE in SAS, and SAEM in NONMEM and MONOLIX with tuned settings, consistently displayed lower RRMSE than the other approaches. For standard dose-response models analyzed through mixed-effects models, differences were identified in the performance of estimation methods available in current software, giving material to modellers to identify suitable approaches based on an accuracy-versus-runtime trade-off.

Modeling longitudinal daily seizure frequency data from pregabalin add-on treatment.

The purpose of this study was to describe longitudinal daily seizure count data with respect to the effects of time and pregabalin add-on therapy. Models were developed in a stepwise manner: base model, time effect model, and time and drug effect (final) model, using a negative binomial distribution with Markovian features. Mean daily seizure count (λ) was estimated to be 0.385 (relative standard error [RSE] 3.09%) and was further increased depending on the seizure count on the previous day. An overdispersion parameter (OVDP), representing extra-Poisson variation, was estimated to be 0.330 (RSE 11.7%). Interindividual variances on λ and OVDP were 84.7% and 210%, respectively. Over time, λ tended to increase exponentially with a rate constant of 0.272 year⁻¹ (RSE 26.8%). A mixture model was applied to classify responders/nonresponders to pregabalin treatment. Within the responders, λ decreased exponentially with respect to dose with a constant of 0.00108 mg⁻¹ (RSE 11.9%). The estimated responder rate was 66% (RSE 27.6%). Simulation-based diagnostics showed the model reasonably reproduced the characteristics of observed data. Highly variable daily seizure frequency was successfully characterized incorporating baseline characteristics, time effect, and the effect of pregabalin with classification of responders/nonresponders, all of which are necessary to adequately assess the efficacy of antiepileptic drugs.