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