问题标签 [beta-distribution]

我有一个数据集，其中包含一个分布在 Beta 分布上的相关右删失变量。请不要提供flexsurv软件包。有必要描述我在 AFT 回归生存模型survreg中使用的 Beta 分布的函数。M. Zhu 一书中的一个例子，我替换了我知道的函数，但我无法描述 Beta 分布的密度函数的一阶和二阶导数的比率：

我相信这是由于基础包rbeta中的函数、gamma函数digamma和函数的导数，但我的数学水平不允许应用它们。trigamma

我需要在r. 我正在寻找可以生成双变量 beta 分布的包/代码。我找不到这个分布的 r 函数。

有谁知道r这个模拟可以使用哪个功能？

我正在上一门机器学习课程，我们得到了第一个统计数据——“编程”练习。

所以练习是这样的：

回想一下“亚马逊的两个卖家价格相同”讲座中的故事。一个有 90 条正面评价和 10 条负面评价。另一个2阳性0阴性。你应该向谁购买？” 写下关于可靠性的后验概率（如讲座中所述）。使用数值积分计算 p(x > y|D1, D2)。您可以使用函数scipy.stats.beta.rvs(a,b,size)生成 Beta 分布的样本。

我们从讲座中了解到以下内容：

应用了两个 Beta 二项式模型：p(x|D1) = Beta(x|91, 11) 和 p(y|D2) = Beta(y|3, 1)

计算卖方 1 比卖方 2 更可靠的概率：

所以我在 Python 中的尝试是这样的：

我收到一条错误消息：

... The maximum number of subdivisions (50) has been achieved....

但最终有一个计算的答案是大约。1.7 . （我们被告知答案是大约 0.7 ）

我的问题是：如何计算 [x > y] 部分，意思是：卖家 1 (x) 比卖家 2 (y) 更可靠的概率？

可以通过单击一个按钮将 Alpha/Beta 版本升级到生产，但由于某种原因，“内部测试”轨道没有“发布到生产”按钮，只有“发布到 alpha/beta”。

我找不到有关该主题的任何文档，并且想知道它是否是故意阻止的，或者我只是遗漏了一些东西。

目前的解决方法是先将发布推广到 Alpha，然后再从 Alpha 推广到生产。但这没有多大意义，因为这些 alpha/beta 阶段目前不用于任何类型的预发布测试。

或者，我每次都可以为生产创建一个单独的版本并只选择相同的 APK，但是当“内部测试”版本中已经指定了一些发行说明和其他元数据并且所有这些信息都应该手动复制时，这也没有意义.

我正在尝试计算可以在此处找到的 A/B 测试数据集的贝叶斯因子。但是，我最终得到一个 NaN，因为 beta 系数的计算结果为零。在计算可能性时，我假设它遵循二项分布。因此，我遵循这个公式：

可能性 = 选择(n,k) * Beta(k+1,n-k+1)

代码可以在下面找到

对于随机模拟，我需要绘制很多随机数，这些随机数是 beta 二项式分布的。

目前我以这种方式实现它（使用python）：

因此对随机数 x 进行采样可以通过以下方式完成：

问题是，对一个随机数进行采样需要 ca。0.5ms，在我的例子中，这主导了整个模拟速度。限制因素是对 beta 函数（或其中的 gamma 函数）的评估。

有谁知道如何加快采样速度？

我正在使用 MATLAB 将 beta 分布拟合到我的数据中。

为此，有几个选项，例如fitdist，它们提供了一组aand并且可以使用或简单地使用and命令b访问统计信息（均值和标准差） 。betastatmeanstd

但是，如果我希望我的分布具有给定的平均值，并且只有 std 是未知的，我该怎么办。

有一个例子。我有一个可变大小的半直方图，其中定义了中心比率及其概率，如下所示：

现在我想为数据拟合一个“beta 分布”，范围在 [0,1] 之间，均值等于vul. 我现在需要的是最好的或std配对。ab

如何对随机效应 beta 分布进行事后检验？

为了构建具有 beta 分布的混合模型，我使用了库 gamlss（我发现在 R 中没有其他方法可以做到这一点）。在这个模型中，有一个具有三个水平的因子。假设 level1、level2 和 level3。摘要显示了 level1 x level2 和 level1 x level3 之间的比较，因为 level1 是基线。

library(gamlss) mymodel <- gamlss(response~ph+dose, family = BE, random=~1|id, data = mydata) 摘要(mymodel)

它显示了 level1 x level2 和 level1 x level3 之间的比较

对于所有三个级别，摘要显示 level1 x level2 和 level1 x level3 之间的比较，因为 level1 是基线。

我想测试level2 x level3之间是否有区别。我怎样才能做到这一点？也许可能有另一个图书馆。谁能帮帮我吗？

在 R 中，以下公式用于计算 beta 分布中点 0 和 0.5 之间的积分，参数为 10 和 20：

结果是：

如何在 Python 中做到这一点？

I have the plot below using

Now, I want to convert the X values into x*100.

Thanks in advance.