An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130915/e3cf9c49/attachment.pl>
[R-es] b
11 messages · Dr. José A. Betancourt B., Jorge I Velez, eric
An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130916/53df4c39/attachment.pl>
An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130915/a39433a4/attachment.pl>
An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130916/cb11e234/attachment.pl>
An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130915/bf5efc9a/attachment.pl>
An embedded and charset-unspecified text was scrubbed... Name: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20130915/47284ac1/attachment.pl>
-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Estimado Jose, quiza podrias adjuntar un archivo con algunos datos para probar nosotros, de modo que tu puedas contrastar tu resultado. Saludos, Eric.
On 09/15/2013 06:14 PM, Dr. José A. Betancourt B. wrote:
LO HICE DE ESTA MANERA Y FUNCIONÓ
a_HTA <- (estudiantes1[,c(1)])
IMC <- (estudiantes1[,c(2)])
estudiantes4 <- data.frame( IMC, a_HTA = factor(c(rep("si", 29),
rep("no",30))))
wt <- wilcox_test(IMC ~ a_HTA, data = estudiantes4,
distribution = "exact", conf.int = TRUE)
print(wt)
De: Jorge I Velez [mailto:jorgeivanvelez en gmail.com]
Enviado el: Sunday, September 15, 2013 4:37 PM
Para: Dr. José A. Betancourt B.
CC: R-help-es
Asunto: Re: [R-es] b
Eso no fue lo que dijiste en tu primer mensaje. Dale una mirada a
?wilcox.test --JIV
2013/9/16 Dr. José A. Betancourt B. <jbetancourt en iscmc.cmw.sld.cu>
Lo que deseamos realizar es una prueba no paramétrica para demostrar que el
IMC encontrado difiere significtivamente entre aquellos que tienen a_HTA
(entecedentas HTA) y los que no lo tienen ( en SPSS de forma fácil utilizo
varios métodos ¿Cuál es la sintaxis para hacerlo posible en R?
a_HTA
IMC
si
26
si
25
si
26
si
21
si
25
si
23
si
28
si
24
si
27
si
22
si
23
si
23
si
27
si
22
si
23
si
25
si
22
si
22
si
22
si
25
si
23
si
23
si
28
si
27
si
27
si
23
si
24
si
26
si
21
no
19
no
18
no
17
no
16
no
18
no
19
no
18
no
17
no
18
no
17
no
18
no
17
no
17
no
18
no
18
no
19
no
17
no
17
no
16
no
18
no
16
no
16
no
16
no
18
no
18
no
19
no
20
no
16
no
18
no
18
De: Jorge I Velez [mailto:jorgeivanvelez en gmail.com]
Enviado el: Sunday, September 15, 2013 10:30 AM
Para: Dr. José A. Betancourt B.
CC: R-help-es
Asunto: Re: [R-es] b
Jose,
Es dificil sin un ejemplo, pero creo que
as.numeric(as.character(lavariable)))
deberia funcionar.
Saludos,
Jorge.-
2013/9/16 Dr. José A. Betancourt B. <jbetancourt en iscmc.cmw.sld.cu>
En varias ocasiones hago una tabla y el tipo de análisis me dice que tiene
que ser numérica la variable, entonces escribo
As.numeric ( la variable) y no me funciona
¿Cuál es la sintaxis adecuada? Y dónde está esa ayuda en R?
--
Este mensaje le ha llegado mediante el servicio de correo electronico que
ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema
Nacional de Salud. La persona que envia este correo asume el compromiso de
usar el servicio a tales fines y cumplir con las regulaciones establecidas
Infomed: http://www.sld.cu/
[[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 -- Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La persona que envia este correo asume el compromiso de usar el servicio a tales fines y cumplir con las regulaciones establecidas Infomed: http://www.sld.cu/ [[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
-----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.12 (GNU/Linux) Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/ iQEcBAEBAgAGBQJSNkQAAAoJEDoUMoeY1RuMEikH/A+apCWC8XwN3gIAZHuyy43i 3MDRWBTzo7d23Ti2pAb+sT6wZVNJiUTZ9kYmihvL8e/MhkJslNtCSjIfjagaVfqk pUTgY3dqrHyYVSuR6g8rYY7juhM6bSbMEw5LrqcrTFzJwOOtSadffOP74G0ApBLg yBKiyMjVBEJVY53aVMPAxaRn4FV+5WbFsSPcWPmyd4YqA891NidaS1KC20Fpn6wq nrrxrHwyxrevfAIV5Y7cmWBROcBB2mwUHeEdf4zaFbQtyqZgl237kGesKc1vSYJO GDyPfTFHVqWZBqWfHadjDt38zxhEv58JCOrhnBt6sKafOS74KbGT6/ZoUhjsJ5w= =uZT3 -----END PGP SIGNATURE-----
Estimado José, por lo que estuve mirando aqui https://en.wikipedia.org/wiki/Wilcoxon_signed-rank_test, en ?wilcox.test y en libros de papel ... los datos deben cumplir ciertas condiciones, quiza la mas restrictiva es que sean pareados ... si no lo son, quiza esta variante del test es mas apropiada: https://en.wikipedia.org/wiki/Mann-Whitney-Wilcoxon_test la que en R se aplica usando el mismo comando wilcox.test(), pero activando la opcion paired = FALSE. Asumiendo que sus datos NO son pareados y que se cumplen las otras condiciones del test, yo lo aplique de la siguiente forma y obtuve el siguiente resultado (usando el archivo que nos envio): tmp <- read.table("ARCHIVO.csv", header=TRUE, sep=",", dec=".", na.strings = "NA") wt <- wilcox.test(IMC ~ a_HTA, data=tmp, paired = FALSE, alternative = c("two.sided"), exact=TRUE) wt 1> wt Wilcoxon rank sum test with continuity correction data: IMC by a_HTA W = 0, p-value = 3.233e-11 alternative hypothesis: true location shift is not equal to 0 de modo que se rechaza ls hipotesis de igualdad de MEDIANAS (NO de medias). Espero que mis comentarios le sirvan, saludos cordiales, Eric Concha M. pd1. dayand, creo que este es tu mismo caso, slds, eric. pd2. los estadisticos que me corrigan si me equivoco por favor, slds.
On 09/15/2013 05:19 PM, Dr. José A. Betancourt B. wrote:
deseamos realizar es una prueba no paramétrica para demostrar que el IMC encontrado difiere significtivamente entre aquellos que tienen a_HTA (entecedentas HTA) y los que no lo tienen ( en SPSS de forma fácil utilizo varios métodos ¿Cuál es la sintaxis para hacerlo posible en R? wilcox.test requiere definir x y y ¿ como es la sintaxis en R? a_HTA IMC si 26 si 25 si 26 si 21 si 25 si 23 si 28 si 24 si 27 si 22 si 23 si 23 si 27 si 22 si 23 si 25 si 22 si 22 si 22 si 25 si 23 si 23 si 28 si 27 si 27 si 23 si 24 si 26 si 21 no 19 no 18 no 17 no 16 no 18 no 19 no 18 no 17 no 18 no 17 no 18 no 17 no 17 no 18 no 18 no 19 no 17 no 17 no 16 no 18 no 16 no 16 no 16 no 18 no 18 no 19 no 20 no 16 no 18 no 18 -- Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La persona que envia este correo asume el compromiso de usar el servicio a tales fines y cumplir con las regulaciones establecidas Infomed: http://www.sld.cu/ [[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
Gracias, tiene mucha lógica lo que dices, en efecto no son pareados y es adecuado especificarlo, lo voy a probar de inmediato gracias -----Mensaje original----- De: neo [mailto:ericconchamunoz en gmail.com] Enviado el: Monday, September 16, 2013 12:14 AM Para: "Dr. José A. Betancourt B." CC: 'R-help-es'; Dayand Toledo Asunto: Re: [R-es] b Estimado José, por lo que estuve mirando aqui https://en.wikipedia.org/wiki/Wilcoxon_signed-rank_test, en ?wilcox.test y en libros de papel ... los datos deben cumplir ciertas condiciones, quiza la mas restrictiva es que sean pareados ... si no lo son, quiza esta variante del test es mas apropiada: https://en.wikipedia.org/wiki/Mann-Whitney-Wilcoxon_test la que en R se aplica usando el mismo comando wilcox.test(), pero activando la opcion paired = FALSE. Asumiendo que sus datos NO son pareados y que se cumplen las otras condiciones del test, yo lo aplique de la siguiente forma y obtuve el siguiente resultado (usando el archivo que nos envio): tmp <- read.table("ARCHIVO.csv", header=TRUE, sep=",", dec=".", na.strings = "NA") wt <- wilcox.test(IMC ~ a_HTA, data=tmp, paired = FALSE, alternative = c("two.sided"), exact=TRUE) wt 1> wt Wilcoxon rank sum test with continuity correction data: IMC by a_HTA W = 0, p-value = 3.233e-11 alternative hypothesis: true location shift is not equal to 0 de modo que se rechaza ls hipotesis de igualdad de MEDIANAS (NO de medias). Espero que mis comentarios le sirvan, saludos cordiales, Eric Concha M. pd1. dayand, creo que este es tu mismo caso, slds, eric. pd2. los estadisticos que me corrigan si me equivoco por favor, slds.
On 09/15/2013 05:19 PM, Dr. José A. Betancourt B. wrote:
deseamos realizar es una prueba no paramétrica para demostrar que el IMC encontrado difiere significtivamente entre aquellos que tienen a_HTA (entecedentas HTA) y los que no lo tienen ( en SPSS de forma fácil
utilizo
varios métodos ¿Cuál es la sintaxis para hacerlo posible en R? wilcox.test requiere definir x y y ¿ como es la sintaxis en R? a_HTA IMC si 26 si 25 si 26 si 21 si 25 si 23 si 28 si 24 si 27 si 22 si 23 si 23 si 27 si 22 si 23 si 25 si 22 si 22 si 22 si 25 si 23 si 23 si 28 si 27 si 27 si 23 si 24 si 26 si 21 no 19 no 18 no 17 no 16 no 18 no 19 no 18 no 17 no 18 no 17 no 18 no 17 no 17 no 18 no 18 no 19 no 17 no 17 no 16 no 18 no 16 no 16 no 16 no 18 no 18 no 19 no 20 no 16 no 18 no 18 -- Este mensaje le ha llegado mediante el servicio de correo electronico que
ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La persona que envia este correo asume el compromiso de usar el servicio a tales fines y cumplir con las regulaciones establecidas
Infomed: http://www.sld.cu/ [[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
-- Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La persona que envia este correo asume el compromiso de usar el servicio a tales fines y cumplir con las regulaciones establecidas Infomed: http://www.sld.cu/
De la manera que me aconsejas no funciona y da este mansaje
Warning message:
In wilcox.test.default(x = c(19L, 18L, 17L, 16L, 18L, 19L, 18L, :
cannot compute exact p-value with ties
como aconsejaste
tmp <- read.table("estudiantes1.csv", header=TRUE, sep=",", dec=".",
na.strings = "NA")
wt <- wilcox.test(IMC ~ a_HTA, data=tmp, paired = FALSE, alternative =
c("two.sided"), exact=TRUE)
wt
como me funciona es asi
estudiantes1<-read.csv("estudiantes1.csv",header=TRUE, sep=",", dec=".")
a_HTA <- (estudiantes1[,c(1)])
IMC <- (estudiantes1[,c(2)])
estudiantes4 <- data.frame( IMC, a_HTA = factor(c(rep("si", 29),
rep("no",30))))
wt <- wilcox_test(IMC ~ a_HTA, data = estudiantes4,
distribution = "exact", conf.int = TRUE)
print(wt)
en la ayuda dice
Independent Two- and K-Sample Location Tests
Description
Testing the equality of the distributions of a numeric response in two or
more independent groups against shift alternatives.
-----Mensaje original-----
De: neo [mailto:ericconchamunoz en gmail.com]
Enviado el: Monday, September 16, 2013 12:14 AM
Para: "Dr. José A. Betancourt B."
CC: 'R-help-es'; Dayand Toledo
Asunto: Re: [R-es] b
Estimado José, por lo que estuve mirando aqui
https://en.wikipedia.org/wiki/Wilcoxon_signed-rank_test, en ?wilcox.test
y en libros de papel ... los datos deben cumplir ciertas condiciones,
quiza la mas restrictiva es que sean pareados ... si no lo son, quiza
esta variante del test es mas apropiada:
https://en.wikipedia.org/wiki/Mann-Whitney-Wilcoxon_test
la que en R se aplica usando el mismo comando wilcox.test(), pero
activando la opcion paired = FALSE. Asumiendo que sus datos NO son
pareados y que se cumplen las otras condiciones del test, yo lo aplique
de la siguiente forma y obtuve el siguiente resultado (usando el archivo
que nos envio):
tmp <- read.table("ARCHIVO.csv", header=TRUE, sep=",", dec=".",
na.strings = "NA")
wt <- wilcox.test(IMC ~ a_HTA, data=tmp, paired = FALSE, alternative =
c("two.sided"), exact=TRUE)
wt
1> wt
Wilcoxon rank sum test with continuity correction
data: IMC by a_HTA
W = 0, p-value = 3.233e-11
alternative hypothesis: true location shift is not equal to 0
de modo que se rechaza ls hipotesis de igualdad de MEDIANAS (NO de medias).
Espero que mis comentarios le sirvan, saludos cordiales,
Eric Concha M.
pd1. dayand, creo que este es tu mismo caso, slds, eric.
pd2. los estadisticos que me corrigan si me equivoco por favor, slds.
On 09/15/2013 05:19 PM, Dr. José A. Betancourt B. wrote:
deseamos realizar es una prueba no paramétrica para demostrar que el IMC encontrado difiere significtivamente entre aquellos que tienen a_HTA (entecedentas HTA) y los que no lo tienen ( en SPSS de forma fácil
utilizo
varios métodos ¿Cuál es la sintaxis para hacerlo posible en R? wilcox.test requiere definir x y y ¿ como es la sintaxis en R? a_HTA IMC si 26 si 25 si 26 si 21 si 25 si 23 si 28 si 24 si 27 si 22 si 23 si 23 si 27 si 22 si 23 si 25 si 22 si 22 si 22 si 25 si 23 si 23 si 28 si 27 si 27 si 23 si 24 si 26 si 21 no 19 no 18 no 17 no 16 no 18 no 19 no 18 no 17 no 18 no 17 no 18 no 17 no 17 no 18 no 18 no 19 no 17 no 17 no 16 no 18 no 16 no 16 no 16 no 18 no 18 no 19 no 20 no 16 no 18 no 18 -- Este mensaje le ha llegado mediante el servicio de correo electronico que
ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La persona que envia este correo asume el compromiso de usar el servicio a tales fines y cumplir con las regulaciones establecidas
Infomed: http://www.sld.cu/ [[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
-- Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La persona que envia este correo asume el compromiso de usar el servicio a tales fines y cumplir con las regulaciones establecidas Infomed: http://www.sld.cu/
Hola Jose, ese mensaje solo dice que no se pudo obtener un valor-p "exacto". Compara el valor que obtienes con tu metodo y pide ver WT usando los pasos que te propongo, aunque obtengas el error que mencionas. Creo que deberian ser el mismo valor-p. Saludos, Eric.
On 09/16/2013 07:24 AM, Dr. José A. Betancourt B. wrote:
De la manera que me aconsejas no funciona y da este mansaje
Warning message:
In wilcox.test.default(x = c(19L, 18L, 17L, 16L, 18L, 19L, 18L, :
cannot compute exact p-value with ties
como aconsejaste
tmp <- read.table("estudiantes1.csv", header=TRUE, sep=",", dec=".",
na.strings = "NA")
wt <- wilcox.test(IMC ~ a_HTA, data=tmp, paired = FALSE, alternative =
c("two.sided"), exact=TRUE)
wt
como me funciona es asi
estudiantes1<-read.csv("estudiantes1.csv",header=TRUE, sep=",", dec=".")
a_HTA <- (estudiantes1[,c(1)])
IMC <- (estudiantes1[,c(2)])
estudiantes4 <- data.frame( IMC, a_HTA = factor(c(rep("si", 29),
rep("no",30))))
wt <- wilcox_test(IMC ~ a_HTA, data = estudiantes4,
distribution = "exact", conf.int = TRUE)
print(wt)
en la ayuda dice
Independent Two- and K-Sample Location Tests
Description
Testing the equality of the distributions of a numeric response in two or
more independent groups against shift alternatives.
-----Mensaje original-----
De: neo [mailto:ericconchamunoz en gmail.com]
Enviado el: Monday, September 16, 2013 12:14 AM
Para: "Dr. José A. Betancourt B."
CC: 'R-help-es'; Dayand Toledo
Asunto: Re: [R-es] b
Estimado José, por lo que estuve mirando aqui
https://en.wikipedia.org/wiki/Wilcoxon_signed-rank_test, en ?wilcox.test
y en libros de papel ... los datos deben cumplir ciertas condiciones,
quiza la mas restrictiva es que sean pareados ... si no lo son, quiza
esta variante del test es mas apropiada:
https://en.wikipedia.org/wiki/Mann-Whitney-Wilcoxon_test
la que en R se aplica usando el mismo comando wilcox.test(), pero
activando la opcion paired = FALSE. Asumiendo que sus datos NO son
pareados y que se cumplen las otras condiciones del test, yo lo aplique
de la siguiente forma y obtuve el siguiente resultado (usando el archivo
que nos envio):
tmp <- read.table("ARCHIVO.csv", header=TRUE, sep=",", dec=".",
na.strings = "NA")
wt <- wilcox.test(IMC ~ a_HTA, data=tmp, paired = FALSE, alternative =
c("two.sided"), exact=TRUE)
wt
1> wt
Wilcoxon rank sum test with continuity correction
data: IMC by a_HTA
W = 0, p-value = 3.233e-11
alternative hypothesis: true location shift is not equal to 0
de modo que se rechaza ls hipotesis de igualdad de MEDIANAS (NO de medias).
Espero que mis comentarios le sirvan, saludos cordiales,
Eric Concha M.
pd1. dayand, creo que este es tu mismo caso, slds, eric.
pd2. los estadisticos que me corrigan si me equivoco por favor, slds.
On 09/15/2013 05:19 PM, Dr. José A. Betancourt B. wrote:
deseamos realizar es una prueba no paramétrica para demostrar que el IMC encontrado difiere significtivamente entre aquellos que tienen a_HTA (entecedentas HTA) y los que no lo tienen ( en SPSS de forma fácil
utilizo
varios métodos ¿Cuál es la sintaxis para hacerlo posible en R? wilcox.test requiere definir x y y ¿ como es la sintaxis en R? a_HTA IMC si 26 si 25 si 26 si 21 si 25 si 23 si 28 si 24 si 27 si 22 si 23 si 23 si 27 si 22 si 23 si 25 si 22 si 22 si 22 si 25 si 23 si 23 si 28 si 27 si 27 si 23 si 24 si 26 si 21 no 19 no 18 no 17 no 16 no 18 no 19 no 18 no 17 no 18 no 17 no 18 no 17 no 17 no 18 no 18 no 19 no 17 no 17 no 16 no 18 no 16 no 16 no 16 no 18 no 18 no 19 no 20 no 16 no 18 no 18 -- Este mensaje le ha llegado mediante el servicio de correo electronico que
ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La persona que envia este correo asume el compromiso de usar el servicio a tales fines y cumplir con las regulaciones establecidas
Infomed: http://www.sld.cu/ [[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
-- Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece Infomed para respaldar el cumplimiento de las misiones del Sistema Nacional de Salud. La persona que envia este correo asume el compromiso de usar el servicio a tales fines y cumplir con las regulaciones establecidas Infomed: http://www.sld.cu/