1
¿Por qué obtengo diferentes predicciones para la expansión polinómica manual y uso la función R `poly`?
¿Por qué obtengo diferentes predicciones para la expansión polinómica manual y uso de la polyfunción R ? set.seed(0) x <- rnorm(10) y <- runif(10) plot(x,y,ylim=c(-0.5,1.5)) grid() # xp is a grid variable for ploting xp <- seq(-3,3,by=0.01) x_exp <- data.frame(f1=x,f2=x^2) fit <- lm(y~.-1,data=x_exp) xp_exp <- data.frame(f1=xp,f2=xp^2) yp <- predict(fit,xp_exp) lines(xp,yp) …