bwplot coding and customise plot

35 views Asked by At

I need some help with a graphic. I want to add the boxplot of the complete data. I've done it on the graph but got two colours on the box plot.

I also need help figuring out how to make the names of the data mechanisms easier to read. For example, increase the width of the ten rectangles.

The code of the plot:

tp <- bwplot(r2~y2|typ, data=d,
         horizontal=TRUE, layout=c(1,10),
         xlab=expression(Y[2]),
         col=c(mdc(2),mdc(1),mdc(3)),strip=FALSE, xlim=c(0,10),
         par.strip.text = list(cex=0.5),
         par.settings = list(layout.heights=list(strip=5)),
         scales = list(y = list(cex = 0.8)),
         strip.left = strip.custom(bg="grey95"))

Here is the dataframe d:

d <- 
structure(list(y1 = c(4.50461095759073, 5.91136212424366, 6.2341726607176, 
4.31517001243106, 4.03468424501696, 4.50461095759073, 5.91136212424366, 
6.2341726607176, 4.31517001243106, 4.03468424501696, 4.50461095759073, 
5.91136212424366, 6.2341726607176, 4.31517001243106, 4.03468424501696, 
4.50461095759073, 5.91136212424366, 6.2341726607176, 4.31517001243106, 
4.03468424501696, 4.50461095759073, 5.91136212424366, 6.2341726607176, 
4.31517001243106, 4.03468424501696, 4.50461095759073, 5.91136212424366, 
6.2341726607176, 4.31517001243106, 4.03468424501696, 4.50461095759073, 
5.91136212424366, 6.2341726607176, 4.31517001243106, 4.03468424501696, 
4.50461095759073, 5.91136212424366, 6.2341726607176, 4.31517001243106, 
4.03468424501696, 4.50461095759073, 5.91136212424366, 6.2341726607176, 
4.31517001243106, 4.03468424501696, 4.50461095759073, 5.91136212424366, 
6.2341726607176, 4.31517001243106, 4.03468424501696, 4.50461095759073, 
5.91136212424366, 6.2341726607176, 4.31517001243106, 4.03468424501696
), y2 = c(4.01008089738578, 4.81674644303593, 5.93752409664717, 
4.09807878631992, 6.13756977661119, 4.01008089738578, 4.81674644303593, 
5.93752409664717, 4.09807878631992, 6.13756977661119, 4.01008089738578, 
4.81674644303593, 5.93752409664717, 4.09807878631992, 6.13756977661119, 
4.01008089738578, 4.81674644303593, 5.93752409664717, 4.09807878631992, 
6.13756977661119, 4.01008089738578, 4.81674644303593, 5.93752409664717, 
4.09807878631992, 6.13756977661119, 4.01008089738578, 4.81674644303593, 
5.93752409664717, 4.09807878631992, 6.13756977661119, 4.01008089738578, 
4.81674644303593, 5.93752409664717, 4.09807878631992, 6.13756977661119, 
4.01008089738578, 4.81674644303593, 5.93752409664717, 4.09807878631992, 
6.13756977661119, 4.01008089738578, 4.81674644303593, 5.93752409664717, 
4.09807878631992, 6.13756977661119, 4.01008089738578, 4.81674644303593, 
5.93752409664717, 4.09807878631992, 6.13756977661119, 4.01008089738578, 
4.81674644303593, 5.93752409664717, 4.09807878631992, 6.13756977661119
), r2 = structure(c(2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 
2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 
2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 
2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), levels = c("Ymis", 
"Yobs"), class = "factor"), typ = c("MARLEFT", "MARLEFT", "MARLEFT", 
"MARLEFT", "MARLEFT", "MARMID", "MARMID", "MARMID", "MARMID", 
"MARMID", "MARRIGHT", "MARRIGHT", "MARRIGHT", "MARRIGHT", "MARRIGHT", 
"MARTAIL", "MARTAIL", "MARTAIL", "MARTAIL", "MARTAIL", "MCAR", 
"MCAR", "MCAR", "MCAR", "MCAR", "MNARLEFT", "MNARLEFT", "MNARLEFT", 
"MNARLEFT", "MNARLEFT", "MNARMID", "MNARMID", "MNARMID", "MNARMID", 
"MNARMID", "MNARRIGHT", "MNARRIGHT", "MNARRIGHT", "MNARRIGHT", 
"MNARRIGHT", "MNARTAIL", "MNARTAIL", "MNARTAIL", "MNARTAIL", 
"MNARTAIL", "Full", "Full", "Full", "Full", "Full", "Full", "Full", 
"Full", "Full", "Full")), row.names = c(NA, -55L), class = "data.frame")
0

There are 0 answers