[R-es] maximo global e intervalo de confianza

PabloEmilio.Verde en uni-duesseldorf.de PabloEmilio.Verde en uni-duesseldorf.de
Vie Jul 17 23:46:48 CEST 2009


Hola Nahuel,

Como esta mi querida Mar del Plata?

Bueno, la pregunta que hiciste es muy complicada y Jorge te ha
dado una buena respuesta.

Sin embargo, el uso de metodos bootstrap con problemas de maximos
u otros estadisticos que estan en el borde del espacio parametrico
es problematica. Uno de los contra-ejemplos clasicos del uso del
bootstrap es justamente construir intervalos de confianza del maximo.

No te mando la respuesta a tu pregunta, pero te mando un par de
ejemplos simples para ilustrarte el problema que podes tener:

# comienzo ...
set.seed(12345)

x1 <- runif(100)*5 # uniform con maximo = 5
x2 <- rnorm(100)   # normal con maximo = Inf

theta.star <- function(dat, Ind)
{
 max(dat[Ind])
}

max(x1)
boot1 <- boot(x1, theta.star, R = 1999)
plot(boot1, jack = T)
boot.ci(boot1)

max(x2)
boot2 <- boot(x2, theta.star, R = 1999)
plot(boot2, jack = T)
boot.ci(boot2)
# fin ...

Saludos,

Pablo


> Buenas tardes Nahuel,
> Una manera de calcularlo es via bootstrap. Solo necesitas una función que
> estime el máximo global de la 2a derivada y luego, utilizando la libreria
> boot, reproducir los resultados usando muestras aleatorias de tus datos
> originales.
>
> Este es un ejemplo:
>
> # Librerias
> require(pspline)
> require(boot)
>
> # Datos
> data(cars)
> attach(cars)
>
> # Función para calcular el máximo de la 2a derivada
> # Los argumentos son d (los datos)
> #                               index (para el boostrap SOLAMENTE)
> #                               k (el orden de la derivada) -- por defecto k
> = 2
>
> foo <- function(d, index, k = 2){
>           ds <- d[index,]
>           fit <- sm.spline(ds$speed, ds$dist)  # debes cambiar esta linea
> con tus datos
>           deriva <- predict(fit, deriv = k)
>           max(deriva$y)
>           }
>
> # Bootstrap
> res <- boot(cars, foo, R = 1000)  # esto toma un poco de tiempo
> plot(res)
> quantile(res$t, probs = c(0.025, 0.975))
> #  2.5%     97.5%
> #  64.71799 103.82926
>
> Para mayor información puedes leer la ayuda en ?boot una vez cargas la
> libreria boot.
>
> Espero sea de utilidad,
>
> Jorge Ivan Velez
>
>
> 2009/7/10 Nahuel farias <>
>
> > Hola a todos
> >
> > Alguien sabe cómo puedo determinar maximos globalespara la derivada
> > 2da de un "sm.spline" y un intervalo de confianza para ese maximo?
> >
> > Muchas gracias
> >
> > --
> > Lic. Nahuel E. Farias
> > Laboratorio de Invertebrados
> > Departamento de Biología
> > Facultad de Ciencias Exactas y Naturales,
> > Universidad Nacional de Mar Del Plata, CC 1245, CP 7600 Mar del Plata,
> > Argentina
> > Tel. + 54 - 0223 - 4752426 (Int. 466)
> >
> > _______________________________________________
> > 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