RG#50: beeswarm box violin plot

# data

set.seed(1234)

bimodal <- c(rnorm(250, -2, 0.6), rnorm(250, 2, 0.6))

uniform <- runif(500, -4, 4)

normal <- rnorm(500, 0, 1.5)

dataf <- data.frame (group = rep(c("bimodal","uniform", "normal"), each = 500), xv = c(bimodal, uniform, normal), cg = rep( c("A","B"), 750))

require(beeswarm)

beeswarm(xv ~ group, data = dataf,method = 'swarm',

pch = 16, pwcol = as.numeric(cg),

xlab = '', ylab = 'xv', labels = c('bimodal', 'uniform', 'normal'))

 

# center 

beeswarm(xv ~ group, data = dataf,

        method = 'center',

        pch = 16, pwcol = as.numeric(cg),

        xlab = '', ylab = 'xv',

        labels = c('bimodal', 'uniform', 'normal'))


# hexagon 

beeswarm(xv ~ group, data = dataf,

        method = 'hex',

        pch = 16, pwcol = as.numeric(cg),

        xlab = '', ylab = 'xv',

        labels = c('bimodal', 'uniform', 'normal'))

# square

dataf$col <- as.numeric(dataf$cg)

dataf$col[dataf$col ==1] = "green4"
dataf$col[dataf$col ==2] = "purple"
beeswarm(xv ~ group, data = dataf,

       method = 'square',

        pch = 16, pwcol = col,

        xlab = '', ylab = 'xv',


        labels = c('bimodal', 'uniform', 'normal'))

RG#18: Violin Box plot

#Violin box plot
#data
 set.seed (45)
xv <- rnorm(500, 50,15)
gv <- c(rep("A",100), rep("B",100), rep("C",100), rep("D",100), rep("E", 100))
dataf <- data.frame(xv, gv)


#plot
require(lattice)
library("latticeExtra")
bwplot(xv, panel = panel.violin, box.ratio = 3, xlab="y(units)")

 bwplot(xv~gv, data=dataf, panel = panel.violin, box.ratio = 3, xlab="y(units)", col = "green4")

# using vioplot 

require(vioplot)

vioplot(xv[1:100], xv[101:200], xv[201:300], xv[301:400], xv[401:500],col="cyan2", border="cyan2", lty=1, lwd=1, rectCol="red",

colMed="white", pchMed=19, add=FALSE, wex=1,drawRect=TRUE  )

vioplot(xv[1:100], xv[101:200], xv[201:300], xv[301:400], xv[401:500],col="white", border="red", lty=1, lwd=1, rectCol="red",

colMed="yellow", pchMed=19, add=TRUE, wex=1,drawRect=TRUE  )

 set.seed (45)
xv <- rnorm(500, 50,15)
pch1 <- rep(1:5, each = 100)

plot(1, type="n", axes=F, xlab="", ylab="")

vioplot(xv[1:100], xv[101:200], xv[201:300], xv[301:400], xv[401:500],col="white", border="red", lty=1, lwd=1, rectCol="red",

colMed="yellow", pchMed=19, add=FALSE, wex=1,drawRect=TRUE )

points(pch1, xv, pch = 19, col = "blue")