slice.design              package:pomp              R Documentation

_D_e_s_i_g_n _m_a_t_r_i_c_e_s _f_o_r _l_i_k_e_l_i_h_o_o_d _s_l_i_c_e_s _a_n_d _p_r_o_f_i_l_e_s

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

     'slice.design' generates a data-frame representing points taken
     along one or more slices through a point in a multidimensional
     space. 'profile.design' generates a data-frame where each row can
     be used as the starting point for a profile likelihood
     calculation.

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

     slice.design(vars, n)
     profile.design(..., vars, n)

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

     ...: Specifies the parameters over which to profile. 

    vars: For the 'slice.design' case: 'vars' is a named list of
          numeric vectors, each of which has length either 1 or 3.
          Variables along which slices are to be taken should have
          length 3, corresponding to the minimum of the range, central
          point, and maximum of the range. For fixed variables, specify
          just the value.

          For the 'profile.design' case: Named list of numeric vectors,
          each of which has length 2, specifying the range over which
          the parameter is to be sampled. 

       n: In the case of 'slice.design', the number of points per
          slice. In the case of 'profile.design', the number of starts
          per profile point. 

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

     'slice.design' returns a data frame with 'n' points per slice. The
     column 'slice' is a factor that tells which slice each point
     belongs to. 

     'profile.design' returns a data frame with 'n' points per profile
     point. The parameters in 'vars' are sampled using 'sobol'.

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

     Aaron A. King kingaa at umich dot edu

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

     'sobol'

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

     ## A single 11-point slice through the point c(A=3,B=8,C=0) along the B direction.
     x <- slice.design(list(A=3,B=c(0,8,10),C=0),n=11)
     dim(x)
     plot(x)
     ## Two slices through the same point along the A and C directions.
     x <- slice.design(list(A=c(0,3,5),B=8,C=c(0,0,5)),n=11)
     dim(x)
     plot(x)

     ## A one-parameter profile design:
     x <- profile.design(p=1:10,vars=list(a=c(0,1),b=c(0,5)),n=20)
     dim(x)
     plot(x)
     ## A two-parameter profile design:
     x <- profile.design(p=1:10,q=3:5,vars=list(a=c(0,1),b=c(0,5)),n=20)
     dim(x)
     plot(x)

