[R-br] simulação de mutinomiais

   Senhores,      estou interessado em simular uma amostra de multinomiais em que suas probabilidades sao geradas por uma regressão logistica, isto é  Y ~ Multinomial(1, Pk0i ,Pk1i ,Pk2i , Pk3i)  em que i = 1, ..., 100 ( tamanho da amostra)  e k= 1,2,3. Assim, k assume 3 valores para cada i. A dúvida é: como faço para informar essa variação das probabilidades na função rmultinom? ou teria outra maneira de fazer isso?     Grato,   Márcio