É muito fácil criar uma matriz enorme e agradável de gráficos de dispersão com histogramas na diagonal para dados multivariados da seguinte maneira:

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

De acordo com a documentação, porém, não parece possível dividir o histograma pelos rótulos dos grupos, como é feito emessa questão. Como você faria isso usando scatterplotMatrix ou uma função semelhante?