ggplot2 - аннотировать текст с надстрочным индексом AND - PullRequest
0 голосов
/ 27 ноября 2018

Я создаю цикл с ggplot2 для большого набора данных.Я посмотрел на эту тему ggplot2 - аннотировать текст с верхним индексом , но объединение этого подхода с моей командой round((... не работает.

вот выдержка из моего фрейма данных full

   country.x    year     emissions      etsemit
   Austria      2005     16194772.5     16539659
   Austria      2006     15039192.4     15275065
   Austria      2007     13757090.8     14124646
   Austria      2008     13582006.8     14572511
   Austria      2009     12526267.6     12767555
   Austria      2010     13852187.5     15506112
   Austria      2011     13666544.9     15131551
   Austria      2012     12256272.5     13121434
   Austria      2013     11224625.0      8074514
   Austria      2014      9499543.9      6426135
   Austria      2015     10623549.8      7514263
   Austria      2016     10448925.8      7142937
   Austria      2017             NA      7795277
   Belgium      2005     29246990.2     25460856
   Belgium      2006     28136794.9     24099282
   Belgium      2007     27435552.7     23706084
   Belgium      2008     25344134.8     23166180
   Belgium      2009     25744709.0     21185552
   Belgium      2010     26341043.0     22073616
   Belgium      2011     22921875.0     18950876
   Belgium      2012     22809482.4     17463388
   Belgium      2013     21242431.6     16728267
   Belgium      2014     20375966.8     15230243
   Belgium      2015     21091058.6     16053800
   Belgium      2016     19792162.1     15027777
   Belgium      2017             NA     15093036

Вот мой код:

ctry <- unique(full$country.x)

for(i in (1:length(ctry))){

#i <- 1  
# Color settings: colorblind-friendly palette
cols <- c("#999999", "#E69F00", "#56B4E9", "#009E73", "#F0E442", "#0072B2", 
"#D55E00", "#CC79A7")

plot.df <- full[full$country.x==ctry[i],]

p <- ggplot() +
geom_line(data=plot.df,aes(x=plot.df$year, y=plot.df$emissions, color='UN 
1.A.1')) +
geom_line(data=plot.df,aes(x=plot.df$year, y=plot.df$etsemit, color='ETS 
20')) +
annotate(geom = 'text', label = round((summary(lm(emissions ~ etsemit, 
data=full))$r.squared) ,3), x = 
Inf, y = Inf, hjust = 1.5, vjust = 2) +
labs(x="Year",y="CO2 Emissions (metric tons)",z="",title=paste("Emissions 
Comparison for",ctry[i])) + 
xlim(2005,2017) +
theme(plot.margin=unit(c(.5,.5,.5,.5),"cm")) +
scale_color_manual(values = cols) +
scale_y_continuous(labels = scales::comma) +
scale_x_continuous(breaks = seq(2005, 2017, by = 5)) +
labs(color="Datasets")
p

ggsave(p,filename=paste("h:/",ctry[i],".png",sep=""),width=6.5, height=6)
}

Все работает гладко, но мне не удаетсячтобы добавить «R ^ 2 =» к функции аннотации:

annotate(geom = 'text', label = round((summary(lm(emissions ~ etsemit, 
  data=full))$r.squared) ,3), x = Inf, y = Inf, hjust = 1.5, vjust = 2)

label делает то, что мне нужно, кроме добавления R^2 = перед значением для R ^ 2.

Вот как это выглядит сейчас:

enter image description here

Я пробовал это:

annotate(geom = 'text', label = bquote("R^2 = "~.(round((summary(lm(emissions ~ etsemit, data=full))$r.squared) ,3))), x = Inf, y = Inf, hjust = 1.5, vjust = 2) , но это приводит к ошибке

Все, что мне нужно, это добавить R ^ 2 -> R^2 = round((... и на этом конкретном графике (это цикл) R^2 = 0.998

Заранее благодарю за любую помощь.

Nordsee

1 Ответ

0 голосов
/ 27 ноября 2018

Полагаю, annotate() на самом деле не принимает выражения типа базовых функций построения графиков, поэтому вместо передачи выражения в качестве выражения вам нужно передать его строковую версию и позволить ggplot проанализировать его для вас.Вы можете сделать это с помощью

annotate(geom = 'text', 
  label = paste0("R^2==", round(summary(lm(emissions ~ etsemit, data=full))$r.squared, 3)), 
  x = Inf, y = Inf, hjust = 1.5, vjust = 2, parse=TRUE)

Большую часть времени, если вы хотите выражение, вы бы использовали bquote, чтобы сделать одно

bquote(R^2==.( round(summary(lm(emissions ~ etsemit, data=full))$r.squared, 3)))

, и вы могли бы deparse(), чтобы получитьстроковая версия, но в этом случае проще использовать paste() Полагаю,

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...