Hola a todos, llevo un par de días intentado usar predict para estimar en un
modelo polinómico. Os mando el fragmento que creo que es muy fácil de entender y
os pongo debajo mi duda.
valores<- c(
73.74,118.96,107.32,114.95,91,122.27,108.28,128.59,111.94,106.31)
dia<-c(1,2,3,4,5,6,7,8,9,10)
# hago un modelo polinomico de grado cuatro
modelo<-lm(valores ~ dia+I(dia^2)+I(dia^3)+I(dia^4))
# ahora quiero predecir el valor en el dia 11 y pongo la siguiente linea pero no
funciona
predict(modelo,11)
¿En qué me equivoco? ¿Me podéis ayudar? Gracias por adelantado,
Fernando
[[alternative HTML version deleted]]
Hola Fernando,
Basicamente necesitas crear un data.frame que incluya el valor de
"dia" que
te interesa y luego pasarlo a predict():
predict(modelo, data.frame(dia = 11))
1
44.52167
Para mas detalles, dale una mirada a ?predict.lm
Saludos,
Jorge Ivan Velez
2011/1/11 Fernando Sanchez <>
> Hola a todos, llevo un par de días intentado usar predict para estimar en
> un
> modelo polinómico. Os mando el fragmento que creo que es muy fácil de
> entender y
> os pongo debajo mi duda.
>
> valores<- c(
> 73.74,118.96,107.32,114.95,91,122.27,108.28,128.59,111.94,106.31)
> dia<-c(1,2,3,4,5,6,7,8,9,10)
> # hago un modelo polinomico de grado cuatro
> modelo<-lm(valores ~ dia+I(dia^2)+I(dia^3)+I(dia^4))
> # ahora quiero predecir el valor en el dia 11 y pongo la siguiente linea
> pero no
> funciona
> predict(modelo,11)
>
> ¿En qué me equivoco? ¿Me podéis ayudar? Gracias por adelantado,
>
> Fernando
>
>
>
>
> [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
[[alternative HTML version deleted]]
No estoy seguro , pero creo que los nuevos valores deben estar en un data.frame. prueba con datos<-data.frame(11) predict(modelo,datos) El 11/01/11 20:12, Fernando Sanchez escribió:> Hola a todos, llevo un par de días intentado usar predict para estimar en un > modelo polinómico. Os mando el fragmento que creo que es muy fácil de entender y > os pongo debajo mi duda. > > valores<- c( 73.74,118.96,107.32,114.95,91,122.27,108.28,128.59,111.94,106.31) > dia<-c(1,2,3,4,5,6,7,8,9,10) > # hago un modelo polinomico de grado cuatro > modelo<-lm(valores ~ dia+I(dia^2)+I(dia^3)+I(dia^4)) > # ahora quiero predecir el valor en el dia 11 y pongo la siguiente linea pero no > funciona > predict(modelo,11) > > ¿En qué me equivoco? ¿Me podéis ayudar? Gracias por adelantado, > > Fernando > > > > > [[alternative HTML version deleted]] > > > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es[[alternative HTML version deleted]]
Hola, gracias a José Luis y Jorge Iván. José Luis: había probado con lo que tú
comentabas del data frame pero no lo conseguía. Con la respuesta de Iván y el
detalle de data.frame(dia=11) Desde luego que nunca se me hubiera ocurrido.
¡¡CASO RESUELTO!! ¡¡ Y NUEVAMENTE GRACIAS A LOS DOS!!
Saludos,
Fernando
Hola Fernando,
Basicamente necesitas crear un data.frame que incluya el valor de
"dia" que te
interesa y luego pasarlo a predict():
predict(modelo, data.frame(dia = 11))
1
44.52167
Para mas detalles, dale una mirada a ?predict.lm
Saludos,
Jorge Ivan Velez
2011/1/11 Fernando Sanchez <>
Hola a todos, llevo un par de días intentado usar predict para estimar en un
modelo polinómico. Os mando el fragmento que creo que es muy fácil de entender y
os pongo debajo mi duda.
valores<- c(
73.74,118.96,107.32,114.95,91,122.27,108.28,128.59,111.94,106.31)
dia<-c(1,2,3,4,5,6,7,8,9,10)
# hago un modelo polinomico de grado cuatro
modelo<-lm(valores ~ dia+I(dia^2)+I(dia^3)+I(dia^4))
# ahora quiero predecir el valor en el dia 11 y pongo la siguiente linea pero no
funciona
predict(modelo,11)
¿En qué me equivoco? ¿Me podéis ayudar? Gracias por adelantado,
Fernando
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
[[alternative HTML version deleted]]