Hola
respecto a la duda de las medidas repetidas yo diría lo siguiente:
el valo de p que observas es casi cero. Si te fijas hay tres asteriscos, que es
lo que indica, o sea, una p mínima, próxima a 0; lo que todos queremos
encontrar. Dos asteristcos y un asterisco también demuestran significación, pero
menor. Se trata de una especie de ayuda visual para localizar mejor el asunto.
Cuando he obtenido ese tipo de números informo de una p de .000
Lo que indicas de subject/error etc se debe, a mi entender, a cuando vas a hacer
medidas repetidas o intrasujetos. El sistema calcula por defecto intersujetos,
de modo que ahí debes especificar qué factor es intrasujetos. En un experiment
puedes tener, por ejemplo tres intersujetos y uno intrasujetos, pues entonces :
resultado=aov(puntuación~factor1*factor2+Error(subject/factor2), data =a)
donde el "factor 2 " es intrasujetos.
Respecto de Tukey, entiendo que para hacer comparaciones planeadas, yo he sudado
tinta en ello. Es sencillo pese a todo, pero por una cosa u otra nunca
encontraba o leía lo que necesitaba. Lo que puedes mirar es lo siguiente.
pairwise.t.test(tr,pat,p.adj="bonferroni") Compara un solo grupo
rt=tiempos respuesta (variable dependiente)
pat=patología (variable independiente)
utilizas pairwise t.test para comparar un solo grupo y lo ajustas por bonferroni
(bastante conservador, como sabrás mejor que yo).
Para más grupos:
> w=TukeyHSD(aov(tr~pat)) pat=patología rt= latencias respuesta
> w
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = tr ~ pat)
$pat
diff lwr upr p adj
tres-dos 53.7000 15.84591 91.55409 0.0026712
uno-dos -106.6917 -144.54576 -68.83758 0.0000000
uno-tres -160.3917 -198.24576 -122.53758 0.0000000
Espero haber ayudado,
saludos
Miguel Lázaro
--- El lun, 21/5/12, r-help-es-request@r-project.org
<r-help-es-request@r-project.org> escribió:
De: r-help-es-request@r-project.org <r-help-es-request@r-project.org>
Asunto: Resumen de R-help-es, Vol 39, Envío 37
Para: r-help-es@r-project.org
Fecha: lunes, 21 de mayo, 2012 17:01
Envíe los mensajes para la lista R-help-es a
r-help-es@r-project.org
Para subscribirse o anular su subscripción a través de la WEB
https://stat.ethz.ch/mailman/listinfo/r-help-es
O por correo electrónico, enviando un mensaje con el texto "help" en
el asunto (subject) o en el cuerpo a:
r-help-es-request@r-project.org
Puede contactar con el responsable de la lista escribiendo a:
r-help-es-owner@r-project.org
Si responde a algún contenido de este mensaje, por favor, edite la
linea del asunto (subject) para que el texto sea mas especifico que:
"Re: Contents of R-help-es digest...". Además, por favor, incluya en
la respuesta sólo aquellas partes del mensaje a las que está
respondiendo.
Asuntos del día:
1. Re: Contactos en Madrid-ayuda con instalación (Juan Jose)
2. Re: Contactos en Madrid-ayuda con instalación (Juan Jose)
3. Anova medidas repetidas (Javier Villacampa González)
4. ayuda con graficos de la distribucion normal (tono ferri vidal)
5. Re: ayuda con graficos de la distribucion normal
(Carlos J. Gil Bellosta )
6. help: Anova medidas repetidas (Javier Villacampa González)
----------------------------------------------------------------------
Message: 1
Date: Mon, 21 May 2012 15:39:10 +0200
From: Juan Jose <jjfgarcia@renta4.es>
To: <miguel.angel.rodriguez.muinos@sergas.es>,
<r-help-es@r-project.org>
Subject: Re: [R-es] Contactos en Madrid-ayuda con instalación
Message-ID: <020701cd3757$19f34c30$4dd9e490$@es>
Content-Type: text/plain
Sí, lo hemos intentado instalar esta mañana como administrador, pero me daba
el error que te comento. He abierto R y en la línea de comandos he tecleado
la línea install.packages("RExcelInstaller", "rcom",
"rsproxy") pero no
dejaba escribir
De: miguel.angel.rodriguez.muinos@sergas.es
[mailto:miguel.angel.rodriguez.muinos@sergas.es]
Enviado el: lunes, 21 de mayo de 2012 15:18
Para: jjfgarcia@renta4.es; r-help-es@r-project.org
Asunto: RE: [R-es] Contactos en Madrid-ayuda con instalación
Es que lo que necesita permisos de administrador es precisamente esa línea.
Debes hacer la instalación de RExcelInstaller como Admin.
Un Saludo,
Miguel.
De: Juan Jose [mailto:jjfgarcia@renta4.es]
Enviado el: lunes, 21 de mayo de 2012 15:15
Para: Rodríguez Muíños, Miguel Ángel; r-help-es@r-project.org
Asunto: RE: [R-es] Contactos en Madrid-ayuda con instalación
Hola Miguel,
Muchas gracias. Sí, me lo han hecho desde informática como administrador,
hemos instalado el R-Project, descargado las librearías, instalado el Rexcel
y Statconn como administrador…pero tengo los problemas que menciono. Incluso
hemos intentado hacer funcionar la macro de Visual desde administrador
también pero nada. Pero está bien que me preguntes esto porque estoy seguro
que alguna cosa básica estoy haciendo mal.
Lo último que intenté es desde RProyect correr estas líneas de código que me
explicaban en una guía, pero me daba error de escritura en la primera
línea….gracias
> install.packages("RExcelInstaller", "rcom",
"rsproxy")
_____
Nota: A información contida nesta mensaxe e os seus posibles documentos
adxuntos é privada e confidencial e está dirixida únicamente ó seu
destinatario/a. Se vostede non é o/a destinatario/a orixinal desta mensaxe,
por favor elimínea. A distribución ou copia desta mensaxe non está
autorizada.
Nota: La información contenida en este mensaje y sus posibles documentos
adjuntos es privada y confidencial y está dirigida únicamente a su
destinatario/a. Si usted no es el/la destinatario/a original de este
mensaje, por favor elimínelo. La distribución o copia de este mensaje no
está autorizada.
See more languages: http://www.sergas.es/aviso_confidencialidad.htm
[[alternative HTML version deleted]]
------------------------------
Message: 2
Date: Mon, 21 May 2012 15:42:47 +0200
From: Juan Jose <jjfgarcia@renta4.es>
To: "''Carlos Ortega''"
<cof@qualityexcellence.es>
Cc: r-help-es@r-project.org
Subject: Re: [R-es] Contactos en Madrid-ayuda con instalación
Message-ID: <020c01cd3757$9b66cb20$d2346160$@es>
Content-Type: text/plain
Hola Carlos
Gracias. Sí, me refería a R tal cual.
Lo hemos instalado bajo Windows XP, versión 32bits y Excel 2003. La versión
de R es la última que está en la página oficial, la .15. Voy a preguntar lo
del sistema operativo pero es el Windows xp, no se si será 32 o 64. Estoy en
ello.
gracias
De: Carlos Ortega [mailto:cof@qualityexcellence.es]
Enviado el: lunes, 21 de mayo de 2012 15:17
Para: Juan Jose
CC: r-help-es@r-project.org
Asunto: Re: [R-es] Contactos en Madrid-ayuda con instalación
Hola,
¿Cuándo dices RPROJECT te refieres al R tal cual?
Es que como tal, no veo que exista ninguna librería que vincule el R con el
Microsoft-Project, como sí existe con la otra que mencionas de RExcel
(aunque realmente se llama RExcelInstaller).
Además del tema que te ha comentado de los permisos de administrador,
especialmente en los directorios de Windows donde se instalarán .dlls, es
importante conocer sobre qué versión de R estás queriendo instalar todo esto
y sobre qué versión de Windows (si tienes un Sistema Operativo de x64 a lo
mejor estas librerías no lo soportan).
Saludos,
Carlos Ortega
www.qualityexcellence.es
El 21 de mayo de 2012 14:52, Juan Jose <jjfgarcia@renta4.es> escribió:
Hola a todos,
Mi nombre es Juanjo. Llevo un tiempo intentado hacer funcionar un programa
que necesita de la instalación de RPROJECT y REXCEL. He mirado varias guías
e instalado los programas pero obtengo varios errores.
El principal es que tras instalar RProject, Rexcel y statcon y al intentar
correr esta macro desde Visual Basic…
Sub Start()
RInterface.StartRServer
RInterface.StopRServer
ActiveWorkbook.Save
End
… obtengo el siguiente error “Error de compilación No puedo encontrar
proyecto o biblioteca”. He revisado los addin de VSB y todo parece correcto.
Revisando la web vi que desde RProject tenía que correr este código:
install.packages("RExcelInstaller", "rcom",
"rsproxy") y pensé que esto
podía ser la solución. Pero resulta que me dice que no los puede instalar
porque el archivo no es writeable.
En definitiva que llevo varias semanas con este tema y por desesperación he
desinstalado todo.
Os quería por favor pedir que o bien alguien me pudiera echar una mano a
través del foro (bien explicando exactamente los pasos s seguir para
instalar todo correctamente) o que si hay un usuario que viva en Madrid se
pudiera acercar por renta4 (paseo de la Habana) e intentar ayudarnos in situ
con el problema. Si este fuera el caso el banco se compromete a pagar el
tiempo dedicado como asesoría externa o invitar a una buena comida en la
zona.
Muchas gracias
Juanjo
Juan José Fernández, CFA | Institutional Sales Trading
T: +34 91 384 85 50 <tel:%2B34%2091%20384%2085%2050> | M: +34 66 987 91
<tel:%2B34%2066%20987%2091%2087> 87
<mailto:dgonzalez@renta4.es> jjfgarcia@renta4.es |
<http://www.renta4.com>
www.renta4.com
Renta4 Equities
Paseo de la Habana, 74
28036 Madrid | Spain
<http://www.renta4.com/> ¡Error! Nombre de archivo no especificado.
La información contenida en este mensaje y los documentos adjuntos es
confidencial y puede ser privilegiada. Si usted no es el destinatario por
favor notifique inmediatamente al remitente, absteniéndose de comunicar el
contenido o reproducirlo por cualquier medio.
The information contained in this e-mail and any attachments are
confidential and may also be privileged. If you are not the named recipient,
please notify the sender inmediatly and do not disclose the contents to
another person, use it for any purpose, or store or copy the informatión in
any way.
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]
------------------------------
Message: 3
Date: Mon, 21 May 2012 15:58:01 +0200
From: Javier Villacampa González
<javier.villacampa.gonzalez@gmail.com>
To: r-help-es@r-project.org
Subject: [R-es] Anova medidas repetidas
Message-ID:
<CAOZ+RxFBj++Vrwm17Z2ur5PkV4vZOe003j95PKaq+DNDVyJ1cg@mail.gmail.com>
Content-Type: text/plain
Hola buenas,
lo primero una vez muchas gracias a todos por todas las veces que me habeis
ayudado.
Llevo días intentando utilizar las tablas ANOVA de medidas repetidas y he
encontrado multiple información, pero la verdad, no me aclaro.
He leido sobre la función ezANOVA o utilizando al aov con Error.
Los problema son:
¿Qué no me termino de aclarar con que significan las tablas. ?Vease
Anova1<- aov(RT ~ Prueba + Error(Subject) , datos) #+
Error(Subject/Prueba)> summary(Anova1)
Error: Subject
Df Sum Sq Mean Sq F value Pr(>F)
Prueba 2 40181770 20090885 18.893 4.114e-07 ***
Residuals 61 64866112 1063379
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 17206 356542956 20722
Que signifca este p-value exactamente, por más que leo solo me lio más
Además no entiendo cuando se emplea Error(Subject) y cuando
Error(Subject/Prueba). En mi caso los sujetos solo hacen una prueba
repeteidas veces así que creo que es Error(Subject)
¿Como se puede calcular las diferencias de Tukey para este tipo de objetos?
Tengo una solución que es programarlo, pero me gustaría que fuese un
paquete. Esté en esta página
http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/
Si teneis alguna lectura que recomendarme estaré gustoso de leerme. Soy
matemático así que las matemáticas no son un problema. Lo malo es
aplicarlas, en la carrera no se entra en este tipo de modelos.
Por último os dejo los datos por si no me he explicado bien.
*Los datos son:*
*RT: *variable continua de tiempos de reacción
*Prueba: *Tres tipos de pruebas neuropsicológicas
*Subject:* Participantes del experimento que hacen durante varios trial la
misma prueba.
head(datos) Subject Prueba card tri cat switch stim isi cod1 resp
Sorp RT Err SorpA SP Aciertos ratioA
1 s00 ngc plus 9 C plus 9 1523 90 0 2.000000
NA 1 0.000000 s00 ngc 354 0.924282
2 s00 ngc bb 10 C trg 1 1543 31 2 1.321928
1116.6 1 2.000000 s00 ngc 354 0.924282
3 s00 ngc rb 11 C trg 2 1534 32 1 1.584963
389.3 1 1.321928 s00 ngc 354 0.924282
4 s00 ngc rb 12 C trg 3 1552 33 1 1.222392
340.4 1 1.584963 s00 ngc 354 0.924282
5 s00 ngc rb 13 C trg 4 1541 34 1 1.000000
279.9 1 1.222392 s00 ngc 354 0.924282
6 s00 ngc bs 14 C trg 5 1545 35 2 2.169925
350.7 1 1.000000 s00 ngc 354 0.924282
tail(datos)
Subject Prueba card tri cat switch stim isi cod1 resp Sorp
RT Err SorpA SP Aciertos ratioA
24507 s74 Swt bs 585 Fo trg 3 1509 13 2 3.902135
374.1 1 3.903184 s74 Swt 330 0.8616188
24508 s74 Swt bs 586 Fo trg 4 1519 14 2 3.901116
339.1 1 3.902135 s74 Swt 330 0.8616188
24509 s74 Swt bb 587 Fo trg 5 1556 15 1 2.888630
870.1 1 3.901116 s74 Swt 330 0.8616188
24510 s74 Swt cros 588 C cros 8 1531 80 0 9.004247
NA 1 2.888630 s74 Swt 330 0.8616188
24511 s74 Swt rs 589 C trg 1 1543 11 1 3.136630
1013.9 1 9.004247 s74 Swt 330 0.8616188
24512 s74 Swt rb 590 C trg 2 1557 12 1 4.134127
983.4 1 3.136630 s74 Swt 330 0.8616188
Muchas gracias por adelantado.Un saludo.
Javier
[[alternative HTML version deleted]]
------------------------------
Message: 4
Date: Mon, 21 May 2012 14:54:53 +0000
From: tono ferri vidal <tonoferri_vidal@hotmail.com>
To: R-help <r-help-es@r-project.org>
Subject: [R-es] ayuda con graficos de la distribucion normal
Message-ID: <SNT134-W280D6B62895A60F4052F0C961D0@phx.gbl>
Content-Type: text/plain
Hola
Alguien me puede ayudar? Necesito graficar 2 funciones de
densidad "normales" superpuestas en un unico grafico.
una densidad es de media 0 y desviacion 1,
la otra es de media 3 y desviacion 1.
saludos
[[alternative HTML version deleted]]
------------------------------
Message: 5
Date: Mon, 21 May 2012 16:56:59 +0200
From: "Carlos J. Gil Bellosta " <cgb@datanalytics.com>
To: tono ferri vidal <tonoferri_vidal@hotmail.com>
Cc: R-help <r-help-es@r-project.org>
Subject: Re: [R-es] ayuda con graficos de la distribucion normal
Message-ID:
<CADg83efU0CVtJjPSVct+sV87DNR5fvm6UOK2HqtaKKN5ZprQYg@mail.gmail.com>
Content-Type: text/plain; charset=ISO-8859-1
Tienes un ejemplo en
http://www.datanalytics.com/blog/2012/04/19/variables-instrumentales-con-r/
Un saludo,
Carlos J. Gil Bellosta
http://www.datanalytics.com
El día 21 de mayo de 2012 16:54, tono ferri vidal
<tonoferri_vidal@hotmail.com> escribió:>
> Hola
>
> Alguien me puede ayudar? Necesito graficar 2 funciones de
> densidad "normales" superpuestas en un unico grafico.
>
> una densidad es de media 0 y desviacion 1,
> la otra es de media 3 y desviacion 1.
>
> saludos
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
------------------------------
Message: 6
Date: Mon, 21 May 2012 17:01:07 +0200
From: Javier Villacampa González
<javier.villacampa.gonzalez@gmail.com>
To: r-help-es@r-project.org
Subject: [R-es] help: Anova medidas repetidas
Message-ID:
<CAOZ+RxFWg6nKRvFhw6jz_uDd5yVeqcYvUBp735aexdMm8hPKBg@mail.gmail.com>
Content-Type: text/plain
Hola buenas,
lo primero una vez muchas gracias a todos por todas las veces que me habeis
ayudado.
Llevo días intentando utilizar las tablas ANOVA de medidas repetidas y he
encontrado multiple información, pero la verdad, no me aclaro.
He leido sobre la función ezANOVA o utilizando al aov con Error.
Los problema son:
¿Qué no me termino de aclarar con que significan las tablas. ?Vease
Anova1<- aov(RT ~ Prueba + Error(Subject) , datos) #+
Error(Subject/Prueba)> summary(Anova1)
Error: Subject
Df Sum Sq Mean Sq F value Pr(>F)
Prueba 2 40181770 20090885 18.893 4.114e-07 ***
Residuals 61 64866112 1063379
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 17206 356542956 20722
Que signifca este p-value exactamente, por más que leo solo me lio más
Además no entiendo cuando se emplea Error(Subject) y cuando
Error(Subject/Prueba). En mi caso los sujetos solo hacen una prueba
repeteidas veces así que creo que es Error(Subject)
¿Como se puede calcular las diferencias de Tukey para este tipo de objetos?
Tengo una solución que es programarlo, pero me gustaría que fuese un
paquete. Esté en esta página
http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/
Si teneis alguna lectura que recomendarme estaré gustoso de leerme. Soy
matemático así que las matemáticas no son un problema. Lo malo es
aplicarlas, en la carrera no se entra en este tipo de modelos.
Por último os dejo los datos por si no me he explicado bien.
*Los datos son:*
*RT: *variable continua de tiempos de reacción
*Prueba: *Tres tipos de pruebas neuropsicológicas
*Subject:* Participantes del experimento que hacen durante varios trial la
misma prueba.
head(datos) Subject Prueba card tri cat switch stim isi cod1 resp
Sorp RT Err SorpA SP Aciertos ratioA
1 s00 ngc plus 9 C plus 9 1523 90 0 2.000000
NA 1 0.000000 s00 ngc 354 0.924282
2 s00 ngc bb 10 C trg 1 1543 31 2 1.321928
1116.6 1 2.000000 s00 ngc 354 0.924282
3 s00 ngc rb 11 C trg 2 1534 32 1 1.584963
389.3 1 1.321928 s00 ngc 354 0.924282
4 s00 ngc rb 12 C trg 3 1552 33 1 1.222392
340.4 1 1.584963 s00 ngc 354 0.924282
5 s00 ngc rb 13 C trg 4 1541 34 1 1.000000
279.9 1 1.222392 s00 ngc 354 0.924282
6 s00 ngc bs 14 C trg 5 1545 35 2 2.169925
350.7 1 1.000000 s00 ngc 354 0.924282
tail(datos)
Subject Prueba card tri cat switch stim isi cod1 resp Sorp
RT Err SorpA SP Aciertos ratioA
24507 s74 Swt bs 585 Fo trg 3 1509 13 2 3.902135
374.1 1 3.903184 s74 Swt 330 0.8616188
24508 s74 Swt bs 586 Fo trg 4 1519 14 2 3.901116
339.1 1 3.902135 s74 Swt 330 0.8616188
24509 s74 Swt bb 587 Fo trg 5 1556 15 1 2.888630
870.1 1 3.901116 s74 Swt 330 0.8616188
24510 s74 Swt cros 588 C cros 8 1531 80 0 9.004247
NA 1 2.888630 s74 Swt 330 0.8616188
24511 s74 Swt rs 589 C trg 1 1543 11 1 3.136630
1013.9 1 9.004247 s74 Swt 330 0.8616188
24512 s74 Swt rb 590 C trg 2 1557 12 1 4.134127
983.4 1 3.136630 s74 Swt 330 0.8616188
Muchas gracias por adelantado.Un saludo.
Javier
[[alternative HTML version deleted]]
------------------------------
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
Fin de Resumen de R-help-es, Vol 39, Envío 37
*********************************************
[[alternative HTML version deleted]]