Barplot avec deux jeux de données

Version imprimable

15/04/2014, 17h23
Francois1435

Barplot avec deux jeux de données

Bonjour, j'ai un souci avec barplot: je fais un calcul pour connaitre le taux d'immigrés dans différentes professions à l'aide de table et cprop (fonction du package questionr qui renvoie simplement les pourcentages en colonnes d'une table). Je fais le calcul une fois sur tout le territoire français et une fois sur la seule région parisienne (à partir d'une base de données nommée d, sur laquelle j'effectue un subset pour isoler les franciliens). J'arrive à avoir un graphique par calcul, avec le code

Code:

barplot(prop[1:2,], ylim = c(0, 100), legend.text = TRUE)

(prop étant le nom du tableau de pourcentage sur le territoire français). J'aimerais avoir les deux graphiques en un seul, c'est à dire avoir pour chaque profession non plus une mais deux colonnes (une pour l'ensemble, une pour l'IDF). Je n'y arrive pas. Quelqu'un a une idée ?

Bonjour,

Vous avez un exemple dans l'aide de la fonction barplot:

Code:

example(barplot)

Vous pouvez y arriver en créant un jeu de données unique à partir de vos deux data.frames (avec rbind par exemple):

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
test.france <- data.frame(
 pourcentage=sample(x=seq(0, 100, 10), 4, replace=FALSE), 
 csp=c("toto", "tata", "titi", "tutu"), 
 region=rep("France", 4), 
 stringsAsFactors=FALSE)
 
test.nord <- data.frame(
 pourcentage=sample(x=seq(0, 100, 10), 4, replace=FALSE), 
 csp=c("toto", "tata", "titi", "tutu"), 
 region=rep("Nord", 4), 
 stringsAsFactors=FALSE)
 
test <- rbind(test.france, test.nord)
 
barplot(tapply(test$pourcentage, list(test$region, test$csp), c), beside=TRUE, legend.text=TRUE)

HTH

Vincent

16/04/2014, 12h36
Francois1435

Je ne comprends pas la totalité du code, donc je dois me tromper quelque part puisque ça ne fonctionne pas chez moi.

Pour commencer, je vais détailler un peu mon cas pour que ce soit clair. J'ai deux jeux de données, l'un, d, avec toute la population, l'autre, didf, avec seulement les parisiens, isolés grâce à subset et une variable région nommée REG.
Sur chaque jeu de donnée je fais un tri croisé entre immigration (oui/non) et profession (a,b,c). Par exemple,

Code:

tab <- table(d$immi,d$prof)

. Sur chacun, je fais ensuite un tableau de pourcentage ainsi:

Code:

cprop(tab)

, puis le barplot. ça fonctionne.
avec

Code:

cbind(tab,tab2)

, j'arrive à obtenir un barplot, mais le tableau comme le diagramme sont mal ordonnés: R se contente de coller les tableaux, quand je voudrais avoir côte à côte les pourcentages d'immigrés dans la profession a pour chaque population, puis la même chose pour la profession b, etc.

Maintenant, concernant le code proposé:
-je ne vois que deux variables, la profession et la région.
-je ne saisis pas le rôle de tapply, de list, et du "c" dans le code (oui ça fait beaucoup...). Je pense avoir bien compris ce que fait une fonction tapply (appliquer une fonction à plusieurs objets ou cellules), mais je ne vois pas ici...

Désolé pour mon ignorance, mais si des explications supplémentaires sont possibles, cela m'aiderait bien.

Merci en tout cas pour la première réponse

J'ai fini par réussir, ça paraitra sans doute anodin pour la plupart des gens, mais je poste quand même comment j'ai fait au cas où.

Concernant le barplot j'ai utilisé le code proposé, bien que j'ai toujours un peu de mal à comprendre ce que signifie la combinaison de tapply et list. Une fois qu'on a un dataframe qui correspond à la forme suggérée par l'exemple de vchouraki ça marche, mais je ne sais pas trop pourquoi.

Concernant le dit dataframe, ça a été plus compliqué. En fait, avec le calcul précédent, j'obtenais un tableau croisé entre les variables CSP et IMMI, avec chaque CSP deux fois (une pour la France, une pour Paris). Or ce qu'il fallait, c'était un tableau avec tous les pourcentages dans une colonne, une colonne CSP, une colonne région (et au départ une colonne IMMI, où l'on ne conserve ensuite que les lignes comprenant les immigrés). Pour constituer le tableau, j'ai donc fait:
Code:

1 2 3 tab.1 <- data.frame(data.frame(cprop(wtd.table(d$IMMI, d$CSP, weights=d$POIDS, normwt=TRUE))), region = rep("France", 15)) tab.2 <- data.frame(data.frame(cprop(wtd.table(didf$IMMI, didf$CSP, weights=didf$POIDS, normwt=TRUE))), region = rep("France", 15)) tab <- rbind(tab.1, tab.2)[tab$Var1 == "immigrés",]
Ce n'est sans doute pas la solution la plus simple, mais je n'ai pas trouvé mieux... s'il y a des idées je prends!

puis le barplot comme proposé:

Code:

barplot(tapply(tab$Freq, list(tab$region, tab$Var2), c), beside = TRUE, ylim = c(0,100), main = "Pourcentage d'immigrés par PCS selon la région", legend.text = TRUE)

voilà, donc je n'ai toujours pas vraiment compris exactement pourquoi, mais ça fonctionne. Merci pour votre aide!

edit: je n'arrive pas à éditer mon premier message pour passer le sujet en résolu, c'est normal?