[R-es] Generacion de binomiales correlacionadas

[Olivier Nuñez]

Es un problema delicado, y como alternativa a lo que JuanMi propone,
puedes también utilizar la metodología expuesta en el articulo  
adjunto.
Consiste básicamente en sumar combinaciones de Bernouilli.
Un saludo. Olivier

--  
____________________________________

Olivier G. Nuñez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es

____________________________________




El 15/12/2010, a las 23:08, Jorge Ivan Velez escribió:

> 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]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es

------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20101216/987debce/attachment-0002.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: Biomial2.pdf
Type: application/pdf
Size: 152811 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20101216/987debce/attachment-0001.pdf>
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20101216/987debce/attachment-0003.html>