So ändern Sie panel-Etiketten und x-Achse sublabels zu einem Gitter bwplot

Bin ich 2 Wochen in RStudio (MacOS) also bitte verzeiht mir, wenn ich ' m mit Blick auf eine offensichtliche Funktion, dass würde mein problem lösen.

Als ein Projekt, ich bin versucht zu reproduzieren, ein box-plot-Diagramm mit 4 plots repräsentieren Net Nutzen, da Krankheit Typ -- "nicht schwere"(0) oder "schwere"(1) -- als x-Achse label, und angesichts der Behandlung -- "Gesprächstherapie"(0) oder "drug therapy"(1) -- als x-Achse sub-label.

Hier ist mein Skript so weit:

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)

Wenn Sie den code ausführen, sehen Sie die box-und-whisker-plot: ich bin fast da Dank der informativen Beiträge in diesem forum über lattice bwplot Funktion.

Allerdings bin ich noch weit davon entfernt, zufrieden mit dem Ergebnis. Ich benutzte die paste0 Funktion zum hinzufügen von string-Deskriptoren, um die X-Achse sublabels für Therapie-Gruppe (ursprünglich mit der Bezeichnung "1,2" und jetzt zeigen, wie "Tx-0, Tx-1"), aber ich würde im Idealfall wie die sublabels zu sagen "Reden-Therapie" und "Medikamentöse Therapie", beziehungsweise. (Ich wusste einfach nur nicht, wie zu tun, Weg mit den vorhandenen numerischen Beschriftungen.)

Ebenso möchte ich die panel-Etiketten um zu sagen, "Nicht Schwerwiegend", wo das label ist derzeit 0 und "Schweren", bei denen das Etikett ist derzeit 1.

InformationsquelleAutor Alex Mayer | 2013-02-05

Den einfachsten Weg soweit ich das sehen kann, ist, konvertieren Sie Ihre tx und dztype Faktoren mit den entsprechenden level-Namen.

tx <- factor(tx, levels=c(0,1), labels=c("Talk Therapy", "Drug Therapy"))
dztype <- factor(dztype, levels=c(0,1), labels=c("Not severe", "Severe"))
bwplot(NBwtp1000 ~ tx | dztype, xlab="Talk Therapy or Drug Therapy", 
       ylab="Net Benefit @ wtp 1000", horizontal=FALSE)

So ändern Sie panel-Etiketten und x-Achse sublabels zu einem Gitter bwplot

Andere Lösung ist die Verwendung von argument scale für die tx Ebenen und argument strip für dztype:

tx <- c(0,0,0,0,1,1,1,1)
dztype <- c(1,0,1,0,0,0,1,1)
bwplot(NBwtp1000 ~ tx | dztype, xlab="Talk Therapy or Drug Therapy", 
       ylab="Net Benefit @ wtp 1000", horizontal=FALSE,
       scales=list(x=list(labels=c("Talk therapy","Drug Therapy"))), 
       strip=strip.custom(var.name="Disease severity", 
                          factor.levels=c(" Not Severe", "Severe"),
                          strip.levels=rep(TRUE,2)))