-1

我想在 R 中绘制一个伽马分布，其中比例参数（alpha）保持不变，但卷积参数发生变化（全部在 1 个图中）。我知道 dgamma 但并不真正知道 x 参数的含义。我想通过提供 alpha 和 beta 来绘制函数。这是可行的吗？

您的解决方案的输出是形状 <- c(2,5,6,12)

> plot(dgamma.wrapper, from=0, to=10)
Error in plot(dgamma.wrapper, from = 0, to = 10) : 
  object 'dgamma.wrapper' not found
> for (i in seq_along(shapes))
+   lines(dgamma(x,shape= shapes[i] , scale = 1) , from=0, to=10, col=i)
Error in dgamma(x, shape = shapes[i], scale = 1) : object 'x' not found
>
4

2 回答 2

1

你可以这样做：

dgamma.wrapper <- function(x)
   dgamma(x,shape= 2 , scale = 1) 
plot(dgamma.wrapper, from=0, to=10)

编辑

对于不止一种形状：

shapes <- c(2,5,6,12)
plot(dgamma.wrapper, from=0, to=10)
for (i in seq_along(shapes))
  curve(dgamma(x,shape= shapes[i] , scale = 1) , from=0, to=10, col=i,add=TRUE)
于 2014-03-10T13:56:30.017 回答
0

我想这就是你要找的。

x <- seq(0,20, .1)
plot(x, dgamma(x, scale=2, shape=1), type="l", ylim=c(0,.5), ylab="y")
for(shape in 2:8){
  lines(x, dgamma(x, scale=2, shape=shape), col=shape)
}

编辑 - 来自维基百科的引述：在概率论和统计学中，累积分布函数 (CDF)，或只是分布函数，描述了具有给定概率分布的实值随机变量 X 在小于或等于 x。在连续分布的情况下，它给出了从负无穷到 x 的概率密度函数下的面积。

这就是本例中 x 所代表的含义。在这里你展示了密度函数。我刚刚定义了一系列值，其中的情节看起来相当有趣。随意将其更改为更适合您的应用程序的内容。

于 2014-03-10T13:57:29.913 回答