グラフィックス参考実例集：数式のプロット のバックアップ(No.1)

• バックアップ一覧
• 差分 を表示
• 現在との差分 を表示
• ソース を表示
• グラフィックス参考実例集：数式のプロット へ行く。
  • 0 (1970-01-01 (木) 18:00:00)
  • 1 (2015-03-01 (日) 01:15:59)

グラフィックス参考実例集：数式のプロット

タイトル等の注釈に expression, substitute 関数の引数として数式表現を与えることにより数式のプロットが可能になる。使える記号や品位は限られており、数式・記号の完全な制御には pictex ドライバーを用い、latex コードに落し、マニュアルで編集する。

例1†

plotmath1 <- function () {
  x <- seq(-4, 4, len = 101)
  y <- cbind(sin(x), cos(x))
  matplot(x, y, type = "l", xaxt = "n",
          main = expression(paste(plain(sin) * phi, "  and  ",
              plain(cos) * phi)),
          ylab = expression("sin" * phi, "cos" * phi), # only 1st is taken
          xlab = expression(paste("Phase Angle ", phi)),
          col.main = "blue")
  axis(1, at = c(-pi, -pi/2, 0, pi/2, pi),
       lab = expression(-pi, -pi/2, 0, pi/2, pi))
}

例2。数式と数値変数の結合†

plotmath2 <- function () {
  plot(1:10, type="n", xlab="", ylab="", main = "plot math & numbers")
  tt <- 1.23 ; mtext(substitute(hat(theta) == that, list(that= tt)))
  for(i in 2:9)
    text(i,i+1, substitute(list(xi,eta) == group("(",list(x,y),")"),
                           list(x=i, y=i+1)))
}

例3†

plotmath3 <- function () {
  plot(1:10, 1:10)
  text(4, 9, expression(hat(beta) == (X^t * X)^{-1} * X^t * y))
  text(4, 8.4, "expression(hat(beta) == (X^t * X)^{-1} * X^t * y)",
       cex = .8) # 二重引用符で囲めば単なる文字列扱いになる
  text(4, 7, expression(bar(x) == sum(frac(x[i], n), i==1, n)))
  text(4, 6.4, "expression(bar(x) == sum(frac(x[i], n), i==1, n))",
       cex = .8)
  text(8, 5, expression(paste(frac(1, sigma*sqrt(2*pi)), " ",
      plain(e)^{frac(-(x-mu)^2, 2*sigma^2)})),
       cex = 1.2)
}

plotmath 機能で扱える数式の一覧†