Graphique à bulles catégoriques pour les études de cartographie
comment créer un graphique à bulles catégorique, en utilisant GNU R, similaire à celui utilisé dans les études de cartographie systématique (voir ci-dessous)?
EDIT: ok, voici ce que j'ai essayé jusqu'à présent. Premièrement, mon ensemble de données (Var1 va à l'axe x, Var2 va à l'axe y):
> grid
Var1 Var2 count
1 Does.Not.apply Does.Not.apply 53
2 Not.specified Does.Not.apply 15
3 Active.Learning..general. Does.Not.apply 1
4 Problem.based.Learning Does.Not.apply 2
5 Project.Method Does.Not.apply 4
6 Case.based.Learning Does.Not.apply 22
7 Peer.Learning Does.Not.apply 6
10 Other Does.Not.apply 1
11 Does.Not.apply Not.specified 15
12 Not.specified Not.specified 15
21 Does.Not.apply Active.Learning..general. 1
23 Active.Learning..general. Active.Learning..general. 1
31 Does.Not.apply Problem.based.Learning 2
34 Problem.based.Learning Problem.based.Learning 2
41 Does.Not.apply Project.Method 4
45 Project.Method Project.Method 4
51 Does.Not.apply Case.based.Learning 22
56 Case.based.Learning Case.based.Learning 22
61 Does.Not.apply Peer.Learning 6
67 Peer.Learning Peer.Learning 6
91 Does.Not.apply Other 1
100 Other Other 1
Ensuite, essayer de tracer les données:
# Based on http://flowingdata.com/2010/11/23/how-to-make-bubble-charts/
grid <- subset(grid, count > 0)
radius <- sqrt( grid$count / pi )
symbols(grid$Var1, grid$Var2, radius, inches=0.30, xlab="Research type", ylab="Research area")
text(grid$Var1, grid$Var2, grid$count, cex=0.5)
Voici le résultat: 
problèmes: les étiquettes des axes sont erronées, les lignes de grille pointillées sont manquant.
3 réponses
Voici la solution de ggplot2. Tout d'abord, radius ajouté comme nouvelle variable à votre base de données.
grid$radius <- sqrt( grid$count / pi )
vous devriez jouer avec la taille des points et les étiquettes de texte à l'intérieur de la parcelle pour l'ajustement parfait.
library(ggplot2)
ggplot(grid,aes(Var1,Var2))+
geom_point(aes(size=radius*7.5),shape=21,fill="white")+
geom_text(aes(label=count),size=4)+
scale_size_identity()+
theme(panel.grid.major=element_line(linetype=2,color="black"),
axis.text.x=element_text(angle=90,hjust=1,vjust=0))
Voici une version utilisant levelplotlatticeExtra.
library(latticeExtra)
levelplot(count~Var1*Var2,data=dat,
panel=function(x,y,z,...)
{
panel.abline(h=x,v=y,lty=2)
cex <- scale(z)*3
panel.levelplot.points(x,y,z,...,cex=5)
panel.text(x,y,label=z,cex=0.8)
},scales=(x=list(abbreviate=TRUE))) ## to get short labels
Pour obtenir la taille de la bulle proprtional le comte , vous pouvez le faire
library(latticeExtra)
levelplot(count~Var1*Var2,data=dat,
panel=function(x,y,z,...)
{
panel.abline(h=x,v=y,lty=2)
cex <- scale(z)*3
panel.levelplot.points(x,y,z,...,cex=5)
panel.text(x,y,label=z,cex=0.8)
})
Je ne l'affiche pas car le rendu n'est pas clair comme dans le cas de la taille du correctif.
Cela vous permettra de commencer en ajoutant les marques de graduation à votre xaxis.
Pour ajouter les lignes, il suffit d'ajouter une ligne à chaque niveau
ggs <- subset(gg, count > 0)
radius <- sqrt( ggs$count / pi )
# ggs$Var1 <- as.character(ggs$Var1)
# set up your tick marks
# (this can all be put into a single line in `axis`, but it's placed separate here to be more readable)
#--------------
# at which values to place the x tick marks
x_at <- seq_along(levels(gg$Var1))
# the string to place at each tick mark
x_labels <- levels(gg$Var1)
# use xaxt="n" to supress the standard axis ticks
symbols(ggs$Var1, ggs$Var2, radius, inches=0.30, xlab="Research type", ylab="Research area", xaxt="n")
axis(side=1, at=x_at, labels=x_labels)
text(ggs$Var1, ggs$Var2, ggs$count, cex=0.5)
notez également qu'au lieu d'appeler l'objet grid je l'ai appelé gg et ggs pour le sous-ensemble. grid est une fonction R. Bien qu'il soit "autorisé" d'écraser la fonction avec un objet, il n'est pas recommandé et peut conduire à des bogues gênants dans la ligne.