1. Introduction
Random effect in the model for association and unobserved heterogeneity is suitably represented by introducing frailty random variable in the model. One way to incorporate frailty in the model is to consider it as an unobserved random proportionality factor that modifies the hazard function of an individual or related individuals. By considering different mortality models, [4] introduced the first univariate frailty model in which he used longevity factor instead of the term frailty. The term frailty was introduced by [5] in the univariate context. [5] and [6] have independently suggested the same model. [7] considered maximum likelihood estimation (M.L.E.) in this model. The estimation is carried out by using E-M algorithm which is suggested by [8] .
In the above frailty model, we observe that the hazard rate is considered to be proportional to baseline hazard rate with constant of proportionality as frailty random variable Z. However, we think that the generalized model should also incorporate the case having no influence of frailty. By considering the hazard rate proportional to the baseline hazard rate with constant of proportionality as
with
being an unknown parameter, one can extend the present frailty model.
Here, we adopt another way as suggested in [9] by incorporating
as a power function modifier of the scale parameter while modeling the conditional location-scale family of survival functions. Let us consider Z as frailty and model the conditional survival function as
(1)
where
is the baseline parameter free survival function. Then the conditional hazard rate is given as
(2)
It can be seen that if
is standard exponential,
(3)
which reduces to the usual frailty model for
.
In this paper, we suggest a new generalized Weibull-exponential frailty model by using the second approach, as used by [9] . To ease the discussion, in Section 2, we first introduce the usual Weibull-exponential frailty model and then we consider the new Weibull-exponential frailty model and study its properties. In Section 3, we discuss maximum likelihood estimation procedure through E-M algorithm. In Section 4, we carry out a simulation study and discuss the results. Fitting of the model on data of the tensile strength of carbon fibers along with likelihood ratio test for
against
is discussed in Section 5. Finally, we provide the conclusion of the fitting of model on data in support of the proposed model.
2. A New Generalized Weibull-Exponential Frailty Model
To overcome or address the problem of heterogeneity in a population resulting from unobserved covariates, [5] [10] and [11] suggested a random effects model for durations. They introduced the frailty and applied it to population data. The classical and commonly applied frailty model assumes a proportional hazards model which is conditional on the random effect (frailty), i.e., the hazard of an individual depends additionally on an unobservable, age-independent random variable Z, which works multiplicatively on the baseline hazard function
.
(4)
In classical frailty model, the conditional distribution of T for an individual with frailty z is
(5)
The base line hazard rate
is obtained as
(6)
Instead of considering modeling of hazard rate, we consider modeling of conditional survival function given the frailty z, which results in the conditional density as
(7)
Mean and variance of the above model can be obtained as below
(8)
(9)
and
(10)
The hazard rate of the model, obtained by replacing
by
in (6), is
(11)
We consider the frailty distribution exponential with mean
,
is known, as
(12)
Here further, we assume that
. This leads to
. Then the joint distribution of T and z is
(13)
Hence, the marginal distribution of T is
(14)
and the conditional distribution of Z given
is
(15)
Estimates of mean for with frailty i.e.
,
and without frailty (i.e.
),
are given by
(16)
(17)
where,
is incomplete gamma function and the above estimates depend on estimates of
and p under the model with frailty and without frailty separately. In (16), above
exists only if
. Using (16), we get
.
3. Maximum Likelihood Estimation
Here, we shall obtain the maximum likelihood estimates of the parameters of gamma-exponential model using E-M algorithm. The likelihood for the gamma-exponential model is given as
(18)
The log-likelihood function of the model is
(19)
To begin with, we presume that ρ is known and obtain the maximum likelihood (m.l.) equations as
(20)
(21)
Then the solution of m.l. equation for θ is
(22)
The m.l. estimator of η is
(23)
Since m.l. equation for p is analytically not solvable, we use Newton-Raphson method. Further we use E-M algorithm as z is unobservable. For that, given the current estimates of p, θ and η, the E-step of the algorithm requires calculations of
and
and they are
(24)
(25)
For different values of ρ in the range of ρ, we obtain the m.l. estimates of p, θ and η using E-M algorithm. The procedure terminates when the likelihood given in (19) is maximized.
4. Simulation Study
In this section, we carry out a simulation study by generating 1000 samples of size 50 each from a generalized Weibull frailty model for
and different values of θ, p and ρ. The samples are used to fit generalized Weibull frailty model under exponential frailty with location and scale parameter equal to 1. The results of this study are given in Table 1. From the table, it can be observed that as the value of θ, p and ρ increases, the frailty model provides accurate estimate of mean whereas estimate of mean without frailty model deviates significantly away from the actual mean. Moreover, the S.E. of the parameter for frailty model is comparatively smaller than without frailty model.
5. Fitting of the Model
In this section, we fit our model on the data of tensile strength of 100 observations of carbon fibers. The data were also used by Nicholas and Padgett in [12] and later discussed in Flaih et al. in [13] . The data were as follows
3.7, 3.11, 4.42, 3.28, 3.75, 2.96, 3.39, 3.31, 3.15, 2.81, 1.41, 2.76, 3.19, 1.59, 2.17, 3.51, 1.84, 1.61, 1.57, 1.89, 2.74, 3.27, 2.41, 3.09, 2.43, 2.53, 2.81, 3.31, 2.35, 2.77, 2.68, 4.91, 1.57, 2.00, 1.17, 2.17, 0.39, 2.79, 1.08, 2.88, 2.73, 2.87, 3.19, 1.87, 2.95, 2.67, 4.20, 2.85, 2.55, 2.17, 2.97, 3.68, 0.81, 1.22, 5.08, 1.69, 3.68, 4.70, 2.03, 2.82, 2.50, 1.47, 3.22, 3.15, 2.97, 2.93, 3.33, 2.56, 2.59, 2.83, 1.36, 1.84, 5.56, 1.12, 2.48, 1.25, 2.48, 2.03, 1.61, 2.05, 3.60, 3.11, 1.69, 4.90, 3.39, 3.22, 2.55, 3.56, 2.38, 1.92, 0.98, 1.59, 1.73, 1.71, 1.18, 4.38, 0.85, 1.80, 2.12, 3.65.
Table 1. Resultant table of simulation study.
We have estimated values of the parameters and their standard error (S.E.) with the help of E-M algorithm implemented by writing a MATLAB program. The estimated value of the parameters and S.E. of the estimates for both models, frailty model and without frailty model, are given in Table 2.
The goodness of fit of the frailty model is done with the help of likelihood ratio test based on marginal likelihood. The likelihood ratio test statistic for testing
is
(26)
The
is rejected if
, where
is upper
percentile of chi-square distribution with one degree of freedom. In the present study on fitting of the model for the data on the tensile strength of carbon fibre,
which is larger than
. This establishes
Table 2. Estimates and S.E. of the parameters.
the goodness of fit of the suggested frailty model.
Finally, it is seen that the estimates of average tensile strength under the model with frailty and without frailty are more or less same but the estimate of S.E. of average tensile strength under frailty model is relatively smaller than that of the model without frailty. Moreover, S.E. of the model parameters for frailty case is significantly smaller than that of model parameters without frailty. Since the sample sizes are same for both models, the extra variability present in the data due to frailty will be extracted by the frailty model and hence, estimate of σ will be smaller as compared to non-frailty model. Now we shall compute the sample size required by the regular model to achieve the same precision (in terms of S.E.) as frailty model. Since,
. Let
and
be the estimates of the standard error under with and without frailty models respectively. Then, from Table 2, we get
Math_79# and
for
. If we wish to achieve
to be same as
, we need
. In terms of percentage increase in sample size, it
amounts to 7% more observations.