Solution (source code)

= 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. \b[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.<On the Use of Trigonometric Parallaxes for the Calibration of Luminosity Systems: Theory by T. E. Lutz and D. H. Kelker>[PASP 85, 573–578; classical uniform-density calibration 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,
$$
p(r\mid\widehat\varpi)\propto
\exp\!\left[-\frac{(\widehat\varpi-1/r)^2}{2\sigma_\varpi^2}\right]p(r),\qquad r>0.
$$
A <proper distance prior for parallax inference> might have $p(r)\propto r^2e^{-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.<Planck 2013 results. XXIX. Planck catalogue of Sunyaev-Zeldovich sources by Planck Collaboration>[abstract and selection-function/completeness analysis, 2013 preprint]{ref} Calibrators, artificial-source tests and a forward model of detection, measurement and follow-up are the practical checks on any correction.