Hola Javier, te faltó poner el sexo dentro del nuevo data.frame
Luciano


El 11 de julio de 2012 14:07, Marcuzzi, Javier Rubén <
javier.ruben.marcuzzi@gmail.com> escribió:

> 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)
>

	[[alternative HTML version deleted]]

