r - Rで「ズームイン」効果を描画する方法

Question

上のスケッチ マップのように、上の図はパラメーター空間のプロット (たとえば、正規分布の平均と分散) であり、下の図は対応する密度プロットであると想像できます。これを行うためのヒントはありますか？ありがとう〜

更新: 拡張機能として、これのインタラクティブ バージョンを作成できますか? たとえば、ポイントの上にマウスを置くと、R はその下に対応するプロットを表示します。

Answer 1

これはインタラクティブ バージョンです。ポイントをクリックすると、対応する密度プロットが表示されます。主に使用さ?identifyれ、@Tylerが提案したように?zoomInPlot。

それがどのように機能するかについての詳細:最初rxlimにrylim定義されているのは、選択したポイントを囲む四角形のサイズであるため、 factor を変更したい場合があり/20ます。複数のクリックの可能性は自明ではidentify()ありません。「最近の」プロットでのみクリックを検出します。

par(mfrow = c(1,2))
plot(1:10) # 1
plot(1:10) # 2
identifyPch(1:10)

は、プロット #2 でのみクリックを検出します (こちらidentifyPch()は から?identify)。この問題では、par(mfg=c(1, 1))次のものが使用されました。

製造

c(i, j) 形式の数値ベクトル。ここで、i と j は、図形の配列内のどの図形が次に描画されるか (設定の場合)、または描画されているか (照会の場合) を示します。配列は、mfcol または mfrow によって既に設定されている必要があります。

zoom <- function (x, y, xlim, ylim, xd, yd) 
{
  rxlim <- x + c(-1, 1) * (diff(range(xd))/20)
  rylim <- y + c(-1, 1) * (diff(range(yd))/20)
  par(mfrow = c(1, 2))
  plot(xd, yd, xlab = "mean", ylab = "sd")
  xext <- yext <- rxext <- ryext <- 0
  if (par("xaxs") == "r") {
    xext <- diff(xlim) * 0.04
    rxext <- diff(rxlim) * 0.04
  }
  if (par("yaxs") == "r") {
    yext <- diff(ylim) * 0.04
    ryext <- diff(rylim) * 0.04
  }
  rect(rxlim[1] - rxext, rylim[1] - ryext, rxlim[2] + rxext, 
       rylim[2] + ryext)
  xylim <- par("usr")
  xypin <- par("pin")
  rxi0 <- xypin[1] * (xylim[2] - (rxlim[1] - rxext))/diff(xylim[1:2])
  rxi1 <- xypin[1] * (xylim[2] - (rxlim[2] + rxext))/diff(xylim[1:2])
  y01i <- xypin[2] * (xylim[4] - (rylim[2] + ryext))/diff(xylim[3:4])
  y02i <- xypin[2] * ((rylim[1] - ryext) - xylim[3])/diff(xylim[3:4])
  mu <- x
  curve(dnorm(x, mean = mu, sd = y), from = -4 * y + mu, to = 4 * y + mu, 
        xlab = paste("mean:", round(mu, 2), ", sd: ", round(y, 2)), ylab = "")
  xypin <- par("pin")
  par(xpd = NA)
  xylim <- par("usr")
  xymai <- par("mai")
  x0 <- xylim[1] - diff(xylim[1:2]) * (xymai[2] + xymai[4] + 
                                         rxi0)/xypin[1]
  x1 <- xylim[1] - diff(xylim[1:2]) * (xymai[2] + xymai[4] + 
                                         rxi1)/xypin[1]
  y01 <- xylim[4] - diff(xylim[3:4]) * y01i/xypin[2]
  y02 <- xylim[3] + diff(xylim[3:4]) * y02i/xypin[2]
  par(xpd = TRUE)
  xend <- xylim[1] - diff(xylim[1:2]) * xymai[2]/(2 * xypin[1])
  xprop0 <- (xylim[1] - xend)/(xylim[1] - x0)
  xprop1 <- (xylim[2] - xend)/(xylim[2] - x1)
  par(xpd = NA)
  segments(c(x0, x0, x1, x1), 
           c(y01, y02, y01, y02), 
           c(xend, xend, xend, xend), 
           c(xylim[4] - (xylim[4] - y01) * xprop0, 
             xylim[3] + (y02 - xylim[3]) * xprop0, 
             xylim[4] - (xylim[4] - y01) * xprop1, 
             xylim[3] + (y02 - xylim[3]) * xprop1))
  par(mfg = c(1, 1))
  plot(xd, yd, xlab = "mean", ylab = "sd")
}

ident <- function(x, y, ...)
{
  ans <- identify(x, y, n = 1, plot = FALSE, ...)
  if(length(ans)) {
    zoom(x[ans], y[ans], range(x), range(y), x, y)
    points(x[ans], y[ans], pch = 19)
    ident(x, y)
  }
}

x <- rnorm(10)
y <- rnorm(10, mean = 5)
par(mfrow = c(1, 2))
plot(x, y, xlab = "mean", ylab = "sd")
ident(x, y)

Answer 2

iplots パッケージはこれにかなり近いですが、技術的には「ズーム」はありません。iplots には、Java GUI で実装されたインタラクティブなリンク プロットがあります。1 つのプロットで点を選択すると、同じ点が他のプロットで強調表示されます。 パッケージのウェブサイトはこちら.

library(iplots)
data(Cars93)
iplot(Cars93$Horsepower, Cars93$MPG.city)
ihist(Cars93$Horsepower)

ここでは、左側の散布図でいくつかのポイントを選択して赤くし、右側のヒストグラムでも強調表示されています (ヒストグラムは密度プロットに最も近いものです)。