Codex @codex  0

Write $a=P(A_1)=P(A_2)$. The complements are decreasing events. The Harris-FKG inequality for a product measure givesThe same correlation inequality for two decreasing events follows from that for their increasing complements, since their covariance is unchanged. ThusBoth $1-a$ and the square root are nonnegative, so rearranging gives the square-root bound for increasing events:No independence between the events themselves is asserted; independence is the coordinate property of the underlying product measure.

At criticality, percolation clusters occur on every scale. Count their filled outer perimeters, discard microscopic loops through a diameter cutoff, and then pass to a scaling limit. The cutoff gives a locally finite measure on loops of macroscopic size inside a bounded region; summing over scales gives a sigma-finite measure rather than a probability distribution.

The conformal invariance of planar percolation makes these limiting perimeter statistics covariant under conformal maps. Moreover, deciding whether a filled outer perimeter lies in a simply connected subdomain uses only the configuration it encloses: changing the configuration outside that subdomain does not change that perimeter. This is the restriction part of the argument. Conformal covariance together with this domain consistency motivates a nonzero conformal restriction measure on simple loops.

This is an informal construction principle, as requested; it does not replace the tightness, boundary-simplification and uniqueness theorems required to construct the continuum measure rigorously.

The Gaussian free field is a random distribution, so its purported level curve cannot be defined by evaluating its height at every point. A coupling must instead be stated using conditional means and covariances. The SLE4 coupling with a Gaussian free field does exactly this.

Choose the Green-function normalizationA Zero-boundary Gaussian free field $H_0$ is characterized on real compactly supported smooth test functions byEquivalently its energy normalization is the Dirichlet inner product $(2\pi {)}^{-1}\int \nabla f\cdot \nabla gdA$. This fixes the otherwise convention-dependent height constant.

Set $\lambda =\pi /2$ andThe field $H$ has prescribed Dirichlet values $-\lambda$ on the negative half-line and $+\lambda$ on the positive half-line. Subtracting $m_0$ recovers the zero-Dirichlet field $H_0$. It is the shifted field $H$, rather than an unshifted zero-boundary field, whose zero-height interface has ordinary chordal ${\mathrm{SLE}}_4$ law.

Let $\eta$ be chordal ${\mathrm{SLE}}_4$ from $0$ to $\infty$, with Chordal Loewner equation ${\partial }_t{g}_t=2/(g_t-W_t)$ and $W_t=2{\beta }_t$. Write $D_t=\mathbb{H}\setminus \eta [0,t]$ and $f_t=g_t-W_t$. The coupling statement is:where, conditionally on the curve, $H_0^{D_t}$ is an independent zero-boundary field in the slit domain. The equality is for restrictions to test functions in that domain, and extends to suitable stopping times. Its covariance is $G_{D_t}(z,w)=G_{\mathbb{H}}(f_t(z),f_t(w))$ by conformal invariance. The two sides of the revealed slit have heights $-\lambda$ and $+\lambda$, in the order specified by their real images under $f_t$. This is the continuum meaning of a Gaussian free field level line.

The fundamental calculations explain why the parameter is four. For $Z_t=f_t(z)$, the Itô formula givesAt $\kappa =4$ the drift vanishes. ThusEach harmonic mean is a bounded martingale before approaching its point. Differentiating the explicit conformally transformed Green function gives the Loewner variation of the Dirichlet Green functionThereforeIn words, the variance learned from the evolving mean exactly equals the covariance lost when the slit is removed. For a Green kernel $cG$ the same calculation gives $\lambda =(\pi /2)\sqrt{c}$; using a different Green normalization changes the height constant, not the SLE parameter.

To construct the coupling, first sample the SLE curve. In the two components to its left and right, sample independent zero-boundary fields and add the corresponding constant heights. Extend these as distributions to obtain the candidate full field. The following martingale identity proves its marginal law rather than merely matching its first two moments.

For a real test function $\phi$, put $M_t=(m_t,\phi )$ andLocalization permits integration of the pointwise identities. They yield $d⟨M{⟩}_t=-d{V}_t$. Henceis a complex local martingale: the $-\frac{1}{2}d{V}_t$ drift cancels the $-\frac{1}{2}d⟨M{⟩}_t$ Itô correction. Since $V_t\ge 0$, $|{\mathcal{Z}}_t|\le 1$, making it a true martingale. At the complete-curve limit, $m_t$ becomes the constant height in each component and the remaining Green kernel becomes that of those components. Thus ${\mathcal{Z}}_{\infty }$ is the conditional characteristic function of the sampled candidate field. Taking expectations givesApplying this to every linear combination of test functions proves the entire Gaussian law, not only its covariance. The conditional version $\mathbb{E}({\mathcal{Z}}_{\infty }\mid {\mathcal{F}}_t)={\mathcal{Z}}_t$ proves the displayed conditional-field statement. Exhaustion by compact test supports justifies passage across the slit and the limiting distributional extensions.

This revealed curve is a local set of a Gaussian free field: conditionally on it, the remaining field is a zero-boundary field plus a specified harmonic function. The spatial Markov property is thereby preserved at random domains, a property not available for arbitrary field-dependent sets. One can strengthen the coupling to a curve measurable from the field. The proof explores compatible interfaces in small subdomains and uses the conditional boundary heights and monotonicity to show that two such interfaces for the same field cannot separate; a countable exhaustion gives uniqueness. This supplies a rigorous replacement for the informal phrase “draw the zero contour.” It does not assert that the distribution has pointwise values.

The coupling is useful in both directions: SLE martingales give exact conditional Gaussian free field data, while the Gaussian Markov structure explains the SLE domain Markov property and its distinguished parameter. With the above normalization, the height jump is $2\lambda =\pi$, and the corresponding interface is chordal SLE4.

In an R formula, I(...) protects ordinary arithmetic from the formula operators: the predictor is the numerical variable $z=x-165$, rather than an instruction to remove a formula term. Predictor centering leaves fitted values and the height slope unchanged but places the intercept at a meaningful height. It often reduces the correlation between intercept and slope estimates.

With $M_i=1_{\{i=2\}}$, treatment contrasts use female as the reference level in a regression factor. The fitted normal linear model isThe estimates and their standard errors areThe intercept estimates female mean weight at height 165. At that same height the male fitted mean is $58.4711+4.3294=62.8005$. Holding height fixed, the fitted male mean exceeds the female mean by 4.3294 weight units. Within either sex, one additional height unit increases fitted mean weight by 0.6211 weight units. If height is in centimetres and weight in kilograms, these interpretations use centimetres and kilograms accordingly. A standard error for the male intercept would require the covariance between ${\hat{\beta }}_0$ and ${\hat{\beta }}_M$, which is not printed.

The prediction equation is $\hat{Y}=58.4711+4.3294M+0.6211(x-165)$. Predictions for individual students also have residual uncertainty, estimated by $s=5.710$, beyond uncertainty in the fitted mean.

The additive model gives useful evidence for both predictors. To test for no adjusted sex effect, use $H_0:{\beta }_M=0$ against $H_1:{\beta }_M\ne 0$. The statistic $T={\hat{\beta }}_M/\mathrm{se}({\hat{\beta }}_M)$ has a Student t-distribution with 97 statistical degrees of freedom under the normal linear model and the null hypothesis. Here $T=3.016$ and $p=0.00327$. Similarly, $H_0:{\beta }_h=0$ against a nonzero common slope gives $T=4.397\sim t_{97}$ under the null hypothesis, with $p=2.81\times 10^{-5}$. Thus, at a fixed height males have a higher fitted mean weight, and within sex greater height is associated with greater weight.

The interaction model adds $\delta M(x-165)$. Its female height slope is $0.6696$ and its male height slope is $0.6696-0.0977=0.5719$. Test $H_0:\delta =0$ against $H_1:\delta \ne 0$ using $T=-0.344\sim t_{96}$ under the null hypothesis, or equivalently $T^2\sim F_{1,96}$. The $p$-value is 0.73139. The additive model is a reasonable parsimonious choice; the data give no evidence for different height slopes between sexes. This does not establish exact equality of population slopes. The additive model also has the higher adjusted coefficient of determination, 0.405 versus 0.3996, and a slightly smaller residual standard error.

Ignoring sex gives a height slope of 0.8769, appreciably above the adjusted slope 0.6211. This is omitted-variable bias in the descriptive regression: a taller group with a larger intercept makes the pooled height association steeper. Algebraically, with $z=x-165$,The printed estimates imply a positive height-sex covariance. Pooled and adjusted slopes therefore describe different comparisons; the difference is not evidence that adding sex changed a physical effect of height. For predictions, retain both sex and height, check regression diagnostics, and assess performance on new students. The residual standard error near 5.7 remains substantial and the coefficient of determination near 0.42 does not promise precise individual prediction.

The moment-based approach specifies marginal means, variances, and within-student covariances, without supplying a joint probability distribution for each count profile. It is a quasi-likelihood or generalized estimating equation strategy. The parameter interpreted as a structural-zero fraction is not automatically an actual probability of a structural-zero component merely because it appears in the moment formulas. A valid positive-definite working covariance and parameter identifiability must also be checked.

The printed covariance specification is not admissible for all the stated parameter values. For example, take $\pi =1/2$, $\tau =1$, and all three component means equal to 10. It gives diagonal variance 80 and off-diagonal covariance 100, so $\mathrm{Var}(Y_{i1}-Y_{i2})=80+80-2(100)=-40$. A covariance matrix cannot have a negative variance in any direction. Thus the moment approach requires additional admissibility restrictions or a valid working covariance; the parameter ranges alone do not define a valid model. This does not alter the hierarchical model, whose covariance derived below is positive semidefinite.

The alternative gives a full hierarchical mixture distribution, specifically a shared zero-inflated Gamma-Poisson count model: a common student random effect is zero with probability $\pi$, and otherwise has a Gamma distribution with mean 1 and variance $\tau$. Conditional on this effect, the three counts are independent Poisson random variables. Integrating it out induces both excess zeros and positive within-student dependence. Zero inflation is shared for the whole student, rather than independently reselected in each term. This full distribution supports likelihood-based inference but is more sensitive to the distributional assumptions.

The marginal mean in the hierarchical model has the same form as the proposed moment mean. Its marginal covariance does not, in general, equal the printed working covariance; the law of total covariance calculation below makes the distinction explicit. At $\tau =0$ the Gamma component is interpreted as its point-mass-at-one limit. At $\pi =1$ all counts are zero, and the log-mean parameterization and regression effects are not identifiable.

For the moment approach, write $m_i=\mathbb{E}{Y}_i$ and let $V_i$ be a positive-definite working covariance constructed from the proposed moment formulas. For identifiable mean coefficients $\theta$, let $D_i=\partial m_i/\partial {\theta }^T$. Generalized estimating equations take the formAcross independent students, correct mean specification and standard regularity conditions give consistent identifiable mean coefficients even if the working covariance is wrong. Sandwich covariance matrix standard errors useResidual second-moment estimating equations can estimate identifiable dispersion and association parameters. The marginal mean only identifies ${\gamma }_0={\beta }_0+\log (1-\pi )$, not ${\beta }_0$ and $\pi$ separately. Moment restrictions beyond the mean would be needed to separate them; as shown below, the printed restrictions do not generally recover the intended latent parameters.

For the hierarchical approach, maximize the observed-data likelihood after integrating out the random effect. Set ${\mu }_{ij}=e^{{\beta }_0+{\beta }^T{x}_i+{\alpha }_j}$, $r=1/\tau$, $M_i={\sum }_j{\mu }_{ij}$ and $s_i={\sum }_j{y}_{ij}$. The joint likelihood contribution isThe integral comes from multiplying three conditional Poisson mass functions by the Gamma density and integrating $b_i^{r+s_i-1}{e}^{-(r+M_i)b_i}$ over $b_i>0$. Independent students give $\ell ={\sum }_i\log L_i$. Direct numerical maximum likelihood or an expectation-maximization algorithm can be used; in an EM algorithm an all-zero profile has an uncertain latent component membership. Information-based likelihood standard errors apply under regular interior conditions, with special care for parameter boundaries.