问题标签 [bernoulli-probability]

我正在练习使用gganimate并想创建类似于这个闪亮应用程序但动画的东西。我取得了一些成功，但遇到了一个非常奇怪的错误，我真的无法解释，希望有人能帮助我。这篇文章的其余部分包含一个完全可重现的示例。

这是我到目前为止所拥有的：

这将创建我正在使用的数据框。这是三个变量，非常简单。我的动画的逻辑是二分结果的概率如何随伯努利分布的p参数而变化。举一个静态图的例子，如果我这样做：

并改变p我的子集，我得到了想要的结果（尽管由于某种原因，有某些数字，看似完全随机，我无法对数据进行子集，即使我可以在数据框中清楚地看到它们。我已经检查了变量是数字等等，所以这不是问题。这是一个附带问题，我不知道它为什么会发生）。

例如，在 p == 0.12 和 p == 0.70 之间切换，并将它们与前面提到的 Shiny 应用程序进行比较，您会看到它匹配，当然轴会有所不同。

当我尝试按如下方式实现动画时：

一切看起来都很棒！轴匹配得很好，动画流畅......但是一旦条达到 p == 0.50，它们都跳到 1，这根本没有反映在数据中。这只是我的计算机/R/ggplot/gganimate 的问题吗？

我正在 STAN（rstan 库）中拟合逻辑模型。我的响应变量没有任何缺失值，但是我的协变量“HB”之一是二进制的并且缺少条目。

因此，目标是在每次迭代时使用伯努利先验（参数为 0.5）来估算二进制向量中缺失的条目。

但是，我遇到了问题：

• 丢失的数据需要在参数或转换参数块中声明为实数或向量；
• 模型块中伯努利分布的实现需要是整数；
• 据我所知，STAN 中没有将实数或向量转换为整数的功能。

我使用了STAN 用户指南第 3.3 节中提供的指南。对于下面的模型，解析器在伯努利赋值行（模型块中的倒数第二行）给我一个错误，说它需要整数。注意：我还尝试在参数块中将 HB_miss 定义为实数并得到相同的错误。

有什么想法可以在STAN中有效地在HB_miss之前分配伯努利吗？

我想从以下df模拟收入场景：price和（估计概率）：est_p

收入是 - s 的总和，price其中：actual_saleTRUE

我创建了一个函数来模拟伯努利试验est_p和price值：

并将其应用于样本df：

看起来不错，实际值 = 中值！但是当我将相同的函数应用于增加的（复制 x 1000 次）样本 -df1000时，实际收入超出了模拟值的范围：

为什么实际收入超出模拟值范围？我在哪里犯了错误，这个问题的正确解决方案是什么？

我有一组伯努利变量，给出具有不同概率的特定值。变量是独立的。我正在尝试为所有可能的结果建立一个简单的离散概率表。我拥有的数据的一个简短示例是：

我试图找到的解决方案是：

其中值是可能值的总和，概率是该值的概率。

我正在使用的示例来自 Excel 工作表。我正在处理的表格是一长串独立测试。每个测试都有一个可能的成功值、一个成功概率和一个不成功值（1 - 失败概率）。概率表是一个计算每个可能结果的概率的表 - 可能值（对该结果的值求和）和该结果的概率。所以第一个可能的结果 3.2 = 1.7 + 1.5 的概率为 0.42 = 0.7 * 0.6。第二个结果是 2.3 = (1.7 + 0.6)，概率为 0.28 = (0.7 * (1 - 0.6)，依此类推。

我有一组伯努利变量，给出具有不同概率的特定值。变量是独立的。我正在尝试为所有可能的结果建立一个简单的离散概率表。我拥有的数据的一个简短示例是：

其中值是可能值的总和，概率是该值的概率。我正在使用的示例来自 Excel 工作表。我正在处理的表格是一长串独立测试。每个测试都有一个可能的成功值、一个成功概率和一个不成功值（概率为（1 - 成功概率））。概率表是一个计算每个可能结果的概率的表 - 可能值（对该结果的值求和）和该结果的概率。所以第一个可能的结果 3.2 = 1.7 + 1.5 的概率为 0.42 = 0.7 * 0.6。第二个结果是 2.3 = (1.7 + 0.6)，概率为 0.28 = (0.7 * (1 - 0.6)，依此类推。

所以我想要得到的解决方案是这样的 (2.29 = 2.3, 0.899 = 0.9:

如何生成带有因子 的张tensorflow量？Bernoulli distributionp

例如：

a = tf.bernoulli(shape=[10,10], p)

生成一个10x100-1 的矩阵，其中矩阵的每个元素为 1 概率p和 0 概率1-p.

我想实现以下功能：

给定两个不同的高斯，我如何创建两个高斯的混合，其中每个数据点来自概率 p = 0.3 的第一个高斯和概率 0.7 的第二个高斯（使用概率的伯努利分布）。

这是我的开始：

为了绘制我的逻辑函数，我使用了：

plot(plogis, from = -10, to = 10) 我的 x 值为：-10 到 10

那么，根据该图的输出，我接下来应该做什么？

编辑：如果我绘制逻辑函数，y = logistic(x) = 1/(1+exp(−x))，我该如何绘制 logit 函数，y = log(x/(1−x)) ，使用逻辑函数的 y 值的输出？

我尝试在 R 中使用 rbern 函数，但我意识到对角线条目并非全为 0。