Dispongo de un archivo de datos con una variable dependiente y, y dos independientes x, z. A continuación un enlace a los datos:
Con ayuda de la comunidad, he implementado el siguiente código para graficar mediante un spineplot la relación entre y ~ x.
#Transforming values to categorical
y <- as.character(y)
y[y>=2]<-"high"
y[y < 2 & y >= 1]<-"medium"
y[y<1]<-"low"
#Sort plot by reorganizing levels in the graph
y <- factor(y, levels = c("low","medium","high"))
#x as a factor
x <- as.factor(x)
#Spineplot
library("RColorBrewer")
print.perc <- function(k, tbl, ndigits=2, str.pct="%", min=0.1, cex=0.8) {
nx <- nrow(tbl)
off <- 0.02
xat <- c(0, cumsum(prop.table(margin.table(tbl, 1)) + off))
posx <- (xat[1L:nx] + xat[2L:(nx + 1L)] - off)/2
ptbl <- prop.table(tbl,1)
lbl <- ifelse(ptbl[k,] >= min,
paste(format(round(100*ptbl[k,], ndigits), nsmall=ndigits), str.pct, "\n",colnames(tbl), sep=""),
"")
text(posx[k], cumsum(ptbl[k,])-ptbl[k,]/2, lbl, cex=cex)
}
paleta <- brewer.pal(nlevels(y),"RdYlBu")
tbl <- spineplot(y ~ as.factor(x),
col=paleta,
xlab = "x",
ylab ="y",
yaxlabels = "")
for (k in 1:nrow(tbl)) print.perc(k, tbl, ndigits=2)
El resultado (utilizando los datos originales) es este:
Me gustaría ahora elaborar un gráfico del mismo tipo pero incluyendo también la variable z. Es decir, se trataría de un gráfico que representase la probabilidad de que se de y, dadas cada una de las combinaciones entre x y z. Para ello, creo que debería tratarse de un gráfico spineplot o heatplot pero en 3D.
Muchas gracias de antemano por los comentarios.