Euler-multinomial models        package:pomp        R Documentation

_E_u_l_e_r-_m_u_l_t_i_n_o_m_i_a_l _m_o_d_e_l_s

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

     Density and random-deviate generation for the Euler-multinomial
     death process with parameters 'size', 'rate', and 'dt'.

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

     reulermultinom(n = 1, size, rate, dt)
     deulermultinom(x, size, rate, dt, log = FALSE)

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

       n: integer; number of random variates to generate.

    size: scalar integer; number of individuals at risk.

    rate: numeric vector of hazard rates.

      dt: numeric scalar; duration of Euler step.

       x: Matrix or vector containing number of individuals that have
          succumbed to each death process.

     log: logical; if TRUE, return logarithm(s) of probabilities.

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

     Direct access to the underlying C routines is available: see the
     header file pomp.h, included with the package.

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

reulermultinom: Returns a 'length(rate)' by 'n' matrix. Each column is
          a different random draw. Each row contains the numbers of
          individuals succumbed to the corresponding process. 

deulermultinom: Returns a vector (of length equal to the number of
          columns of 'x') containing the probabilities of observing
          each column of 'x' given the specified parameters ('size',
          'rate', 'dt'). 

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

     Aaron A. King kingaa at umich dot edu

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

     'euler'

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

     print(x <- reulermultinom(5,size=100,rate=c(a=1,b=2,c=3),dt=0.1))
     deulermultinom(x,size=100,rate=c(1,2,3),dt=0.1)

