Es bastante fácil construir una matriz de gran dispersión agradable con histogramas en la diagonal para datos multivariados de la siguiente manera:

scatterplotMatrix(somedata[1:points.count,],groups=somedata[1:points.count,class],
                by.groups=TRUE,diagonal="histogram")

Sin embargo, según la documentación, no parece posible dividir el histograma por las etiquetas de grupo como se hace enesta pregunta. ¿Cómo harías eso usando scatterplotMatrix o una función similar?