José Trujillo Carmona
2010-Nov-24 14:22 UTC
[R-es] Límites de confianza de la mediana en distribuciones simétricas
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]]
Carlos Ortega
2010-Nov-24 14:39 UTC
[R-es] Límites de confianza de la mediana en distribuciones simétricas
Hola,
No aparece porque no has especificado a la función wilcox.test que lo
calcule (mirar los parámetros de la función conf.int y conf.level):
---------------------------------------------
wilcox.test(iris$Sepal.Width, iris$MEDIANA,
+ alternative=''two.sided'',paired=TRUE, conf.int=T)
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
Incluso el test permite modificar el intervalo (conf.level):
wilcox.test(iris$Sepal.Width, iris$MEDIANA,
+ alternative=''two.sided'',paired=TRUE, conf.int=T,
conf.level=0.80)
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
80 percent confidence interval:
0.04994611 0.14996317
sample estimates:
(pseudo)median
0.09994806
Saludos,
Carlos Ortega
2010/11/24 José Trujillo Carmona <trujillo@unex.es>
> 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]]
Jorge Ivan Velez
2010-Nov-24 14:40 UTC
[R-es] Límites de confianza de la mediana en distribuciones simétricas
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]]
David Méndez
2010-Nov-25 10:56 UTC
[R-es] Límites de confianza de la mediana en distribuciones simétricas
Enviado desde mi HTC ----- Mensaje original ----- De: Oscar PerpiñanLamigueiro <oscar.perpinan en upm.es> Enviado: Miércoles, 24 de Noviembre de 2010 04:09 p.m. Para: trujillo en unex.es CC: r-help-es en r-project.org Asunto: Re: [R-es] Límites de confianza de la mediana en distribuciones simétricas 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[No se incluye el mensaje original entero]