R中的伽玛分布

时间:2014-03-10 13:49:33

标签: r gamma

我想在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
> 

2 个答案:

答案 0 :(得分: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)

答案 1 :(得分: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在这个例子中所代表的含义。在这里显示密度函数。我刚刚定义了一系列值,其中的情节看起来相当有趣。随意将其更改为更适合您的应用程序的东西。