Page 127 - Contributed Paper Session (CPS) - Volume 7
P. 127

CPS2044 Mohd Asrul Affendi A. et al.
                                0,     ℎ      
                            = {
                            
                                1,                ℎ      
               The corresponding log‐likelihood is;
                                                              −1             
                                                              
                                                                                 
                (, , ) = ∑(1 −  ) [  (( [ ] ) ( [ ] )  ) − [( [ ] ) ]]
                                   
                                       
                                                                ′
                                                                                   ′
                                                    ′
                           =1                                               
                                                                     −1
                                                                    
                              + ∑  [  (( [  + ] ) ( [  + ] )  )
                                       
                                    
                                                                    ′
                                                    ′
                                =1                             
                                                               
                              + [(         )  + (             ) − (              ) ]]
                                                        ′
                                        ′
                                  [ ]   [  + ]  [  + ]
                                                                          ′
                                        
                                                                          
                                                       
               The  model  parameters  a,β,γ  can  be  obtained  using  the  Newton‐Raphson
               algorithm  to  find  the  solution  of  the  likelihood  equations  given  above.
               However, we have used the SPLUS function nlminb to obtain the results.

               Simulation  of  Parametric  Weibull  Time-Varying  Covariate  Model:
               Suppose W is a random variable with cumulative distribution function F, it
               suffices  that  F(w)  follows  uniform  distribution,  U~UNIF[0,1].  Similarly,  the
               survival function 1‐F(w) follows U~UNIF[0,1]. Therefore, the survival function
               in the case of parametric Weibull time covariate model can be define as;
                                     
                      [− (      ) ] ,                                                                                        <  
                                   ′
                              [ ]
                  =
                                                            
                       [− [(  ′  ) + (     ′      ) − (      ′     ) ]] ,    ≥  
                    {         [ ]  [  + ]  [  + ]
               Thus, by inverting the above piece‐wise Weibull function, the survival time T
               can be obtained as;
                                         1                                          
                              ′
                          [ ][−log ()],                                                                    <  [− (  ′  ) ]
                                                                              [ ]
                  =                                                                   .
                                                                 1                 
                                                                  
                      [  + ] [− log() − (     ) + (     ) ] ,   ≥  [− (     ) ]
                          ′
                                                          ′
                                                                                   ′
                                               ′
                    {                     [ ]  [  + ]    [ ]
               Simulation  Studies:  To  illustrate  the  simulation  strategy  and  estimation
               method, we used the following parameters; α=1, β=2 and  γ=1. The covariate
               x_i~Binomial(1,0.5), tci~Weibull(shape=1,scale=0.04). Censoring rates 20% and
               40% were used to check the effect of censoring rate on parameter estimates.
               Also, sample sizes n=100, and n=200 were used to study the effect of sample
               size  on  parameter  estimates.  Performance  metrics  used  to  assess  the
               estimating methods are standard error (SE), bias and Mean square error (MSE).


                                                                  114 | I S I   W S C   2 0 1 9
   122   123   124   125   126   127   128   129   130   131   132