Let the survival time t for individuals in a population have the density
function as follows:
f(t) = αt^α-1 exp(x^T β - t^α exp{x^T β})
where xi^T β = β0 + δiβ1, δi
is 1 if the ith patient was in the treatment group and
0 otherwise. Code a Newton-Raphson algorithm and find the maximum likelihood
estimators of α, β0 and β1.
The table below has the length of remission (in weeks) [t] for acute leukemia patients
in the treatment and control group for a clinical trial.
Treatment 6 13 16 22 23 20 32 32
Control 1 2 2 3 4 5 8 11

"R programming"