simmixstock0            package:mixstock            R Documentation

_S_i_m_u_l_a_t_e _a _s_a_m_p_l_e _f_r_o_m _s_o_u_r_c_e_s _a_n_d _m_i_x_e_d _s_t_o_c_k

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

     Given information on marker frequencies in sources, the
     contributions of each source to a mixed stock, and sample sizes,
     generates multinomial samples of markers from sources and the
     mixed stock.

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

     simmixstock0(sourcefreq, input.freq, sourcesampsize, mixsampsize, rseed=NULL)

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

sourcefreq: Matrix (markers=rows, source=cols) of true marker
          frequencies in sources

input.freq: True contributions from each source to the mixed stock

sourcesampsize: Sample size from (each) source

mixsampsize: Sample size from mixed stock

   rseed: Random-number seed

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

     Multinomial sample. 

 mixsamp: Sample from mixed stock (vector)

sourcesamp: Sample from sources (matrix, markers=rows, source=cols)

     ...

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

     Ben Bolker

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

     true.freq <- matrix(c(0.65,0.33,0.01,0.01,
                           0.33,0.65,0.01,0.01),ncol=2)
     true.contrib <- c(0.9,0.1)
     x <- simmixstock0(true.freq,true.contrib,sourcesampsize=50,mixsampsize=100,rseed=1004)

