[R-es] Generacion de binomiales correlacionadas
J. Miguel Marin
jmmarin en est-econ.uc3m.es
Mie Dic 15 23:29:42 CET 2010
Hola,
yo creo que deberías hacerlo vía cópulas.
Mira en http://cran.r-project.org/web/packages/copula/index.html
Es la manera habitual de generar distribuciones multivariantes, a
partir de unas distribuciones marginales fijadas asignando una medida
de correlación entre ellas.
> 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
>
jm~
_______________________________
J. Miguel Marin
http://www.est.uc3m.es/jmmarin
Dep. of Statistics
University Carlos III of Madrid
Spain (E.U.)
Más información sobre la lista de distribución R-help-es