[R-es] Weighted Likelihood

[Mauricio Monsalvo]

Hola Olivier,
Tengo un set de datos (N = 309) con un "ponderador", que puedo interpretar
como la representatividad del dato en la población.  Tal que si
datos$Ponderador =  3,5556, implica que mi observación "vale" por 3,5556.
Estuve buscando una forma de utilizar el WEIGHT BY Ponderador . del SPSS
pero, "Unlike SPSS, weights in R are applied at the operation level rather
than globally." https://stackoverflow.com/questions/
36838581/easy-way-like-in-spss-to-weight-my-variables-
by-another-variable-in-r-studio
La discusión se extiende un poco más acá: https://stackoverflow.
com/questions/7026549/weight-data-with-r-part-ii Sinceramente, no me da la
nafta para discutir la cuestión de fondo.
Lo cierto es que necesito "ponderar" mis casos en varias funciones, la
mayoría de las cuales son distribuciones de frecuencia. Pueden sintetizarse
en esta forma:

*Capacitación al personal docente, según Región y Dependencia*

*Recibieron Capacitación*

*N*

*Porc*

*PorcAc*

*NEA*

*NOA*

*No EPI*

*Semi EPI*

*EPI*

Sí, regularmente

191

61,4

61,4

69,8

57,7

36

67,8

81,1

Sí, ocasionalmente

66

21,2

82,6

21,9

20,9

25

24,8

12,2

Nunca

53

17

99,6

7,3

21,4

38

7,4

6,7

NS/NC

1

0,3

99,9

1
(Creo que no es necesario el código).​
Pondero mi matriz:
>
​ library(
​​
survey
​)
​
​
> datos <- svydesign(id = ~1 #  is a formula for no clusters
+                      , weights = ~datos$Ponderador
+                      , data = datos
​ ​
#
data may be a character string giving the name of a table​
+                      # , fpc = ~fpc
+                      )
> class(datos)
[1] "survey.design2" "survey.design"
​
Ahora bien, cuando intento llevar a cabo el análisis, resulta:
​
> datos[, .N, by =
​ ​
.(B2a)]
Independent Sampling design (with replacement)
svydesign(id = ~1, weights = ~datos$Ponderador, data = datos)
​
¿Solo puedo correr funciones contenidas en la librería "​survey", como por
ejemplo svytable? En esta entrada parece confirmarlo
https://www.r-bloggers.com/social-science-goes-r-weighted-survey-data/,
pero necesito que me lo digan en Castellano :(
​ ​
Le di una mirada a:
library(weights) #
https://cran.r-project.org/web/packages/weights/weights.pdf
Pero no me resultó. Creo que tiene un alcance limitado y no

​¿No existe una forma alternativa para ponderar datos[, .N, by =
.(B2a)]  utilizando
la variable Ponderador?

​Muchas gracias.​
Saludos



El 16 de abril de 2015, 6:40, Olivier Nuñez <onunez en unex.es> escribió:

> No das muchos detalles sobre porque quieres ponderar los datos, pero
> intuyo que es por el diseño.
> Si es así, la función svydesign del paquete te permite especificar bien la
> probabilidad de selección del dato (probs) en la muestra o de manera
> equivalente, o el "sampling weight" que puedes interpretar como la
> representatividad del dato en la población. Una vez especificado el diseño
> mediante esta función, puedes llevar cualquier análisis estandár de tus
> datos. Un saludo. Olivier
>
> ----- Mensaje original -----
>
> De: "Víctor Nalda Castellet" <victor.nalda.castellet en gmail.com>
> Para: "Olivier Nuñez" <onunez en unex.es>
> CC: "r-help-es" <r-help-es en r-project.org>
> Enviados: Jueves, 16 de Abril 2015 11:25:10
> Asunto: Re: [R-es] Weighted Likelihood
>
> Hola,
>
> He buscado por el pdf de la documentación qué función es la que me puede
> servir pero no encuentro ninguna.
> ¿Tenías alguna función especifica de ese paquete pensada?
>
>
> Un saludo, y muchas gracias por la ayuda.
>
> El 16 de abril de 2015, 10:45, Víctor Nalda Castellet <
> victor.nalda.castellet en gmail.com > escribió:
>
>
>
> ¡Muchas gracias Olivier!
>
>
> Un saludo.
>
> El 16 de abril de 2015, 10:44, Olivier Nuñez < onunez en unex.es > escribió:
>
> <blockquote>
> Mira el paquete survey.
> Un saludo. Olivier
>
> ----- Mensaje original -----
> De: "Víctor Nalda Castellet" < victor.nalda.castellet en gmail.com >
> Para: "r-help-es" < r-help-es en r-project.org >
> Enviados: Miércoles, 15 de Abril 2015 14:21:33
> Asunto: [R-es] Weighted Likelihood
>
> Buenas tardes,
>
> Estoy intentando ajustar distribuciones utilizando un vector de ponderación
> en los datos (Weighted Likelihood). ¿Existen paquetes en R que resuelven
> esto? He mirado ya el paquete "wle" pero no me permite introducir los pesos
> mediante los cuales ponderar los datos.
>
> En un primer momento, se me ha ocurrido realizar lo siguiente: repetir
> cada elemento del vector datos tantas veces como indique el vector pesos
> “modificado” (multiplicado por una potencia de 10, para que sean todos
> números enteros). Y con este vector de datos repetidos realizar el ajuste.
>
> Pero me gustaría conocer si existen otras formas "más elegantes" de
> hacerlo.
>
>
> Un saludo, y muchas gracias por adelantado.
>
> --
> Víctor Nalda Castellet
>
> [[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
>
>
>
>
>
>
> --
> Víctor Nalda Castellet
>
> </blockquote>
>
>
>
>
> --
> Víctor Nalda Castellet
>
>
>         [[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
>



-- 
Mauricio

	[[alternative HTML version deleted]]