Solution (source code)

= Solution

When the event comes from the one-covariate subject, the <Schoenfeld function for a three-person binary risk set> is
$$
\boxed{s(\beta)=1-\frac{e^\beta}{2+e^\beta}=\frac2{2+e^\beta}.}
$$
Hence
$$
\boxed{s(-\infty)=1,\qquad s(0)=\frac23,\qquad s(\infty)=0.}
$$
These are limits where necessary. At a very negative coefficient the model assigns almost no event <probability> to the observed one-covariate subject, producing the largest positive discrepancy. At zero coefficient all three subjects have equal <hazard functions>, so the expected event covariate is $1/3$ and the discrepancy is $2/3$. At a very positive coefficient the observed subject is predicted to have the event almost surely, so the discrepancy tends to zero. The function decreases strictly and stays positive at every finite coefficient: this one event alone favors increasing $\beta$, while the other events in the complete <Cox partial likelihood> determine its overall estimate.