Buenas tardes,

Estoy interesado en generar observaciones de una distribucion binomial
bivariada en la que hay _cierto_ grado de correlacion (denotemoslo rho).
Podria por favor alguien indicarme como hacerlo en R?

Este es el contexto. Supongamos que se tienen dos experimentos en los que la
variable respuesta _sigue_ una distribucion binomial, i.e., X_i
~Binomial(n_i, p_i), i=1,2 y que, por ahora, n_i y p_i son conocidos. El
interes principal es calcular T = P(X_1=1, X_2=1).  Si las X_i's fueran
independientes (caso 1), seria suficiente generar binomiales con los
parametros correspondientes a cada tipo de experimento (via rbinom) y
calcular lo que se necesita usando table(). Ahora, si X_i's *no* son
independientes (caso 2), el calculo de T no puede hacerse de la misma forma.
Observe que independiente del modelo (caso) bajo el cual estemos trabajando,
T _podria_ calcularse como una funcion de rho, n_i y p_i. Desafortunadamente
no tengo claro como _modelar_ el segundo caso utilizando R.

Alguna sugerencia?

Muchas gracias,
      Jorge Ivan Velez

	[[alternative HTML version deleted]]

