2017-06-14 15:35 GMT-04:00 Gerardo Gold <[email protected]>:
> 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
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
