Biomedical Engineering Reference
In-Depth Information
difference in the properties of the estimators, due to these conventions. Thus
for convenience, we make use the convention L < T R in this chapter. Then
8
<
(L;R]
if L < R
the resulting interval I is of the form I =
where (a;1) is
:
[T;T]
if L = R;
treated as (a;1]. The common likelihood function for interval-censored data
is given by Peto (1973) as
Y
(f(Ti))
if
))
if
(F(R
if
) F(L
if
))
1
if
L(F) =
(2.1)
if
where
if
is the indicator function of the event fL
if
= R
if
g, F is the cdf of T, f
the density function of T under the parametric setup but f(t) = F(t)F(t)
under the nonparametric setup. No further assumptions are made in Peto
(1973). The derivation of the nonparametric maximum likelihood estimator
(NPMLE) under L does not need any further assumptions. Peto supposes
that the MLE or NPMLE under L is consistent and asymptotically ecient
without further assumptions. However, this is not true (see a counterexample
as follows).
Example 1 Suppose that a random sample of IC data are of the forms:
(L
if
;R
if
) = (1; 2) or (2; 3) or (1;1), each having n
if
replications. Assume
a parametric setup that fT
T
(3) = p and fT
T
(2) = 1 p. By Equation (2.1),
the likelihood based on (L
if
;R
if
)'s is L = p
n
2
(1 p)
n
1
(1)
n
3
. The MLE of p is
p = n
2
=(n
1
+ n
2
).
We shall show that without proper assumptions, the MLE is inconsistent.
Let (U;V ) be a random censoring vector, V = U + 1,
8
<
3=4
if t = 2
if
UjT
(0jt) =
andf
UjT
(2jt) = 1 f
UjT
(0jt)
(2.2)
:
1=2
if t = 3
8
<
(1;U)
if T U
and the observable random vector (L;R)=
(V;1)
if T > V
:
(U;V )
if U < T V:
Search WWH ::
Custom Search