0 环境与数据准备

#---环境建设---#
rm(list = ls());#清空变量空间
library(datasets)
library(export)

#---数据读入---#
dataname <- "iris"
rawdata <- datasets::iris

1、GGally_ggpairs

#散点图矩阵——Iris(GGally_ggpairs)
library(GGally)
ggpairs(rawdata, columns=1:5, aes(color=Species)) + 
  ggtitle("散点图矩阵——Iris(GGally_ggpairs)")+
  theme_bw() 
graph2png(file = paste("散点图矩阵(GGally_ggpairs) of", dataname))

生成结果:

2、pairs

#散点图矩阵——Iris(pairs)
panel.hist <- function(x, ...)
{
  usr <- par("usr"); on.exit(par(usr))
  par(usr = c(usr[1:2], 0, 1.5) )
  h <- hist(x, plot = FALSE)
  breaks <- h$breaks; nB <- length(breaks)
  y <- h$counts; y <- y/max(y)
  rect(breaks[-nB], 0, breaks[-1], y, col = "cyan", ...)
}
pairs(rawdata[1:5], main = "散点图矩阵——Iris(pairs)", pch = 21,
      panel = panel.smooth,
      diag.panel = panel.hist,
      bg = c("#1b9e77", "#d95f02", "#7570b3")[unclass(rawdata$Species)])
graph2png(file = paste("散点图矩阵(pairs) of", dataname))

生成结果: 

3、lattice_splom

#散点图矩阵——Iris(lattice_splom)
library(lattice)
splom(rawdata[1:5], 
      groups=rawdata$Species, 
      main="散点图矩阵——Iris(lattice_splom)")
graph2png(file = paste("散点图矩阵(lattice_splom) of", dataname))

生成结果:

4、WVPlots_PairPlot

#散点图矩阵——Iris(WVPlots_PairPlot)
library(WVPlots) 
PairPlot(rawdata, 
         colnames(rawdata)[1:5], 
         "散点图矩阵——Iris(WVPlots_PairPlot)", 
         group_var = "Species")
graph2png(file = paste("散点图矩阵(WVPlots_PairPlot) of", dataname))

生成结果:

5、YaleToolkit_gpairs

#散点图矩阵——Iris(YaleToolkit_gpairs)
library(YaleToolkit)
library(gpairs)
gpairs(iris, upper.pars = list(scatter = 'stats'),
       scatter.pars = list(pch = 1:3,
                           col = as.numeric(iris$Species)),
       stat.pars = list(verbose = FALSE))
graph2png(file = paste("散点图矩阵(YaleToolkit_gpairs) of", dataname))

致谢:

学习过程中,参考了网络上各位大侠的美文,在此一并表示感谢!

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐