Estoy buscando una curva método de ajuste que proporciona los límites superior e inferior de los intervalos de confianza para el conjunto de la curva, que se ilustra en la figura a continuación. Suponga que el rojo sólido de la curva es la equipada con la curva de la que podría ser obtenido por el ordinario de la curva de herramientas de montaje. Sin embargo, la curva de alguna manera predice los valores esperados. Estoy buscando la discontinua curvas que son los límites superior e inferior de los armarios de la curva. Tenga en cuenta que mis datos han heterocedasticidad propiedades. Cualquier ayuda es muy apreciada.
Respuesta
¿Demasiados anuncios?
Usted puede utilizar el geom_smooth() la función de la ggplot2 biblioteca en R.
He aquí un ejemplo:
x = seq(0, 100, by=0.5)
y = sqrt(x)
y = y + rnorm(n = length(y),mean = 0,sd = 3)
df = data.frame(cbind(x,y))
require(ggplot2)
ggplot(data = df, aes(x = x,y = y)) + geom_smooth()
Este es el resultado:
geom_smooth llama a un ajuste de curvas método llamado de loess, que hace de local de regresión. Con las opciones por defecto, para cada x, se considera un barrio que contiene 75% de los puntos de datos y se ajusta a una ecuación cuadrática utilizando mínimos cuadrados ponderados. Se asume que los errores están distribuidos normalmente y calcula los intervalos de confianza como se describe en las páginas 44-46 de http://www.netlib.org/a/cloess.pdf
Aquí está la documentación de loess: https://stat.ethz.ch/R-manual/R-devel/library/stats/html/loess.html
