Pith. sign in

REVIEW 3 major objections 6 minor 18 references

A mixture distribution approach for assessing genetic impact from twin study

T0 review · 3 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Pooling monozygotic and dizygotic twins into one likelihood with a shared variance makes dizygotic correlation estimation root-n consistent even for homogeneous traits, so the genetic impact $\delta=\rho_M-\rho_D$ can be inferred with…

desk verdict Root-n result for rho_D is interesting and likely right, but the equal-variance assumption is load-bearing and the proof lives in an unavailable supplement. read the letter →

arxiv 2507.13605 v1 pith:VQ22N32G submitted 2025-07-18 stat.ME stat.AP

Zonghui Hu , Pengfei Li , Dean Follmann , Jing Qin This is my paper · ORCID
classification stat.MEstat.AP MSC 62F1262F0362H2062P10
keywords correlationcoefficientmaximumlikelihoodestimationmixturedistributionunorderedpairstwinstudygeneticimpactroot-nconsistencyhomogeneitytest
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 addresses a classic blind spot in twin studies: the two measurements in a twin pair arrive without a meaningful label, and standard correlation measures either ignore the order or are biased by it. The authors model the unordered pair as a two-component mixture of bivariate normals and fit it jointly to monozygotic and dizygotic twins, assuming the two zygosity groups share a variance. The central result is that this combined likelihood estimates the dizygotic correlation $\rho_D$ at root-n speed even in the homogeneous case where the dizygotic means are equal—the case that makes the usual separate mixture estimator slow down to $n^{-1/4}$. Consequently the genetic impact $\delta=\rho_M-\rho_D$, Falconer's heritability-style contrast, is consistently estimable with bootstrap confidence intervals. If true, this gives a practical route to assessing genetic control of immune traits and any other unordered paired data.

What carries the argument

The load-bearing object is the mixture density $f(y_1,y_2;\mu_1,\mu_2,\rho,\sigma^2)=0.5\phi(y_1,y_2;\mu_1,\mu_2,\rho,\sigma^2)+0.5\phi(y_1,y_2;\mu_2,\mu_1,\rho,\sigma^2)$, which says the observed pair is the labeled pair in either of two equally likely orders. The combined log-likelihood is $l_C=l_M+l_D$, where $l_M$ is the ordinary bivariate normal log-likelihood for MZ pairs with $\mu_1=\mu_2$ and $l_D$ is the mixture log-likelihood for DZ pairs, with $\sigma^2$ shared across zygosity. Sharing the variance is what does the work: the regular MZ part estimates $\sigma^2$ and $\rho_M$ at the usual $n^{-1/2}$ rate, and because $\sigma^2$ is then effectively known inside the DZ mixture, the poorly identified mean labels $\mu_{D1}$ and $\mu_{D2}$ may drift at $n^{-1/4}$ without dragging $\rho_D$ and $\sigma^2$ below $n^{-1/2}$. The likelihood-ratio statistic for homogeneity then inherits the boundary-type mixture limit $0.5\chi_0^2+0.5\chi_1^2$, with a small-sample adjustment $a_n=0.5+6.828/n_M$.

What would settle it

Simulate homogeneous DZ data with $\mu_{D1}=\mu_{D2}$ but $\sigma_D^2\neq\sigma_M^2$, fit the combined likelihood, and track the empirical standard deviation of $\hat\rho_D$ as $n$ grows; if it does not shrink like $n^{-1/2}$, the equal-variance condition is essential. Equivalently, compute a likelihood-ratio test of $\sigma_M^2=\sigma_D^2$ on the immune-trait data and check whether bootstrap coverage of $\delta$ degrades when the test rejects.

Watch

Extended reading notes

Core claim

The paper's central claim is Theorem 2: under the combined likelihood $l_C$ built from MZ and DZ data with a common variance $\sigma^2$, the maximum likelihood estimators satisfy $\hat\rho_M-\rho_{M0}=O_p(n^{-1/2})$ and $\hat\rho_D-\rho_{D0}=O_p(n^{-1/2})$ even when the DZ means coincide, i.e. $\mu_{D1}=\mu_{D2}=\mu_{D0}$. In that homogeneous case the separate DZ mixture alone only gives $\hat\rho_D-\rho_{D0}=O_p(n_D^{-1/4})$ (Theorem 1). The combined estimator therefore keeps the genetic impact $\delta=\rho_M-\rho_D$ at root-n consistency no matter whether the trait is homogeneous or heterogeneous, and simulations show $\hat\delta$ is nearly unbiased and symmetrically distributed, with bootstrap intervals near nominal coverage. The same likelihood yields a homogeneity test for $\mu_{D1}=\mu_{D2}$ whose null limiting distribution is $0.5\chi_0^2+0.5\chi_1^2$, with a simulation-calibrated small-sample adjustment.

Load-bearing premise

The load-bearing premise is that MZ and DZ twins have the same variance, because the MZ data anchor the variance in the DZ mixture and the claimed root-n rate for $\rho_D$ collapses if that equality fails; the paper checks this assumption only informally in Section 5, with no formal test.

Editorial extensions

If this is right

  • If the combined-mixture result holds, genetic impact $\delta=\rho_M-\rho_D$ can be estimated with ordinary bootstrap inference for every trait, because $\hat\delta$ is root-n consistent whether or not the trait is homogeneous.
  • Pearson's correlation should not be used for DZ twins when a trait may be heterogeneous: the paper's formula shows it underestimates $\rho_D$ by an amount growing with $(\mu_{D1}-\mu_{D2})^2$, so it overstates $\delta$, while the mixture estimator removes that bias.
  • The homogeneity test with limiting distribution $0.5\chi_0^2+0.5\chi_1^2$ and the small-sample adjustment $a_n=0.5+6.828/n_M$ provides a practical check for whether DZ means differ.
  • The method transfers to any unordered paired data, not only twins, because the mixture only assumes random labeling and equal variance for the groups being compared.

Reading between the lines

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

  • The equal-variance assumption $\sigma_M^2=\sigma_D^2$ is checked only informally in Section 5; a natural extension would be a sensitivity analysis or a formal likelihood-ratio test, and the root-n claim for $\rho_D$ would need re-examination if that equality fails.
  • The empirical correction $a_n=0.5+6.828/n_M$ is fit by simulation with $n_M$ as the only covariate; an analytic higher-order expansion could replace it and would clarify how the slow $n^{-1/4}$ mean estimators affect small-sample calibration.
  • In the immune-trait application, traits flagged as possibly heterogeneous (B cells, CD4, CD8) are exactly those where Pearson's $\delta$ exceeds the mixture $\delta$; this suggests a diagnostic rule that a large gap between Pearson and mixture estimates is itself evidence of order-labeling bias.
  • For other latent-label bivariate data, such as allele-specific expression or two-channel assays without strand assignment, the same combined-likelihood trick may recover root-n correlation estimation even when component means are equal, though the paper does not test these settings.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 6 minor

Summary. The manuscript proposes a bivariate mixture model for twin data in which the ordering within each pair is unknown. For monozygotic (MZ) twins the model reduces to a bivariate normal with equal means; for dizygotic (DZ) twins it is a 50-50 mixture of two bivariate normals with swapped means. The authors study maximum likelihood estimation from the separate MZ/DZ likelihoods and from a combined likelihood that imposes a common variance σ². Theorem 2 claims that under the combined likelihood, ρD and ρM are root-n consistent even when μD1 = μD2, while the separate mixture has slower rates (Theorem 1). This motivates inference on the genetic impact δ = ρM − ρD. A likelihood ratio homogeneity test for μD1 = μD2 is proposed (Theorem 3). Simulations and an application to six immune traits are presented.

Significance. If the central claim is correct, the combined-likelihood approach is a practical advance: it gives standard-rates inference for the genetic impact without genotype data, in a setting (unordered pairs with potential heterogeneity) where Pearson correlations are biased. The paper is honest about the difficulty of the homogeneity test and supplies a data application. The simulations support the claimed rates for ρD and δ under the assumed equal-variance model. However, the significance is currently conditional: the main theorems are not verifiable from the submitted text (proofs are in a missing supplement), Theorem 1's rates appear inconsistent with standard local expansions, and the equal-variance assumption on which the root-n result for ρD rests is only informally checked. These points must be resolved before the practical claims are accepted.

major comments (3)
  1. [Section 3, Theorems 1–3] All three theorems are stated without proof in the main text; the reader is referred to a supplementary document that is not part of the submitted material. Since Theorem 2 is the central claim and Theorem 1 contains non-standard rates, the manuscript cannot be evaluated for correctness as presented. Please provide the complete proofs (ideally in the main text or an accessible supplement) and state the regularity conditions and the parameterization used.
  2. [Section 2, Eq. (3) and Section 5] The combined likelihood (3) is constructed under the equal-variance assumption σ_M² = σ_D² = σ², and the heuristic in Section 3.2 explicitly uses the MZ likelihood to anchor σ² for the DZ mixture. If the true variances differ, the DZ likelihood is misspecified and the claimed root-n consistency of ρD, and hence of δ, is not guaranteed; the intervals in Table 4 would then be unreliable. Section 5 checks this assumption only by noting that separate variance estimates are "quite close", with no test or numerical evidence. Please add a formal or simulation-based assessment of sensitivity to variance inequality (e.g., coverage of bootstrap intervals for δ under σ_D² ≠ σ_M²), or modify the method to be robust to this assumption.
  3. [Section 3.1, Theorem 1] The rates in Theorem 1 are surprising and need justification. With α = 0.5 fixed, the DZ mixture density at μD1 = μD2 = μD0 depends on (μD1 − μD2) only through quadratic and higher powers, so the standard local expansion gives bμD1 − bμD2 = Op(n^{−1/4}) for the difference and root-n for the mean (μD1 + μD2)/2, while ρD and σ_D² should be root-n since the model reduces to a single bivariate normal. The stated rates n^{−1/8} for the individual means and n^{−1/4} for ρD and σ_D² are not reconciled with this heuristic or with the separate-likelihood simulations in Table 1 (which do show sd ratios consistent with n^{−1/4} for ρD). Please clarify whether α is actually estimated, or provide the proof and explain the n^{−1/8} rate.
minor comments (6)
  1. [Introduction] The study is described as including 324 twins, with 96 MZ and 288 DZ twins; these numbers sum to 384. Please correct the total and verify the sample sizes used in Table 4.
  2. [Section 2] The displayed formula for Pearson's correlation coefficient is garbled: the second factor in the numerator appears to be missing Y_{i2}. Please correct the formula so that the convergence result (2) is connected to it.
  3. [Section 3.1] In Theorem 1, "bσ2 D2" appears to be a typo for "bσ2_D".
  4. [Figure 1] The caption contains misspellings ("qunatile") and the phrase "both both"; please fix these presentation issues.
  5. [Section 3.3] The adjusted limiting distribution relies on the empirically fitted constant a = 6.828 in a_n = 0.5 + 6.828/nM. Since this is a fitted quantity, please report its precision and show type I error simulations of the homogeneity test at the adjusted level, not only QQ plots.
  6. [Table 4] Please state the sample sizes used in the application and describe the bootstrap procedure (number of resamples, method for the CI).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: root-n consistency is a derived theorem and the fitted correction constant is an explicit empirical adjustment, not a prediction.

full rationale

The paper's central claim (Theorem 2) is a stated asymptotic result with proof deferred to the supplementary document; the main text offers only a heuristic explanation (Section 3.2) that the MZ likelihood makes sigma^2 'effectively known' for the DZ mixture. This depends on the equal-variance assumption sigma_M^2 = sigma_D^2 in equation (3), but that is a modeling condition, not a definitional equivalence. The combined likelihood is constructed from separate MZ and DZ likelihood pieces, and the claimed root-n rate for rho_D is not imposed by construction; it is asserted as a theorem. The only fitted quantity in the paper is the Bartlett-type correction constant a = 6.828 in Section 4, explicitly labeled an 'empirical formula' for the homogeneity test's finite-sample adjustment; it is not used to establish the central asymptotic claim and is not presented as a prediction. No self-citations are load-bearing; the cited prior work on mixture testing is external and contextual. The paper itself flags the equal-variance check as informal in Section 5 ('the limiting distribution of the test statistic is difficult to derive and beyond the scope of this work'), which is a robustness limitation rather than circularity. The deferred proofs and the empirical correction are acknowledged limitations, but no equation or fitted parameter is recycled as its own output. Therefore no significant circularity is present, and the correct score is 0.

Assumptions & free parameters 1 free parameters · 4 assumptions · 0 invented entities

The central claim rests on normality of the underlying ordered pairs, equal MZ/DZ variance, and random labeling. The first is standard for twin trait modeling; the second is load-bearing for the combined likelihood but only informally checked; the third is plausible by exchangeability. No invented entities. The only fitted quantity is the homogeneity test correction constant.

free parameters (1)
  • a (empirical correction constant) = 6.828
    Empirical constant in the adjustment a_n = 0.5 + a/n_M for the homogeneity test's null distribution; estimated from 10,000 simulated datasets in Table 3. This is a finite-sample correction, not part of the central root-n result.
assumptions (4)
  • domain assumption Underlying ordered pair (Y1*,Y2*) follows a bivariate normal distribution
    Used throughout Section 2 to justify the mixture model and formula (2).
  • domain assumption Equal variance in MZ and DZ twins, sigma^2_M = sigma^2_D
    Stated in Section 2 before equation (3); the combined likelihood relies on it. Informal check in Section 5.
  • domain assumption Labeling is random, so mixture weight alpha = 0.5
    Section 2: 'it is safe to fix alpha = 0.5 as labeling is random.'
  • domain assumption MZ twins share identical mean, mu_1 = mu_2
    Section 2: 'For MZ twins, mu1 = mu2.'

how reviews work

0 comments
Cite this review

Pith. "Pith review of A mixture distribution approach for assessing genetic impact from twin study." pith.science (2026). https://pith.science/paper/VQ22N32G

@misc{pith2026250713605,
  author       = {Pith},
  title        = {Pith review of: A mixture distribution approach for assessing genetic impact from twin study},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VQ22N32G}},
  note         = {Machine review of arXiv:2507.13605}
}
read the original abstract

This work was motivated by a twin study with the goal of assessing the genetic control of immune traits. We propose a mixture bivariate distribution to model twin data where the underlying order within a pair is unclear. Though estimation from mixture distribution is usually subject to low convergence rate, the combined likelihood, which is constructed over monozygotic and dizygotic twins combined, reaches root-n consistency and allows effective statistical inference on the genetic impact. The method is applicable to general unordered pairs.

Figures

Figures reproduced from arXiv: 2507.13605 by the authors.

Figure 1
Figure 1. Quantile-quantile plot of RC n under H0: the unadjusted limiting distribution on the top, the adjusted limiting distribution on the bottom, nM = nD = 100 on the left and 400 on the right, at ρM = 0.9 and ρD = 0.3. with Φ the cumulative distribution function of the standard normal distribution, works well. For MZ twin data (Yi1, Yi2), the order within a pair does not matter and it has the same distribution as the und… view at source ↗
Figure 2
Figure 2. Homogeneity test results: the observed log-likelihood ratios for the six immune [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

18 extracted references · 16 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence := #2 '...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

  3. [3]

    , Joosten, I

    Aguirre-Gamboa, R. , Joosten, I. , Urbano, P. C. , van der Molen, R. G. , van Rijssen, E. , van Cranenbroek, B. , Oosting, M. , Smeekens, S. , Jaeger, M. , Zorro, M. et al. (2016). Differential effects of environmental and genetic factors on t and b cell immune traits. Cell reports 17, 2474--2487

    work page 2016

  4. [4]

    , Hall, M

    Ahmadi, K. , Hall, M. , Norman, P. , Vaughan, R. , Snieder, H. , Spector, T. & Lanchbury, J. (2001). Genetic determinism in the relationship between human cd4+ and cd8+ t lymphocyte populations? Genes and Immunity 2, 381--387

    work page 2001

  5. [5]

    Bartlett, M. S. (1937). Properties of sufficiency and statistical tests. Proceedings of The Royal Society A 160, 268--282

    work page 1937

  6. [6]

    , Jojic, V

    Brodin, P. , Jojic, V. , Gao, T. , Bhattacharya, S. , Angel, C. J. L. , Furman, D. , Shen-Orr, S. , Dekker, C. L. , Swan, G. E. , Butte, A. J. et al. (2015). Variation in the human immune system is largely driven by non-heritable influences. Cell 160, 37--47

    work page 2015

  7. [7]

    & Sun, J

    Charnigo, R. & Sun, J. (2004). Testing homogeneity in a mixture distribution via the L_2 -distance between competing models. Journal of the American Statistical Association 99, 488--498

    work page 2004

  8. [8]

    & Chen, J

    Chen, H. & Chen, J. (2001). The likelihood ratio test for homogeneity in finite mixture models. The Canadian Journal of Statistics 29, 201--215

    work page 2001

Show all 18 references
  1. [9]

    & Gassiat, E

    Dacunha-Castelle, D. & Gassiat, E. (1999). Testing the order of a model using locally conic parametrization: Population mixtures and stationary ARMA processes. The Annals of Statistics 27, 1178--1209

    1999

  2. [10]

    & Phillips, A

    Davies, P. & Phillips, A. (1988). Nonparametric tests of population differences and estimation of the probability of misidentification with unidentified paired data. Biometrika 75, 753--760

    1988

  3. [11]

    Ernst, M. D. , Guerra, R. & Schucany, W. R. (1996). Scatterplots for unordered pairs. The American Statistician 50, 260--265

    1996

  4. [12]

    Frankham, R. (1996). Introduction to quantitative genetics (4th edn): by D ouglas S . F alconer and T rudy FC M ackay L ongman. Trends in Genetics 12, 280

    1996

  5. [13]

    Hinkley, D. (1973). Two-sample tests with unordered pairs. Journal of the Royal Statistical Society. Series B (Methodological) , 337--346

    1973

  6. [14]

    , Chen, J

    Li, S. , Chen, J. , Guo, J. , Jing, B.-Y. , Tsang, S.-Y. & Xue, H. (2015). Likelihood ratio test for multi-sample mixture model and its application to genetic imprinting. Journal of the American Statistical Association 110, 867--877

    2015

  7. [15]

    & Shao, Y

    Liu, X. & Shao, Y. (2003). Asymptotics for likelihood ratio tests under loss of identifiability. The Annals of Statistics 31, 807--832

    2003

  8. [16]

    , Quaye, L

    Roederer, M. , Quaye, L. , Mangino, M. , Beddall, M. H. , Mahnke, Y. , Chattopadhyay, P. , Tosi, I. , Napolitano, L. , Barberio, M. T. , Menni, C. et al. (2015). The genetic architecture of the human immune system: a bioresource for autoimmunity and disease pathogenesis. Cell ...

    2015

  9. [17]

    Sedgwick, P. (2013). Intraclass correlation coefficient. BMJ: British Medical Journal 346

    2013

  10. [18]

    Serfling, R. J. (2000). Approximation Theorems of Mathematical Statistics. Wiley, New York

    2000

Pith tools

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