Pith. sign in

REVIEW 2 major objections 4 minor 102 references

Learning High-dimensional Gaussians from Censored Data

T0 review · 2 major / 4 minor · reviewed 2026-08-16 · deepseek-v4-flash

Pith's one-line read Two censored-Gaussian models are learnable with polynomial samples

desk verdict The self-censoring result is a solid extension of DGTZ18, but the linear-thresholding half has a real gap in the Langevin gradient-sampling analysis that needs to be fixed before the runtime claim is credible. read the letter →

arxiv 2504.19446 v1 pith:YAPI73K6 submitted 2025-04-28 cs.LG cs.CCmath.STstat.MLstat.TH

classification cs.LGcs.CCmath.STstat.MLstat.TH MSC 62H1262F1062G07
keywords missingnotatrandomGaussiandistributionlearningcensoreddatatruncatedsamplesself-censoringlinearthresholdingprojectedstochasticgradientdescentLangevinMonteCarlo
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper asks when a high-dimensional Gaussian can be recovered from samples in which coordinates are missing not at random, with the missingness rule known. It answers for two mechanisms. Under self-censoring, where coordinate $i$ is seen iff $y_i$ satisfies a known predicate, if every pair of coordinates is jointly observed with probability at least $\alpha$, both mean and covariance can be estimated to small total-variation error using $\tilde{O}(d^2(\lambda_{\max}/\lambda_{\min})^2/(\alpha\varepsilon^2))$ samples. Under linear-thresholding missingness, where coordinate $i$ is seen iff $v_i^\top y \le b_i$, the mean can be estimated in polynomial time with known covariance, provided any small subset of coordinates is observed with probability at least $\alpha$ and a $\gamma$-anchoring set of always-observed coordinates exists. If these results hold, MNAR corruption is not inherently an obstacle to efficient Gaussian estimation.

What carries the argument

For self-censoring, the machinery is coordinate-pair decomposition plus the truncated-Gaussian SGD used in dimensions 1 and 2, with Assumption 1.1 ensuring enough complete pairs. For linear thresholding, the machinery is the negative log-likelihood $\ell(\mu)$ whose gradient requires sampling from the conditional Gaussian $N(\mu,\Sigma)$ restricted to $\{y: S(y)=A,\, y_A=x\}$; because linear thresholds define convex sets, conditioning a Gaussian on a convex set can only reduce variance, which yields strong convexity of $\ell$. The gradient samples are obtained by projected Langevin Monte Carlo rather than rejection sampling, since rejection costs would grow exponentially when $\mu$ drifts from $\mu^*$; the anchoring assumption bounds the conditional pattern probabilities and keeps the sampler polynomial. A projection step keeps iterates in a Mahalanobis ball so the stochastic gradients have bounded centralized second moment.

What would settle it

Instantiate the linear-thresholding model with $d=2$, $\Sigma=I$, coordinate 2 always observed, coordinate 1 observed when $y_1 \le t$, and run Algorithm 3 as $t$ sweeps from 0 toward infinity. The theorem predicts the sample bound grows polynomially in $1/\alpha$ and $1/\gamma$; if the observed degradation in reaching $\|\hat\mu - \mu^*\|_\Sigma \le \varepsilon$ is exponential in $1/\gamma$, the assumptions as stated are not sufficient.

Watch

Extended reading notes

Core claim

The central claim is that estimation under MNAR Gaussian data is computationally tractable exactly when the missingness rule has a decomposable or anchored structure. For self-censoring, the paper shows the $d$-dimensional problem reduces to one- and two-dimensional truncated-Gaussian subproblems, one per coordinate or coordinate pair; each subproblem is solved with the known truncated-sample SGD, and the pieces are assembled into a consistent covariance matrix. For linear-thresholding missingness with known $\Sigma$, the paper defines the population negative log-likelihood as a function of $\mu$ and proves it is $\lambda$-strongly convex with $\lambda = \alpha\beta/\lambda_{\max}(\Sigma)$, so projected SGD converges. The key algorithmic step is a projected Langevin Monte Carlo sampler that draws approximately from the conditional Gaussian over unseen coordinates given the observed pattern, with bias and second-moment bounds that make the PSGD convergence theorem applicable. The result is Theorem 1.5: with $\mathrm{poly}(d, 1/\alpha, 1/\beta, 1/\gamma, \lambda_{\max}/\lambda_{\min}, 1/\varepsilon, \log(1/\delta))$ samples, the estimator meets $\|\hat\mu - \mu^*\|_\Sigma \le \varepsilon$.

Load-bearing premise

The load-bearing premise is Assumption 1.4: there must be an anchoring set of coordinates that is observed in every sample and whose values make every missingness pattern occur with probability either 0 or at least $\gamma$; without such a set, the gradient sampler's cost becomes exponential and the strong-convexity estimate depending on pattern probabilities collapses.

Editorial extensions

If this is right

  • Under self-censoring with every pair observed with probability at least $\alpha$, both mean and covariance are recoverable to TV error $O(\varepsilon)$ with $\tilde{O}(d^2(\lambda_{\max}/\lambda_{\min})^2/(\alpha\varepsilon^2))$ samples, so heavy coordinate-wise missingness is not a barrier.
  • Under linear-thresholding missingness with known covariance, the mean can be found in polynomial time whenever small subsets are observed with probability at least $\alpha$ and a $\gamma$-anchor set exists, giving a finite-sample MNAR guarantee without an all-or-nothing observation model.
  • The self-censoring result needs only pairwise observation, not a constant fraction of fully observed samples, so it strictly extends the earlier truncated-Gaussian setting.
  • The $\Omega(1/\sqrt{\lambda_{\min}})$ lower bound shows some eigenvalue dependence is unavoidable for self-censoring, even when uncensored distributions are far apart in TV distance.
  • Because the likelihood is strongly convex only for linear-thresholding missingness, the same optimization route will not directly work for general self-censoring rules.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural next step, left open by the paper, is unknown $\Sigma$ under linear-thresholding missingness; the Hessian then involves a fourth moment, but the anchoring condition may still make a second-order method or a moment-based correction work.
  • The anchoring assumption is reminiscent of anchor words in topic modeling; if anchor coordinates are domain-known, the algorithm applies directly, and the dependence on $1/\gamma$ makes it worth measuring $\gamma$ in practice before deployment.
  • The coordinate-pair reduction for self-censoring suggests a general principle: when a missingness model factorizes over small subsets, estimation can be assembled from small truncated subproblems, potentially extending to other exponential families with pairwise-sufficient statistics.
  • For practical deployment, the $\mathrm{poly}(1/\gamma)$ and $\mathrm{poly}(1/\alpha)$ factors mean the algorithm's usefulness depends on how large these probabilities are; identifying real missingness mechanisms where $\gamma$ is not tiny would be a valuable empirical study.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper studies high-dimensional Gaussian estimation from missing-not-at-random (MNAR) data under two missingness models. In the self-censoring model, coordinate i is observed exactly when the underlying value lies in a known set S_i; the authors reduce estimation to one- and two-dimensional truncated-Gaussian subroutines from DGTZ18 and prove poly(d,1/epsilon) sample complexity under a pairwise-observation assumption. In the linear-thresholding model, coordinate i is observed when v_i^T y <= b_i, the covariance is known, and the goal is mean estimation; the paper proposes a projected stochastic gradient descent on the negative log-likelihood, with gradients sampled by a projected Langevin algorithm, and claims poly(d,1/alpha,1/beta,1/gamma,lambda_max/lambda_min,1/epsilon,log(1/delta)) samples and time under Assumptions 1.3 and 1.4.

Significance. If the linear-thresholding result were fully established, it would be a substantial advance: it would give the first polynomial-time algorithm for high-dimensional MNAR mean estimation under explicit, non-asymptotic conditions, going beyond the all-or-nothing censoring of DGTZ18. The self-censoring part is a clean and mostly convincing reduction that leverages published external subroutines, so it is not circular. The paper is commendably explicit about its assumptions and gives concrete pseudocode. However, the proof of the central runtime claim for the Langevin gradient sampler has a serious gap, so the main advertised theorem (Theorem 1.5) is not established as written.

major comments (2)
  1. [C.3, Lemma C.2 and Theorem 4.6] The proof that Algorithm 5 mixes in polynomial time is not sound. In the notation of Algorithm 5, K={z in R^{d-|A|} : S(x∘z)=A} is an affine slice of the missingness pattern; hence the probability Pr_{w∼N(0,Sigma)}[w in K] used in the proof of Lemma C.2 is zero, and if K is instead interpreted in the conditional space N(mu_cond,Sigma_cond), Assumption 1.4 gives no lower bound on its probability, since it conditions on the anchor set C rather than on the full observed vector y_A=x. For example, take d=3, C={1}, S(y)={1,2} iff y_3>0 (and {1} otherwise), and Corr(y_2,y_3)=-0.9. Then Assumption 1.4 holds with gamma=1/2, but for an observed sample with A={1,2} and x_2=t, the measure of K under the conditional Gaussian is Phi(-0.9t/sqrt(0.19)), which decays as exp(-Omega(t^2)). The proof does not show that such t cannot occur with non-negligible probability, nor does it adapt the Langevin radius to the actual measure of K. Consequently the polynomial bound on the number of Langevin steps in Theorem 4.6 is not established, and the running time claimed in Theorem 1.5 has no term that could absorb this exponential degradation.
  2. [4.3.1, Lemma 4.5] The strong-convexity proof is incomplete as printed. The text contains the garbled string 'enecccfckgvkkjlfnlclcggguvjldtdurlrlltnecnid' and skips the block-matrix manipulations needed to pass from Eq. (14) to the final inequality v^T nabla^2 ell(mu) v >= alpha v_H^T Sigma_HH^{-1} v_H. The displayed formulas suggest the identity may be correct, but a referee cannot verify the lemma without the missing algebra, and strong convexity is load-bearing for the SGD convergence argument. The proof needs to be rewritten with the omitted steps supplied.
minor comments (4)
  1. [Appendix C.3] Lemma C.2 is stated twice with the same label; the second occurrence should be renumbered or removed to avoid confusion.
  2. [Algorithm 5] Line 4 uses B_Sigma(W^{-1}mu_cond,R), but B_Sigma is a ball in R^d while the argument W^{-1}mu_cond lies in R^{d-|A|}; the intended metric for the projected Langevin dynamics should be clarified.
  3. [Algorithm 2] Line 7 says to remove rows containing at least one '*' but a missing entry is not previously defined; the intended definition of a missing-value indicator should be stated.
  4. [Lemma 3.4] The lower-bound proof concludes that o(max{1/epsilon, sqrt(lambda)/epsilon}) samples are insufficient, but the lemma statement is phrased as o(1/sqrt(lambda_min)) for epsilon = Omega(sqrt(lambda_min)); the relationship between these parameters should be spelled out for the reader.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central reductions rely on published external subroutines, and no fitted parameter or prior conclusion is relabeled as a prediction.

full rationale

The paper's derivation chain is a reduction, not a restatement of its assumptions. For self-censoring, Assumption 1.1 is used only to guarantee that enough samples contain each pair of coordinates; the per-pair estimates are delegated to the published DGTZ18 truncated-Gaussian algorithm (Theorem 3.1), whose stated assumptions do not include the self-censoring model. For linear thresholding, Assumptions 1.3 and 1.4 are explicit hypotheses of Theorem 1.5, not consequences of the algorithm, and the proof uses the generic SGD convergence lemma (Lemma 6 of CDIZ22) and the projected-Langevin guarantee of BEL18; both are external, parameter-free lemmas with stated assumptions that do not contain the target mean-estimation theorem. The authors' earlier DGTZ18 and CDIZ22 works are cited heavily, but in each cited instance they supply an independent published result rather than a conclusion that is being derived in this paper. No fitted parameter is renamed as a prediction, and no uniqueness theorem from the authors' prior work is used to force the model choice. There is a separate correctness concern that is not circularity: in the proof of Lemma C.1 the paper states 'Setting δ = γ/2, we get that: Pr[v∈BΣ(0,R1)∩K]≥γ/2', which implicitly requires a lower bound on Pr[v∈K]; Assumption 1.4 bounds pattern probabilities only after conditioning on the anchor set C, not on the full observed set A, so this step is unsupported. That gap affects the validity of the Langevin runtime guarantee but does not make the claim reduce to its own inputs by definition. Therefore no structural circularity is present, and the appropriate circularity score is 0.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The paper does not fit any constant to data; alpha, beta, gamma are lower-bound assumptions, not fitted values. The algorithmic parameters (eta, R, M) are prescribed in terms of d, lambda_max, lambda_min, S, gamma, epsilon. The main burden falls on the modeling assumptions (pairwise visibility, small-subset visibility, anchoring) and on imported results from prior work by the same research group (DGTZ18, CDIZ22), which are published and not re-derived here.

assumptions (6)
  • domain assumption Assumption 1.1: for every pair i,j, Pr[y_i in S_i, y_j in S_j] >= alpha.
    Invoked in Section 3 to guarantee enough samples for 1D and 2D truncated subproblems; if it fails, the decomposition into pairwise problems has no sample support.
  • domain assumption Assumption 1.3: any subset A of size at most beta d is observed with probability at least alpha.
    Used in Lemma 4.5 for strong convexity and in the initialization (Algorithm 2) to ensure every block of coordinates is present with non-negligible probability.
  • domain assumption Assumption 1.4: existence of a gamma-anchoring set C as in Definition 4.1.
    Needed for the Langevin gradient estimator: conditioned on y_C the pattern probabilities are bounded below, which controls the cost of sampling from K and the bias bound in Corollary 4.7.
  • standard math DGTZ18 Theorem 3.1 (truncated Gaussian estimation) can be applied to 1D and 2D subproblems.
    Imported as a black box in Section 3; the paper does not re-prove it, so the central bounds for self-censoring inherit the correctness of this published subroutine.
  • standard math KP77 Corollary 2.1: conditioning a Gaussian on a convex set does not increase its variance in the Loewner order.
    Used in Lemma 4.2 and Lemma 4.5 to lower-bound the Hessian of the negative log-likelihood; this is a classical result cited and not proved.
  • standard math BEL18 Theorem 1: convergence of projected Langevin Monte Carlo to the truncated Gaussian.
    Basis for Theorem 4.6, which states that Algorithm 5's output is within TV epsilon of the conditional Gaussian; assumed to hold with the stated polynomial parameters.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Learning High-dimensional Gaussians from Censored Data." pith.science (2026). https://pith.science/paper/YAPI73K6

@misc{pith2026250419446,
  author       = {Pith},
  title        = {Pith review of: Learning High-dimensional Gaussians from Censored Data},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YAPI73K6}},
  note         = {Machine review of arXiv:2504.19446}
}
abstract

We provide efficient algorithms for the problem of distribution learning from high-dimensional Gaussian data where in each sample, some of the variable values are missing. We suppose that the variables are missing not at random (MNAR). The missingness model, denoted by $S(y)$, is the function that maps any point $y$ in $R^d$ to the subsets of its coordinates that are seen. In this work, we assume that it is known. We study the following two settings: (i) Self-censoring: An observation $x$ is generated by first sampling the true value $y$ from a $d$-dimensional Gaussian $N(\mu*, \Sigma*)$ with unknown $\mu*$ and $\Sigma*$. For each coordinate $i$, there exists a set $S_i$ subseteq $R^d$ such that $x_i = y_i$ if and only if $y_i$ in $S_i$. Otherwise, $x_i$ is missing and takes a generic value (e.g., "?"). We design an algorithm that learns $N(\mu*, \Sigma*)$ up to total variation (TV) distance epsilon, using $poly(d, 1/\epsilon)$ samples, assuming only that each pair of coordinates is observed with sufficiently high probability. (ii) Linear thresholding: An observation $x$ is generated by first sampling $y$ from a $d$-dimensional Gaussian $N(\mu*, \Sigma)$ with unknown $\mu*$ and known $\Sigma$, and then applying the missingness model $S$ where $S(y) = {i in [d] : v_i^T y <= b_i}$ for some $v_1, ..., v_d$ in $R^d$ and $b_1, ..., b_d$ in $R$. We design an efficient mean estimation algorithm, assuming that none of the possible missingness patterns is very rare conditioned on the values of the observed coordinates and that any small subset of coordinates is observed with sufficiently high probability.

Figures

Figures reproduced from arXiv: 2504.19446 by the authors.

Figure 1
Figure 1. In this example, the self-censoring missingness mechanism is as follows: Each coordinate [PITH_FULL_IMAGE:figures/full_fig_p020_1.png] view at source ↗
Figure 2
Figure 2. An illustration of convex sets in Section 4.3.2. Proof of Theorem 4.6. With R2 as in Lemma C.2, set R = R2 + S, so that L = BΣ(µ, R) ∩ K. Since R2 > R1, Lemma C.1 implies that L contains a ball of radius r. On the other hand, by Lemma C.2, Pr[w ∈ L/ ] ≤ ε 4 Pr[w ∈ K] for w ∼ N (µ, Σ), which implies that the truncation of N (µ, Σ) to K and to L are at most ε/2 far from each other in TV distance. We can now use the ma… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

102 extracted references · 64 canonical work pages

  1. [1]

    Statistical methods in medical research

    Peter Armitage, Geoffrey Berry, and John Nigel Scott Matthews. Statistical methods in medical research . John Wiley & Sons, 2008

  2. [2]

    Classification of alcohols obtained by qcm sensors with different characteristics using abc based neural network

    M Fatih Adak, Peter Lieberzeit, Purim Jarujamrus, and Nejat Yumusak. Classification of alcohols obtained by qcm sensors with different characteristics using abc based neural network. Engineering Science and Technology, an International Journal , 23(3):463--469, 2020

  3. [3]

    Multiple imputation for missing data: A cautionary tale

    Paul D Allison. Multiple imputation for missing data: A cautionary tale. Sociological methods & research , 28(3):301--309, 2000

  4. [4]

    Regression analysis when the dependent variable is truncated normal

    Takeshi Amemiya. Regression analysis when the dependent variable is truncated normal. Econometrica: Journal of the Econometric Society , pages 997--1016, 1973

  5. [5]

    Regression models: Censored, sample selected, or truncated data , volume 111

    Richard Breen et al. Regression models: Censored, sample selected, or truncated data , volume 111. Sage, 1996

  6. [6]

    Multiple regression and estimation of the mean of a multivariate normal distribution

    Alvin J Baranchik. Multiple regression and estimation of the mean of a multivariate normal distribution. Technical report, STANFORD UNIV CALIF, 1964

  7. [7]

    The art of progressive censoring

    N Balakrishnan and Erhard Cramer. The art of progressive censoring. Statistics for industry and technology , 2014

  8. [8]

    Robustly learning mixtures of k arbitrary gaussians

    Ainesh Bakshi, Ilias Diakonikolas, He Jia, Daniel M Kane, Pravesh K Kothari, and Santosh S Vempala. Robustly learning mixtures of k arbitrary gaussians. In Proceedings of the 54th Annual ACM SIGACT Symposium on Theory of Computing , pages 1234--1247, 2022

Show all 102 references
  1. [9]

    Sampling from a log-concave distribution with projected langevin monte carlo

    S \'e bastien Bubeck, Ronen Eldan, and Joseph Lehec. Sampling from a log-concave distribution with projected langevin monte carlo. Discrete & Computational Geometry , 59(4):757--783, 2018

  2. [10]

    Handling missing data in survey research

    J Michael Brick and Graham Kalton. Handling missing data in survey research. Statistical methods in medical research , 5(3):215--238, 1996

  3. [11]

    Maximum likelihood estimation of the multivariate normal mixture model

    Otilia Boldea and Jan R Magnus. Maximum likelihood estimation of the multivariate normal mixture model. Journal of the American Statistical Association , 104(488):1539--1549, 2009

  4. [12]

    Identification in missing data models represented by directed acyclic graphs

    Rohit Bhattacharya, Razieh Nabi, Ilya Shpitser, and James M Robins. Identification in missing data models represented by directed acyclic graphs. In Uncertainty in Artificial Intelligence , pages 1149--1158. PMLR, 2020

  5. [13]

    Mean and variance of truncated normal distributions

    Donald R Barr and E Todd Sherrill. Mean and variance of truncated normal distributions. The American Statistician , 53(4):357--361, 1999

  6. [14]

    Polynomial learning of distribution families

    Mikhail Belkin and Kaushik Sinha. Polynomial learning of distribution families. In 2010 IEEE 51st Annual Symposium on Foundations of Computer Science , pages 103--112. IEEE, 2010

  7. [15]

    Double/debiased machine learning for treatment and structural parameters, 2018

    Victor Chernozhukov, Denis Chetverikov, Mert Demirer, Esther Duflo, Christian Hansen, Whitney Newey, and James Robins. Double/debiased machine learning for treatment and structural parameters, 2018

  8. [16]

    Rethinking the truncated normal distribution

    Jinho Cha, Byung Rae Cho, and Julia L Sharp. Rethinking the truncated normal distribution. International Journal of Experimental Design and Process Optimisation , 3(4):327--363, 2013

  9. [17]

    High-dimensional robust mean estimation via gradient descent

    Yu Cheng, Ilias Diakonikolas, Rong Ge, and Mahdi Soltanolkotabi. High-dimensional robust mean estimation via gradient descent. In International Conference on Machine Learning , pages 1768--1778. PMLR, 2020

  10. [18]

    What makes a good fisherman? linear regression under self-selection bias

    Yeshwanth Cherapanamjeri, Constantinos Daskalakis, Andrew Ilyas, and Manolis Zampetakis. What makes a good fisherman? linear regression under self-selection bias. arXiv preprint arXiv:2205.03246 , 2022

  11. [19]

    What makes a good fisherman? linear regression under self-selection bias

    Yeshwanth Cherapanamjeri, Constantinos Daskalakis, Andrew Ilyas, and Manolis Zampetakis. What makes a good fisherman? linear regression under self-selection bias. In Proceedings of the 55th Annual ACM Symposium on Theory of Computing , pages 1699--1712, 2023

  12. [20]

    Learning mixtures of structured distributions over discrete domains

    Siu-On Chan, Ilias Diakonikolas, Xiaorui Sun, and Rocco A Servedio. Learning mixtures of structured distributions over discrete domains. In Proceedings of the twenty-fourth annual ACM-SIAM symposium on Discrete algorithms , pages 1380--1394. SIAM, 2013

  13. [21]

    Missing not at random in end of life care studies: multiple imputation and sensitivity analysis on data from the action study

    Giulia Carreras, Guido Miccinesi, Andrew Wilcock, Nancy Preston, Daan Nieboer, Luc Deliens, Mogensm Groenvold, Urska Lunder, Agnes van der Heide, and Michela Baccini. Missing not at random in end of life care studies: multiple imputation and sensitivity analysis on data from t...

  14. [22]

    List decodable mean estimation in nearly linear time

    Yeshwanth Cherapanamjeri, Sidhanth Mohanty, and Morris Yau. List decodable mean estimation in nearly linear time. In 2020 IEEE 61st Annual Symposium on Foundations of Computer Science (FOCS) , pages 141--148. IEEE, 2020

  15. [23]

    De-biased machine learning of global and local parameters using regularized riesz representers

    Victor Chernozhukov, Whitney Newey, and Rahul Singh. De-biased machine learning of global and local parameters using regularized riesz representers. arXiv preprint arXiv:1802.08667 , 2018

  16. [24]

    A simple and general debiased machine learning theorem with finite sample guarantees

    Victor Chernozhukov, Whitney K Newey, and Rahul Singh. A simple and general debiased machine learning theorem with finite sample guarantees. arXiv preprint arXiv:2105.15197 , 2021

  17. [25]

    On the solution of estimating equations for truncated and censored samples from normal populations

    A Clifford Cohen. On the solution of estimating equations for truncated and censored samples from normal populations. Biometrika , 44(1/2):225--236, 1957

  18. [26]

    Truncated and censored samples: theory and applications

    A Clifford Cohen. Truncated and censored samples: theory and applications . CRC press, 1991

  19. [27]

    Learning from untrusted data

    Moses Charikar, Jacob Steinhardt, and Gregory Valiant. Learning from untrusted data. In Proceedings of the 49th Annual ACM SIGACT Symposium on Theory of Computing , pages 47--60, 2017

  20. [28]

    Learning mixtures of gaussians

    Sanjoy Dasgupta. Learning mixtures of gaussians. In 40th Annual Symposium on Foundations of Computer Science (Cat. No. 99CB37039) , pages 634--644. IEEE, 1999

  21. [29]

    Efficient statistics, in high dimensions, from truncated samples

    Constantinos Daskalakis, Themis Gouleakis, Chistos Tzamos, and Manolis Zampetakis. Efficient statistics, in high dimensions, from truncated samples. In 2018 IEEE 59th Annual Symposium on Foundations of Computer Science (FOCS) , pages 639--649. IEEE, 2018

  22. [30]

    Computationally and statistically efficient truncated regression

    Constantinos Daskalakis, Themis Gouleakis, Christos Tzamos, and Manolis Zampetakis. Computationally and statistically efficient truncated regression. In Conference on Learning Theory , pages 955--960. PMLR, 2019

  23. [31]

    Robustly learning any clusterable mixture of gaussians

    Ilias Diakonikolas, Samuel B Hopkins, Daniel Kane, and Sushrut Karmalkar. Robustly learning any clusterable mixture of gaussians. arXiv preprint arXiv:2005.06417 , 2020

  24. [32]

    Simplified estimation from censored normal samples

    Wilfrid J Dixon. Simplified estimation from censored normal samples. The Annals of Mathematical Statistics , pages 385--391, 1960

  25. [33]

    Estimation of parameters of truncated or censored exponential distributions

    Walter L Deemer Jr and David F Votaw Jr. Estimation of parameters of truncated or censored exponential distributions. The Annals of Mathematical Statistics , 26(3):498--504, 1955

  26. [34]

    Recent advances in algorithmic high-dimensional robust statistics

    Ilias Diakonikolas and Daniel M Kane. Recent advances in algorithmic high-dimensional robust statistics. arXiv preprint arXiv:1911.05911 , 2019

  27. [35]

    Robustly learning a gaussian: Getting optimal error, efficiently

    Ilias Diakonikolas, Gautam Kamath, Daniel M Kane, Jerry Li, Ankur Moitra, and Alistair Stewart. Robustly learning a gaussian: Getting optimal error, efficiently. In Proceedings of the Twenty-Ninth Annual ACM-SIAM Symposium on Discrete Algorithms , pages 2683--2702. SIAM, 2018

  28. [36]

    Robust estimators in high-dimensions without the computational intractability

    Ilias Diakonikolas, Gautam Kamath, Daniel Kane, Jerry Li, Ankur Moitra, and Alistair Stewart. Robust estimators in high-dimensions without the computational intractability. SIAM Journal on Computing , 48(2):742--864, 2019

  29. [37]

    Statistical query lower bounds for learning truncated gaussians

    Ilias Diakonikolas, Daniel M Kane, Thanasis Pittas, and Nikos Zarifis. Statistical query lower bounds for learning truncated gaussians. In The Thirty Seventh Annual Conference on Learning Theory , pages 1336--1363. PMLR, 2024

  30. [38]

    A statistical taylor theorem and extrapolation of truncated densities

    Constantinos Daskalakis, Vasilis Kontonis, Christos Tzamos, and Emmanouil Zampetakis. A statistical taylor theorem and extrapolation of truncated densities. In Conference on Learning Theory , pages 1395--1398. PMLR, 2021

  31. [39]

    Detecting low-degree truncation

    Anindya De, Huan Li, Shivam Nadimpalli, and Rocco A Servedio. Detecting low-degree truncation. In Proceedings of the 56th Annual ACM Symposium on Theory of Computing , pages 1027--1038, 2024

  32. [40]

    Maximum likelihood from incomplete data via the em algorithm

    Arthur P Dempster, Nan M Laird, and Donald B Rubin. Maximum likelihood from incomplete data via the em algorithm. Journal of the Royal Statistical Society: Series B (Methodological) , 39(1):1--22, 1977

  33. [41]

    Testing convex truncation

    Anindya De, Shivam Nadimpalli, and Rocco A Servedio. Testing convex truncation. In Proceedings of the 2023 Annual ACM-SIAM Symposium on Discrete Algorithms (SODA) , pages 4050--4082. SIAM, 2023

  34. [42]

    Truncated linear regression in high dimensions

    Constantinos Daskalakis, Dhruv Rohatgi, and Emmanouil Zampetakis. Truncated linear regression in high dimensions. Advances in Neural Information Processing Systems , 33:10338--10347, 2020

  35. [43]

    Efficient truncated linear regression with unknown noise variance

    Constantinos Daskalakis, Patroklos Stefanou, Rui Yao, and Emmanouil Zampetakis. Efficient truncated linear regression with unknown noise variance. Advances in Neural Information Processing Systems , 34, 2021

  36. [44]

    The relative performance of full information maximum likelihood estimation for missing data in structural equation models

    Craig K Enders and Deborah L Bandalos. The relative performance of full information maximum likelihood estimation for missing data in structural equation models. Structural equation modeling , 8(3):430--457, 2001

  37. [45]

    Properties and applications of hh functions

    RA Fisher. Properties and applications of hh functions. Mathematical tables , 1:815--852, 1931

  38. [46]

    Efficient parameter estimation of truncated boolean product distributions

    Dimitris Fotakis, Alkis Kalavasis, and Christos Tzamos. Efficient parameter estimation of truncated boolean product distributions. In Conference on Learning Theory , 2020

  39. [47]

    An examination into the registered speeds of american trotting horses, with remarks on their value as hereditary data

    Francis Galton. An examination into the registered speeds of american trotting horses, with remarks on their value as hereditary data. Proceedings of the Royal Society of London , 62(379-387):310--315, 1898

  40. [48]

    Stochastic relaxation, gibbs distributions, and the bayesian restoration of images

    Stuart Geman and Donald Geman. Stochastic relaxation, gibbs distributions, and the bayesian restoration of images. IEEE Transactions on pattern analysis and machine intelligence , (6):721--741, 1984

  41. [49]

    Learning mixtures of gaussians in high dimensions

    Rong Ge, Qingqing Huang, and Sham M Kakade. Learning mixtures of gaussians in high dimensions. In Proceedings of the forty-seventh annual ACM symposium on Theory of computing , pages 761--770, 2015

  42. [50]

    Sample selection bias as a specification error

    James J Heckman. Sample selection bias as a specification error. Econometrica: Journal of the econometric society , pages 153--161, 1979

  43. [51]

    What to do about missing values in time-series cross-section data

    James Honaker and Gary King. What to do about missing values in time-series cross-section data. American journal of political science , 54(2):561--581, 2010

  44. [52]

    Robust and heavy-tailed mean estimation made simple, via regret minimization

    Sam Hopkins, Jerry Li, and Fred Zhang. Robust and heavy-tailed mean estimation made simple, via regret minimization. Advances in Neural Information Processing Systems , 33:11902--11912, 2020

  45. [53]

    The method of simulated scores for the estimation of ldv models

    Vassilis A Hajivassiliou and Daniel L McFadden. The method of simulated scores for the estimation of ldv models. Econometrica , pages 863--896, 1998

  46. [54]

    Robust mean estimation on highly incomplete data with arbitrary outliers

    Lunjia Hu and Omer Reingold. Robust mean estimation on highly incomplete data with arbitrary outliers. In International Conference on Artificial Intelligence and Statistics , pages 1558--1566. PMLR, 2021

  47. [55]

    Estimation of averages in truncated samples

    Charles N Haas and Peter A Scheff. Estimation of averages in truncated samples. Environmental science & technology , 24(6):912--919, 1990

  48. [56]

    Robust estimation of a location parameter

    Peter J Huber. Robust estimation of a location parameter. In Breakthroughs in statistics , pages 492--518. Springer, 1992

  49. [57]

    Social experimentation, truncated distributions, and efficient estimation

    Jerry A Hausman and David A Wise. Social experimentation, truncated distributions, and efficient estimation. Econometrica: Journal of the Econometric Society , pages 919--938, 1977

  50. [58]

    Robust learning of mixtures of gaussians

    Daniel M Kane. Robust learning of mixtures of gaussians. In Proceedings of the 2021 ACM-SIAM Symposium on Discrete Algorithms (SODA) , pages 1246--1258. SIAM, 2021

  51. [59]

    What to expect of classifiers? reasoning about logistic regression with missing features

    Pasha Khosravi, Yitao Liang, YooJung Choi, and Guy Van den Broeck. What to expect of classifiers? reasoning about logistic regression with missing features. arXiv preprint arXiv:1903.01620 , 2019

  52. [60]

    On the learnability of discrete distributions

    Michael Kearns, Yishay Mansour, Dana Ron, Ronitt Rubinfeld, Robert E Schapire, and Linda Sellie. On the learnability of discrete distributions. In Proceedings of the twenty-sixth annual ACM symposium on Theory of computing , pages 273--282, 1994

  53. [61]

    Reduction of variance for gaussian densities via restriction to convex sets

    Marek Kanter and Harold Proppe. Reduction of variance for gaussian densities via restriction to convex sets. Journal of Multivariate Analysis , 7(1):74--81, 1977

  54. [62]

    Efficient truncated statistics with unknown truncation

    Vasilis Kontonis, Christos Tzamos, and Manolis Zampetakis. Efficient truncated statistics with unknown truncation. In 2019 IEEE 60th Annual Symposium on Foundations of Computer Science (FOCS) , pages 1578--1595. IEEE, 2019

  55. [63]

    The prevention and treatment of missing data in clinical trials

    Roderick J Little, Ralph D'Agostino, Michael L Cohen, Kay Dickersin, Scott S Emerson, John T Farrar, Constantine Frangakis, Joseph W Hogan, Geert Molenberghs, Susan A Murphy, et al. The prevention and treatment of missing data in clinical trials. New England Journal of Medicin...

  56. [64]

    tail" functions; when the

    Alice Lee. Table of the gaussian" tail" functions; when the" tail" is larger than the body. Biometrika , 10(2/3):208--214, 1914

  57. [65]

    A fast spectral algorithm for mean estimation with sub-gaussian rates

    Zhixian Lei, Kyle Luh, Prayaag Venkat, and Fred Zhang. A fast spectral algorithm for mean estimation with sub-gaussian rates. In Conference on Learning Theory , pages 2598--2612. PMLR, 2020

  58. [66]

    On robust mean estimation under coordinate-level corruption

    Zifan Liu, Jong Ho Park, Theodoros Rekatsinas, and Christos Tzamos. On robust mean estimation under coordinate-level corruption. In International Conference on Machine Learning , pages 6914--6924. PMLR, 2021

  59. [67]

    Statistical analysis with missing data , volume 793

    Roderick JA Little and Donald B Rubin. Statistical analysis with missing data , volume 793. John Wiley & Sons, 2019

  60. [68]

    Agnostic estimation of mean and covariance

    Kevin A Lai, Anup B Rao, and Santosh Vempala. Agnostic estimation of mean and covariance. In 2016 IEEE 57th Annual Symposium on Foundations of Computer Science (FOCS) , pages 665--674. IEEE, 2016

  61. [69]

    Limited-dependent and qualitative variables in econometrics

    Gangadharrao S Maddala. Limited-dependent and qualitative variables in econometrics . Number 3. Cambridge university press, 1986

  62. [70]

    Graphical models for inference with missing data

    Karthika Mohan, Judea Pearl, and Jin Tian. Graphical models for inference with missing data. Advances in neural information processing systems , 26, 2013

  63. [71]

    Semiparametric inference for nonmonotone missing-not-at-random data: the no self-censoring model

    Daniel Malinsky, Ilya Shpitser, and Eric J Tchetgen Tchetgen. Semiparametric inference for nonmonotone missing-not-at-random data: the no self-censoring model. Journal of the American Statistical Association , pages 1--9, 2021

  64. [72]

    Full law identification in graphical models of missing data: Completeness results

    Razieh Nabi, Rohit Bhattacharya, and Ilya Shpitser. Full law identification in graphical models of missing data: Completeness results. In International Conference on Machine Learning , pages 7153--7163. PMLR, 2020

  65. [73]

    Parameter estimation for multivariate generalized gaussian distributions

    Fr \'e d \'e ric Pascal, Lionel Bombrun, Jean-Yves Tourneret, and Yannick Berthoumieu. Parameter estimation for multivariate generalized gaussian distributions. IEEE Transactions on Signal Processing , 61(23):5960--5971, 2013

  66. [74]

    On the systematic fitting of curves to observations and measurements

    Karl Pearson. On the systematic fitting of curves to observations and measurements. Biometrika , 1(3):265--303, 1902

  67. [75]

    On the generalised probable error in multiple normal correlation

    Karl Pearson and Alice Lee. On the generalised probable error in multiple normal correlation. Biometrika , 6(1):59--68, 1908

  68. [76]

    Learning from censored and dependent data: The case of linear dynamics

    Orestis Plevrakis. Learning from censored and dependent data: The case of linear dynamics. In Conference on Learning Theory , pages 3771--3787. PMLR, 2021

  69. [77]

    Finite-sample guarantees for high-dimensional dml

    Victor Quintas-Martinez. Finite-sample guarantees for high-dimensional dml. arXiv preprint arXiv:2206.07386 , 2022

  70. [78]

    Holoclean: Holistic data repairs with probabilistic inference

    Theodoros Rekatsinas, Xu Chu, Ihab F Ilyas, and Christopher R \'e . Holoclean: Holistic data repairs with probabilistic inference. arXiv preprint arXiv:1702.00820 , 2017

  71. [79]

    Non-response models for the analysis of non-monotone ignorable missing data

    James M Robins and Richard D Gill. Non-response models for the analysis of non-monotone ignorable missing data. Statistics in medicine , 16(1):39--56, 1997

  72. [80]

    Some thoughts on the distribution of earnings

    Andrew Donald Roy. Some thoughts on the distribution of earnings. Oxford economic papers , 3(2):135--146, 1951

  73. [81]

    Semiparametric regression estimation in the presence of dependent censoring

    Andrea Rotnitzky and James M Robins. Semiparametric regression estimation in the presence of dependent censoring. Biometrika , 82(4):805--820, 1995

  74. [82]

    Analysis of semi-parametric regression models with non-ignorable non-response

    Andrea Rotnitzky and James Robins. Analysis of semi-parametric regression models with non-ignorable non-response. Statistics in medicine , 16(1):81--102, 1997

  75. [83]

    Semiparametric regression for repeated outcomes with nonignorable nonresponse

    Andrea Rotnitzky, James M Robins, and Daniel O Scharfstein. Semiparametric regression for repeated outcomes with nonignorable nonresponse. Journal of the american statistical association , 93(444):1321--1339, 1998

  76. [84]

    Sensitivity analysis for selection bias and unmeasured confounding in missing data and causal inference models

    James M Robins, Andrea Rotnitzky, and Daniel O Scharfstein. Sensitivity analysis for selection bias and unmeasured confounding in missing data and causal inference models. IMA VOLUMES IN MATHEMATICS AND ITS APPLICATIONS , 116:1--94, 2000

  77. [85]

    Robust learning with missing data

    Marco Ramoni and Paola Sebastiani. Robust learning with missing data. Machine Learning , 45:147--170, 2001

  78. [86]

    Inference and missing data

    Donald B Rubin. Inference and missing data. Biometrika , 63(3):581--592, 1976

  79. [87]

    Learning mixtures of arbitrary gaussians

    Arora Sanjeev and Ravi Kannan. Learning mixtures of arbitrary gaussians. In Proceedings of the thirty-third annual ACM symposium on Theory of computing , pages 247--257, 2001

  80. [88]

    Missing data as a causal and probabilistic problem

    Ilya Shpitser, Karthika Mohan, and Judea Pearl. Missing data as a causal and probabilistic problem. Technical report, CALIFORNIA UNIV LOS ANGELES DEPT OF COMPUTER SCIENCE, 2015

  81. [89]

    Adjusting for nonignorable drop-out using semiparametric nonresponse models

    Daniel O Scharfstein, Andrea Rotnitzky, and James M Robins. Adjusting for nonignorable drop-out using semiparametric nonresponse models. Journal of the American Statistical Association , 94(448):1096--1120, 1999

  82. [90]

    Estimation of the mean of a multivariate normal distribution

    Charles M Stein. Estimation of the mean of a multivariate normal distribution. The annals of Statistics , pages 1135--1151, 1981

  83. [91]

    Review of inverse probability weighting for dealing with missing data

    Shaun R Seaman and Ian R White. Review of inverse probability weighting for dealing with missing data. Statistical methods in medical research , 22(3):278--295, 2013

  84. [92]

    Necessary and sufficient conditions for explicit solutions in the multivariate normal estimation problem for patterned means and covariances

    Ted H Szatrowski. Necessary and sufficient conditions for explicit solutions in the multivariate normal estimation problem for patterned means and covariances. The Annals of Statistics , pages 802--810, 1980

  85. [93]

    Statistical methods for robust inference in causal and missing data models

    Eric Joel Tchetgen Tchetgen. Statistical methods for robust inference in causal and missing data models . Harvard University, 2006

  86. [94]

    Missing value estimation methods for dna microarrays

    Olga Troyanskaya, Michael Cantor, Gavin Sherlock, Pat Brown, Trevor Hastie, Robert Tibshirani, David Botstein, and Russ B Altman. Missing value estimation methods for dna microarrays. Bioinformatics , 17(6):520--525, 2001

  87. [95]

    Estimation of relationships for limited dependent variables

    James Tobin. Estimation of relationships for limited dependent variables. Econometrica: journal of the Econometric Society , pages 24--36, 1958

  88. [96]

    Semiparametric theory and missing data

    Anastasios A Tsiatis. Semiparametric theory and missing data. 2006

  89. [97]

    A survey of sampling from contaminated distributions

    John W Tukey. A survey of sampling from contaminated distributions. Contributions to probability and statistics , pages 448--485, 1960

  90. [98]

    Discrete choice models for nonmonotone nonignorable missing data: Identification and inference

    Eric J Tchetgen Tchetgen, Linbo Wang, and BaoLuo Sun. Discrete choice models for nonmonotone nonignorable missing data: Identification and inference. Statistica Sinica , 28(4):2069, 2018

  91. [99]

    Flexible imputation of missing data

    Stef Van Buuren. Flexible imputation of missing data . CRC press, 2018

  92. [100]

    Bond returns, liquidity, and missing data

    Arthur Warga. Bond returns, liquidity, and missing data. Journal of Financial and Quantitative Analysis , 27(4):605--617, 1992

  93. [101]

    Learning distributions generated by one-layer relu networks

    Shanshan Wu, Alexandros G Dimakis, and Sujay Sanghavi. Learning distributions generated by one-layer relu networks. Advances in neural information processing systems , 32, 2019

  94. [102]

    Inverse probability weighted estimation for general missing data problems

    Jeffrey M Wooldridge. Inverse probability weighted estimation for general missing data problems. Journal of econometrics , 141(2):1281--1301, 2007

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.