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@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

	[[alternative HTML version deleted]]

