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 <[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 [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
