 Research
 Open Access
 Published:
Distributional regression in clinical trials: treatment effects on parameters other than the mean
BMC Medical Research Methodology volume 22, Article number: 56 (2022)
Abstract
Background
The classical linear model is widely used in the analysis of clinical trials with continuous outcomes. However, required model assumptions are frequently not met, resulting in estimates of treatment effect that can be inefficient and biased. In addition, traditional models assess treatment effect only on the mean response, and not on other aspects of the response, such as the variance. Distributional regression modelling overcomes these limitations. The purpose of this paper is to demonstrate its usefulness for the analysis of clinical trials, and superior performance to that of traditional models.
Methods
Distributional regression models are demonstrated, and contrasted with normal linear models, on data from the LIPID randomized controlled trial, which compared the effects of pravastatin with placebo in patients with coronary heart disease. Systolic blood pressure (SBP) and the biomarker midregional proadrenomedullin (MRproADM) were analysed. Treatment effect was estimated in models that used response distributions more appropriate than the normal (BoxCoxt and Johnson’s S_{u} for MRproADM and SBP, respectively), applied censoring below the detection limit of MRproADM, estimated treatment effect on distributional parameters other than the mean, and included random effects for longitudinal observations. A simulation study was conducted to compare the performance of distributional regression models with normal linear regression, under conditions mimicking the LIPID study. The R package gamlss (Generalized Additive Models for Location, Scale and Shape), which implements maximum likelihood estimation for distributional regression modelling, was used throughout.
Results
In all cases the distributional regression models fit the data well, in contrast to poor fits obtained for traditional models; for MRproADM a small but significant treatment effect on the mean was detected by the distributional regression model and not the normal model; and for SBP a beneficial treatment effect on the variance was demonstrated. In the simulation study distributional models strongly outperformed normal models when the response variable was nonnormal and heterogeneous; and there was no disadvantage introduced by the use of distributional regression modelling when the response satisfied the normal linear model assumptions.
Conclusions
Distributional regression models are a rich framework, largely untapped in the clinical trials world. We have demonstrated a sample of the capabilities of these models for the analysis of trials. If interest lies in accurate estimation of treatment effect on the mean, or other distributional features such as variance, the use of distributional regression modelling will yield superior estimates to traditional normal models, and is strongly recommended.
Trial registration
The LIPID trial was retrospectively registered on ANZCTR on 27/04/2016, registration number ACTRN12616000535471.
Background
The classical linear model is widely used in the analysis of clinical trials with continuous outcomes. However when the required model assumptions are not met, estimates of treatment effect can be inefficient and biased. Since Nelder and Wedderburn’s [1] seminal paper introducing generalized linear models (GLMs), there has been a revolution in the development of regression methodology, enabled by an accompanying exponential increase in computing power over the same period. The current stateoftheart “distributional regression” framework, first proposed by [2] as Generalized Additive Models for Location, Scale and Shape (GAMLSS), accommodates any computable parametric response distribution; the facility for modelling any distributional parameter (not just the mean), smooth terms for continuous covariates; and random effects for modelling clustered (e.g. longitudinal) observations. Within this rich environment, the analyst has the ability to model trials outcomes using appropriate assumptions, resulting in estimates of treatment effect which should be efficient and asymptotically unbiased.
The use of distributional regression in clinical trials has broader motivations than efficient and unbiased statistical modelling. The mean of a biomarker may not be the only feature of the biomarker distribution that governs prognosis, and the beneficial effects of treatment may be mediated through other distributional characteristics. An example is blood pressure variability. Studies have shown that blood pressure variability is prognostic for cardiovascular outcomes and clinical trials assessing interventions that decrease blood pressure variability have been recommended [3]. Such trials would assess the effect of treatment on the variance of a biomarker distribution, for which distributional regression is well suited.
Development of regression modelling
In the case of a twoarmed trial with continuous outcome y, the classical linear model specifies
where t _{i} is an indicator variable for treatment, β _{t} is the treatment effect and x _{i} = (x _{i1}, …, x _{ip})^{⊤} is a covariate vector whose elements may be continuous or binary. While taking logarithms or other transformations of y and/or some of the covariates is sometimes helpful for satisfying the normality and linearity assumptions, the restrictive model assumptions of (conditional) normality, homoscedasticity and linearity are in practice frequently not met. Nevertheless this model remains popular for analysis.
Relaxation of the classical linear model assumptions started with the GLM: “Theoretical and applied statistics were both convulsed by the publication of the GLM paper by Nelder & Wedderburn (1972).” [4]. The revolutionary aspects of the GLM were the extension of the choice of response distribution to any member of the exponential family of distributions, which includes the normal, Poisson, binomial, Gamma, inverse Gaussian and Tweedie distributions; inclusion of a link function g(·) in the model specification for μ _{i}:
where g(·) is any monotonic differentiable function; and an algorithm for the computation of the maximum likelihood estimates which was computationally feasible within computing constraints at the time. Since then development of the regression framework has included broadening the linear predictor to include smooth terms (defined below); the extension of allowed response distributions to any parametric distribution which is twice differentiable and computable; and the ability to specify model equations similar to (1) for any distribution parameter, not just the mean. For example, a normal linear model which displays variance which is dependent on the treatment and/or covariate(s) (i.e. heteroscedasticity) is accommodated:
The superscripts on the β's denote the parameter to which the coefficient corresponds. So, for example, \({\beta}_t^{\mu }\) is the treament effect on μ, while \({\beta}_t^{\sigma }\) is the treatment effect on σ ^{2}. The logarithmic link for σ ^{2} guarantees positivity of parameter estimates \({\hat{\sigma}}_i^2\); other link functions may also be used [5]. Note that not all of x _{1}, …, x _{p} and t need be present in the model equations for both μ and σ ^{2}.
GAMLSS modelling is implemented in the R package gamlss [6, 7], in which over one hundred response distributions are currently available. These enable modelling of the following types of outcome variable, without the need for data transformation:

continuous outcomes (positive or negative);

continuous, nonnegative outcomes;

counts (bounded or unbounded);

proportions (continuous, on the interval zero to one);

any of the above with inflated probabilities at zero, and/or one in the case of proportions.
Parameter estimation in gamlss is performed using the method of maximum likelihood, resulting in the coefficient estimates having the usual asymptotic properties of maximum likelihood estimates, viz. normality, unbiasedness, consistency and efficiency. Bayesian estimation for GAMLSS models, which we do not cover in this paper, is available in the R package bamlss [8].
In this paper we focus on the modelling of continuous outcomes.
Model terms
The model equations given above are of the following (linear) form, for generic parameter θ, and which we denote as \({\eta}_i^{\theta }\):
Should the effect of a covariate x _{j} on \({\eta}_i^{\theta }\) be nonlinear, we can, in order of preference:

find a nonlinear transformation of x _{j} that captures the relationship; or

use a “smooth function” of x _{j}, discussed below; or

categorize x _{j} and enter it into the model as a factor.
A nonlinear transformation, if justified by the data (and possibly the underlying science), is arguably the best approach to a nonlinear relationship. However, should an appropriate transformation not be obvious (as is frequently the case), smooth functions are an excellent alternative. They were introduced into regression modelling as spline terms in Generalized Additive Models, by [9]. (A good recent reference is [10].) The linear predictor becomes
where the \({s}_j^{\theta}\left({x}_{ij}\right)\) can be linear terms (\({\beta}_j^{\theta }{x}_{ij}\)) or smooth functions. There are a few alternatives available for smooth functions, most notably splines, lowess and fractional polynomials. Mathematically and computationally, spline terms are fairly easily handled as they are composed of a series of linear terms (or basis functions), so the model retains its linear structure. There is a tradeoff between complexity (curve “too wiggly,” too sensitive to local variations) and simplicity (curve “too smooth,” not sensitive to important variations), which is resolved in model estimation by using penalties for curve complexity.
Random effects
Clinical trials often involve dependence structures induced by design features such as clustered randomization, longitudinal and repeated measurement on individuals, multiple endpoints, crossover of treatments and many others. Distributional regression using the GAMLSS model has the capacity to introduce random effects into the linear predictor for any or all of the parameters being modelled. Such random effects models can capture the effect of these dependence structures on parameters other than the mean, using mixed GAMLSS models that are analogous to standard mixed models for correlated data. The software and computational tools described below have the capacity to incorporate these random effects in a straightforward manner, making the flexible distributional regression framework available in contexts where dependence must be incorporated into the model.
The GAMLSS model
The full GAMLSS model for a twoarmed trial is
where \(\mathcal{D}\left(\cdot \right)\) is a parametric distribution with computable first and second derivatives and having up to four parameters μ, σ, ν and τ; g _{k}(·) is the link function for the kth parameter, which is monotonic and differentiable; and the functions s _{j}(·) represent a variety of different effects: linear or nonlinear effects of continuous covariates, smooth terms (usually implemented mathematically as splines), spatial effects, or random effects.
Note that the response distribution \(\mathcal{D}\) may have one, two, three or four distribution parameters; of these, not all need to be modelled with covariates; and the sets of covariates in the models for distribution parameters may be the same, disjoint or overlapping.
Methods
Software
Estimation of GAMLSS models is implemented using penalized maximum likelihood, or Bayesian estimation, in the R packages gamlss [7] and bamlss [8] respectively. We have used gamlss for estimation throughout.
LIPID trial
The LongTerm Intervention with Pravastatin in Ischaemic Disease (LIPID) trial is a doubleblind, randomized controlled trial of patients with stable coronary heart disease and a broad range of cholesterol levels. The trial compared the effects of pravastatin with those of placebo in 9014 patients over a mean followup of 6.1 years. Both treatment arms received advice on a cholesterollowering diet. The primary outcome of the trial was mortality from coronary heart disease, and they found a relative risk reduction of 24% (95% CI, 12–35, p < 0.001). The trial also found lower overall mortality and lower incidence of all cardiovascular outcomes in patients treated with pravastatin [11].
Description of biomarker MRproADM
The LIPID trial has also published data on the long term effects of treatment with pravastatin [12] and the analysis of eight biomarkers used in coronary heart disease [13]. One of the eight biomarkers measured in the LIPID study is plasma midregional proadrenomedullin (MRproADM), which is a surrogate marker for adrenomedullin release. Adrenomedullin acts as a vasodilator and has important roles in microcirculation and endothelial dysfunction. Lower baseline and reductions in the change of MRproADM have been associated with decreased risk of major clinical events, even after adjustment for other important biomarkers such as Btype natriuretic peptide (BNP).
Results
MRproADM
The detection limit of MRproADM is 0.05 nmol/L, so it is assumed to be leftcensored at this point. Figure 1 shows density plots of MRproADM at 12 months; there appears to be a small beneficial treatment effect.
GAMLSS analysis
GAMLSS has the functionality for the creation of censored versions of any of its distributions; in this case, we found the BoxCox t (BCT) distribution [14], leftcensored at 0.05, to provide the best fit, according to the AIC criterion. (This was found using the gamlss function fitDist, which performs a search of all appropriate distributions (in this case, continuous distributions on the positive real line), and ranks them in order of Akaike information criterion (AIC)). The BCT is a fourparameter, nonnegative continuous distribution, with μ being the median (approximately), σ the (approximate) coefficient of variation, ν a skewness parameter and τ a kurtosis parameter.
We contrast the leftcensored normal and leftcensored BCT regression models. For comparability with the normal model, we initially only fit the treatment and covariate effects to the parameter μ of the BCT (“reduced model”).

1.
Reduced leftcensored normal model:

2.
Reduced leftcensored BCT model:
where y _{i} is the level of MRproADM activity at month 12 and, for both models, the last argument (0.05) signifies the leftcensoring cutoff. The model equation for μ is
where t _{i} is an indicator for treatment and x _{i} is baseline MRproADM activity. The logarithmic link for μ ensures nonnegativity; the assumption of proportionality implicit in (3) is shown to be satisfied in scatterplots of the logarithms of month 12 and baseline biomarker activity, which are strongly linear (Supplementary material).
We use the AIC as criterion for selection of the covariates in the “extended models” for both distributions. Covariates available for analysis were sex and baseline MRproADM activity. Parameter estimates for the four models are shown in Table 1, and a plot of the treatment effect on the parameter μ as estimated in the reduced models, with 95% confidence intervals, in Fig. 2. (Because of the log link, we display \(\exp \left({\hat{\beta}}_{\mathrm{t}}^{\mu}\right)\), with the value 1 being consistent with no treatment effect.) Note that the BCT model concludes a significant beneficial treatment effect on the median, whereas the normal model does not. The extended models both include sex as a predictor for σ, and the BCT extended model also includes a treatment effect for σ and ν, and baseline MRproADM in the models for σ, ν and τ.
Quantilequantile plots of the normalized quantile residuals are given in Fig. 3; the BCT models provide good fits, whereas the normal model fits are poor. A comparison of the four models using the AIC is given in Table 2. Both BCT models are preferred to the normal models, and of these the extended BCT model is preferred.
Figure 4 shows fitted probability density functions (PDFs) of the reduced normal and extended BCT models, with baseline MRproADM activity at its median and male gender, superimposed on the histogram of the response data in the neighbourhood of the median of baseline MRproADM and male gender. It can be seen that the BCT model successfully captures the sharp peak in the observed data, whereas the normal model is unable to model this feature and compensates with greater variance than is observed.
Systolic blood pressure
SBP is the pressure in arteries during the contraction of the heart muscle, and is a major risk factor for cardiovascular disease in people over 50 years old. SBP generally increases with age, as blood vessels stiffen and plaque builds up over time. SBP was measured 12monthly from baseline in the LIPID trial. We initially analyse change from baseline at 72 months, and subsequently change from baseline at 48 and 72 months, as a longitudinal model. Density plots of SBP at 72 months and SBP change from baseline at 72 months, by treatment group, are shown in Fig. 5.
Models for change in SBP at 72 months
The Johnson’s S _{u} (JSU) distribution [14] was found to provide the best fit to change in SBP at 72 months. The JSU is a fourparameter continuous distribution, with parameters μ (the mean), σ (the standard deviation), ν (skewness parameter) and τ (kurtosis parameter). The JSU looks somewhat similar to the normal distribution when ν = 0, but is capable of skewness and more kurtosis (‘more peaked’) than the normal. Plots of the JSU pdf are shown in the Supplementary Material.
Model selection using the AIC, with the JSU as response distribution, resulted in the extended JSU model given in Table 3. Reduced JSU and normal models, and the extended normal model, are also shown.
In the extended JSU model, the treatment effect on the mean μ is nonsignificant (p = 0.952). However treatment is significant for the standard deviation σ (p = 0.030), in the direction of reduced variability for the active treatment. The extended normal model also finds a significant treatment effect on σ, and no significant treatment effect on μ. Parameter estimates for μ and σ are similar for the extended JSU and normal models; however, for the reduced models they are quite different.
Normal plots of the normalized quantile residuals of the four models are shown in Fig. 6. As for the previous example, the normal models display a lack of fit, particularly in the upper tail, whereas the JSU models fit particularly well.
Longitudinal analysis of change in SBP at months 48 and 72
It may be of interest to model SBP trajectories over the followup period. In this case we are dealing with longitudinal observations, in which each participant has repeated observations which are assumed to be correlated. Modelling of correlated outcomes is achieved using one of two approaches: random effects, or generalized estimating equations (GEE). In the distributional regression framework, random effects methodology is well developed and we illustrate their application in this example. (Although in principle GEE models could also be incorporated in the distributional regression framework, this methodology does not appear to have been developed.)
To facilitate comparison with the extended JSU model for change in SBP at 72 months, we have implemented a random effects model for change in SBP at 48 and 72 months, using the same set of covariates as the former model as well as an effect for month (72 vs 48) for μ and σ. A random intercept for subject, for the parameter μ, has been included in order to model withinsubject correlation. Estimates are given in Table 4: as for the model for month 72, there is a significant treatment effect for σ (p = 0.001) but not for μ (p = 0.381). The estimated multiplicative effect of treatment on σ is exp(− 0.048) = 0.953, i.e. a decrease in standard deviation of 4.7%. There is also a significant effect of month (72 vs 48) on both μ (p = 0.006) and σ (p = 0.012), in the direction of increased μ and σ at month 72 compared with month 48.
The normal plot of the normalized quantile residuals of the longitudinal model (Fig. 7) shows a wellfitting model.
Simulation study
We have demonstrated the usefulness of distributional regression in the analysis of clinical trials. In order to provide evidence that it delivers reliable estimates of treatment effect that are superior to those based on the normal linear model, we have conducted a simulation study, based on the MRproADM biomarker in the LIPID study. The simulation used the treatment allocations and baseline MRproADM of the n = 6, 539 participants who had no missing values for MRproADM at month 12. Baseline MRproADM was included in all models for μ, and in models for σ in the case of extended models.
In the first simulation, we used the extended normal model for 12month MRproADM (see Table 1) as generating model. Using fitted means \({\hat{\mu}}_i\) and variances \({\hat{\sigma}}_i^2\), we simulated 500 samples with normally distributed responses from \(\mathcal{N}\left({\hat{\mu}}_i,{\hat{\sigma}}_i^2\right)\), for i = 1, …, n. For each sample (j = 1, …, 500), we

chose the best distribution using the AIC as selection criterion. We chose from the following distributions: normal, JSU, BCT and skew normal;

using the chosen distribution and a treatment effect specified on the mean only in the case of the normal (i.e. the reduced normal model), and μ and σ for distributions other than the normal, estimated the treatment effect on μ and its standard error.
In the second simulation, we used the extended BCT model for 12month MRproADM (see Table 1) as generating model. Using fitted values \({\hat{\mu}}_i\), \({\hat{\sigma}}_i\), \(\hat{\nu}\) and \(\hat{\tau}\), we simulated 500 samples with responses from \(\mathrm{BCT}\left({\hat{\mu}}_i,{\hat{\sigma}}_i,\hat{\nu},\hat{\tau}\right)\), for i = 1, …, n. For each sample (j = 1, …, 500), we

chose the best distribution as above,

using the chosen distribution and a treatment effect specified on μ and σ, estimated the treatment effect on μ and its standard error, and

compared the estimates with those obtained from a reduced normal model.
Table 5 shows the results of the response distribution selection for both simulations. In the case of the normal generating model, the normal was chosen in 86% of samples, and the JSU in the remaining 14%. For the BCT generating model, the BCT was chosen in 100% of samples.
Table 6 shows the results of the first simulation, in which the generating model is the normal and the chosen distributions were the normal and JSU. As the parameter μ is the mean for both the normal and JSU distributions, parameter estimates for the model for μ are comparable. For the samples for which the JSU distribution was chosen, all of the 95% confidence intervals for the treatment effect on the mean included the true (simulated) effect. For the samples for which the normal distribution was chosen, 96% of the 95% confidence intervals cover the true (simulated) effect. The average standard error of the treatment effect on the mean was, in both cases, 0.004.
These results are reassuring that, when the underlying model is normal, distributional regression does not introduce any disadvantage. In the 14% of cases in which the normal was not the chosen distribution, the JSU was chosen. This distribution is similar to the normal, with more flexibility provided by its four parameters. The coverage rate for \({\beta}_t^{\mu }\) under JSU estimation was better than the normal; and standard errors provided by the two estimating models were the same (to three decimal places).
Results of the second simulation, in which the generating model is the BCT, are shown in Table 7. In all cases the BCT was the chosen distribution; we compare this with estimation under the normal linear (reduced) model. The parameter μ is the mean in case of the normal estimating model, and the (approximate) median in the case of the BCT estimating model. While these are both location parameters, they are not the same and as a result their parameter estimates are not directly comparable. Instead we compare the rejection rates of the hypothesis of no treatment effect on μ:
for the normal and BCT estimating models, using the Wald test. Standard errors for coefficient estimates, based on the inverse of the observed information matrix, are available from gamlss output and were used to generate the Wald statistics. For the BCT estimating model we also tested the hypothesis of no overall treatment effect on μ and σ:
using the likelihood ratio test statistic with null distribution the chisquare with two degrees of freedom. Results are shown in Table 7.
While estimation under the extended BCT model produces acceptably high rejection rates, the rejection rate for \({\beta}_t^{\mu }\) of the reduced normal model is poor. This confirms the result observed in Table 1 and Fig. 2, in which treatment effect on the parameter μ of the BCT model was significant, but the corresponding estimate under the normal model failed to reach significance.
Discussion and conclusions
We have demonstrated some of the richness of the distributional regression modelling framework, and its capability of accurately modelling continuous outcomes with features different from the normal distribution.^{Footnote 1} Since its first publication, GAMLSS has gained increasing popularity in a wide variety of fields of application [15]. Distributional regression offers not only a very wide choice of response distribution, but also the ability to model parameters other than the mean; random effects to accommodate clustered observations; smooth terms for modelling nonlinear effects; and zero inflation and censoring. Particularly noteworthy in the context of clinical trials, is the ability to model a treatment effect on a feature of the response distribution other than the mean, for example the variance. Traditionally the analysis of trials is focused on estimation of the treatment effect, implicitly assumed to be on the mean. In fact a treatment effect on the variance of the response may be of importance, for example for the outcome blood pressure, and conventional normal modelling does not have the capability to detect this. This points to a need for a rethink of the interpretation of the notion of “treatment effect” when treatment is included in distribution parameters other than the mean (or median).
In our two examples, under the normality assumption the lack of fit was not improved by the addition of a model for σ; however the use of more appropriate response distributions resulted in wellfitting models. Where the normal linear model is used in situations in which model assumptions are not met, estimates of treatment effect based on the misspecified model will be biased [16] show that, under fairly general conditions, a large class of regression models (Gaussian, binomial, Poisson) yield asymptotically correct Type I errors for hypotheses on treatment effect, even when the models are incorrectly specified. This robustness to violation of model assumptions suggests that, if the hypothesis concerning treatment effect is their only goal, analysts need not be too concerned about lack of model fit as long as the conditions are met. This has not been our experience; we have found in our applications and via simulation, that use of the normal model when the underlying generating model is nonnormal and heteroscedastic, can lead to poor power for the detection of a treatment effect. Consequently if interest lies in accurate detection and estimation of treatment effect on the mean, or distributional features other than the mean, use of distributional regression modelling will yield superior estimates and is strongly recommended.
Availability of data and materials
The LIPID data set is held by the National Health and Medical Research Council Clinical Trials Centre, University of Sydney. Complete individual patient data have been provided to other research groups for metaanalysis. Proposals for analyses or collaborative studies by researchers are welcome and should be submitted to the corresponding author. The R code used to analyse the data is available from the corresponding author on request.
Notes
Discrete outcomes may also be modelled in this framework; we have not discussed these in this paper.
Abbreviations
 AIC:

Akaike information criterion
 BCT:

BoxCox t distribution
 GAMLSS:

Generalized Additive Models for Location, Scale and Shape
 GEE:

generalized estimating equations
 GLM:

generalized linear model
 JSU:

Johnson’s Su distribution
 LIPID:

LongTerm Intervention with Pravastatin in Ischaemic Disease
 MRproADM:

midregional proadrenomedullin
 PDF:

probability density function
 SBP:

systolic blood pressure
 SE:

standard error
References
Nelder JA, Wedderburn RWM. Generalized linear models. J R Stat Soc Ser A. 1972;135:370–84.
Rigby RA, Stasinopoulos DM. Generalized additive models for location, scale and shape. J R Stat Soc: Ser C: Appl Stat. 2005;54:507–54.
Muntner P, Whittle J, Lynch AI, Colantonio LD, Simpson LM, Einhorn PT, et al. Visittovisit variability of blood pressure and coronary heart disease, stroke, heart failure, and mortality. Ann Intern Med. 2015;163:329–38.
Aitkin M. A history of the GLIM statistical package. Int Stat Rev. 2018;86:275–99.
Robledo KP, Marschner IC. A new algorithm for fitting semiparametric variance regression models. Comput Stat. 2021;36:1–23.
R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2020.
Stasinopoulos DM, Rigby RA, Heller GZ, Voudouris V, De Bastiani F. Flexible regression and smoothing: Using GAMLSS in R. Chapman & Hall/CRC: Boca Raton; 2017.
Umlauf N, Klein N, Zeileis A. BAMLSS: Bayesian additive models for location, scale and shape (and beyond). J Comput Graph Stat. 2018;27:612–27.
Hastie TJ, Tibshirani RJ. Generalized additive models. London: Chapman & Hall; 1990.
Wood SN. Generalized additive models. An introduction with R. 2nd ed. Chapman & Hall/CRC: Boca Raton; 2017.
LongTerm Intervention with Pravastatin in Ischaemic Disease (LIPID) Study Group. Prevention of cardiovascular events and death with pravastatin in patients with coronary heart disease and a broad range of initial cholesterol levels. N Engl J Med. 1998;339:1349–57.
Hague WE, Simes J, Kirby A, Keech AC, White HD, Hunt D, et al. Longterm effectiveness and safety of pravastatin in patients with coronary heart disease: Sixteen years of followup of the LIPID study. Circulation. 2016;133:1851–60.
Tonkin AM, Blankenberg S, Kirby A, Zeller T, Colquhoun DM, FunkeKaiser A, et al. Biomarkers in stable coronary heart disease, their modulation and cardiovascular risk: The LIPID biomarker study. Int J Cardiol. 2015;201:499–507.
Rigby RA, Stasinopoulos MD, Heller GZ, De Bastiani F. Distributions for modeling location, scale, and shape: Using GAMLSS in R. Chapman & Hall/CRC: Boca Raton; 2019.
Stasinopoulos MD, Rigby RA, De Bastiani F. Stat Model. 2018;18:248–73. https://doi.org/10.1177/1471082X18759144.
Rosenblum M, Van Der Laan MJ. Using regression models to analyze randomized trials: Asymptotically valid hypothesis tests despite incorrectly specified models. Biometrics. 2009;65:937–45.
Acknowledgements
Not applicable.
Funding
Not applicable.
Author information
Affiliations
Contributions
LIPID trial data was available to the authors at NHMRC Clinical Trials Centre, University of Sydney. All authors participated in discussions and concept development; the data was analysed and first version of the manuscript was written by GH; all authors revised the manuscript; all authors approved the final version.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The LIPID study was approved by the ethics committee at each of the 87 participating centers. A full list is given in the Supplementary Material. All patients gave written informed consent for the trial and separately for further clinic or remote followup. The study was performed in accordance with the Declaration of Helsinki.
Consent for publication
Not applicable.
Competing interests
Not applicable.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Supplementary material: MRproADM scatterplot.
Additional file 2.
LIPID site and Ethics committee.
Additional file 3.
Supplementary material: Plot of JSU distribution.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Heller, G.Z., Robledo, K.P. & Marschner, I.C. Distributional regression in clinical trials: treatment effects on parameters other than the mean. BMC Med Res Methodol 22, 56 (2022). https://doi.org/10.1186/s12874022015348
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12874022015348
Keywords
 Distributional regression
 Clinical trials
 Treatment effect
 Gamlss