Appliquer fonction cor.test à un subset

[jeaneudes trihan]

Bonjourà tous,

J'essaie de tester une corrélation de pearson à un sous groupe de données. J'ai essayé d'utiliser le code suivant, mais R ne semble pas tenir compte du subset, car il m'affiche le même résultat de corrélation que pour cor.test(column_1, column_2)

Code : Tout sélectionner

> df <- data.frame(column_1=1:10,column_2=c(1:5,6,4,3,11,9),column_3=rep(c("group_A","group_B"),each=5))

> with(subset(df, column_3==c("group_A")), cor.test(column_1, column_2, alternative = ("two.sided"), method = c("pearson")))

J'ai également essayé cet autre code, mais sans succès

Code : Tout sélectionner

> df <- data.frame(column_1=1:10,column_2=c(1:5,6,4,3,11,9),column_3=rep(c("group_A","group_B"),each=5))

> cor.test(column_1, column_2, subset = column_3==c("Group_A))

Quelqu'un peut-il me dire ou je me trompe?

Merci beaucoup d'avance.
Amicalement
Jean-Eudes TRIHAN

[Eric Wajnberg]

Il y a plusieurs solutions possibles. En voici une :

Code : Tout sélectionner

lapply(split(df,df$column_3),function(x) cor.test(x$column_1,x$column_2)) 

HTH, Eric.

[jeaneudes trihan]

Merci beaucoup. C'est parfait.
Je profite juste de cette conversation (même si ce n'est pas vraiment le sujet initial), mais existe-t-il un moyen de représenter graphiquement cette formule par un plot qui donnerait par exemple les données selon le subset "Column_3" avec une couleur pour le "Group_A" et une autre couleur pour le "Group_B" et donc les droites de régression dans les couleurs correspondantes?

Amicalement, J-Eudes

[Eric Wajnberg]

Encore une fois, il existe plusieurs solutions. Avec l'argument "function()" de la fonction lapply() que j'ai utilisée, vous pouvez utiliser n'importe quelle fonction de votre choix, y compris une fonction que vous écrivez vous même et qui fait le graphique que vous souhaitez..

HTH, Eric.