 Solution

The observed population is the parent population filtered by an astronomical selection function and measurement errors. Malmquist bias arises because a flux-limited galaxy selection admits luminous objects from larger volumes: in Euclidean space $d_{\max }\propto L^{1/2}$ and $V_{\max }\propto L^{3/2}$. Thus observed luminosity counts are weighted by $V_{\max }$. Even a nominal standard candle with intrinsic scatter is preferentially observed on its bright side near the limit; ignoring this can underestimate distances. A volume-limited sample or a justified $1/V_{\max }$ weighting helps, subject to completeness, geometry and evolution assumptions. Separately, Eddington bias occurs when noise scatters a steeply rising population across a detection threshold: more faint objects enter than bright objects leave.

Lutz-Kelker bias concerns noisy stellar parallaxes and population/selection weighting. Locally uniform spatial density gives $p(r)\propto r^2$ and, by changing variables, $p(\varpi )\propto {\varpi }^{-4}$ over the model's valid finite range. More distant stars are available to scatter into a parallax-limited sample than nearby stars to scatter out. Inverting an uncertain parallax also introduces a distinct nonlinear-estimation problem. There is no universal magnitude correction for every star: spatial density, luminosity, extinction and how the sample was selected all matter. The classical treatment motivates this issue; it should not replace an explicit generative model.^[ref]

Bayes' theorem combines the measurement likelihood function with a justified prior distribution. For example, with a Gaussian parallax error and consistent distance/parallax units,A proper distance prior for parallax inference might have $p(r)\propto r^2{e}^{-r/L}$ with a justified scale $L$; extending $r^2$ to infinity gives an improper posterior here. Negative measured parallaxes remain valid noisy data. For selection conditioning in Bayesian inference, distinguish individual inference from population inference: if selection is determined entirely by the exact data already conditioned on, it adds no further information about that individual's parameters. For population parameters $\Lambda$, the selected-data likelihood must instead include the detection normalization $\alpha (\Lambda )=\int S(d)p(d\mid \Lambda )dd$. Applying a cut twice is also an error.

Selection can depend on surface brightness, angular resolution, dust, sky coverage, colour and observing cadence as well as total flux. Crowding and unresolved binaries bias distance calibrators. Transit surveys favour large planets and short periods, with an additional geometric transit probability approximately $R_{\star }/a$ for small planets on circular orbits; injection/recovery tests quantify the remaining transit survey completeness. Cluster catalogues select noisy mass proxies from a steep abundance distribution. Planck's 2013 SZ catalogue explicitly characterized detection completeness and reliability rather than identifying every detection with an unbiased mass census.^[ref] Calibrators, artificial-source tests and a forward model of detection, measurement and follow-up are the practical checks on any correction.