Discussion :

[KaguraSohma]

Bonjour à tous,

J'aimerais pouvoir juxtaposer deux graphiques. Je m'explique: j'ai un graphique qui affiche les barycentres de plusieurs matrices et un autre graphique qui affiche le barycentre de ces barycentres. J'aimerais donc savoir s'il est possible d'avoir sur un même graphique les barycentres + le barycentre de ces barycentres (avec par exemple un point différent pour le différencier).

Si vous avez une idée, n'hésitez pas! Merci!

[pitipoisson]

Bonjour,

Ça dépend des fonctions que tu utilises les représenter.
Si tu souhaites une aide précise, je t'encourage à poster du code.

Si c'est juste pour ajouter des points, tu peux utiliser la fonction points().

[KaguraSohma]

Bon, je vais essayer d'expliquer sans que ce soit trop compliqué. Le sujet est dans un contexte biologique/médical. J'ai en fait plusieurs nappes de pression (ce sont en fait des images .jpeg qui montrent les pressions d'un certain dispositif. Le patient doit en fait s'asseoir sur ce dispositif, qui est une sorte de coussin, et le pc est relié à ce dispositif pour calculer les pressions exercées sur ce coussin).
J'ai donc par exemple 750 nappes de pression pour un patient. Je dois donc calculer le barycentre de chaque nappe de pression.
Voici le code de cette fonction:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
# Calcul du barycentre pour une nappe de pression
CmdCalculBarycentreNappeDePression <- function(NappePression) 
{
	#print("#### $x ####")
	LaDataMoyenneX <- 0
	LaDataMoyenneY <- 0
	LaSomme <- 0
	for (i in 1:32)
	{
		for (j in 1:32)
		{
			val <- NappePression[i,j]
			LaDataMoyenneX <- LaDataMoyenneX + val * i 
			LaDataMoyenneY <- LaDataMoyenneY + val * j
			LaSomme <- LaSomme + val
		}
	}
	LaBarycentre$x <<- LaDataMoyenneX/LaSomme
	LaBarycentre$y <<- LaDataMoyenneY/LaSomme
	#print(LaBarycentre$x)
	return(LaBarycentre)
	#print("#### $y ####")
	#print(LaBarycentre$y)
}

J'ai ensuite fait une fonction qui calcule le barycentre de toutes les nappes. La voici:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
CmdBarycentreNappe <- function()
{
	vec <- (CmdCalculBarycentreNappeDePression(AmgDataPressure$p[[1]]))
 
	for (i in 2:AmgNappeMax)
	{
		vec <- c(vec, CmdCalculBarycentreNappeDePression(AmgDataPressure$p[[i]]))
		mat <- matrix((vec), ncol=2, byrow=T)
	}
	print(mat)
	plot((mat), xlim =c(0,32), ylim=c(0,32),xlab="Coordonnée x des barycentres", ylab="Coordonnée y des barycentres", main= "Nuage de points des barycentres")	
}

Cette fonction permet aussi d'afficher le nuage de points de ces barycentres.

Pour finir, voici la fonction qui calcule le barycentre de tous les barycentres:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
#Calcul du barycentre moyen (moyenne des barycentres)
CmdBarycentreMoyen <- function()
{
	LaBarycentreMoyenX <- 0
	LaBarycentreMoyenY <- 0
	for (i in 1:AmgNappeMax)
	{
		LaBarycentre <- CmdCalculBarycentreNappeDePression(AmgDataPressure$p[[i]])
		LaBarycentreMoyenX <- LaBarycentreMoyenX + LaBarycentre$x
		LaBarycentreMoyenY <- LaBarycentreMoyenY + LaBarycentre$y
	}
	LaBarycentreMoyen$x <<- LaBarycentreMoyenX/AmgNappeMax
	LaBarycentreMoyen$y <<- LaBarycentreMoyenY/AmgNappeMax
#	print(LaBarycentreMoyen)
	return(LaBarycentreMoyen)
#	plot(LaBarycentreMoyen)
}

J'aimerais donc, à partir de ça, afficher dans un seul graphique, les barycentres de toutes les nappes de pression ainsi que le barycentre des barycentres.

Si quelqu'un a le courage de lire tout ça et d'essayer de comprendre, ce serait vraiment sympa!

Merci!

[pitipoisson]

Ah bin là je pense que tu peux simplement utiliser la fonction points()

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
BM <- CmdBarycentreMoyen()
 
points(BM$x, BM$y, <arguments optionnels>)

[KaguraSohma]

Super, ça marche!! C'était vraiment pas compliqué en fait...

Merci beaucoup!