[R-es] Límites de confianza de la mediana en distribuciones simétricas

[Oscar Perpiñan Lamigueiro]

Hola,

Por si te interesa una generalización robusta de este test, el método de Mee (por ejemplo) (http://books.google.com/books?id=_tAJr4ooOM8C&lpg=PA174&ots=mfFzPdT7Pk&dq=mee%20wilcoxon&hl=es&pg=PA174#v=onepage&q=mee%20wilcoxon&f=false) está incluido dentro del paquete WRS (por ahora sólo en R-Forge): http://r-forge.r-project.org/projects/wrs/

Saludos.

Oscar.


-- 
Oscar Perpiñán Lamigueiro
Profesor Ayudante Doctor
Dpto. de Ingeniería Eléctrica
EUITI-UPM
On Wed, 24 Nov 2010 09:40:22 -0500
Jorge Ivan Velez <jorgeivanvelez en gmail.com> wrote:

> Buenos dias Jose,
> 
> La function wilcox.test si proporciona el IC que necesitas. De ?wilcox.test:
> 
> conf.int:   a confidence interval for the location parameter. (Only present
> if argument conf.int = TRUE)
> 
> En tu caso, seria algo como
> 
> 
>   wilcox.test(iris$Sepal.Width, iris$MEDIANA, alternative = 'two.sided',
> paired = TRUE, conf.int = TRUE)
> 
> Wilcoxon signed rank test with continuity correction
> 
> data:  iris$Sepal.Width and iris$MEDIANA
> V = 6817, p-value = 0.03015
> alternative hypothesis: true location shift is not equal to 0
> 95 percent confidence interval:
>  0.04999865 0.15002894
> sample estimates:
> (pseudo)median
>     0.09994806
> 
> Otra alternativa es usar bootstrap:
> 
>    require(boot)
> 
> b <- boot(iris$Sepal.Width, function(x, d) median(x[d]), R = 1000)
> 
> quantile(b$t - 2.95, probs = c(0.025, 0.975))
> #   2.5% 97.5%
> #   0.05  0.15
> 
> Un saludo,
> Jorge
> 
> 
> 2010/11/24 José Trujillo Carmona <>
> 
> > Por si alguno pudiera ayudarme.
> >
> > Al realizar el t.test para una muestra, junto con el valor de t y el
> > p-valor, la función proporciona la estimación de la media y su INTERVALO
> > DE CONFIANZA.
> >
> > Desde el punto de vista de la estadística de rangos esto se puede hacer
> > mediante:
> >
> >  > iris$MEDIANA <- with(iris, 2.95)
> >
> >  > median(iris$Sepal.Width - iris$MEDIANA, na.rm=TRUE) # median difference
> > [1] 0.05
> >
> >  > wilcox.test(iris$Sepal.Width, iris$MEDIANA,
> > alternative='two.sided',paired=TRUE)
> >
> >    Wilcoxon signed rank test with continuity correction
> >
> > data:  iris$Sepal.Width and iris$MEDIANA
> > V = 6817, p-value = 0.03015
> > alternative hypothesis: true location shift is not equal to 0
> >
> >
> > Ahora bien, como se puede ver la función wilcox.test no proporciona los
> > límites de confianza de la mediana que se puede obtener a partir del
> > test de Wilcoxon (v.g. en "Applied Nonparametric Statistics" W. W.
> > Daniel, pp.49-56)
> >
> > ¿Alguien conoce algún paquete que contenga la función que busco?
> >
> > Muchas gracias.
> >
> >
> > --
> > _____---^---_____
> >
> > Univ. de Extremadura
> > Dept. Matemáticas.
> > Despacho B29
> > Tf: + 34 924 289 300
> > Ext. 86823
> >
> >
> >        [[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]]
>