[R-es] Logit y ponderaciones

[J. Miguel Marin]

Hola,

En la librería "survey" tienes el comando "svyglm" (Survey-weighted 
generalised linear models)

Creo que es lo que buscas.

Un saludo

> Hola,
>
> Hace un tiempo que estoy intentando hacer un logit con R utilizando
> ponderaciones, pero creo que no doy con la solución y agradecería si
> alguien pudiera ayudarme. Os cuento:
>
> Los datos que utilizo son "dato por individuo" y provienen de encuesta. Las
> ponderaciones se corresponden con proporciones de población para solventar
> desvíos en las cuotas.
>
> El resumen de los datos que utilizo es el siguiente:
>
> dep       indep            pond
> 0:25494   in :19874   Min.   :0.2678
> 1: 2811   out: 8426   1st Qu.:0.4718
>           NR :    5   Median :1.0000
>                       Mean   :0.9982
>                       3rd Qu.:1.0000
>                       Max.   :2.3048
>
> La expresión que utilizo es la siguiente:
>
> mod1<-glm(data=prova,dep~indep,family=binomial(link="logit"),weights=pond)
>
> Pero me avisa con el siguiente warning:
>
> "In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!"
>
> Cuando vi el aviso la primera vez me sorprendió y busqué un poco. Creo que
> el problema viene porque las ponderaciones en glm no son exactamente lo que
> yo estoy utilizando ni están pensadas para mi tipo de datos. Pero no
> encuentro cómo puedo hacerlo y no me parece normal obviar el warning.
>
> De momento estoy barajando la posibilidad de replicar los datos de tal
> forma que no sean necesarias las ponderaciones, o quizás generar recuentos
> y trabajar a partir de allí. Pero no tengo demasiado claro ni si es la
> manera ni si acabaré "cochineando" demasiado. Supongo que no ser
> estadístico también influye en no encontrar la solución...
>
> Tenéis alguna idea?
>
> Muchas gracias,
>
> Manel.
>
> 	[[alternative HTML version deleted]]
>
>




jm~

_______________________________

        J. Miguel Marin

http://www.est.uc3m.es/jmmarin

    Dep. of Statistics
University Carlos III of Madrid
        Spain (E.U.)