4 min read

図のラベルを楽に付けたい (R)

## Warning in library(package, lib.loc = lib.loc, character.only = TRUE,
## logical.return = TRUE, : there is no package called 'agricolae'
## Warning in library(package, lib.loc = lib.loc, character.only = TRUE,
## logical.return = TRUE, : there is no package called 'slackr'

expression関数内に変数を動的に組み込みたかった話

研究生活では、図を書くことが多い
例えば以下のような図

plot1 <-
  data_frame(x = seq(from = 0, to = 100, length.out = 500), y = log(x) + rnorm(500, 0, 0.1)) %>%
  ggplot(aes(x = x, y = y)) +
  theme_bw(base_size = 20, base_family = "serif") +
  geom_point(alpha = .5)

plot1 +
  xlab(expression(paste("Light intensity ["  * mu * "mol" ~ m^-2 ~ s^-1 * "]"))) +
  ylab(expression(paste("Net photosynthetic rate [" * mu * "mol" ~ m^-2 ~ s^-1 * "]")))

光の強さに応じて、葉の光合成速度が大きくなるなど
似たような単位で似たような図をたくさん書く
単位によく現れるギリシャ文字、上・下付き、イタリックといった書式を
サポートするのは、expression関数だ
毎度expression(paste("Light intensity [" * mu * "mol" ~ m^-2 ~ s^-1 * "]"))なんて書くのは大変だし、expression(paste("Net photosynthetic rate [" * mu * "mol" ~ m^-2 ~ s^-1 * "]"))に変えるときにコピペするのもキレイではない
なので、軸のテキスト部分を変数化して楽をしようとした

expression関数の中では文字列を評価してくれない

labeler1 <- function(text){
  expression(paste(text ~ "["  * mu * "mol" ~ m^-2 ~ s^-1 * "]"))
}
plot1 + xlab(labeler1("value2"))

pasteを先に実行すると、muがないことを怒られる

paste(text, "[", mu, "g", m^-2, "]") %>% expression
## Error in paste(text, "[", mu, "g", m^-2, "]"): object 'mu' not found

遅延評価 (lazy evaluation)、変数の作る環境あたりが原因だということはわかった

どうすればいいのか

base::bquote()でいけた

labeler2 <- function(text){
  bquote(.(text) ~ "["  * mu * "mol" ~ m^-2 ~ s^-1 * "]")
}
plot1 + xlab(labeler2("Light intensity")) + ylab(labeler2("Net photosynthetic rate"))   

bquote()内では、.(変数)とすることで、変数を受け取ることができる
それ以外はexpressionと同じ使い方
文字間を詰める場合は*で、スペースを開ける場合には~で結合

plot(0, 0, type = "n", xlim = c(0, 1), ylim = c(0, 4), xlab = "", ylab = "")

moji <- letters

par(family = "Times New Roman")
text(.5, 3, labels = bquote(.(moji[1]) ~ "[" * mu * "g" ~ m^-2 * "]"), cex = 2)
text(.5, 2, labels = bquote(.(moji[2]) ~ "[" ~ mu ~ "g" ~ m^-2 ~ "]"), cex = 2)
text(.5, 1, labels = bquote(.(moji[3]) ~ "[" ~~ mu ~~ "g" ~ m^-2 ~~ "]"), cex = 2)

bquoteの使い方の曖昧和訳

bquote {base}

expression関数に部分的に変数をいれる

概要

LISPのbackquote macroに似たやつ
.()で囲まれた部分以外を引数whereで指定された環境の下で評価する

使用法

bquote(expr, where = parent.frame())

引数

expr: 言語オブジェクト (name, call, expression) where: 言語オブジェクトを評価する環境

expression関数で同じことできないのか

要勉強
言語オブジェクトと環境のあたりをもう少し理解できればうまく行く気がする

参考

bquote (R Documentation)
言語オブジェクト (どんな鳥も@Hatena::Diary)

session_info()
##  setting  value                       
##  version  R version 3.5.0 (2018-04-23)
##  system   x86_64, darwin15.6.0        
##  ui       X11                         
##  language (EN)                        
##  collate  en_US.UTF-8                 
##  tz       Australia/Brisbane          
##  date     2018-05-16                  
## 
##  package      * version   date       source                          
##  assertthat     0.2.0     2017-04-11 CRAN (R 3.5.0)                  
##  backports      1.1.2     2017-12-13 CRAN (R 3.5.0)                  
##  base         * 3.5.0     2018-04-24 local                           
##  bindr          0.1.1     2018-03-13 CRAN (R 3.5.0)                  
##  bindrcpp       0.2.2     2018-03-29 CRAN (R 3.5.0)                  
##  bitops       * 1.0-6     2013-08-17 CRAN (R 3.5.0)                  
##  blogdown       0.6       2018-04-18 CRAN (R 3.5.0)                  
##  bookdown       0.7       2018-02-18 CRAN (R 3.5.0)                  
##  codetools      0.2-15    2016-10-05 CRAN (R 3.5.0)                  
##  colorspace     1.3-2     2016-12-14 CRAN (R 3.5.0)                  
##  compiler       3.5.0     2018-04-24 local                           
##  data.table   * 1.10.4-3  2017-10-27 CRAN (R 3.5.0)                  
##  datasets     * 3.5.0     2018-04-24 local                           
##  devtools     * 1.13.5    2018-02-18 CRAN (R 3.5.0)                  
##  digest         0.6.15    2018-01-28 CRAN (R 3.5.0)                  
##  doParallel     1.0.11    2017-09-28 CRAN (R 3.5.0)                  
##  doRNG          1.6.6     2017-04-10 CRAN (R 3.5.0)                  
##  dplyr        * 0.7.4     2017-09-28 CRAN (R 3.5.0)                  
##  evaluate       0.10.1    2017-06-24 CRAN (R 3.5.0)                  
##  foreach      * 1.4.4     2017-12-12 CRAN (R 3.5.0)                  
##  ggplot2      * 2.2.1     2016-12-30 CRAN (R 3.5.0)                  
##  glue           1.2.0     2017-10-29 CRAN (R 3.5.0)                  
##  graphics     * 3.5.0     2018-04-24 local                           
##  grDevices    * 3.5.0     2018-04-24 local                           
##  grid         * 3.5.0     2018-04-24 local                           
##  gridExtra    * 2.3       2017-09-09 CRAN (R 3.5.0)                  
##  gtable       * 0.2.0     2016-02-26 CRAN (R 3.5.0)                  
##  htmltools      0.3.6     2017-04-28 CRAN (R 3.5.0)                  
##  iterators      1.0.9     2017-12-12 CRAN (R 3.5.0)                  
##  knitr        * 1.20      2018-02-20 CRAN (R 3.5.0)                  
##  labeling       0.3       2014-08-23 CRAN (R 3.5.0)                  
##  lazyeval       0.2.1     2017-10-29 CRAN (R 3.5.0)                  
##  lubridate    * 1.7.4     2018-04-11 CRAN (R 3.5.0)                  
##  magrittr     * 1.5       2014-11-22 CRAN (R 3.5.0)                  
##  MASS         * 7.3-49    2018-02-23 CRAN (R 3.5.0)                  
##  memoise        1.1.0     2017-04-21 CRAN (R 3.5.0)                  
##  methods      * 3.5.0     2018-04-24 local                           
##  munsell        0.4.3     2016-02-13 CRAN (R 3.5.0)                  
##  parallel       3.5.0     2018-04-24 local                           
##  pforeach     * 1.3       2018-04-25 Github (hoxo-m/pforeach@2c44f3b)
##  pillar         1.2.1     2018-02-27 CRAN (R 3.5.0)                  
##  pkgconfig      2.0.1     2017-03-21 CRAN (R 3.5.0)                  
##  pkgmaker       0.22      2014-05-14 CRAN (R 3.5.0)                  
##  plyr         * 1.8.4     2016-06-08 CRAN (R 3.5.0)                  
##  purrr          0.2.4     2017-10-18 CRAN (R 3.5.0)                  
##  R6             2.2.2     2017-06-17 CRAN (R 3.5.0)                  
##  RColorBrewer * 1.1-2     2014-12-07 CRAN (R 3.5.0)                  
##  Rcpp           0.12.16   2018-03-13 CRAN (R 3.5.0)                  
##  RCurl        * 1.95-4.10 2018-01-04 CRAN (R 3.5.0)                  
##  registry       0.5       2017-12-03 CRAN (R 3.5.0)                  
##  reshape2     * 1.4.3     2017-12-11 CRAN (R 3.5.0)                  
##  rlang          0.2.0     2018-02-20 CRAN (R 3.5.0)                  
##  rmarkdown      1.9       2018-03-01 CRAN (R 3.5.0)                  
##  rngtools       1.2.4     2014-03-06 CRAN (R 3.5.0)                  
##  rprojroot      1.3-2     2018-01-03 CRAN (R 3.5.0)                  
##  scales       * 0.5.0     2017-08-24 CRAN (R 3.5.0)                  
##  stats        * 3.5.0     2018-04-24 local                           
##  stringi        1.2.2     2018-05-02 cran (@1.2.2)                   
##  stringr      * 1.3.1     2018-05-10 cran (@1.3.1)                   
##  tibble         1.4.2     2018-01-22 CRAN (R 3.5.0)                  
##  tidyr        * 0.8.0     2018-01-29 CRAN (R 3.5.0)                  
##  tools          3.5.0     2018-04-24 local                           
##  utils        * 3.5.0     2018-04-24 local                           
##  withr          2.1.2     2018-03-15 CRAN (R 3.5.0)                  
##  xfun           0.1       2018-01-22 CRAN (R 3.5.0)                  
##  xtable         1.8-2     2016-02-05 CRAN (R 3.5.0)                  
##  yaml           2.1.18    2018-03-08 CRAN (R 3.5.0)