Suppose two samples of right-censored failure time data: $$\boldsymbol{X_1} = \{(X_{i1}, \delta_{i1}) = (\min(T_{i1}, C_{i1}), 1_{T_{i1} < C_{i1}}): i = 1, ..., n_1$$ $$\boldsymbol{X_2} = \{(X_{i2}, \delta_{i2}) = (\min(T_{i2}, C_{i2}), 1_{T_{i2} < C_{i2}}): i = 1, ..., n_2$$ where $T_{ij}$ are the underlying failure times and $C_{ij}$ are the censoring times. Suppose for each sample, the failure times are distributed according to survival functions $S_1$ and $S_2$, respectively. We are interested in testing whether $S_{T_2}(t) = S_{T_1}^\beta(t)$ and so we propose the hypotheses: $$H_0: \beta = 1, H_1: \beta \neq 1$$ Deriving the log-likelihood function for $\beta$ and differentiating with respect to $\beta$, we obtain the score function under the null hypothesis: $$S(\beta=1) = \sum_{i=1}^{n_{2}} \left\{ \delta_{i2} + \log(S_{T_2}(X_{i2})) \right\}$$
Questions:
- Is the expected value of the score function equal to 0? How can one show this?
- How does this likelihood based test reduce to the log-rank test? Is there any connection?
- Are there any references to formal likelihood based testing procedures for survival analysis?