Model-based analyses of bioequivalence crossover trials using the stochastic approximation expectation maximisation algorithm

Development and comparison of model-integrated evidence approaches for bioequivalence studies with pharmacokinetic end points

By applying nonlinear mixed-effect (NLME) models, model-integrated evidence (MIE) approaches are able to analyze bioequivalence (BE) data with pharmacokinetic end points that have sparse sampling, which is problematic for non-compartmental analysis (NCA). However, MIE approaches may suffer from inflation of type I error due to underestimation of parameter uncertainty and to the assumption of asymptotic normality. In this study, we developed a MIE BE analysis method that is based on a pre-defined model and consists of several steps including model fitting, uncertainty assessment, simulation, and BE determination. The presented MIE approach has several improvements compared with the previously reported model-integrated methods: (1) treatment, sequence, and period effects are only added to absorption parameters (such as relative bioavailability and rate of absorption) instead of all PK parameters; (2) a simulation step is performed to generate confidence intervals of the pharmacokinetic metrics for BE assessment; and (3) in an effort to maintain type I error, two more advanced parameter uncertainty evaluation approaches are explored, a nonparametric (case resampling) bootstrap, and sampling importance resampling (SIR). To evaluate the developed method and compare the uncertainty assessment methods, simulation experiments were performed for BE studies using a two-way crossover design with different amounts of information (sparse to rich designs) and levels of variability. Based on the simulation results, the method using SIR for parameter uncertainty quantification controls type I error at the nominal level of 0.05 (i.e., the significance level set for BE evaluation) even for studies with small sample size and/or sparse sampling. As expected, our MIE approach for BE assessment exhibited higher power than the NCA-based method, especially as the data becomes sparser and/or more variable.

Model-based bioequivalence approach for sparse pharmacokinetic bioequivalence studies: Model selection or model averaging?

Conventional pharmacokinetic (PK) bioequivalence (BE) studies aim to compare the rate and extent of drug absorption from a test (T) and reference (R) product using non-compartmental analysis (NCA) and the two one-sided test (TOST). Recently published regulatory guidance recommends alternative model-based (MB) approaches for BE assessment when NCA is challenging, as for long-acting injectables and products which require sparse PK sampling. However, our previous research on MB-TOST approaches showed that model misspecification can lead to inflated type I error. The objective of this research was to compare the performance of model selection (MS) on R product arm data and model averaging (MA) from a pool of candidate structural PK models in MBBE studies with sparse sampling. Our simulation study was inspired by a real case BE study using a two-way crossover design. PK data were simulated using three structural models under the null hypothesis and one model under the alternative hypothesis. MB-TOST was applied either using each of the five candidate models or following MS and MA with or without the simulated model in the pool. Assuming T and R have the same PK model, our simulation shows that following MS and MA, MB-TOST controls type I error rates at or below 0.05 and attains similar or even higher power than when using the simulated model. Thus, we propose to use MS prior to MB-TOST for BE studies with sparse PK sampling and to consider MA when candidate models have similar Akaike information criterion.

Uncertainty Computation at Finite Distance in Nonlinear Mixed Effects Models-a New Method Based on Metropolis-Hastings Algorithm

The standard errors (SE) of the maximum likelihood estimates (MLE) of the population parameter vector in nonlinear mixed effect models (NLMEM) are usually estimated using the inverse of the Fisher information matrix (FIM). However, at a finite distance, i.e. far from the asymptotic, the FIM can underestimate the SE of NLMEM parameters. Alternatively, the standard deviation of the posterior distribution, obtained in Stan via the Hamiltonian Monte Carlo algorithm, has been shown to be a proxy for the SE, since, under some regularity conditions on the prior, the limiting distributions of the MLE and of the maximum a posterior estimator in a Bayesian framework are equivalent. In this work, we develop a similar method using the Metropolis-Hastings (MH) algorithm in parallel to the stochastic approximation expectation maximisation (SAEM) algorithm, implemented in the saemix R package. We assess this method on different simulation scenarios and data from a real case study, comparing it to other SE computation methods. The simulation study shows that our method improves the results obtained with frequentist methods at finite distance. However, it performed poorly in a scenario with the high variability and correlations observed in the real case study, stressing the need for calibration.

Evaluation of model-based bioequivalence approach for single sample pharmacokinetic studies

In a traditional pharmacokinetic (PK) bioequivalence (BE) study, a two-way crossover study is conducted, PK parameters (namely the area under the time-concentration curve [AUC] and the maximal concentration [C max ]) are obtained by noncompartmental analysis (NCA), and the BE analysis is performed using the two one-sided test (TOST) method. For ophthalmic drugs, however, only one sample of aqueous humor, in one eye, per eye can be obtained in each patient, which precludes the traditional BE analysis. To circumvent this issue, the U.S. Food and Drug Administration (FDA) has proposed an approach coupling NCA with either parametric or nonparametric bootstrap (NCA bootstrap). The model-based TOST (MB-TOST) has previously been proposed and evaluated successfully for various settings of sparse PK BE studies. In this paper, we evaluate, via simulations, MB-TOST in the specific setting of single sample PK BE study and compare its performance to NCA bootstrap. We performed BE study simulations using a published PK model and parameter values and evaluated multiple scenarios, including study design (parallel or crossover), sampling times (5 or 10 spread across the dosing interval), and geometric mean ratio (of 0.8, 0.9, 1, and 1.25). Using the simulated structural PK model, MB-TOST performed similarly to NCA bootstrap for AUC. ForC max , the latter tended to be conservative and less powerful. Our research suggests that MB-TOST may be considered as an alternative BE approach for single sample PK studies, provided that the PK model is correctly specified and the test drug has the same structural model as the reference drug.

Semiparametric generalized estimating equations for repeated measurements in cross-over designs

A model for cross-over designs with repeated measures within each period was developed. It was obtained using an extension of generalized estimating equations that includes a parametric component to model treatment effects and a non-parametric component to model time and carry-over effects; the estimation approach for the non-parametric component is based on splines. A simulation study was carried out to explore the model properties. Thus, when there is a carry-over effect or a functional temporal effect, the proposed model presents better results than the standard models. Among the theoretical properties, the solution is found to be analogous to weighted least squares. Therefore, model diagnostics can be made by adapting the results from a multiple regression. The proposed methodology was implemented in the data sets of the cross-over experiments that motivated the approach of this work: systolic blood pressure and insulin in rabbits.

Relative bioavailability of delamanid 50 mg tablets dispersed in water in healthy adult volunteers

AIM

Delamanid is a novel drug for the treatment of drug-resistant tuberculosis, manufactured as 50-mg solid and 25-mg dispersible tablets. We evaluated the effects of dispersing the 50-mg tablet, focusing on the relative bioavailability.

METHODS

Delamanid, 50-mg tablets administered dispersed vs swallowed whole, was investigated in a phase I, four-period, crossover study. Two of three dose strengths of delamanid (25, 50 or 100 mg) were given to healthy adult participants, in both whole and dispersed forms, with a 7-day washout period. Blood samples were collected over 168 h after each dose. Delamanid and its metabolite DM-6705 were analysed with a validated liquid chromatography tandem mass spectrometry assay. The pharmacokinetics of both analytes were analysed using nonlinear mixed-effect modelling. Palatability and acceptability were determined using a standardized questionnaire.

RESULTS

Twenty-four participants completed the study. The bioavailability of dispersed tablets was estimated to be 107% of whole tablets, with a 90% confidence interval of 99.7-114%, fulfilling bioequivalence criteria. The two formulations were not significantly different regarding either bioavailability or its variability. Bioavailability increased at lower doses, by 34% (26-42%) at 50 mg and by 74% (64-86%) at 25 mg, relative to 100 mg. The majority of participants (93%) found the dispersed formulation acceptable in palatability across all delamanid doses.

CONCLUSIONS

Dispersed 50-mg delamanid tablets have similar bioavailability to tablets swallowed whole in adult volunteers. This can be an option for children and other patients who cannot swallow whole tablets, improving access to treatment.

Impact of model misspecification on model-based tests in PK studies with parallel design: real case and simulation studies

This article evaluates the performance of pharmacokinetic (PK) equivalence testing between two formulations of a drug through the Two-One Sided Tests (TOST) by a model-based approach (MB-TOST), as an alternative to the classical non-compartmental approach (NCA-TOST), for a sparse design with a few time points per subject. We focused on the impact of model misspecification and the relevance of model selection for the reference data. We first analysed PK data from phase I studies of gantenerumab, a monoclonal antibody for the treatment of Alzheimer’s disease. Using the original rich sample data, we compared MB-TOST to NCA-TOST for validation. Then, the analysis was repeated on a sparse subset of the original data with MB-TOST. This analysis inspired a simulation study with rich and sparse designs. With rich designs, we compared NCA-TOST and MB-TOST in terms of type I error and study power. With both designs, we explored the impact of misspecifying the model on the performance of MB-TOST and adding a model selection step. Using the observed data, the results of both approaches were in general concordance. MB-TOST results were robust with sparse designs when the underlying PK structural model was correctly specified. Using the simulated data with a rich design, the type I error of NCA-TOST was close to the nominal level. When using the simulated model, the type I error of MB-TOST was controlled on rich and sparse designs, but using a misspecified model led to inflated type I errors. Adding a model selection step on the reference data reduced the inflation. MB-TOST appears as a robust alternative to NCA-TOST, provided that the PK model is correctly specified and the test drug has the same PK structural model as the reference drug.

On the curved exponential family in the Stochastic Approximation Expectation Maximization Algorithm

The Expectation-Maximization Algorithm (EM) is a widely used method allowing to estimate the maximum likelihood of models involving latent variables. When the Expectation step cannot be computed easily, one can use stochastic versions of the EM such as the Stochastic Approximation EM. This algorithm, however, has the drawback to require the joint likelihood to belong to the curved exponential family. To overcome this problem, [16] introduced a rewriting of the model which “exponentializes” it by considering the parameter as an additional latent variable following a Normal distribution centered on the newly defined parameters and with fixed variance. The likelihood of this new exponentialized model now belongs to the curved exponential family. Although often used, there is no guarantee that the estimated mean is close to the maximum likelihood estimate of the initial model. In this paper, we quantify the error done in this estimation while considering the exponentialized model instead of the initial one. By verifying those results on an example, we see that a trade-off must be made between the speed of convergence and the tolerated error. Finally, we propose a new algorithm allowing a better estimation of the parameter in a reasonable computation time to reduce the bias.

New Model-Based Bioequivalence Statistical Approaches for Pharmacokinetic Studies with Sparse Sampling

In traditional pharmacokinetic (PK) bioequivalence analysis, two one-sided tests (TOST) are conducted on the area under the concentration-time curve and the maximal concentration derived using a non-compartmental approach. When rich sampling is unfeasible, a model-based (MB) approach, using nonlinear mixed effect models (NLMEM) is possible. However, MB-TOST using asymptotic standard errors (SE) presents increased type I error when asymptotic conditions do not hold. In this work, we propose three alternative calculations of the SE based on (i) an adaptation to NLMEM of the correction proposed by Gallant, (ii) the a posteriori distribution of the treatment coefficient using the Hamiltonian Monte Carlo algorithm, and (iii) parametric random effects and residual errors bootstrap. We evaluate these approaches by simulations, for two-arms parallel and two-period, two-sequence cross-over design with rich (n = 10) and sparse (n = 3) sampling under the null and the alternative hypotheses, with MB-TOST. All new approaches correct for the inflation of MB-TOST type I error in PK studies with sparse designs. The approach based on the a posteriori distribution appears to be the best compromise between controlled type I errors and computing times. MB-TOST using non-asymptotic SE controls type I error rate better than when using asymptotic SE estimates for bioequivalence on PK studies with sparse sampling.

Efficient model-based bioequivalence testing

The classical approach to analyze pharmacokinetic (PK) data in bioequivalence studies aiming to compare two different formulations is to perform noncompartmental analysis (NCA) followed by two one-sided tests (TOST). In this regard, the PK parameters area under the curve (AUC) and $C_{\max}$ are obtained for both treatment groups and their geometric mean ratios are considered. According to current guidelines by the U.S. Food and Drug Administration and the European Medicines Agency, the formulations are declared to be sufficiently similar if the $90\%$ confidence interval for these ratios falls between $0.8$ and $1.25 $. As NCA is not a reliable approach in case of sparse designs, a model-based alternative has already been proposed for the estimation of $\rm AUC$ and $C_{\max}$ using nonlinear mixed effects models. Here we propose another, more powerful test than the TOST and demonstrate its superiority through a simulation study both for NCA and model-based approaches. For products with high variability on PK parameters, this method appears to have closer type I errors to the conventionally accepted significance level of $0.05$, suggesting its potential use in situations where conventional bioequivalence analysis is not applicable.

Robust designs accounting for model uncertainty in longitudinal studies with binary outcomes

To optimize designs for longitudinal studies analyzed by mixed-effect models with binary outcomes, the Fisher information matrix can be used. Optimal design approaches, however, require a priori knowledge of the model. We aim to propose, for the first time, a robust design approach accounting for model uncertainty in longitudinal trials with two treatment groups, assuming mixed-effect logistic models. To optimize designs given one model, we compute several optimality criteria based on Fisher information matrix evaluated by the new approach based on Monte-Carlo/Hamiltonian Monte-Carlo. We propose to use the DDS-optimality criterion, as it ensures a compromise between the precision of estimation of the parameters, and hence the Wald test power, and the overall precision of parameter estimation. To account for model uncertainty, we assume candidate models with their respective weights. We compute robust design across these models using compound DDS-optimality. Using the Fisher information matrix, we propose to predict the average power over these models. Evaluating this approach by clinical trial simulations, we show that the robust design is efficient across all models, allowing one to achieve good power of test. The proposed design strategy is a new and relevant approach to design longitudinal studies with binary outcomes, accounting for model uncertainty.

Design and inference for 3-stage bioequivalence testing with serial sampling data

A bioequivalence test is to compare bioavailability parameters, such as the maximum observed concentration (Cmax ) or the area under the concentration-time curve, for a test drug and a reference drug. During the planning of a bioequivalence test, it requires an assumption about the variance of Cmax or area under the concentration-time curve for the estimation of sample size. Since the variance is unknown, current 2-stage designs use variance estimated from stage 1 data to determine the sample size for stage 2. However, the estimation of variance with the stage 1 data is unstable and may result in too large or too small sample size for stage 2. This problem is magnified in bioequivalence tests with a serial sampling schedule, by which only one sample is collected from each individual and thus the correct assumption of variance becomes even more difficult. To solve this problem, we propose 3-stage designs. Our designs increase sample sizes over stages gradually, so that extremely large sample sizes will not happen. With one more stage of data, the power is increased. Moreover, the variance estimated using data from both stages 1 and 2 is more stable than that using data from stage 1 only in a 2-stage design. These features of the proposed designs are demonstrated by simulations. Testing significance levels are adjusted to control the overall type I errors at the same level for all the multistage designs.

Flexible multivariate nonlinear models for bioequivalence problems

Modelling the concentration of a drug in the bloodstream over time is usually done using compartment models. In pharmacokinetic data, they turn into highly nonlinear mixed-effects models (NLMEMs) when we take the heterogeneity between subjects into account. Fitting of NLMEMs can be difficult and may involve complex algorithms, with convergence critically depending on the initial values and maybe requiring data transformations. In this article, we propose a flexible alternative to the usual parametric compartment models, inspired by the Multivariate SuperImposition by Translation and Rotation (MSITAR) model but adapted to be applicable in this new field. A fully parametric one-compartment NLMEM is considered for comparison. We make use of a Bayesian approach and illustrate the method on a real dataset where the interest lies in contrasting the average and individual bioequivalence of a test and reference formulation of an anti-hypertensive drug.

Population Pharmacokinetic Analysis of the Oral Absorption Process and Explaining Intra-Subject Variability in Plasma Exposures of Imatinib in Healthy Volunteers

BACKGROUND AND OBJECTIVE

Imatinib mesylate is presently the first-line treatment for chronic myeloid leukemia (CML). The aim of this study was to investigate the absorption and distribution kinetics of imatinib in healthy Iranian volunteers using nonlinear mixed effects modeling (NLMEM) to assess the overall, intra- and inter-subject variabilities in pharmacokinetic parameters after oral administration.

METHODS

This analysis was based on data from 24 healthy subjects who participated in a bioequivalence study after administering a single dose of 200 mg of each formulation. Imatinib concentrations were quantified using a validated liquid chromatography method. To simultaneously describe the imatinib pharmacokinetic profiles obtained with both formulations, a population pharmacokinetic model was applied to data using SAEM algorithm implemented in MONOLIX, whilst simulations were used by numerical solving of ordinary differential equations to calculate secondary parameters in individuals for bioequivalence studies.

RESULTS

According to goodness-of-fit criteria, a two-compartment open model with sequential zero- then first-order absorption and first-order elimination was used as the structural pharmacokinetic model. Inter-individual variability (IIV) was considered for all parameters. Typical population estimates (% IIV) were fraction of the drug absorbed with a zero-order kinetic (Fr) of 0.153 (47.9 %) in period (Tk0) of 0.714 h (47.4 %), first-order absorption rate constant (k a) of 0.94 h(-1)(31.2 %), oral clearance of 19 L/h (27.9 %), central volume of distribution (V c/F) of 139 L (21.5 %), apparent peripheral volume of distribution (V p/F) of 130 L (29.7 %) and the apparent inter-compartment clearance (Q/F) of 29.6 L/h (41.8 %). Body mass index (BMI) was the only covariate found to significantly affect V p /F. The coefficient of variation for intra-individual plasma exposure (AUC0-∞) was 27.8 %.

CONCLUSIONS

Analyses using NLMEM for imatinib exhibited absorption complexities such as two input rates and medium to high intra-individual variability in drug exposure.

Approximate testing in two-stage nonlinear mixed models

We investigate here small sample properties of approximate F-tests about fixed effects parameters in nonlinear mixed models. For estimation of population fixed effects parameters as well as variance components, we apply the two-stage approach. This method is useful and popular when the number of observations per sampling unit is large enough. The approximate F-test is constructed based on large sample approximation to the distribution of nonlinear least squares estimates of subject-specific parameters. We recommend a modified test statistic that takes into consideration approximation to the large sample Fisher information matrix (See [1]). Our main focus is on comparing finite sample properties of broadly used approximate tests (Wald test and likelihood ratio test) and the modified F-test under the null hypothesis, especially accuracy of p-values (See [2]). For that purpose two extensive simulation studies are conducted based on pharmacokinetic models (See [3, 4]).

Statistical power calculations for mixed pharmacokinetic study designs using a population approach

Simultaneous modelling of dense and sparse pharmacokinetic data is possible with a population approach. To determine the number of individuals required to detect the effect of a covariate, simulation-based power calculation methodologies can be employed. The Monte Carlo Mapped Power method (a simulation-based power calculation methodology using the likelihood ratio test) was extended in the current study to perform sample size calculations for mixed pharmacokinetic studies (i.e. both sparse and dense data collection). A workflow guiding an easy and straightforward pharmacokinetic study design, considering also the cost-effectiveness of alternative study designs, was used in this analysis. Initially, data were simulated for a hypothetical drug and then for the anti-malarial drug, dihydroartemisinin. Two datasets (sampling design A: dense; sampling design B: sparse) were simulated using a pharmacokinetic model that included a binary covariate effect and subsequently re-estimated using (1) the same model and (2) a model not including the covariate effect in NONMEM 7.2. Power calculations were performed for varying numbers of patients with sampling designs A and B. Study designs with statistical power >80% were selected and further evaluated for cost-effectiveness. The simulation studies of the hypothetical drug and the anti-malarial drug dihydroartemisinin demonstrated that the simulation-based power calculation methodology, based on the Monte Carlo Mapped Power method, can be utilised to evaluate and determine the sample size of mixed (part sparsely and part densely sampled) study designs. The developed method can contribute to the design of robust and efficient pharmacokinetic studies.

New approach to assess bioequivalence parameters using generalized gamma mixed-effect model (model-based asymptotic bioequivalence test)

In the pharmacokinetic (PK) study under a 2x2 crossover design that involves both the test and reference drugs, we propose a mixed-effects model for the drug concentration-time profiles obtained from subjects who receive different drugs at different periods. In the proposed model, the drug concentrations repeatedly measured from the same subject at different time points are distributed according to a multivariate generalized gamma distribution, and the drug concentration-time profiles are described by a compartmental PK model with between-subject and within-subject variations. We then suggest a bioequivalence test based on the estimated bioavailability parameters in the proposed mixed-effects model. The results of a Monte Carlo study further show that the proposed model-based bioequivalence test is not only better on maintaining its level but also more powerful for detecting the bioequivalence of the two drugs than the conventional bioequivalence test based on a non-compartmental analysis or the one based on a mixed-effects model with a normal error variable. The application of the proposed model and test is finally illustrated by using data sets in two PK studies.

The use of pharmacometrics to optimize biosimilar development

Pharmacometric approaches can assist in biosimilar development by leveraging quantitative knowledge of the originator product characteristics such as dose-exposure and exposure-response information to support a targeted approach to clinical studies. The degree to which these approaches can be applied relies on the level of information known about the originator and information that supports application of the originator model to the biosimilar. A model-based approach testing the hypothesis that the biosimilar PK and/or PK/PD profile is similar to the originator in the target patient population is aligned with the central comparability exercise required for the biosimilar approval. This Commentary details the key opportunities in study design and study analysis where pharmacometrics approaches can aid biosimilar development.

Clinical trial simulation to evaluate power to compare the antiviral effectiveness of two hepatitis C protease inhibitors using nonlinear mixed effect models: a viral kinetic approach

BACKGROUND

Models of hepatitis C virus (HCV) kinetics are increasingly used to estimate and to compare in vivo drug's antiviral effectiveness of new potent anti-HCV agents. Viral kinetic parameters can be estimated using non-linear mixed effect models (NLMEM). Here we aimed to evaluate the performance of this approach to precisely estimate the parameters and to evaluate the type I errors and the power of the Wald test to compare the antiviral effectiveness between two treatment groups when data are sparse and/or a large proportion of viral load (VL) are below the limit of detection (BLD).

METHODS

We performed a clinical trial simulation assuming two treatment groups with different levels of antiviral effectiveness. We evaluated the precision and the accuracy of parameter estimates obtained on 500 replication of this trial using the stochastic approximation expectation-approximation algorithm which appropriately handles BLD data. Next we evaluated the type I error and the power of the Wald test to assess a difference of antiviral effectiveness between the two groups. Standard error of the parameters and Wald test property were evaluated according to the number of patients, the number of samples per patient and the expected difference in antiviral effectiveness.

RESULTS

NLMEM provided precise and accurate estimates for both the fixed effects and the inter-individual variance parameters even with sparse data and large proportion of BLD data. However Wald test with small number of patients and lack of information due to BLD resulted in an inflation of the type I error as compared to the results obtained when no limit of detection of VL was considered. The corrected power of the test was very high and largely outperformed what can be obtained with empirical comparison of the mean VL decline using Wilcoxon test.

CONCLUSION

This simulation study shows the benefit of viral kinetic models analyzed with NLMEM over empirical approaches used in most clinical studies. When designing a viral kinetic study, our results indicate that the enrollment of a large number of patients is to be preferred to small population sample with frequent assessments of VL.

Pharmacokinetic similarity of biologics: analysis using nonlinear mixed-effects modeling

Our objective was to show, using two examples, that a pharmacokinetic (PK) similarity analysis can be performed using nonlinear mixed-effects models (NLMEM). We used two studies that compared different biosimilars: a three-way crossover trial with somatropin and a parallel-group trial with epoetin-α. For both data sets, the results of NLMEM-based analysis were compared with those of noncompartmental analysis (NCA). For the latter analysis, we performed an NLMEM-based equivalence Wald test on secondary parameters of the model: the area under the curve and the maximal concentration. Somatropin PK was described by a one-compartment model and epoetin-α PK by a two-compartment model with linear and Michaelis-Menten elimination. For both studies, similarity of PK was demonstrated by means of both NCA and NLMEM, and both methods led to similar results. Therefore, for establishing similarity, PK data can be analyzed by either of the methods. NCA is an easier approach because it does not require data modeling; however, NLMEM leads to a better understanding of the underlying biological system.

Design evaluation and optimisation in crossover pharmacokinetic studies analysed by nonlinear mixed effects models

Bioequivalence or interaction trials are commonly studied in crossover design and can be analysed by nonlinear mixed effects models as an alternative to noncompartmental approach. We propose an extension of the population Fisher information matrix in nonlinear mixed effects models to design crossover pharmacokinetic trials, using a linearisation of the model around the random effect expectation, including within-subject variability and discrete covariates fixed or changing between periods. We use the expected standard errors of treatment effect to compute the power for the Wald test of comparison or equivalence and the number of subjects needed for a given power. We perform various simulations mimicking crossover two-period trials to show the relevance of these developments. We then apply these developments to design a crossover pharmacokinetic study of amoxicillin in piglets and implement them in the new version 3.2 of the r function PFIM.