Me faltan 2 semanas para usar RStudio (MacOS), así que perdóneme si me olvido de una función obvia que resolvería mi problema
Como proyecto, estoy intentando reproducir un gráfico de diagrama de caja con 4 gráficos que representan el Beneficio neto, dado el tipo de enfermedad "no grave" (0) o "grave" (1), como una etiqueta del eje x, y tratamiento dado - "terapia de conversación" (0) o "terapia de medicamentos" (1) - como una subetiqueta del eje x.

Aquí está mi guión hasta ahora:

tx <- c(0,0,0,0,1,1,1,1)
dztype <- c(1,0,1,0,0,0,1,1)
NBwtp1000 <- c(-5500,-4000,-5000,-1000,-5000,-5000,-2800,-2000)
require(lattice)
bwplot(NBwtp1000 ~ paste0("Tx ", tx) | paste0("Disease Severity ", dztype),
       xlab="Talk Therapy (Tx 0) or Drug Therapy (Tx 1)", 
       ylab="Net Benefit @ wtp 1000", horizontal=FALSE)

Si ejecuta el código, verá el diagrama de caja y bigotes: Casi estoy allí gracias a algunas publicaciones informativas en este foro sobre celosías.bwplot función.

Sin embargo, todavía estoy lejos de estar contento con el resultado. Usé elpaste0 función para agregar descriptores de cadena a las subetiquetas del eje X para el grupo de tratamiento (originalmente etiquetado "1,2" y ahora se muestra como "Tx 0, Tx 1"), pero lo ideal sería que esos sublabeles dijeran "Terapia de conversación" y " Terapia de drogas ", respectivamente. (Simplemente no sabía cómo eliminar las etiquetas numéricas existentes).
De manera similar, me gustaría que las etiquetas del panel digan "No grave", donde la etiqueta es actualmente 0 y "Grave", donde la etiqueta es actualmente 1.