A Bayesian approach to linear regression in astronomy
Date Issued
2016
Author(s)
Abstract
Linear regression is common in astronomical analyses. I discuss a Bayesian hierarchical modelling of data with heteroscedastic and possibly correlated measurement errors and intrinsic scatter. The method fully accounts for time evolution. The slope, the normalization, and the intrinsic scatter of the relation can evolve with the redshift. The intrinsic distribution of the independent variable is approximated using a mixture of Gaussian distributions whose means and standard deviations depend on time. The method can address scatter in the measured independent variable (a kind of Eddington bias), selection effects in the response variable (Malmquist bias), and departure from linearity in form of a knee. I tested the method with toy models and simulations and quantified the effect of biases and inefficient modelling. The R-package LIRA (LInear Regression in Astronomy) is made available to perform the regression.
Volume
455
Issue
2
Start page
2149
Issn Identifier
0035-8711
Ads BibCode
2016MNRAS.455.2149S
Rights
open.access
File(s)![Thumbnail Image]()
Loading...
Name
82_A Bayesian approach to linear regression in astronomy_MNRAS-2016-Sereno-2149-62.pdf
Description
PDF editoriale
Size
939.68 KB
Format
Adobe PDF
Checksum (MD5)
32b915b6d89ed397386fd380068b63bf
