[R-es] Generacion de binomiales correlacionadas

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]]