simul                 package:rhosp                 R Documentation

_s_i_m_u_l_a_t_e _t_h_e _f_i_r_s_t _m_o_d_e_l _o_f _t_h_e _h_o_s_p_i_t_a_l _r_i_s_k

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

     simul simulate the first model of the hospital that is to say 
     simulating the sequence of variables Xi (which follow the disxi
     distribution passed as an argument) and the sequence of Zi (where
     the probability of succes p is passed as an argument) simul
     returns T the vector of durations between two declared side
     effects

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

     simul(nbPatient, disXi, disP, toplot = TRUE)

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

nbPatient: the number of patient of the simulation 

   disXi: the distribution of the variable Xi : disXi is a 3 elements
          list : rangen <=>  a random positive variable generator ;
          nbparam <=> number of parameter of this distribution and
          param <=> a list of parameters 

    disP: the side effect probability (success  probability of Zi) p :
          disP is a 3 elements list : disfun <=> a distribution
          function ; nbparam <=> number of parameter of this
          distribution and param <=> a list of parameters 

  toplot: a logical variable to plot the variable Zi 

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

     T the vector of durations between two declared side effects

_N_o_t_e:

     further notes

_A_u_t_h_o_r(_s):

     CJ

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

     simul2

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

     arg1Exp<-list(rangen=rexp,nbparam=1,param=list(1/3));
     arg1Bin<-list(rangen=rbinom,nbparam=2,param=list(1,1/20));
     arg1Unif<-list(rangen=runif,nbparam=2,param=list(0,20));
     arg1Lnorm<-list(rangen=rlnorm,nbparam=2,param=list(1/4,1));

     arg2Exp<-list(disfun=pexp,nbparam=1,param=list(1/5));
     arg2Cst<-list(disfun=pcst<-function(x,p) p ,nbparam=1,param=list(1/3));
     arg2Comp<-list(disfun=pcomp<-function(x,mu1,mu2,mu3){1-1/3*exp(-mu1* x)-1/2*exp(-mu2 *x)-1/6*exp(-mu3 *x)}
     ,nbparam=3,param=list(1/3,1/5,1/10));
     arg2Gamma<-list(disfun=pgamma,nbparam=2,param=list(3,1/3));
     arg2Lnorm<-list(disfun=plnorm,nbparam=2,param=list(1/20,2));

     T<-simul(100,arg1Exp,arg2Exp)
     T<-simul(100,arg1Bin,arg2Cst)
     T<-simul(100,arg1Unif,arg2Comp)
     T<-simul(100,arg1Lnorm,arg2Gamma)

