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 best solutions below