r - x軸の特定の部分にラベルを付けて正規分布をプロットするにはどうすればよいですか？

Question

次のコードを使用して、Rで標準正規分布を作成しています。

x <- seq(-4, 4, length=200)
y <- dnorm(x, mean=0, sd=1)
plot(x, y, type="l", lwd=2)

x軸は、平均と、平均の上下の3つの標準偏差のポイントでラベル付けする必要があります。これらのラベルを追加するにはどうすればよいですか？

Answer 1

最も簡単な (しかし一般的ではない) 方法は、x 軸の制限を制限することです。+/- 1:3 シグマはそのようにラベル付けされ、平均は 0 としてラベル付けされます - 平均からの偏差が 0 であることを示します。

plot(x,y, type = "l", lwd = 2, xlim = c(-3.5,3.5))

もう 1 つのオプションは、より具体的なラベルを使用することです。

plot(x,y, type = "l", lwd = 2, axes = FALSE, xlab = "", ylab = "")
axis(1, at = -3:3, labels = c("-3s", "-2s", "-1s", "mean", "1s", "2s", "3s"))

Answer 2

この回答のコードを使用すると、作成をスキップして関数xで使用することができます。curve()dnorm

curve(dnorm, -3.5, 3.5, lwd=2, axes = FALSE, xlab = "", ylab = "")
axis(1, at = -3:3, labels = c("-3s", "-2s", "-1s", "mean", "1s", "2s", "3s"))

しかし、これは指定されたコードをもう使用しません。

Answer 3

R 組み込み関数を使用せずに何かを行う難しい方法が好きな場合、または R 外でこれを行いたい場合は、次の式を使用できます。

x<-seq(-4,4,length=200)
s = 1
mu = 0
y <- (1/(s * sqrt(2*pi))) * exp(-((x-mu)^2)/(2*s^2))
plot(x,y, type="l", lwd=2, col = "blue", xlim = c(-3.5,3.5))

Answer 4

一般的な例: Normal(2, 1)

f <- function(x) dnorm(x, 2, 1)
plot(f, -1, 5)

これは非常に一般的なもので、f は任意のパラメーターを使用して自由に定義できます。たとえば、次のようになります。

f <- function(x) dbeta(x, 0.1, 0.1)
plot(f, 0, 1)

Answer 5

この目標のために、私は特にラティスが大好きです。find P(a < X < b) などの確率の問題を処理するときに通常必要となる、曲線の下の特定の領域などのグラフィック情報を簡単に実装します。

library(lattice)

e4a <- seq(-4, 4, length = 10000)            # Data to set up out normal
e4b <- dnorm(e4a, 0, 1)

         xyplot(e4b ~ e4a,                   # Lattice xyplot
               type = "l",
               main = "Plot 2",
               panel = function(x,y, ...){
                   panel.xyplot(x,y, ...)
                   panel.abline( v = c(0, 1, 1.5), lty = 2)  #set z and lines

                   xx <- c(1, x[x>=1 & x<=1.5], 1.5)         #Color area
                   yy <- c(0,   y[x>=1 & x<=1.5], 0) 
                   panel.polygon(xx,yy, ..., col='red')
               })

この例では、 と の間の領域z = 1をz = 1.5目立たせます。問題に応じて、このパラメーターを簡単に移動できます。

軸ラベルは自動です。