Muy buenas, quiero calcular correlaciones de Pearson entre dos variables (a,b) teniendo en cuenta una tercera (c). Para ello estoy usando una función llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce una función similar en alguna librería de R? Por otro lado, para ver si me cuadraban los resultados, he hecho una regresión lineal entre las variables "a" y "c" y luego he hecho otra regresión entre los residuos del modelo resultante y la variable "b". Creo que esto sería más o menos lo mismo que hacer la correlación parcial entre a y b tendiendo en cuenta c, y de hecho los resultados son similares, pero en el caso de la regresión el valor de P no sale significativo. ¿Sería válido este método? Me interesa también representar gráficamente la correlación parcial, por eso he usado también este otro método. Cualquier ayuda es bienvenida y agradecida. Saludos, Javier
Hola Javier, John Fox tiene esta funcionalidad integrada en Rmcdr. Quizas https://stat.ethz.ch/pipermail/r-help/2009-March/193523.html pueda servirte para lo que necesitas. Seria interesante comparar los resultados que obtienes utilizando la funcion pcor.test() y con partial.cor() Saludos, Jorge Ivan Velez 2011/4/26 Javier Martinez <>> Muy buenas, > > quiero calcular correlaciones de Pearson entre dos variables (a,b) > teniendo en cuenta una tercera (c). Para ello estoy usando una función > llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en > realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce > una función similar en alguna librería de R? Por otro lado, para ver > si me cuadraban los resultados, he hecho una regresión lineal entre > las variables "a" y "c" y luego he hecho otra regresión entre los > residuos del modelo resultante y la variable "b". Creo que esto sería > más o menos lo mismo que hacer la correlación parcial entre a y b > tendiendo en cuenta c, y de hecho los resultados son similares, pero > en el caso de la regresión el valor de P no sale significativo. ¿Sería > válido este método? Me interesa también representar gráficamente la > correlación parcial, por eso he usado también este otro método. > > Cualquier ayuda es bienvenida y agradecida. > > Saludos, > > Javier > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >[[alternative HTML version deleted]]
miguel.angel.rodriguez.muinos en sergas.es
2011-Apr-26 20:30 UTC
[R-es] Correlaciones parciales
Hola Javier. Prueba con la función cor.test() del paquete "stats". Un Saludo, Miguel. ________________________________________ De: r-help-es-bounces en r-project.org [r-help-es-bounces en r-project.org] En nombre de Javier Martinez [javi.martinez.lopez en gmail.com] Enviado el: martes, 26 de abril de 2011 21:47 Para: r-help-es Asunto: [R-es] Correlaciones parciales Muy buenas, quiero calcular correlaciones de Pearson entre dos variables (a,b) teniendo en cuenta una tercera (c). Para ello estoy usando una función llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce una función similar en alguna librería de R? Por otro lado, para ver si me cuadraban los resultados, he hecho una regresión lineal entre las variables "a" y "c" y luego he hecho otra regresión entre los residuos del modelo resultante y la variable "b". Creo que esto sería más o menos lo mismo que hacer la correlación parcial entre a y b tendiendo en cuenta c, y de hecho los resultados son similares, pero en el caso de la regresión el valor de P no sale significativo. ¿Sería válido este método? Me interesa también representar gráficamente la correlación parcial, por eso he usado también este otro método. Cualquier ayuda es bienvenida y agradecida. Saludos, Javier _______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es 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
miguel.angel.rodriguez.muinos en sergas.es
2011-Apr-26 20:52 UTC
[R-es] Correlaciones parciales
Ups, olvídalo... no es esto lo que buscas!!! Sorry. ________________________________________ De: r-help-es-bounces en r-project.org [r-help-es-bounces en r-project.org] En nombre de Rodríguez Muíños, Miguel Ángel Enviado el: martes, 26 de abril de 2011 22:30 Para: javi.martinez.lopez en gmail.com; r-help-es en r-project.org Asunto: Re: [R-es] Correlaciones parciales Hola Javier. Prueba con la función cor.test() del paquete "stats". Un Saludo, Miguel. ________________________________________ De: r-help-es-bounces en r-project.org [r-help-es-bounces en r-project.org] En nombre de Javier Martinez [javi.martinez.lopez en gmail.com] Enviado el: martes, 26 de abril de 2011 21:47 Para: r-help-es Asunto: [R-es] Correlaciones parciales Muy buenas, quiero calcular correlaciones de Pearson entre dos variables (a,b) teniendo en cuenta una tercera (c). Para ello estoy usando una función llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce una función similar en alguna librería de R? Por otro lado, para ver si me cuadraban los resultados, he hecho una regresión lineal entre las variables "a" y "c" y luego he hecho otra regresión entre los residuos del modelo resultante y la variable "b". Creo que esto sería más o menos lo mismo que hacer la correlación parcial entre a y b tendiendo en cuenta c, y de hecho los resultados son similares, pero en el caso de la regresión el valor de P no sale significativo. ¿Sería válido este método? Me interesa también representar gráficamente la correlación parcial, por eso he usado también este otro método. Cualquier ayuda es bienvenida y agradecida. Saludos, Javier _______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es 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 _______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es 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
He probado partial.cor() y no sé por qué pero me da un error del tipo: "sistema es computacionalmente singular: número de condición recíproco = 0" No entiendo bien qué hace esa función. Gracias a ambos. Saludos, Javier 2011/4/26 Jorge Ivan Velez <jorgeivanvelez en gmail.com>:> Hola Javier, > John Fox tiene esta funcionalidad integrada en Rmcdr. Quizas > https://stat.ethz.ch/pipermail/r-help/2009-March/193523.html pueda servirte > para lo que necesitas. > Seria interesante comparar los resultados que obtienes utilizando la funcion > pcor.test() y con partial.cor() > Saludos, > Jorge Ivan Velez > > 2011/4/26 Javier Martinez <> >> >> Muy buenas, >> >> quiero calcular correlaciones de Pearson entre dos variables (a,b) >> teniendo en cuenta una tercera (c). Para ello estoy usando una función >> llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en >> realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce >> una función similar en alguna librería de R? Por otro lado, para ver >> si me cuadraban los resultados, he hecho una regresión lineal entre >> las variables "a" y "c" y luego he hecho otra regresión entre los >> residuos del modelo resultante y la variable "b". Creo que esto sería >> más o menos lo mismo que hacer la correlación parcial entre a y b >> tendiendo en cuenta c, y de hecho los resultados son similares, pero >> en el caso de la regresión el valor de P no sale significativo. ¿Sería >> válido este método? Me interesa también representar gráficamente la >> correlación parcial, por eso he usado también este otro método. >> >> Cualquier ayuda es bienvenida y agradecida. >> >> Saludos, >> >> Javier >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es en r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- http://oikos.inf.um.es
¿Podrías escribir las sentencia completa? "computacionalmente singular" hace pensar en que la matriz de varianza y covarianzas sea singular. Entre otras posibilidades de error en la sintaxis, lo pregunta más obvia es ¿alguna de las variables es combinación lineal de otras? Saludos. Javier Martinez escribió:> He probado partial.cor() y no sé por qué pero me da un error del tipo: > > "sistema es computacionalmente singular: número de condición recíproco = 0" > > No entiendo bien qué hace esa función. > > Gracias a ambos. > > Saludos, > > Javier > > 2011/4/26 Jorge Ivan Velez <jorgeivanvelez en gmail.com>: > >> Hola Javier, >> John Fox tiene esta funcionalidad integrada en Rmcdr. Quizas >> https://stat.ethz.ch/pipermail/r-help/2009-March/193523.html pueda servirte >> para lo que necesitas. >> Seria interesante comparar los resultados que obtienes utilizando la funcion >> pcor.test() y con partial.cor() >> Saludos, >> Jorge Ivan Velez >> >> 2011/4/26 Javier Martinez <> >> >>> Muy buenas, >>> >>> quiero calcular correlaciones de Pearson entre dos variables (a,b) >>> teniendo en cuenta una tercera (c). Para ello estoy usando una función >>> llamada "pcor.test" (http://www.yilab.gatech.edu/pcor.html), que en >>> realidad no está en ningún paquete de R, que yo sepa. ¿Alguien conoce >>> una función similar en alguna librería de R? Por otro lado, para ver >>> si me cuadraban los resultados, he hecho una regresión lineal entre >>> las variables "a" y "c" y luego he hecho otra regresión entre los >>> residuos del modelo resultante y la variable "b". Creo que esto sería >>> más o menos lo mismo que hacer la correlación parcial entre a y b >>> tendiendo en cuenta c, y de hecho los resultados son similares, pero >>> en el caso de la regresión el valor de P no sale significativo. ¿Sería >>> válido este método? Me interesa también representar gráficamente la >>> correlación parcial, por eso he usado también este otro método. >>> >>> Cualquier ayuda es bienvenida y agradecida. >>> >>> Saludos, >>> >>> Javier >>> >>> _______________________________________________ >>> R-help-es mailing list >>> R-help-es en r-project.org >>> https://stat.ethz.ch/mailman/listinfo/r-help-es >>> >> > > > >-- _____---^---_____ Univ. de Extremadura Dept. Matemáticas. Despacho B29 Tf: + 34 924 289 300 Ext. 86823