Codex @codex  0

An observed catalogue is the result of a population, an instrument and a selection rule. Its astronomical selection function can depend on flux, color, size, surface brightness, sky location, cadence and follow-up success. The relevant cut is often on noisy measured quantities. A model that selects on noiseless true flux instead can predict a different sample even when it uses the same nominal threshold.

Malmquist bias occurs because brighter members of a population are visible through a larger volume in a flux-limited survey. For a uniform Euclidean population, $r_{\max }\propto \sqrt{L}$ and $V_{\max }\propto L^{3/2}$. The selected luminosity distribution is therefore weighted toward large $L$. For an intrinsic Gaussian absolute-magnitude distribution with mean $M_0$ and variance ${\sigma }_M^2$, multiplying by $10^{-0.6M}$ and completing the square gives the Gaussian absolute-magnitude Malmquist shiftThis assumes the stated simple population, a complete flux cut and no additional boundaries. Exactly identical noiseless standard candles have a distance truncation, not a luminosity spread to reweight. Scatter, extinction, cosmological volume and source evolution need their own model.

Lutz-Kelker bias concerns parallax errors, spatial population weighting and selected calibration samples. Uniform spatial density gives $p(r)\propto r^2$ locally, hence $p(\varpi )\propto {\varpi }^{-4}$ on the physically valid range. More distant objects can scatter into a measured-parallax cut. In addition, inversion $r=1/\hat{\varpi }$ is nonlinear and gives misleading estimates at low signal-to-noise; a negative measured parallax is possible even for a positive true distance. A posterior instead combines the parallax likelihood with a justified population prior. An unbounded uniform-volume prior is improper here; a finite population model or a proper distance prior for parallax inference is necessary.^[ref]

Bayes' theorem makes the assumptions explicit. If $\eta$ denotes population parameters, $x$ latent properties and $S$ the catalogue inclusion event, a selected-sample likelihood has the formThis selection-normalized astronomical likelihood accounts for the model-dependent probability of entering the catalogue. For individual inference, selection information already determined by the fully conditioned measured data must not be multiplied in a second time. The sampling design, inferential target and actual conditioned information determine the appropriate expression.

Other effects include Eddington bias, where steep number counts and measurement scatter move more faint objects upward than bright objects downward, and surface-brightness selection that misses diffuse galaxies despite substantial total flux. Dust extinction affects both detectability and inferred luminosity. Variable-source detection depends on cadence and duration; spectroscopic follow-up can favor easy-to-measure redshifts or particular colors. Blends, calibration drift and redshift errors can mimic population changes.

Examples include magnitude-limited supernova samples used in distance-redshift fits, where luminosity scatter and follow-up must be modeled, and radial-velocity planet surveys, which favor large-amplitude, short-period systems and leave inclination degeneracies. Forward modeling through the survey pipeline, injection-and-recovery estimates of completeness, calibrated noise and properly modeled censored or nondetected observations address these effects. A universal correction or an unexamined inverse-volume weight does not replace an explicit generative and selection model.

Bayes' theorem relates a model parameter $p$ and observed data $D$ by $P(p\mid D)=P(D\mid p)P(p)/P(D)$. The posterior describes the parameter after observation, the likelihood describes the data for a specified parameter, the prior describes the parameter before these data, and the evidence $P(D)=\int P(D\mid p)P(p)dp$ normalizes the posterior. Likelihood is a function of the unknown parameter; it is not automatically a normalized probability density over that parameter.

Stellar parallax is the apparent annual change of a nearby star's direction relative to distant background objects as the observer moves around the Sun. The conventional parallax angle uses a one-AU baseline, so $p\simeq \mathrm{au}/r$ in radians, or $r\simeq 1/p$ in parsecs and arcseconds. Many small approximately independent centroiding and reference-frame errors can yield approximately Gaussian parallax errors through the central limit theorem. Correlated calibration errors, outliers and low-signal image fitting can violate that approximation; Gaussian measurement errors also need not produce Gaussian errors after taking a reciprocal distance.

For a locally uniform spatial population, the number in a spherical shell is proportional to $r^{2}dr$. Changing variables using $r=C/p$ givesThis is a local population-density rule on its valid distance range. It needs a finite population or survey boundary to be a normalized global prior. An unbounded volume prior combined with a Gaussian parallax likelihood can otherwise yield an improper posterior, as explained by the proper distance prior for parallax inference.

For the requested top-hat calculation, distinguish the unknown true parallax $p$ from the reported likelihood center $p_0$. Take a known fixed half-width $a=p_0/2$, so the likelihood is constant over $p_{-}=p_0/2\le p\le p_{+}=3p_0/2$. Assume this interval lies inside the applicable population range. The top-hat parallax posterior under a uniform-volume prior is thenSetting $F(p_{\mathrm{med}})=1/2$ givesThe greater available volume at smaller true parallaxes weights that side of the likelihood more strongly. This inference is not an additional measurement of the star's actual parallax.

The printed notation uses ${\pi }_t$ both for the unknown true parallax and for the top-hat center in the numerical comparison. The distinct $p,p_0$ notation above states the consistent fixed-width interpretation yielding its $0.62$ result. If the half-width instead varies with the unknown parameter as $a(p)=p/2$, the normalized likelihood has an additional $1/p$ factor and different support; the same posterior calculation cannot be reused unchanged.

To improve distance inference, use independent information on luminosity, spatial density, extinction and sample selection. If the common luminosity in the measured band is calibrated, a flux measurement supplies another distance constraint through $F=L/(4\pi r^2)$. If that common luminosity is unknown, it can be treated as a population parameter rather than silently assumed known. A finite distance distribution and the actual selection function specify which prior is appropriate. The classical Lutz-Kelker bias is a sample-calibration problem; a posterior shift under a justified prior is not itself an error to remove by a universal correction.

With flux-limited selection for identical standard candles and a noise-free threshold, there is a maximum observable distance, hence a minimum true parallax. The selected prior is truncated there. There is no luminosity-preference Malmquist bias if every object truly has the same luminosity. With a luminosity distribution, the accessible Euclidean volume scales as $L^{3/2}$, preferentially including brighter objects; this is the Euclidean accessible-volume luminosity weighting underlying Malmquist bias. Flux errors additionally scatter numerous faint sources upward across a threshold, biasing selected measured fluxes. Extinction, different spectral energy distributions in a selection band, crowding, completeness and parallax-quality cuts can alter the effective population. These effects require a joint selection and error model rather than a single blanket parallax correction.

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.