Hola a todos:

Estaba pensando como explicar algo para unas personas que no son estadísticos, necesitan el número para tomar decisiones, dicen ¿cuánto sería en porcentaje? pero no conocen de modelos mixtos, no por ser ignorantes o incapaces, todo lo contrario, simplemente nadie les dijo nada al respecto y son buenos profesionales en su profesión.

Por lo cuál escribí algo muy sencillo, seguramente al procesarlo se darán cuenta, no interesa el resultado, simplemente es mostrar algo que no es porcentaje pero que es útil, pero me encontré con un problema, es con el predict, seguramente al ejecutar el siguiente código en R se darán cuenta (está escrito en este).

ano_nacimiento <- c(2000, 2000, 2001, 2001, 2001, 2001, 2002)
ano_nacimiento <- as.factor(ano_nacimiento)
sexo <- c('macho', 'hembra', 'macho', 'hembra', 'macho', 'hembra', 'macho')
clima <- c('bueno', 'regular', 'malo','bueno', 'regular', 'malo', 'bueno')
peso <- c(354, 251, 327, 328, 301, 270, 330)
dataset <- data.frame(ano_nacimiento, sexo, clima, peso)
library(nlme)
modelo <- lme(peso ~ sexo-1 +clima, random=~1|ano_nacimiento, dataset, na.action="na.exclude")
summary(modelo)
modelo$coefficients$fixed
modelo$coefficients$random
f <-fitted(modelo)
r <-resid(modelo)
dataset2 <- data.frame(ano_nacimiento, sexo, clima, peso, f, r)
dataset2
plot(modelo)
predict(modelo)
# predecir por ejemplo, ano_nacimiento 2002, sexo hembra (no está en los datos), clima bueno y malo
predecir <- data.frame(ano_nacimiento=2002, clima=c('bueno','malo'))
predict(modelo, predecir, level=0:1)

No me doy cuenta el error, ¿alguna sugerencia o material para leer?

Desde ya muchas gracias.
Javier Marcuzzi


	[[alternative HTML version deleted]]

