[R-es] Regresión ponderada

[Javier Marcuzzi]

Estimado Freddy

Esa parte me gusta de nlme vs lme4, no se en cuánto es técnicamente fundamentado a mi gusto o es solo una subjetividad desde mi desconocimiento al no venir de una carrera estadística o matemática.

Javier Rubén Marcuzzi

De: Freddy Omar López Quintero
Enviado: miércoles, 14 de junio de 2017 22:01
Para: Gerardo Gold
CC: Lista R
Asunto: Re: [R-es] Regresión ponderada

2017-06-14 15:35 GMT-04:00 Gerardo Gold <gerardo.gold.bou en gmail.com>:

> Para el caso de datos agrupados, donde ya se tiene la desviación estándar o
> varianza calculada ajustar un modelo ponderado es fácil, pero no se como
> hacerlo con datos sin agrupar.
>

​Si no entiendo mal, el núcleo del problema son las varianzas de la
relación. Una opción sería incluir entonces una estructura para la
varianza. Tomando arbitrariamente el conjunto de datos swiss como ejemplo,
podría ser algo como:

library(nlme)
> gls(Fertility~Education, data=swiss, weights=varIdent(form=~1|Education))​


​y/o con el argumento correlation de la misma función.​

​Quizás sirva.

¡Salud!​



-- 
«Pídeles sus títulos a los que te persiguen, pregúntales
cuándo nacieron, diles que te demuestren su existencia.»

Rafael Cadenas

	[[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]]