[R-es] Calcular Error en modelo lineal
Jorge I Velez
jorgeivanvelez en gmail.com
Jue Abr 21 16:09:03 CEST 2016
Buenos dias Jesus,
Esos valores son _aproximados_. Las estimaciones, de acuerdo con teoria de
regresion, podrian obtenerse de manera puntual y construir intervalos de
_confianza_ y _prediccion_ alrededor de estos. Ten en cuenta que el 2do
tipo de intervalos de calcula para observaciones _futuras_.
En R puedes calcularlos de la siguiente manera:
## IC de confianza
## ver ?predict.lm para mas detalles
R> data.frame(y, predict(modelo, interval = "confidence"))
y fit lwr upr
#1 8.35 9.938571 6.580445 13.29670
#2 12.42 11.804286 9.134239 14.47433
#3 18.00 15.535714 13.664949 17.40648
#4 17.58 17.401429 15.396872 19.40599
#5 17.97 19.267143 16.798908 21.73538
#6 20.76 21.132857 18.014915 24.25080
## intervalos de prediccion para x = 25
R> predict(modelo, newdata = data.frame(x = 25), interval = "prediction")
# fit lwr upr
#1 23.93143 17.69035 30.17251
Lo anterior significa que E[y|x=25] = 23.93 y el intervalo de prediccion
del 95% es (17.69, 30.17).
Espero sea de utilidad.
Saludos cordiales,
Jorge.-
2016-04-21 8:56 GMT-05:00 Jesús Para Fernández <j.para.fernandez en hotmail.com
>:
> Buenas, una pregunta.
>
> Si yo estoy calculando un modelo lineal, el caso más simple, 1 variable
> respuesta y una variable explicativa y creo un modelo, me da un R2 del 80%
> y quiero ver como es esa relacion entre las variables, para calcular el
> error de predicción del modelo, basta con ver el intervalo de confianza del
> modelo e irme a los extremos?
>
> Por si no me he expresado bien, un ejemplo tonto:
>
> y=c(8.35,12.42,18.00,17.58,17.97,20.76)
> x=c(10,12,16,18,20,22)
>
> modelo<-lm(y~x)
> summary(modelo)
> library(MASS)
> Confint(modelo, level=0.95)
>
> Con esto tengo un modelo :
> > Confint(RegModel.1, level=0.95)
> Estimate 2.5 % 97.5 %
> (Intercept) 0.6100000 -6.8296312 8.049631
> x 0.9328571 0.4919135 1.373801
>
> Es decir, el intervalo de confianza de la respuesta y en función de los
> valores x, sería de la forma (ymin,ymax), siendo:
> ymin = -6,82 + 0,49*x
> ymax = 8,05 + 1,37*x
>
> Es esto correcto???
>
> Gracias
> Jesús
>
>
>
>
> [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
[[alternative HTML version deleted]]
Más información sobre la lista de distribución R-help-es