Este taller está basado en los ejemplos del libro Data Analysis Using Regression and Multilevel/Hierarchical Models. Los autores del libro dispusieron los datos y los scripts de R en el siguiente enlace http://www.stat.columbia.edu/~gelman/arm/. Visite el enlace para descargar el material necesario.

Preguntas

Además del paquete lme4 que se ha utilizado hasta ahora en el curso, existe otro paquete que se puede utilizar para modelos mixtos, ese paquete se llama nlme.

  1. Visite la página de documentación del paquete nlme.

  2. Descargue el manual de referencia del paquete nlme.

  3. Instale el paquete nlme.

install.packages("nlme")
  1. Cargue el paquete nlme.
library(nlme)
  1. Consulte la ayuda de la principal función lme.
help(lme)

  1. Consulte la ayuda de la base de datos Orthodont, escriba en la consola help("Orthodont").

  2. Dibuje los datos de Orthodont para que entienda de que se tratan los datos.

plot(Orthodont)

  1. En la ayuda de la función lme hay dos ejemplos que usan la base de datos Orthodont, replique los ejemplos y escriba en su cuaderno las expresiones de los dos modelos.

  2. Vuelva a repetir todo el Taller capítulo 12 parte 1 de usando el paquete nlme y la función lme.

  3. Construya en su cuaderno una para comparar la forma de ajustar un modelo mixto con los paquetes nlme y lme4.

Nota: ambos paquetes son los más utilizados para modelos mixtos. Cada uno tiene sus puntos a favor y por esa razón es obligatorio y necesario manejarlos.