Wyznaczanie etykiet za pomocą R

Question

Próbuję wykreślić podzestawy danych za pomocą wykresu aspektu. Jest to 2 kolumny x4 rząd działek. Pierwsza kolumna z wartościami od 120 do 150 i druga kolumna z wartościami od 180 do 250, a zmienne są przypisywane jako mniejsze lub większe w pliku danych. Mam problem z dodaniem konkretnych etykiet ('120 < = symbol grecki (alfa) < = 150', '180 < = symbol grecki (alfa) < = 250) na zmienne kolumny ("Lesser", "Greater") . Próbowałem to robić jak:

require(graphics) 
library(ggplot2) 
hp <- ggplot2(data) + .... 
xlow <- paste("120 <", expression(alpha), " < 150") 
xhi <- paste("180 <", expression(alpha), " < 250") 
mf_labeller <- function(var, value){ 
    value <- as.character(value) 
    if (var=="regime") { 
    value[value=="Lesser"] <- xlow 
    value[value=="Greater"] <- xhi 
} 
    return(value) 
} 
(hp %+% data) + facet_grid(param~regime, scales="free", labeller=mf_labeller) 

dostaję „Błąd wartości ==«Wielkopolska»: porównanie jest niedozwolone dla wyrażenia”. Próbowałem także wymuszać dane. Jeśli robię etykietowania jako w jakikolwiek sposób:

data2 <- transform(data, regime = factor(regime, levels=c("Lesser", "Greater"), labels=c(expression(paste("120 <", alpha," < 150")), expression("180 < alpha < 250")))) 

dostaję całe wyrażenie w brzmieniu dosłownym, a nie pożądanych działek. Może czegoś brakuje. Każda pomoc jest naprawdę doceniana !! góry dziękuję ..

Fragmenty mojego próbce wynosi:

regime,param,line,XX,Var,sner 
Lesser,Rise,VII,AA,4.968624,0.1275248 
Lesser,Rise,VII,BB,3.719405,0.08470305 
Lesser,Rise,VII,CC,7.608773,0.177848 
Lesser,Rise,VII,DD,9.874395,0.1367159 
Lesser,Text,VII,AA,4.968624,0.1275248 
Lesser,Text,VII,BB,3.719405,0.08470305 
Lesser,Text,VII,CC,7.608773,0.177848 
Lesser,Text,VII,DD,9.874395,0.1367159 
Lesser,Chant,VII,AA,0.1771826,0.186758 
Lesser,Chant,VII,BB,0.3611497,0.5484656 
Lesser,Chant,VII,CC,0.7719002,0.8864444 
Lesser,Chant,VIII,DD,1.829022,0.2639881 
Greater,Rise,VII,AA,4.968624,0.1275248 
Greater,Rise,VII,BB,3.719405,0.08470305 
Greater,Rise,VII,CC,7.608773,0.177848 
Greater,Rise,VII,DD,9.874395,0.1367159 
Greater,Text,VII,AA,4.968624,0.1275248 
Greater,Text,VII,BB,3.719405,0.08470305 
Greater,Text,VII,CC,7.608773,0.177848 
Greater,Text,VII,DD,9.874395,0.1367159 
Greater,Chant,VII,AA,0.1771826,0.186758 
Greater,Chant,VII,BB,0.3611497,0.5484656 
Greater,Chant,VII,CC,0.7719002,0.8864444 
Greater,Chant,VIII,DD,1.829022,0.2639881 

i kod użyłem: Pisałem wcześniej wyraz w niewłaściwy sposób tutaj w poście, ale ja używać odpowiedniego jednego w kodzie .

x <- read.table("sample.csv", header=T, sep=',') 
require(graphics) 
library(ggplot2) 
ppi <- 300 
png("figure.png", width=6*ppi, height=6*ppi, res=ppi) 
hp <- ggplot(data=x,aes(x=XX, y=Var, colour=line, group = line)) + geom_errorbar(aes(ymin=Var-sner, ymax=Var+sner, colour=line), width=.3) + geom_line(aes(ymax=Var+sner), size=0.7) + geom_point(aes(ymax=Var+sner), shape=21, size=2,fill="white") + theme_bw() + theme(axis.text.x = element_text(angle=90,vjust=0.25), panel.grid.minor=element_blank(), panel.grid.major=element_blank(), panel.background =element_blank(), legend.position="none") + scale_colour_hue(l=40) 
xlow <- expression(paste("120 <",alpha," < 150")) 
xhi <- expression(paste("180 <", alpha," < 250")) 
.. earlier code block as alternative here... 
data2 <- transform(x, regime = factor(regime, levels=c("Lesser", "Greater"), labels=c(bquote(120<.(alpha)~phantom()<150), bquote(180<.(alpha)~phantom()<250)) 
(hp %+% data2) + facet_grid(param~regime, scales="free", labeller=label_bquote) 

Próbowałem jak wcześniej z mf_labeller zbyt .. góry dziękuję

Answer 1

Nie mamy danych do testu (ani nie mamy kod, którego używasz), ale wcześniejsze doświadczenie z siatki spisek nauczył mnie, że może być potrzebna bquote lub jakaś kombinacja sapply i as.expression. Spróbuj:

..., labels=c(bquote(120<.(alpha)~phantom()<150), bquote(180<.(alpha)~phantom()<250)) 

W pierwszej sekcji kodu bym nie myślał należy używać expression wewnątrz paste ale raczej konstruowania etykiet, a potem buduje wektor ekspresyjny.

Jeśli poczujesz się zdesperowany, możesz najpierw zbudować fabułę, a potem wymusić, by etykiety wyglądały inaczej. How can I use grid to edit a ggplot2 object to add math expressions to facet labels?

Answer 2

Nie mam prawdziwej odpowiedzi, ale wyrażam to w nadziei, że jest to pomocne. Moje rozwiązanie powoduje błąd z dwoma mniejszymi niż znaki, ale działa z jednym, co uderza mnie jako zachowanie błędne. (Edit patrz komentarz poniżej lub Brian Diggs odpowiedzieć przypuszczalnie/mam nadzieję, że powyższe dlaczego to nie jest wadliwy.)

przykładowe dane, które działa, trzeba pamiętać, jestem zastępując plus dla jednego mniej niż w twolabs:

twolabs <- c("120<alpha+150","150<alpha+250") 
df<- data.frame(x=rnorm(mean=4,40), 
    y=rnorm(mean=4,40), 
    labs=sample(twolabs,40,replace=T)) 

Ta fabuła działa pięknie.

ggplot(df,aes(x=x,y=y))+geom_point()+facet_grid(.~labs,labeller=label_parsed) 

Jednak dane z dwoma mniej niż znaków zgłasza błąd:

twolabs <- c("120<alpha<150","150<alpha<250") 
df2<- data.frame(x=rnorm(mean=4,40), 
    y=rnorm(mean=4,40), 
    labs=sample(twolabs,40,replace=T)) 

> ggplot(df2,aes(x=x,y=y))+geom_point()+facet_grid(.~labs,labeller=label_parsed) 
Error in parse(text = x) : <text>:1:10: unexpected '<' 
1: 120<alpha< 

za pomocą label_parsed, myślę nazywamy ostatecznie nazwać gridDevices :: plotmath, na co warto. Czy ktoś może to zrobić dalej?

Answer 3

Formularz chcesz dla ekspresji plotmath jest

expression(120 <= alpha~phantom() <= 150) 
expression(180 <= alpha~phantom() <= 250) 

które można uzyskać w ten sposób lub pośrednio przez

parse(text="120 <= alpha~phantom() <= 150") 
parse(text="180 <= alpha~phantom() <= 250") 

Najłatwiej jest, aby ustawić poziom czynnika do tych ciągów i użyj funkcji label_parsed.

Mocking to się z danymi mtcars, ponieważ nie zapewniają powtarzalną przykład:

dat <- mtcars 
dat$regime <- factor(dat$am) 
levels(dat$regime) <- list("120 <= alpha~phantom() <= 150" = "0", 
          "180 <= alpha~phantom() <= 250" = "1") 

ggplot(aes(x=wt, y=mpg), data=dat) + 
    geom_point() + 
    facet_grid(cyl ~ regime, labeller=label_parsed) 

Jeśli naprawdę chcesz użyć własnej drukarki etykiet, a nie chcą zmień poziomy czynników twoich podstawowych danych, możesz to zrobić. Start z mock danych z „mniejszego” i poziomy „większy”:

dat <- mtcars 
dat$regime <- factor(dat$am, levels=c(0,1), labels=c("Lesser", "Greater")) 

Następnie funkcja etykietująca i kod ggplot byłoby:

mf_labeller <- function(var, value){ 
    value <- as.character(value) 
    if (var=="regime") { 
    value[value=="Lesser"] <- "120 <= alpha~phantom() <= 150" 
    value[value=="Greater"] <- "180 <= alpha~phantom() <= 250" 
    value <- lapply(value, function(x) parse(text=x)) 
    } 
    return(value) 
} 

ggplot(aes(x=wt, y=mpg), data=dat) + 
    geom_point() + 
    facet_grid(cyl ~ regime, labeller=mf_labeller) 

wykres wygląda tak samo.