Mark F. J. Steel
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Mark F. J. Steel has a total of 43 co-authors. Here are the top co-authors by number of shared publications:
- Jim E. Griffin
Department of Statistical Science University College London London UK
16 shared publications - Jacek Osiewalski
Cracow University of Economics
15 shared publications Carmen Fernández
11 shared publications- Francisco Javier Rubio
Department of Statistical Science, University College London, London, UK
9 shared publications Eduardo Ley
7 shared publicationsGary Koop
University of Edinburgh
6 shared publications
Latest Publications
- 1 December 2025
- Published by Institute of Mathematical Statistics in The Annals of Applied Statistics
- Vol. 19 (4) , 3416-3437
Media bias and polarization through the lens of a Markov switching latent space network model
News outlets are now more than ever incentivized to provide their audience with slanted news, while the intrinsic homophilic nature of online social media may exacerbate polarized opinions. Here we propose a new dynamic latent space model for time-varying online audience-duplication networks, which exploits social media content to conduct inference on media bias and polarization of news outlets. We contribute to the literature in several directions: (1) Our model provides a novel measure of media bias that combines information from both network data and text-based indicators; (2) we endow our model with Markov-switching dynamics to capture polarization regimes while maintaining a parsimonious specification; (3) we contribute to the literature on the statistical properties of latent space network models. The proposed model is applied to a set of data on the online activity of national and local news outlets from four European countries in the years 2015 and 2016. We find evidence of a strong positive correlation between our media slant measure and a well-grounded external source of media bias. In addition, we provide insight into the polarization regimes across the four countries considered.
- 1 December 2025
- Published by Institute of Mathematical Statistics in The Annals of Applied Statistics
- Vol. 19 (4) , 3416-3437
Media bias and polarization through the lens of a Markov switching latent space network model
News outlets are now more than ever incentivized to provide their audience with slanted news, while the intrinsic homophilic nature of online social media may exacerbate polarized opinions. Here we propose a new dynamic latent space model for time-varying online audience-duplication networks, which exploits social media content to conduct inference on media bias and polarization of news outlets. We contribute to the literature in several directions: (1) Our model provides a novel measure of media bias that combines information from both network data and text-based indicators; (2) we endow our model with Markov-switching dynamics to capture polarization regimes while maintaining a parsimonious specification; (3) we contribute to the literature on the statistical properties of latent space network models. The proposed model is applied to a set of data on the online activity of national and local news outlets from four European countries in the years 2015 and 2016. We find evidence of a strong positive correlation between our media slant measure and a well-grounded external source of media bias. In addition, we provide insight into the polarization regimes across the four countries considered.
- Open Access
- 1 January 2025
- Published by Institute of Mathematical Statistics in Bayesian Analysis
- Vol. -1 (-1) , 1-26
Model Uncertainty in Latent Gaussian Models with Univariate Link Function
We consider a class of latent Gaussian models with a univariate link function (ULLGMs). These are based on standard likelihood specifications (such as Poisson, Binomial, Bernoulli, Erlang, etc.), but incorporate a latent normal linear regression framework on a transformation of a key scalar parameter. We allow for model uncertainty regarding the covariates included in the regression. The ULLGM class typically accommodates extra dispersion in the data and has clear advantages for deriving theoretical properties and designing computational procedures. We formally characterize posterior existence under a convenient and popular improper prior and show that ULLGMs inherit the consistency properties from the latent Gaussian model. We propose a simple and general Markov chain Monte Carlo algorithm for Bayesian model averaging in ULLGMs. Simulation results suggest that the framework provides accurate results that are robust to some degree of misspecification. The methodology is successfully applied to measles vaccination coverage data from Ethiopia and to data on bilateral migration flows between 38 countries.
- 1 January 2025
- Published by Institute of Mathematical Statistics in Bayesian Analysis
- Vol. -1 (-1) , 1-26
Open AccessModel Uncertainty in Latent Gaussian Models with Univariate Link Function
We consider a class of latent Gaussian models with a univariate link function (ULLGMs). These are based on standard likelihood specifications (such as Poisson, Binomial, Bernoulli, Erlang, etc.), but incorporate a latent normal linear regression framework on a transformation of a key scalar parameter. We allow for model uncertainty regarding the covariates included in the regression. The ULLGM class typically accommodates extra dispersion in the data and has clear advantages for deriving theoretical properties and designing computational procedures. We formally characterize posterior existence under a convenient and popular improper prior and show that ULLGMs inherit the consistency properties from the latent Gaussian model. We propose a simple and general Markov chain Monte Carlo algorithm for Bayesian model averaging in ULLGMs. Simulation results suggest that the framework provides accurate results that are robust to some degree of misspecification. The methodology is successfully applied to measles vaccination coverage data from Ethiopia and to data on bilateral migration flows between 38 countries.
- Open Access
- 16 February 2024
- Published by Oxford University Press (OUP) in Journal of the Royal Statistical Society Series B: Statistical Methodology
- Vol. 86 (3) , 580-582
Gregor Steiner and Mark Steel’s contribution to the Discussion of ‘Parameterizing and simulating from causal models’ by Evans and Didelez
Gregor Steiner, Mark Steel; Gregor Steiner and Mark Steel’s contribution to the Discussion of ‘Parameterizing and Simulating from Causal Models’ by Evans and Di
- 16 February 2024
- Published by Oxford University Press (OUP) in Journal of the Royal Statistical Society Series B: Statistical Methodology
- Vol. 86 (3) , 580-582
Open AccessGregor Steiner and Mark Steel’s contribution to the Discussion of ‘Parameterizing and simulating from causal models’ by Evans and Didelez
Gregor Steiner, Mark Steel; Gregor Steiner and Mark Steel’s contribution to the Discussion of ‘Parameterizing and Simulating from Causal Models’ by Evans and Di
- Open Access
In Search of Lost (Mixing) Time: Adaptive Markov chain Monte Carlo schemes for Bayesian variable selection with very large p
The availability of data sets with large numbers of variables is rapidly increasing. The effective application of Bayesian variable selection methods for regression with these data sets has proved difficult since available Markov chain Monte Carlo methods do not perform well in typical problem sizes of interest. The current paper proposes new adaptive Markov chain Monte Carlo algorithms to address this shortcoming. The adaptive design of these algorithms exploits the observation that in large $p$ small $n$ settings, the majority of the $p$ variables will be approximately uncorrelated a posteriori. The algorithms adaptively build suitable non-local proposals that result in moves with squared jumping distance significantly larger than standard methods. Their performance is studied empirically in high-dimensional problems (with both simulated and actual data) and speedups of up to 4 orders of magnitude are observed. The proposed algorithms are easily implementable on multi-core architectures and are well suited for parallel tempering or sequential Monte Carlo implementations.Open AccessIn Search of Lost (Mixing) Time: Adaptive Markov chain Monte Carlo schemes for Bayesian variable selection with very large p
The availability of data sets with large numbers of variables is rapidly increasing. The effective application of Bayesian variable selection methods for regression with these data sets has proved difficult since available Markov chain Monte Carlo methods do not perform well in typical problem sizes of interest. The current paper proposes new adaptive Markov chain Monte Carlo algorithms to address this shortcoming. The adaptive design of these algorithms exploits the observation that in large $p$ small $n$ settings, the majority of the $p$ variables will be approximately uncorrelated a posteriori. The algorithms adaptively build suitable non-local proposals that result in moves with squared jumping distance significantly larger than standard methods. Their performance is studied empirically in high-dimensional problems (with both simulated and actual data) and speedups of up to 4 orders of magnitude are observed. The proposed algorithms are easily implementable on multi-core architectures and are well suited for parallel tempering or sequential Monte Carlo implementations. - Open Access
Flexible linear mixed models with improper priors for longitudinal and survival data
We propose a Bayesian approach using improper priors for hierarchical linear mixed models with flexible random effects and residual error distributions. The error distribution is modelled using scale mixtures of normals, which can capture tails heavier than those of the normal distribution. This generalisation is useful to produce models that are robust to the presence of outliers. The case of asymmetric residual errors is also studied. We present general results for the propriety of the posterior that also cover cases with censored observations, allowing for the use of these models in the contexts of popular longitudinal and survival analyses. We consider the use of copulas with flexible marginals for modelling the dependence between the random effects, but our results cover the use of any random effects distribution. Thus, our paper provides a formal justification for Bayesian inference in a very wide class of models (covering virtually all of the literature) under attractive prior structures that limit the amount of required user elicitation.
Open AccessFlexible linear mixed models with improper priors for longitudinal and survival data
We propose a Bayesian approach using improper priors for hierarchical linear mixed models with flexible random effects and residual error distributions. The error distribution is modelled using scale mixtures of normals, which can capture tails heavier than those of the normal distribution. This generalisation is useful to produce models that are robust to the presence of outliers. The case of asymmetric residual errors is also studied. We present general results for the propriety of the posterior that also cover cases with censored observations, allowing for the use of these models in the contexts of popular longitudinal and survival analyses. We consider the use of copulas with flexible marginals for modelling the dependence between the random effects, but our results cover the use of any random effects distribution. Thus, our paper provides a formal justification for Bayesian inference in a very wide class of models (covering virtually all of the literature) under attractive prior structures that limit the amount of required user elicitation.
- Open Access
- 10 July 2017
- Published by Springer Nature in Statistical Methods & Applications
- Vol. 27 (2) , 207-218
Discussion of “Nonparametric Bayesian Inference in Applications”: Bayesian nonparametric methods in econometrics
- 10 July 2017
- Published by Springer Nature in Statistical Methods & Applications
- Vol. 27 (2) , 207-218
Open AccessDiscussion of “Nonparametric Bayesian Inference in Applications”: Bayesian nonparametric methods in econometrics
- Open Access
- 1 July 2017
- Published by Elsevier in Econometrics and Statistics
- Vol. 3, 73-88
Incorporating unobserved heterogeneity in Weibull survival models: A Bayesian approach
- 1 July 2017
- Published by Elsevier in Econometrics and Statistics
- Vol. 3, 73-88
Open AccessIncorporating unobserved heterogeneity in Weibull survival models: A Bayesian approach
- Open Access
Methods and Tools for Bayesian Variable Selection and Model Averaging in Univariate Linear Regression
In this paper we briefly review the main methodological aspects concerned with the application of the Bayesian approach to model choice and model averaging in the context of variable selection in regression models. This includes prior elicitation, summaries of the posterior distribution and computational strategies. We then examine and compare various publicly available {\tt R}-packages for its practical implementation summarizing and explaining the differences between packages and giving recommendations for applied users. We find that all packages reviewed lead to very similar results, but there are potentially important differences in flexibility and efficiency of the packages.
Open AccessMethods and Tools for Bayesian Variable Selection and Model Averaging in Univariate Linear Regression
In this paper we briefly review the main methodological aspects concerned with the application of the Bayesian approach to model choice and model averaging in the context of variable selection in regression models. This includes prior elicitation, summaries of the posterior distribution and computational strategies. We then examine and compare various publicly available {\tt R}-packages for its practical implementation summarizing and explaining the differences between packages and giving recommendations for applied users. We find that all packages reviewed lead to very similar results, but there are potentially important differences in flexibility and efficiency of the packages.
- 14 July 2016
- Published by Oxford University Press (OUP) in Journal of the Royal Statistical Society Series A: Statistics in Society
- Vol. 180 (2) , 613-631
Bayesian Survival Modelling of University Outcomes
Summary: Dropouts and delayed graduations are critical issues in higher education systems world wide. A key task in this context is to identify risk factors associated with these events, providing potential targets for mitigating policies. For this, we employ a discrete time competing risks survival model, dealing simultaneously with university outcomes and its associated temporal component. We define survival times as the duration of the student's enrolment at university and possible outcomes as graduation or two types of dropout (voluntary and involuntary), exploring the information recorded at admission time (e.g. educational level of the parents) as potential predictors. Although similar strategies have been previously implemented, we extend the previous methods by handling covariate selection within a Bayesian variable selection framework, where model uncertainty is formally addressed through Bayesian model averaging. Our methodology is general; however, here we focus on undergraduate students enrolled in three selected degree programmes of the Pontificia Universidad Católica de Chile during the period 2000–2011. Our analysis reveals interesting insights, highlighting the main covariates that influence students’ risk of dropout and delayed graduation.
- 14 July 2016
- Published by Oxford University Press (OUP) in Journal of the Royal Statistical Society Series A: Statistics in Society
- Vol. 180 (2) , 613-631
Bayesian Survival Modelling of University Outcomes
Summary: Dropouts and delayed graduations are critical issues in higher education systems world wide. A key task in this context is to identify risk factors associated with these events, providing potential targets for mitigating policies. For this, we employ a discrete time competing risks survival model, dealing simultaneously with university outcomes and its associated temporal component. We define survival times as the duration of the student's enrolment at university and possible outcomes as graduation or two types of dropout (voluntary and involuntary), exploring the information recorded at admission time (e.g. educational level of the parents) as potential predictors. Although similar strategies have been previously implemented, we extend the previous methods by handling covariate selection within a Bayesian variable selection framework, where model uncertainty is formally addressed through Bayesian model averaging. Our methodology is general; however, here we focus on undergraduate students enrolled in three selected degree programmes of the Pontificia Universidad Católica de Chile during the period 2000–2011. Our analysis reveals interesting insights, highlighting the main covariates that influence students’ risk of dropout and delayed graduation.
- Open Access
Bayesian modelling of skewness and kurtosis with two-piece scale and shape distributions
We formalise and generalise the definition of the family of univariate double two--piece distributions, obtained by using a density--based transformation of unimodal symmetric continuous distributions with a shape parameter. The resulting distributions contain five interpretable parameters that control the mode, as well as the scale and shape in each direction. Four-parameter subfamilies of this class of distributions that capture different types of asymmetry are discussed. We propose interpretable scale and location-invariant benchmark priors and derive conditions for the propriety of the corresponding posterior distribution. The prior structures used allow for meaningful comparisons through Bayes factors within flexible families of distributions. These distributions are applied to data from finance, internet traffic and medicine, comparing them with appropriate competitors.
Open AccessBayesian modelling of skewness and kurtosis with two-piece scale and shape distributions
We formalise and generalise the definition of the family of univariate double two--piece distributions, obtained by using a density--based transformation of unimodal symmetric continuous distributions with a shape parameter. The resulting distributions contain five interpretable parameters that control the mode, as well as the scale and shape in each direction. Four-parameter subfamilies of this class of distributions that capture different types of asymmetry are discussed. We propose interpretable scale and location-invariant benchmark priors and derive conditions for the propriety of the corresponding posterior distribution. The prior structures used allow for meaningful comparisons through Bayes factors within flexible families of distributions. These distributions are applied to data from finance, internet traffic and medicine, comparing them with appropriate competitors.