R语言之可视化（28）蜜蜂图

stripchart(decrease ~ treatment, data = OrchardSprays, 
  vertical = TRUE, log = "y", method = 'jitter', jitter = 0.2, cex = 1,
  pch = 16, col = rainbow(8),
  main = 'stripchart')
beeswarm(decrease ~ treatment, data = OrchardSprays, 
  log = TRUE, pch = 16, col = rainbow(8),
  main = 'beeswarm')

data(breast)
  beeswarm(time_survival ~ ER, data = breast,
    pch = 16, pwcol = 1 + as.numeric(event_survival),
    xlab = "", ylab = "Follow-up time (months)",
    labels = c("ER neg", "ER pos"))
  legend("topright", legend = c("Yes", "No"),
    title = "Censored", pch = 16, col = 1:2)

## Generate some random data
set.seed(123)
distro <- list(runif = runif(100, min = -3, max = 3), 
               rnorm = rnorm(100))
for (m in c("swarm", "center", "hex", "square")) {
  beeswarm(distro, 
    col = 2:3, pch = 16,
    method = m, 
    main = paste('method = "', m, '"', sep = ''))
}

for (ii in c("none", "gutter", "wrap", "random", "omit")) {
  beeswarm(distributions, 
    pch = 21, col = 2:4, bg = "#00000050",
    corral = ii, 
    main = paste('corral = "', ii, '"', sep = ''))
}

 boxplot(len ~ dose, data = ToothGrowth, 
    outline = FALSE,     ## avoid double-plotting outliers, if any
    main = 'boxplot + beeswarm')
  beeswarm(len ~ dose, data = ToothGrowth, 
    col = 4, pch = 16, add = TRUE)

  beeswarm(len ~ dose, data = ToothGrowth, 
    col = 4, pch = 16,
    main = 'beeswarm + bxplot')
  bxplot(len ~ dose, data = ToothGrowth, add = TRUE)

 beeswarm(distributions, col = 2:4, 
    main = 'Default')
  beeswarm(distributions, col = 2:4, side = -1, 
    main = 'side = -1')
  beeswarm(distributions, col = 2:4, side = 1, 
    main = 'side = 1')
  beeswarm(distributions, col = 2:4, priority = "descending", 
    main = 'priority = "descending"')
  beeswarm(distributions, col = 2:4, priority = "random", 
    main = 'priority = "random"')  
  beeswarm(distributions, col = 2:4, priority = "density", 
    main = 'priority = "density"')  

if(!require(beeswarm)) install.packages(beeswarm)
data(breast)
beeswarm(time_survival ~ event_survival, data = breast,
  method = 'swarm',
  pch = 16, pwcol = as.numeric(ER),
  xlab = '', ylab = 'Follow-up time (months)',
  labels = c('Censored', 'Metastasis'))
legend('topright', legend = levels(breast$ER),
  title = 'ER', pch = 16, col = 1:2)</pre>

beeswarm(time_survival ~ event_survival, data = breast,
method = 'swarm',
pch = 16, pwcol = as.numeric(ER),
xlab = '', ylab = 'Follow-up time (months)',
labels = c('Censored', 'Metastasis'))

boxplot(time_survival ~ event_survival, data = breast, add = T, names = c("",""), col="#0000ff22")