msBIC                 package:mixPHM                 R Documentation

_P_H_M _m_o_d_e_l _s_e_l_e_c_t_i_o_n _w_i_t_h _B_I_C

_D_e_s_c_r_i_p_t_i_o_n:

     This function fits models for different proportionality
     restrictions.

_U_s_a_g_e:

     msBIC(x, K, method = "all", Sdist = "weibull", EMoption = "classification", EMstop = 0.0001, maxiter = 1000)

_A_r_g_u_m_e_n_t_s:

       x: Data frame or matrix of dimension n*p with survival times
          ('NA''s allowed).

       K: A vector with number of mixture components.

  method: A vector with the methods provided in 'phmclust':  With
          '"separate"' no restrictions are imposed, '"main.g"' relates
          to a group main effect,  '"main.p"' to the variables main
          effects. '"main.gp"' reflects the proportionality assumption
          over groups and variables. '"int.gp"' allows for interactions
          between groups and variables. If method is '"all"', each
          model is fitted.

   Sdist: Various survival distrubtions such as '"weibull"',
          '"exponential"', and '"rayleigh"'.

EMoption: '"classification"' is based on deterministic cluster
          assignment,  '"maximization"' on deterministic assignment,
          and '"randomization"'  provides a posterior-based randomized
          cluster assignement.

  EMstop: Stopping criterion for EM-iteration.

 maxiter: Maximum number of iterations.

_D_e_t_a_i_l_s:

     Based on the output BIC matrix, model selection can be performed
     in terms of the number of mixture components and imposed
     proportionality restrictions.

_V_a_l_u_e:

     Returns an object of class 'BICmat' with the following values: 

  BICmat: Matrix with BIC values

       K: Vector with different components

  method: Vector with proportional hazard methods

   Sdist: Survival distribution

_S_e_e _A_l_s_o:

     'screeBIC'

_E_x_a_m_p_l_e_s:

     ##Fitting 3 Weibull proportional hazard models (over groups, pages) for K=2,3 components
     data(webshop)
     res <- msBIC(webshop, K = c(2,3), method = c("main.p","main.g"), maxiter = 10)
     res 

