r - 绘制伽马概率密度函数
问题描述
我试图在 R 中绘制 gamma 概率密度函数,其中 y∈(0,10) for (k = 1,μ = 1), (k = 2, μ = 1), (k = 2, μ = 2 )。在 R 中,
在 R 中,pgamma 函数接受:
pgamma(q, shape, rate = 1, scale = 1/rate, alpha = shape, beta = scale, lower.tail = TRUE, log.p = FALSE)
在 R 中,我尝试过:
pgamma(1,1,rate=1,scale = 1/rate, alpha = shape, beta = scale, lower.tail = True, log.p = False)
但我得到了消息
Error in pgamma(1, 1, rate = 1, scale = 1/rate, lower.tail = TRUE, log.p = FALSE) :
object 'rate' not found
这是我第一次绘制伽马分布,希望能得到一些帮助。
解决方案
下面使用基本 R 图形绘制三个密度。
首先,你想要的参数值。我假设您mu
的定义在Gamma 分布的维基百科页面中。
k <- c(1, 2, 2)
mu <- c(1, 1, 2)
theta <- mu/k
现在,情节。
plot(0, 0, xlim = c(0, 10), ylim = c(0, 1), type = "n")
for(i in seq_along(k))
curve(dgamma(x, shape = k[i], scale = theta[i]), from = 0, to = 10, col = i, add = TRUE)
推荐阅读
- mysql - 如何通过控制台 SQL 仅选择数字字段和仅文本字段?
- kubernetes - 我可以在 minikube 上禁用 Tls 验证吗?
- tensorflow - 如何部署 TensorFlow 概率模型
- android - 无法使用 post 方法从 API 获取所有数据
- powershell - 将字符串行放入“数组”
- linux - 为什么在其中创建新文件时目录的更改时间(ctime)会更改?
- c# - C# 在自定义对象中一起添加唯一项
- ag-grid - 在枢轴模式和服务器端渲染中显示图例的问题
- python - 将同一个 excel 中的多张工作表导入 pandas 的一个数据框中
- javascript - 页面刷新后如何保值?