Isa García Barón
2014-Nov-14 15:37 UTC
[R-es] Cómo aplicar weights a las observaciones en un GLM binomial
Gracias por la ayuda Jose Luis. pero o no te he entendido bien o mi duda es tan sencilla que no me he explicado. SI yo tampoco he entendido mal tu explicación, mi problema es cómo obtengo ese "tus.pesos" para introducir, por ejemplo, en la función: library(survey) # objeto del diseño muestral ddatos <- svydesign(id=~1, weights =~ tus.pesos, data = tus.datos) # en caso de una reg logística modelo <- svyglm(respuesta~ var1 + var2, family = binomial, design = ddatos) Después si se proseguir con la realización del modelo, pero lo que no sé es cómo pondero la muestra, es decir si ese "tus.pesos" es una vecotr, osea una nueva variable en mi dataset, qué valores debe contener. Perdonad las molestias, Un saludo y gracias Isa [[alternative HTML version deleted]]
Carlos Ortega
2014-Nov-14 15:57 UTC
[R-es] Cómo aplicar weights a las observaciones en un GLM binomial
Hola, Hubo un hilo en el R-Help que trató de este tema. ---------------------------------------- Weights in binomial glm http://r.789695.n4.nabble.com/Weights-in-binomial-glm-td1991249.html ---------------------------------------- Saludos, Carlos Ortega www.qualityexcellence.es El 14 de noviembre de 2014, 16:37, Isa García Barón <xanamb en gmail.com> escribió:> Gracias por la ayuda Jose Luis. pero o no te he entendido bien o mi duda es > tan sencilla que no me he explicado. > > SI yo tampoco he entendido mal tu explicación, mi problema es cómo obtengo > ese "tus.pesos" para introducir, por ejemplo, en la función: > > library(survey) > > # objeto del diseño muestral > > ddatos <- svydesign(id=~1, weights =~ tus.pesos, data = tus.datos) > > # en caso de una reg logística > > modelo <- svyglm(respuesta~ var1 + var2, family = binomial, design > ddatos) > > > Después si se proseguir con la realización del modelo, pero lo que no sé es > cómo pondero la muestra, es decir si ese "tus.pesos" es una vecotr, osea > una nueva variable en mi dataset, qué valores debe contener. > > Perdonad las molestias, > > Un saludo y gracias > > Isa > > [[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 >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
José Luis Cañadas
2014-Nov-14 20:02 UTC
[R-es] Cómo aplicar weights a las observaciones en un GLM binomial
El vie, 14 de noviembre de 2014 17:09, Carlos Ortega < cof en qualityexcellence.es> escribió:> Hola, > > Hubo un hilo en el R-Help que trató de este tema. > > ---------------------------------------- > Weights in binomial glm > http://r.789695.n4.nabble.com/Weights-in-binomial-glm-td1991249.html > ---------------------------------------- > > Saludos, > Carlos Ortega > www.qualityexcellence.es > > El 14 de noviembre de 2014, 16:37, Isa García Barón <xanamb en gmail.com> > escribió: > > > Gracias por la ayuda Jose Luis. pero o no te he entendido bien o mi duda > es > > tan sencilla que no me he explicado. > > > > SI yo tampoco he entendido mal tu explicación, mi problema es cómo > obtengo > > ese "tus.pesos" para introducir, por ejemplo, en la función: > > > > library(survey) > > > > # objeto del diseño muestral > > > > ddatos <- svydesign(id=~1, weights =~ tus.pesos, data = tus.datos) > > > > # en caso de una reg logística > > > > modelo <- svyglm(respuesta~ var1 + var2, family = binomial, design > > ddatos) > > > > > > Después si se proseguir con la realización del modelo, pero lo que no sé > es > > cómo pondero la muestra, es decir si ese "tus.pesos" es una vecotr, osea > > una nueva variable en mi dataset, qué valores debe contener. > > > > Perdonad las molestias, > > > > Un saludo y gracias > > > > Isa > > > > [[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 > > > > > > -- > Saludos, > Carlos Ortega > www.qualityexcellence.es > > [[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 >[[alternative HTML version deleted]]
José Luis Cañadas
2014-Nov-14 20:12 UTC
[R-es] Cómo aplicar weights a las observaciones en un GLM binomial
Hola. Efectivamente, no leí bien tu pregunta. Quizá este post de Raúl Vaquerizo te pueda ayudar. Compara el sobremuestreo en reg logística con asignar pesos a las observaciones (e ilustra como calcular las ponderaciones) http://analisisydecision.es/sobremuestreo-y-pesos-a-las-observaciones-ahora-con-r/ Espero que esta vez si te sea útil. Saludos El 14/11/2014 17:09, "Carlos Ortega" <cof en qualityexcellence.es> escribió:> Hola, > > Hubo un hilo en el R-Help que trató de este tema. > > ---------------------------------------- > Weights in binomial glm > http://r.789695.n4.nabble.com/Weights-in-binomial-glm-td1991249.html > ---------------------------------------- > > Saludos, > Carlos Ortega > www.qualityexcellence.es > > El 14 de noviembre de 2014, 16:37, Isa García Barón <xanamb en gmail.com> > escribió: > > > Gracias por la ayuda Jose Luis. pero o no te he entendido bien o mi duda > es > > tan sencilla que no me he explicado. > > > > SI yo tampoco he entendido mal tu explicación, mi problema es cómo > obtengo > > ese "tus.pesos" para introducir, por ejemplo, en la función: > > > > library(survey) > > > > # objeto del diseño muestral > > > > ddatos <- svydesign(id=~1, weights =~ tus.pesos, data = tus.datos) > > > > # en caso de una reg logística > > > > modelo <- svyglm(respuesta~ var1 + var2, family = binomial, design > > ddatos) > > > > > > Después si se proseguir con la realización del modelo, pero lo que no sé > es > > cómo pondero la muestra, es decir si ese "tus.pesos" es una vecotr, osea > > una nueva variable en mi dataset, qué valores debe contener. > > > > Perdonad las molestias, > > > > Un saludo y gracias > > > > Isa > > > > [[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 > > > > > > -- > Saludos, > Carlos Ortega > www.qualityexcellence.es > > [[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 >[[alternative HTML version deleted]]