Case–control study 2026-10-14
A case–control study selects subjects according to an outcome: cases have the outcome and controls represent an appropriate comparison population. It compares their prior exposures or genetic variants to study association. Its retrospective sampling must be accounted for when interpreting risks; selection bias and confounding can distort an unadjusted association. In genetics, matched family pseudo-control genotypes can replace unrelated controls for a conditional transmission disequilibrium test.
Convenience sampling 2026-10-25
Recruitment based mainly on accessibility rather than known probability selection from a defined target population. A larger convenience sample can reduce some sampling variability without eliminating selection bias or making respondents relevant to the intended construct.
Hypercalcaemia 2026-10-05
Hypercalcaemia is an abnormally elevated serum calcium concentration. A diagnosis in medical records depends on testing and diagnostic practice, which can introduce selection bias into an observational study.
Longitudinal attrition 2026-10-25
Loss of participants or outcome observations across follow-up. If remaining under observation depends on exposure and determinants of the outcome, analyzing only retained participants can induce selection bias.
Multiplicity bias from a detector trigger 2026-10-18
A trigger changes which particle populations are common in the recorded sample. Let be the number of particles in an event, and stipulate an independent Bernoulli distribution for each particle's ability to produce an accepted telescope signal, with success probability . If any successful particle suffices, the event acceptance isThis is a deliberately simplified model. Its includes the joint geometrical acceptance and response of the whole telescope; it is not an assumed universal single-counter efficiency. Actual shower directions, spatial correlations and shared electronics can violate the independence assumption.
The conditional probability formula givesFor chosen underlying fractions , and , the two acceptances are and . Five-particle events consequently form about of recorded events. A population that is only of incoming events can occupy nearly of the photographs. This is selection bias, not evidence that the source population changed.
A coincidence window trades timing acceptance against background, while a trigger enriches events with more particles
. Original calculations with the explicitly chosen rates, Gaussian timing spread and independent-particle acceptance model in the text. No historical event counts or measured timing curves are reproduced. Details. Observational study Created 2026-09-28 Updated 2026-10-25
An observational study measures exposures and outcomes without assigning the exposure by an intervention. Its causal interpretation therefore requires assumptions that address confounding and selection bias.
Conditional on the printed two-group test output, the male-coded group has higher mean efficacy by about 1.034 scale points, with a 95% confidence interval of about 0.785 to 1.282 points. This is strong evidence against equal population means under the test's assumptions. It does not establish a causal effect, and the accompanying overall mean is inconsistent with the group summaries.
Here the coding identifies group 0 as male and group 1 as female in this dataset. The counts are 1215 and 1251, respectively, or approximately 49.27% and 50.73% of 2466 observations. These are sample compositions, not evidence that the sample represents every population of interest. The continuous efficacy scale runs from 0 to 16, with higher scores indicating stronger efficacy. Its printed sample median is 10 and its quartiles are 8 and 12, so its interquartile range is 4 points. The printed overall sample mean, 9.348, requires the separate consistency qualification below.
Let and denote the two population means. The null hypothesis is , and the alternative hypothesis is the two-sided test . The two sample means giveThe positive sign means higher mean efficacy in group 0, given this coding; it does not say that every member of that group has higher efficacy than every member of group 1.
Welch's t-test estimates uncertainty from the two groups separately rather than assuming equal population variances:The printed implies a standard error of approximately points. The noninteger 2277.4 statistical degrees of freedom comes from the Welch approximation; it is not a count of fractional observations. Using the corresponding Student t-distribution gives the displayed p-value, approximately . The test and interval calculations agree to the precision of the printed numbers.[ref]
Student's t-Test in the R documentation. Arguments and Details: two-sided default and unequal-variance Welch calculation; Value: estimates, statistic, degrees of freedom and interval
A p-value is the probability, under the null hypothesis and the assumed sampling model, of a test result at least as incompatible with equal means as the observed result. It is not the probability that the null hypothesis is true, the probability that the finding is a chance mistake, or the size of the mean difference. At a prespecified significance level such as , this result rejects the equal-mean null hypothesis; it also rejects at 1%. A Type I error is rejection when that hypothesis is true, and a valid level- procedure controls its long-run probability. A Type II error is failure to reject a false hypothesis; statistical power is the probability of rejection under a specified alternative. Neither error probability for this one realized claim is supplied by its p-value.
The 95% confidence interval comes from , with critical value about 1.961. Its positive endpoints support the same direction and exclude zero. The frequentist 95% statement concerns coverage over repeated samples, not a 95% posterior probability assigned to this fixed interval. Reporting the difference in scale points supplies an effect size; its practical importance also depends on what a point means substantively. A very small p-value alone does not answer that question.
The inferential calculation requires independent observations and an appropriate sampling model. With these large groups, approximately normal sampling distributions of the means can be more plausible than normality of individual bounded scores, but a large sample size does not correct dependence, clustering, selection bias or measurement problems. The output does not establish random assignment or remove possible confounding, so the estimated difference is an association.
There is a material source inconsistency. For exhaustive groups with complete outcome data, the sample mean must be the sample-size-weighted mean of the group means:which cannot round to 9.348. Figure 7 independently reports all 2466 efficacy values present and again gives 9.348; Question 6 refers back to the same continuous scale. The pooled-mean compatibility under rounding rules out explaining the discrepancy by the displayed decimal precision. Thus the two-group test supports the stated conclusion if its own group summaries and test fields are trusted, while the overall mean cannot simultaneously be accepted as a coherent summary of those groups. No unprinted exclusion or uniquely correct replacement mean is established.
The fitted age slope is very small and negative, about efficacy points per year, with no conventional evidence against a zero linear slope: , . The model explains only about 0.0367% of the observed efficacy variation. These results concern a linear association, not proof of no possible age relationship or a causal effect of ageing.
The descriptive statistics report 2466 observed ages and efficacy scores, with no missing values for either variable. Age ranges from 14 to 97 years, with sample mean 52.689 years and sample standard deviation 19.064 years. Efficacy ranges from 0 to 16, with sample mean 9.348 and sample standard deviation 3.176 scale points. The standard deviations use different units and should not be compared as though their numerical ratio measured relative substantive variability. The efficacy mean retains the conflict with Question 1's group summaries; that conflict is not resolved by this separate regression output.
A simple linear regression models , where is efficacy and is age in years. Ordinary least squares chooses the coefficients to minimize the residual sum of squares. From the printed coefficients the fitted line isHolding to this model, ten additional years correspond to fitted scale points. For example, the fitted mean is about 9.453 at age 20 and 9.261 at age 80, a difference of about points. These are fitted means, not exact predictions of individual scores. The line does not imply that every person's efficacy falls at that rate, nor that cross-sectional age differences identify within-person changes over time.
The intercept is the predicted mean at age zero, about 9.517. Since zero is outside the observed 14–97 age range, its direct interpretation requires extrapolation. The printed standard error of 0.188043 and yield the extremely small displayed p-value for a zero intercept. That test does not test whether age is a useful predictor. Recentring age at its observed sample mean would put the intercept near , the fitted mean at a substantively relevant reference age. This agrees with the rounded overall efficacy mean.
The age standard error is 0.003356 points per year. The coefficient divided by this standard error gives , agreeing with the printed t value. Estimating an intercept and one regression slope from 2466 observations leaves residual statistical degrees of freedom. Under a correctly specified normal linear model, the Student t-distribution provides the stated two-sided p-value. At 5%, or 1%, it does not reject . The p-value is not a 34.2% probability that age has no effect, and nonrejection does not establish a practically negligible effect size.[ref]
Summarizing Linear Model Fits in the R documentation. Value: coefficient estimate, standard error, two-sided t and p columns; residual degrees of freedom
A 95% confidence interval calculated from the printed estimate and standard error is approximatelyor about points per decade. It includes zero and gives the scale of slopes compatible with these summaries and assumptions. Assessing a practically negligible effect would require an explicit substantive tolerance and an appropriate equivalence test, not simply the fact that . A nonlinear association could also coexist with a near-zero overall linear slope.
The coefficient of determination is , meaning that the fitted line explains about of the sample's centred efficacy sum of squares compared with an intercept-only fit. This is not 36.7%, a fraction of people whose efficacy is predicted correctly, or the proportion of efficacy caused by age. The adjusted coefficient of determination, approximately , is slightly negative because its penalty for estimating the extra coefficient exceeds the tiny improvement in fit. It need not be constrained to be nonnegative and does not mean a negative amount of explained variation.
The overall on 1 and 2464 statistical degrees of freedom, with , tests the same zero-slope hypothesis as the age t test because there is only one predictor. The simple-regression output identities give and ; the printed numbers agree allowing their displayed precision. The signed Pearson correlation coefficient is consequently about . These are related summaries of the same weak linear association, not independent corroborating tests.
The residual standard error of 3.177 points estimates the typical unexplained error scale under the fitted model. It is close to the efficacy sample standard deviation, as expected when the line explains almost none of the variation. It is not a standard error of the slope or a guarantee that each individual's prediction error lies within 3.177 points. The residual summary spans roughly to , with sample median about 0.573 and central quartiles about and . An ordinary least-squares fit with an intercept has zero mean residual, not necessarily zero median residual. These few quantiles cannot establish normality, constant conditional spread, absence of outliers or linearity. The apparent last-decimal discrepancies among rounded error scales and other summary fields do not force an additional contradiction: their rounding intervals admit compatible regression summaries.
The conventional uncertainty calculations assume independent, mean-zero errors conditional on age, an appropriate linear mean, and homoscedasticity; normal errors give the usual exact finite-sample t and F calibration. Large samples can relax some distributional requirements but do not remove confounding, selection bias or dependence. A residual-versus-fitted plot could reveal curvature or changing spread, and a normal probability plot could assess error shape; neither plot can be reconstructed from the printed residual quantiles. A prediction interval for an individual includes residual variation and is wider than a confidence interval for the mean response. The supplied output supports a very weak, imprecisely signed linear association with little predictive gain from age alone, while leaving nonlinear, causal and design-based questions unresolved.
Taken at face value, the printed two-sided p-value rejects equality of average political efficacy at the 5% level, with the non-employed group scoring higher. The effect magnitude and strength of evidence cannot be reported as one coherent result because the estimates, interval, test statistic and p-value contradict one another.
Let and be the expected political efficacy scores in the two employment groups in the target population. The null hypothesis is and the two-sided alternative hypothesis is . A Welch's t-test compares the difference of sample means with its standard error, estimating the two group variances separately:The reference Student t-distribution has approximate, generally noninteger statistical degrees of freedom; it does not require equality of the two variances. Observations should be independent, with finite group variances; normal distributions give the usual small-sample rationale, while large samples permit an approximate normal sampling distribution under suitable conditions. A large sample cannot correct a biased sampling procedure or dependence between observations.[ref]
Two-Sample t-Test for Equal Means by National Institute of Standards and Technology. §1.3.5.3, unequal-variance statistic and Welch–Satterthwaite degrees of freedom
A p-value is a probability calculated under the null hypothesis and the sampling assumptions: for this two-sided test it measures the probability of an absolute test statistic at least as large as the observed one. It is not the probability that the null hypothesis is true, the probability that chance caused the observed difference, or a measure of the effect size. If the significance level was chosen in advance, the printed falls just below it. This would reject at 5%, but not at 1%. The significance level describes the long-run probability of a type I error, rejecting a true null hypothesis; a type II error instead fails to reject a false null hypothesis. Statistical power is the rejection probability under a specified alternative hypothesis. None of these errors for this particular realised study can be read as a posterior probability from .
Both the printed sample means and the positive confidence interval point toward greater political efficacy in group 0, the non-employed group. Their magnitudes disagree. Subtracting the printed sample means gives scale units. Yet the printed 95% confidence interval has midpoint , not . A two-sided Welch interval is centred on the difference of sample means. This violates the consistency identities for a two-sided t-test; the interval and estimates cannot describe the same reported comparison. If the interval alone were trusted, it would exclude zero and indicate an increase between about and units, but that is explicitly a conditional interpretation. A 95% confidence interval refers to repeated-sample coverage of a fixed statistical parameter, not a 95% posterior probability for that realised interval.[ref]
There are additional arithmetic warnings. The printed at statistical degrees of freedom gives a two-sided p-value about , not . Thus taking the test statistic as authoritative would imply much stronger evidence, not merely a marginal 5% result. The employment counts sum to , but their percentages are and , not and . If those counts also apply to the two political efficacy sample means, their weighted arithmetic mean would be about , not the reported overall ; separate complete-case group sizes are not supplied, so that last cross-output comparison is conditional on a common sample. No uniquely corrected difference, standard error or interval follows from these contradictory fields.
The political efficacy summary describes a bounded observed scale from to , with median and middle quartiles and . The interquartile range is ; the overall sample mean being below the median does not, by itself, establish a particular distributional shape. These descriptive statistics cannot supply the missing coherent group variances or a standardized effect size. Even a corrected small p-value would establish an association under a statistical model, rather than a causal effect of employment. Age, other background covariates, selection bias and the quality of the political efficacy measurement remain relevant. A defensible report gives the conditional 5% decision and direction, identifies the contradictions, and withholds a single numerical effect-and-uncertainty claim until the output is checked.
Let and be the population means of the newspaper trust score for the unvaccinated and vaccinated groups. A two-sided hypothesis test comparesThe observed difference of sample means is : the vaccinated group has a one-unit higher average trust score. The negative test statistic agrees with subtraction in this order. A Student's t-test scales the estimated difference by its estimated standard error; means that the observed difference is 3.6 estimated standard errors below the equal-mean null hypothesis.
Using the reported reference Student t-distribution with 747 statistical degrees of freedom,Thus the equal-mean null hypothesis is rejected at a prespecified 5%, 1% or 0.1% significance level, subject to the test's assumptions. A p-value is a probability of a result at least this extreme under the specified null hypothesis and sampling model. It is neither the probability that the null hypothesis is true nor the probability that this particular conclusion is mistaken. A chosen significance level controls the procedure's type I error probability under that model.
The estimated difference is more informative than statistical significance alone. If the reported statistic uses this one-unit contrast, its estimated standard error is . The corresponding approximate 95% confidence interval isThis interval concerns a population mean difference under repeated sampling; it is not a range containing 95% of individual trust scores. Its practical interpretation also depends on the trust scale and the importance of a one-unit difference.
The groups are disjoint, so an independent-samples comparison is appropriate if the sampled observations are independent. The exact pooled Student's t-test additionally assumes common population variance and normal distributions; Welch's t-test permits unequal variances and uses an approximate reference distribution. A bounded or ordinal rating need not be exactly normal, so treating it numerically and using a large-sample approximation also requires judgment about group sizes and the score's meaning.
There is a bookkeeping qualification. The later outputs use 748 observations. If these are also the complete observations for this two-group comparison, the pooled test has statistical degrees of freedom, and the Welch degrees of freedom upper bound is also 746. The printed 747 therefore needs the original test call or a different sample count to reconcile it. The reported-calibration calculation above answers the stated test; it does not invent group sizes or silently change the reference distribution.
The evidence supports an association between higher newspaper trust and recorded vaccination, rather than a causal effect of trust. An observational study does not rule out confounding, selection bias or reverse causation. A statistical hypothesis test cannot supply those missing causal assumptions.
Several mechanisms could explain an observed association without a causal effect of power on friendliness. No exposure was assigned, so this is an observational study. Age, organizational role, wealth, gender, industry, public prominence or other common causes could produce confounding. Reverse direction is also possible: interpersonal behavior might influence advancement and hence power. Adjusting for age alone would not remove every alternative explanation.
The measured outcome is perceived friendliness from a photograph. Facial expression, pose, lighting, clothing, attractiveness and familiarity could affect those ratings independently of everyday behavior. If raters recognize an executive or know their status, expectations could alter the outcome; financial experts' power ratings may themselves involve judgment and measurement error. These are construct validity issues, since a photographic impression need not measure friendliness toward colleagues.
The preselected thirty executives are not a probability sample of all workers or even necessarily of all executives. Selection bias, narrow representation and restriction of range could alter associations and limit external validity. Multiple ratings of the same executive are also not independent executive observations; ignoring that structure would be pseudoreplication. These possibilities apply whether the reported association is positive, negative or inconclusive; the unavailable numerical findings are not assumed.
Bootstrapping estimates sampling uncertainty by repeating an analysis on resampled data. In the ordinary nonparametric bootstrap, draw observations with replacement from the empirical sample, retaining its size, and recalculate the statistic many times. The resulting distribution can estimate a standard error, assess bias or form a justified bootstrap confidence interval. It is useful when an analytic sampling distribution is difficult to obtain.
The resampling unit must match the design: keep paired measurements together, and resample independent clusters when observations are clustered. Resampling individual ratings as if they were independent lecturers would misrepresent uncertainty. The ordinary bootstrap assumes the sample and chosen resampling scheme adequately represent the relevant data-generating variation; it does not repair selection bias, guarantee accurate inference from a tiny sample, or succeed for every statistic. Specify the interval method and computational settings, since a simple percentile interval is not universally appropriate.
Past exam of the mathematics course of the University of Cambridge 2001 iii Paper 29 3 b ii Solution Created 2026-10-03 Updated 2026-10-14
Two reasons are selection bias and unequal ascertainment.
- The slaughter cohort need not be representative of the adult dairy herd. Culling can select animals with poor health, while its age distribution may also differ. Eligibility from 24 months includes many cattle below the usual age at clinical onset, and cases before 30 months are rare; differing age composition can strongly alter the underlying risk. These effects can work in different directions, so age and reason for slaughter should be examined rather than assuming a direction.
- A focused veterinary examination at slaughter and repeated routine observations by farmers have different opportunities to detect subtle clinical signs. Early disease can remain unreported in the field, while inspection can reveal previously accumulated undetected cases. Thus one is a single-time screen for existing disease and the other a stream of newly recognized cases; even after converting denominators to comparable time units, prevalence, incidence and detection sensitivity are not identical quantities.
Standardize the sampled population and the ascertainment process before interpreting the numerical detection-rate ratio causally.
Past exam of the mathematics course of the University of Cambridge 2001 iii Paper 41 1 iii Solution Created 2026-10-03 Updated 2026-10-14
A Type Ia supernova can be standardized as a standard candle by correcting its peak brightness using the light curve shape and colour index. Its host spectrum gives the redshift; the corrected apparent magnitude and calibrated absolute magnitude give the distance modulusThe resulting Hubble diagram measures the shape of the luminosity distance as a function of redshift. A spatially flat matter-and-vacuum model, for example, predictsNonzero cosmological spatial curvature changes both the expansion rate and the conversion from radial to transverse comoving distance. Comparing the observed Hubble diagram with these predictions constrains the matter density parameter, cosmological constant and cosmological spatial curvature; a single low-redshift observation cannot separate them. The nearby sample calibrates the common intercept, while distant Type Ia supernovae supply the redshift-dependent shape. Relative faintness at a specified redshift means a greater inferred luminosity distance.
Type Ia supernova cosmology requires correction for interstellar extinction, possible evolution of the standardized luminosity, and selection bias. The absolute magnitude–Hubble constant degeneracy means that an uncalibrated sample constrains the shape rather than independently determining .
Past exam of the mathematics course of the University of Cambridge 2006 iii Paper 44 1 a Solution Created 2026-10-03 Updated 2026-10-14
Use a longitudinal study with repeated measurements of the same individuals. Recruit people across several starting ages and measure the same ear, with a standardized anatomical definition and instrument, at baseline and subsequent visits over several years. Repeat measurements at each visit to estimate measurement error; where practical, conceal previous measurements from the observer. Record sex, body size, birth cohort and observer, and document reasons for loss to follow-up.
A simple random-intercept linear mixed model for ear length at elapsed follow-up time isThe individual intercept absorbs persistent differences between people. Alternatively, differences remove that intercept directly. Estimate the within-person regression coefficient and its confidence interval, accounting for the correlation of repeated measurements. Compare it with the cross-sectional slope rather than treating the two slopes as automatically identical.
This longitudinal study directly checks whether an individual's ears tend to lengthen as time passes. Several entry-age cohorts help assess whether the change varies with starting age. Calendar-period changes and selective dropout can still affect interpretation, so repeated observation alone does not eliminate every possible source of confounding or selection bias.
Past exam of the mathematics course of the University of Cambridge 2006 iii Paper 44 1 c Solution Created 2026-10-03 Updated 2026-10-14
A cross-sectional age association need not equal within-person growth. The following alternatives each suggest a concrete extension of the study.
- Birth-cohort differences. Older participants were born in a different period. Persistent differences established early in life could produce a positive age slope even if every person's adult ear length were constant. Repeat measurements in a longitudinal study of several birth cohorts, separating within-person changes from differences in their initial measurements. Calendar age, birth year and observation year are algebraically related, so all three unrestricted linear effects cannot be identified without additional assumptions.
- Confounding by sex, body size or ancestry. Suppose a characteristic is associated both with ear length and with age among the sampled patients. Omitting it then changes the apparent age slope. Measure these characteristics, examine their distributions across ages, and fit a multiple linear regression or compare age slopes within appropriately defined groups. A guessed overall sex proportion supplies none of the required age-by-sex information.
- Selection bias from attendance or consent. Consulting patients and people willing to be measured may differ from the general population. For example, attendance could select different combinations of age and body characteristics. Sample from a population register, recruit people regardless of consultation, record refusals and compare age-specific recruitment patterns; compare the resulting slope with the practice-based slope.
- Survivorship bias. Ear length or a correlated characteristic might be associated with survival. Then people with larger ears could be overrepresented among older survivors even without growth. Follow a population cohort, measure ears before subsequent deaths or dropout, and investigate whether survival or retention depends on initial ear length after accounting for age and other covariates. This tests the proposed selection mechanism rather than assuming it.
- Age-related measurement error. Different practitioners may see different age mixes and have different measurement offsets; observers might also measure older ears differently. Use one protocol, calibrate instruments, randomly allocate observers where possible, conceal previous readings, obtain duplicate independent measurements and record observer identity. Analyze observer effects and the repeat-measurement differences before interpreting an age slope biologically.
These explanations are hypotheses to be checked, not established facts about the participants. Neither a small p-value nor an approximately linear scatter plot excludes them.
Past exam of the mathematics course of the University of Cambridge 2010 iii Paper 34 1 g Solution Created 2026-10-03 Updated 2026-10-14
The groups were not formed by randomization. Their background health, access to care, clinic, country and reasons for selecting a delivery method can predict both the method and the outcome, creating confounding and selection bias. Moreover, the label describing the reference delivery may be determined after labour: patients whose attempted nonoperative delivery develops a complication can move into an operative category, leaving a selected uncomplicated reference group. Comparing intended strategies would avoid treating this post-baseline classification as a baseline assignment.
A useful comparison would therefore align eligibility and follow-up time, compare intended delivery strategies among patients genuinely eligible for both, and account for sufficient pretreatment predictors. Logistic regression cannot create comparable patients where the positivity assumption fails or remove unmeasured confounding.
Past exam of the mathematics course of the University of Cambridge 2015 iii Paper 35 1 c Solution Created 2026-10-03 Updated 2026-10-14
A complete-case analysis changes the mix of people contributing to the logistic regression. Let mean that both required variables are observed, and write . Within , Bayes' theorem givesIf completeness depends differently on the two responses, and that selection also varies with partner count or the other predictors, the extra term changes the fitted regression coefficients. For example, people with high counts and poor condom use can be under-represented, weakening or otherwise distorting the apparent association. Dropping people with missing also excludes some whose membership of the target group is unknown.
The fitted conditional association among complete responders need not equal the target-population association. No direction of selection bias follows without specifying the response mechanism. Selection on predictors alone need not bias a correctly specified conditional logistic regression; it is the outcome-related selection described here that removes that reassurance.
Past exam of the mathematics course of the University of Cambridge 2015 iii Paper 35 1 f Solution Created 2026-10-03 Updated 2026-10-06
Treat the extra measurements as auxiliary variables for multiple imputation. Include alcohol consumption and disease history as predictors in both the and conditional distributions, together with age, sex and the other incomplete response. If these auxiliary variables themselves have missing entries, give them appropriate models and update them within the same multiple imputation by chained equations sweeps. Preserve observed entries throughout.
Useful auxiliary measurements can improve predictions, precision and the plausibility of missing at random. Their association with the missing responses reduces uncertainty in the imputed values; association with nonresponse may explain some selection that otherwise depends on unobserved values. This can reduce selection bias under the enlarged missing at random model and reduce between-imputation variance. It does not guarantee elimination of missing not at random dependence. They need not be added to the final substantive logistic regression merely because they help imputation: the final analysis should still estimate its intended association.
Past exam of the mathematics course of the University of Cambridge 2018 iii Paper 207 3 e Solution Created 2026-10-03 Updated 2026-10-05
For the electronic-record observational study, three strengths are:
- The large dataset improves precision of the estimated statistical association and allows detection of associations too small for a modest study.
- Adjusting for age, sex, and ancestry addresses specified measured sources of confounding.
- Routinely collected hypercalcaemia and migraine diagnoses give evidence from clinical practice, rather than depending solely on recalled calcium intake.
Three weaknesses are:
- Residual confounding remains possible: for example, other diseases, medicines, or health behaviours may affect both serum calcium and migraine.
- Co-occurrence does not establish temporal order, leaving reverse causality possible. Migraine, its management, or the decision to investigate it may alter recorded serum calcium.
- Selection bias and diagnostic misclassification can arise because patients enter health records and receive testing nonrandomly. Moreover, hypercalcaemia is not a direct measure of dietary calcium intake, limiting the claimed intervention's interpretation.
For two-sample Mendelian randomization, three strengths are:
- Genetic variants are inherited before adult migraine develops, making ordinary disease-to-genotype reverse causality implausible.
- Under instrumental-variable independence, genetic assignment reduces the environmental confounding that undermines a conventional observational study.
- The large outcome sample and combined genetic risk score improve precision, while concordant weighted-median Mendelian randomization estimator and MR-Egger regression results offer checks under different assumptions about horizontal pleiotropy.
Three weaknesses are:
- Horizontal pleiotropy can violate the exclusion restriction, while population stratification or inherited family effects can violate instrumental-variable independence. In the supplied figure, one variant has a much larger serum calcium association than the others, motivating a check for dependence on that single variant.
- The score explains only 1.25% of serum calcium variance, so instrument strength must be assessed rather than assumed. Different populations in the two samples can also undermine transport of genetic associations. A small variance fraction is not itself proof of a weak instrument in a large sample.
- A lifelong genetically influenced shift in serum calcium is not equivalent to changing dietary calcium intake. The robustness methods also need assumptions: the weighted-median Mendelian randomization estimator needs more than half of the weight from valid instruments, and MR-Egger regression needs the InSIDE assumption. The nonsignificant MR-Egger regression intercept is compatible with zero average directional horizontal pleiotropy, but with eight variants it cannot prove all instruments valid.
Thus the agreement supports the hypothesis, but does not establish the effect of a dietary intervention.
