我试图建立一个矩阵A (200,200)，其中A[i][j] = A[j][i]遵循概率p的伯努利分布，如果元素在同一个指数群中(即0-49，50-99,100-149,150-199)，如果不是，则是概率q。
到目前为止，我试图像下面这样直截了当地实现，
import numpy as np
from scipy.stats import bernoulli
def link_weight(p=0.05, q=0.04):
A = np.zeros((200, 200))
for i in range(len(A)):
for j in range(len(A
浏览 3提问于2021-09-07得票数 2
回答已采纳
我想用haskell编写一个bernoulli函数bernoulli:: Integer -> Rational，使用下面的算法计算给定整数的bernoulli数。
定义中使用"frac“和"binom”函数来计算二项式。到目前为止，这就是我所拥有的：
fact :: Integer -> Integer
fact i = foldr (*) 1 [1..i]
binom :: Integer -> Integer -> Integer
binom n k = (fact n) `div` (fact k* fact (n-k))
bern
浏览 4提问于2021-12-01得票数 2
回答已采纳
更新2:我不知道为什么投票仍在进行(2014年3月)。自从很多年前我提出这个问题以来，这个问题似乎已经解决了。确保您使用的是最新版本的boost.
更新:为了格式化数字，可能需要初始化C++流，并且当用加载共享库时，初始化不会发生。
我打电话给你
cout << 1 << "!" << endl;
在通过boost.python导出到共享库的方法中。它不会打印任何东西，但如果我打印出来
cout << "%" << "!" << endl;
它起作用了。
这一点很重要，因
浏览 9提问于2010-12-06得票数 15
回答已采纳
在比较bernoulli_distribution的默认构造函数(50/50的真/假概率)和uniform_int_distribution{0, 1} (一致的概率为0或1)时，我发现bernoulli_distributions至少比uniform_int_distribution慢2倍和6倍以上，尽管它们给出了相同的结果。
我希望bernoulii_distribition能够更好地执行，因为它是专门为两个结果的概率而设计的，是正确的还是错误的；然而，它没有。
鉴于以上和下面的性能指标，伯努利发行版在uniform_int_distributions?上有实际用途吗？
结果超过5次运行(
浏览 6提问于2021-04-14得票数 1
回答已采纳
我是这个论坛的新手，也是新手，所以请原谅我:)。
我想写一个bash脚本来进行CPU压力测试，同时观察时钟的速度。因此，我做了以下工作：
1.)对于CPU压力测试，我使用以下代码创建了一个名为"bernoulli“的脚本：
#!/bin/bash
# argument 1: n
function bernoulli()
{
if (( $1 < 3 ))
then
echo 1
else
echo $(( $(bernoulli $(( $1 - 1 ))) + $(bernoulli $(( $1 - 2 ))) )
浏览 7提问于2016-04-20得票数 1
回答已采纳
4回答
我有原始R码如下所示。
Bernoulli <- rbinom(1000, 1, 0.5)
mix.sample <- rep(0, 1000) #reserve storage
for (i in 1:1000) { #for each Bernoulli realization
if (Bernoulli[i] == 1){ #sample corresponding normal component
mix.sample[i] <- rnorm(1, mean=10, sd=1)
}
else {
mix.sample[i]
浏览 1提问于2018-11-24得票数 2
回答已采纳
我开始使用C++，我正在努力解决一个编译问题
我有一个源文件"binomial.cpp“，其中定义了类的方法：
#include "binomial.hpp"
using namespace std;
int Bernoulli::operator()(){
return (rand() < p*RAND_MAX) ? a : b;
};
int Binomial::operator()(){
int result(0);
for(int i(0);i<n;++i){
in
浏览 1提问于2020-09-10得票数 0
回答已采纳
2回答
我试图实现一个返回n:th bernoulli数的方法，如下所示：
Object subclass: #Bernoulli.
Bernoulli class extend [
"******************************************************
* Psuedo code for bernoulli method I'm working from:
*
* function B(n)
浏览 6提问于2016-05-04得票数 1
回答已采纳
我正在练习使用gganimate，并想要创建类似于但动画的东西。我已经取得了一些成功，但我得到了一个非常奇怪的错误，我实在无法解释，并希望有人能帮助我。这篇文章的其余部分包含了一个完全可复制的例子。
以下是我到目前为止所拥有的：
bernoulli <- cbind.data.frame(c(rep(1,101), rep(0,101)), c(seq(1,0,by=-0.01), seq(0,1,by=0.01)), c(seq(1,0,by=-0.01), seq(1,0,by=-0.01)))
names(bernoulli) <- c("success",
浏览 0提问于2019-07-23得票数 3
回答已采纳
在Pytorch中，您可以执行以下操作：
x = torch.bernoulli(my_data)
在tensorflow中有类似的功能吗？输入可以是二维张量，如(批，连)？
我试过tensorflow.contrib.distributions.Bernoulli：
import numpy as np
tmp_x1 = np.random.rand(20,5)
new_data_2 = tf.convert_to_tensor(tmp_x1)
from tensorflow.contrib.distributions import Bernoulli
tmp2_x1 = Bernoul
浏览 0提问于2018-10-31得票数 0
回答已采纳
4回答
我正在尝试计算不同n的bernoulli多项式。我得到一个错误，我不能迭代整数，我看过不同的帖子，但仍然不能正确。我读到我必须将整数更改为字符串，以便迭代它们，但它仍然不起作用。你能帮上忙吗？
import scipy.special
def cub(x):
return x*x*x
def bernoulli(k,n,x):
a=scipy.special.binom(n, k)
return cub(k/n)*a*x**k*(1-x)**(n-k)
def suma(n,x):
s=0
for str(k) in range (0,n)
浏览 2提问于2019-11-30得票数 0
问题摘要：
为什么我的密度从我的样本如此不同的pmf和我如何进行这个模拟，使pmf和样本的估计是相似的。
问题：
我用scipy模拟了一个独立的伯努利试验的样本。现在，我试图获取我创建的样本的密度直方图，并将其与pmf (概率质量函数)进行比较。我希望密度直方图能显示每个在pmf附近徘徊的两个桶，但相反，我在pmf值之上有2个桶，在5。有人能告诉我如何创建一个密度直方图，而不是针对Bernoulli？我尝试过与其他几个发行版进行类似的模拟，它似乎运行得很好。我在这里缺少了什么，你能告诉我如何操作我的代码来完成这个任务吗？
import numpy as np
import matplotlib
浏览 1提问于2020-08-28得票数 0
回答已采纳
我试着计算第n项，但它给了我错误的答案
import math
def bernoulli(m):
if m == 0:
return 1
else:
t = 0
for k in range(0, m):
t += math.comb(m, k) * bernoulli(k) / (m - k + 1)
return 1 - t
def pn(n, x):
sum = 0
for i in range(n):
sum += ((bernoulli(2 *
浏览 4提问于2021-12-01得票数 1
回答已采纳
当第一行是1，1/2，1/3....这里有一张图片来支持这个问题。
是否存在比朴素O(n^2)方法更有效的方法？
我在研究Bernoulli数时遇到了这一点，随后就达到了“Akiyama-Tanigawa算法”。
其中一种方法是简单地预先计算结果并将其存储在表中。由于Bernoulli数增长非常快，对于大多数实际目的，对于更大的n，我们不需要Bernoulli数。考虑Bernoulli(400)-它大约-(10^550)。
但是仅从算法上看，有比O(n^2)更好的方法吗？
浏览 3提问于2013-04-15得票数 13
如何在不需要转换到特定发行版的情况下从MathDotNet中的发行版中抽取一个样本？
我有一个分布d，它可以是任意随机变量，作为一个IDistribution传递。现在，我想试一试。我想要做到这一点，因为必须尽可能少地对实际发行版本身进行强制转换(我不希望对Bernoulli、Normal等真正特定的发行类型进行大量强制转换的大型case语句。
我已经尝试了以下代码，用于IDistribution d，它的类型为Bernoulli，其平均值为Bernoulli
Console.WriteLine("Mean is " + ((Bernoulli)d).Mean);
Consol
浏览 0提问于2015-07-27得票数 1
我有一个用Django Python框架编写的rest API，我想将进入API的流量以固定的比率(例如。流A的30%和流B的70% ).What是在不使用数据库的情况下实现上述目标的最佳方式。
浏览 19提问于2020-04-15得票数 1
问题
谁知道为什么我可以将groovy闭包保存到JMeter公开的哈希表中，而不是保存一个groovy函数呢？
极小例子
在测试计划中，在线程组中，在脚本语言Groovy3.0.3的JSR223采样器中，我有以下最小的失败示例(当粘贴到JSR223中的" script“框中或从文件名”test.groovy“引用时都会失败)：
def bern = { p -> Math.random() < p }
props.put("bern", bern)
closure = props.get("bern")
log.warn("pro
浏览 0提问于2020-10-12得票数 1
回答已采纳
1回答
我在用具有的中编辑现有的Java实现时遇到了困难。最初的实现将Bernoulli数的计算放在主要方法中。我做了一个新的类来返回单个Bernoulli数的计算。我做错了什么？
import java.math.BigInteger;
public class Bernoulli {
public static void main(String[] args) {
int N = 20;
System.out.println(bern(N));
}
public static BigRational bern(int N) {
浏览 3提问于2015-09-12得票数 0
回答已采纳
1回答
C语言中的随机图
、、、
我有一个关于随机图的作业。我不能理解这个问题。有没有人能向我说明一下我应该做什么？
设N为正整数，p为0- 1之间的数。(N, p)随机图是通过以下步骤生成的图：
绘制N个顶点，分别用1, 2, . . . , N表示；对于每对不同顶点的(u, v)，用概率p将两个顶点连接到一条边。如果在任意两个顶点之间有一条路，则称图是连通的。
在本实验中，您将编写代码来生成大型随机图，并调查此类图的连通性。
我们将N修改为500,000，但让p在{0.05, 0.10, 0.15, ..., 0.95}中变化。对于每个p值，您需要创建100个(N, p)随机图。您需要开发一个方法(当然，在您的程序中实现它
浏览 2提问于2012-11-30得票数 1
3回答
我有一个变量A，它是Bernoulli分布式的，A = pymc.Bernoulli('A', p_A)，但是对于p_A我没有硬值，我想对它进行采样。我知道它应该很小，所以我想使用指数分布的p_A = pymc.Exponential('p_A', 10)。
但是，指数分布可以返回大于1的值，这将放弃A。有没有一种不需要在我自己的p_A -decorated函数中重新实现Bernoulli或指数分布的方法来包围@pymc.stochastic的输出？
浏览 0提问于2015-12-04得票数 1
回答已采纳
3回答
单元测试两位数失败
、、、、
我正在测试方差法(单位)，但出现了双重问题。
public double Variance()
{
return Probablity * (1 - Probablity);
}
[Test]
public void Variance_WithSuccessProbablity02_Returns016()
{
var bernoulli = new BernoulliDistribution(0.2);
Assert.AreEqual(bernoulli.Variance(), 0.16);
浏览 5提问于2012-08-09得票数 1
回答已采纳
试图通过转移“做贝叶斯数据分析”(Kruschke)一书中的一些模型来学习PyMC。一个基本的例子(来自Ch。9)假设一组硬币是按照p~Bern(theta)分配的，其中θ来自具有固定参数的Beta分布(“薄荷”)。下面是我如何对其进行编码(在PyMC2中)：
import pymc as pm
import numpy as np
from matplotlib import pyplot as plt
import seaborn as sbn
from pymc.Matplot import plot as mcplot
from pymc import Bernoulli, Beta
浏览 0提问于2015-03-09得票数 2
回答已采纳
2回答
Python或Numpy中是否有内置库可根据各种常见发行版生成随机数，如：
正常
泊松
指数型
伯努利
还有其他人呢？
有这样的多变量分布库吗？
浏览 8提问于2009-10-31得票数 13
回答已采纳
以下是我对Bernoulli和多项朴素Bayes之间的区别的理解：
Bernoulli显式地建模了一个特性的存在/缺失，而多项式则没有。
我遗漏了什么吗？
有人能直观地解释为什么这种差异很重要吗？也许用一个例子，很明显你会用其中一个而不是另一个
在维基百科的页面上，他们说：
注意，带有Bernoulli事件模型的朴素贝叶斯分类器与频率被截断为1的多项式NB分类器不同。
为什么会这样呢？
浏览 0提问于2018-02-08得票数 10
回答已采纳
1回答
我正在尝试从雅典娜那里获取N行的随机样本。但是由于我想要从中抽取样本的表格很大，所以
SELECT
id
FROM mytable
ORDER BY RANDOM()
LIMIT 100
可能是因为ORDER BY要求将所有数据发送到单个节点，然后该节点对数据进行打乱和排序。
我知道，但它允许对一定百分比的行进行采样，而不是对某些行进行采样。有没有更好的方法来做这件事？
浏览 0提问于2017-06-13得票数 20
回答已采纳
我完全是SAS的新手。我想在我现有的表中添加两列随机生成的0和1。 有人能给我一些反馈吗？ 我知道如何在R中做到这一点，但绝对不是在SAS中。 提前谢谢你
浏览 6提问于2019-04-30得票数 0
回答已采纳
我们认为Bernoulli朴素贝叶斯假设所有连续特征的高斯分布。如果我在数据集中也有分类特性，会发生什么？
在将Bernoulli朴素转换应用于数据集中之前，数据中哪种类型的先验转换是合适的？
浏览 0提问于2018-05-25得票数 2
1回答
指数伯努利抽样
、、、、
Bernoulli是一种概率分布。我需要从指数bernoulli中取样，并返回一个二进制值(即0或1)。我找到了这个算法，我想实现它，但是我不理解算法的第3步：r1 = r1 & (2^h - 1 )。有人能帮忙吗？
浏览 27提问于2022-01-14得票数 0
2回答
我正在尝试使用Microsoft.ML.Probabilistic，示例在C#中。
例如：
static void Main(string[] args)
{
var firstHeads = Variable.Bernoulli(0.5);
var secondHeads = Variable.Bernoulli(0.5);
var bothHeads = firstHeads & secondHeads;
}
现在有了一点魔力，&操作符被C#编译器解析为
Variable<bool>.operator &(Variable<
浏览 2提问于2020-07-09得票数 2
回答已采纳
我试图实现一个返回Bernoulli数的函数：
(defn binom [N K]
(if (or (= K 0) (= N K))
(int 1)
(int (+ (binom (- N 1) (- K 1))
(binom (- N 1) K)))))
(defn bernoulli [M]
(if (= M 0)
(int 1)
(for [K (range 0 (- M 1))]
(* (binom M K)
(be
浏览 4提问于2016-05-08得票数 3
回答已采纳
如何在Bernoulli distribution的tensorflow中生成一个因子为p的张量？ 例如： a = tf.bernoulli(shape=[10,10], p) 生成0-1的矩阵10x10，其中矩阵的每个元素都是概率为p的1，概率为1-p的0。
浏览 63提问于2020-01-16得票数 0
我想使用TABLESAMPLE从满足特定条件的行中随机抽取PostgreSQL。
它运行得很好：
select * from customers tablesample bernoulli (1);
但我不知道如何在脚本中嵌入条件。举个例子
select * from customers where last_name = 'powell' tablesample bernoulli (1);
抛出此错误：
SQL 42601：ERROR:语法错误位于或接近“表样”位置: 71
浏览 0提问于2020-01-27得票数 4
1回答
我试图用不同内部参数的各种实现来创建一个特性：
pub trait ERP {
fn new() -> Self;
fn sample(&self) -> f64;
}
pub struct Bernoulli {
p: f64
}
impl ERP for Bernoulli {
fun new(p: f64) -> Bernoulli {
Bernoulli { p: p }
}
fun sample(&self) -> f64 { self.p } // Filler code
浏览 2提问于2015-07-02得票数 10
回答已采纳
1回答
目标是从已知参数的分布中获得样本。 例如，自定义分布是p(X|θ)，其中θ是K维的参数向量，X是N维的随机向量。 现在我们知道(1)θ是已知的；(2) p(X|θ)未知，但我知道p(X|θ)∝f(X，θ)，f是一个已知函数。 pymc3可以从p(X|theta)中进行这样的采样吗？ 目的不是从参数的后验分布中抽样，而是想从自定义的分布中抽样。 从一个简单的从伯努利分布中抽样的例子开始。我做了以下工作： import pymc3 as pm
import numpy as np
import scipy.stats as stats
import pandas as pd
import the
浏览 34提问于2019-07-01得票数 2
回答已采纳
我创建了两个具有特定属性的随机变量(x和y)。现在，我想从这两个变量中从头开始创建一个数据帧。不幸的是，我输入的内容似乎是错误的。我怎么才能正确地做这件事？
# creating variable x with Bernoulli distribution
from scipy.stats import bernoulli, binom
x = bernoulli.rvs(size=100,p=0.6)
# form a column vector (n, 1)
x = x.reshape(-100, 1)
print(x)
# creating variable y with normal
浏览 4提问于2022-10-22得票数 0
回答已采纳
1回答
我正在制作一个2000次迭代的伯努利实验的直方图。我得到了很多0和1，所以我绘制了结果，得到了一个很好的直方图，频率从0到1000。现在我必须将频率除以2000，然后再次绘制，但我该如何做呢？ 我试着得到值的向量和频率的个别化，但最后我不知道如何把值放在x轴上，把新的频率放在y轴上。 Bernoulli = function(p){
x = runif(1)
if(x <= p){
0
}
else{
1
}
}
Experiment = function(){
vec = c(1,2000)
for(i in 1:2000){
v
浏览 24提问于2019-08-31得票数 0
回答已采纳
1回答
Python函数
、、、
背景
在Julia编程中，我们有一个Rand函数，它是：
兰德(rng=GLOBAL_RNG，S，dims.)
它从S指定的一组值中选择一个随机元素或随机元素数组，并给出第三arg的维数数组。
问题
我想用同样的逻辑来模拟一个保证政策问题：
m = 60000
n = len(policy.AGE)
k=1/10
@time begin
Death = zeros(Bool, m, n)
Accident = rand(Bernoulli(k), m, n)
这是朱莉娅代码，但是当把它带到python时，我找不到rand(Bernoulli(k)，m，n)的正确形式。
我知道
浏览 2提问于2020-11-01得票数 2
回答已采纳
我用PyMC编写了以下程序：
import pymc
from pymc.Matplot import plot as mcplot
def testit( passed, test_p = 0.8, alpha = 5, beta = 2):
Pi = pymc.Beta( 'Pi', alpha=alpha, beta=beta)
Tj = pymc.Bernoulli( 'Tj', p=test_p)
@pymc.deterministic
def flipper( Pi=Pi, Tj=Tj):
浏览 2提问于2014-08-14得票数 1
回答已采纳
一个Bernoulli分布等价于一个二项分布，只有一个试用版，即Apache库中的BinomialDistribution(1,p)。一个Bernoulli分布显然比一般的二项分布简单得多。
如果性能很重要，我是应该通过子类BernoulliDistribution(p)来滚动自己的AbstractIntegerDistribution，还是只使用BinomialDistribution(1,p)就可以获得几乎相同的性能？
浏览 1提问于2015-02-26得票数 1
回答已采纳
#include "Functions.hpp"
#include <random>
#include <algorithm>
#include <iostream>
using std::numeric_limits;
using std::cout;
using std::endl;
std::deque<bool> retBernoulli(double p, int seed, int size)
{
std::bernoulli_distribution dist(p);
std::mt19
浏览 5提问于2022-11-03得票数 0
2回答
cancer <- read.csv('breast-cancer-wisconsin.data', header = FALSE, na.strings="?")
cancer <- cancer[complete.cases(cancer),]
names(cancer)[11] <- "class"
cancer[, 11] <- factor(cancer[, 11], labels = c("benign", "malignant"))
library(gbm)
浏览 6提问于2014-06-02得票数 3
回答已采纳
2回答
我目前正在寻找如何从表中获取随机数据样本的最佳方法(例如，在单元格中)。我知道PRESTO提供了随机()函数或TABLESAMPLE BERNOULLI/SYSTEM。问题是当查询表中有大量记录时，需要花费大量的时间，这不适合与JayDeBeApi合作，后者可能会关闭等待响应时间过长的连接。
我更喜欢使用TABLESAMPLE BERNOULLI/SYSTEM，它以要获取的记录的百分比作为参数。与允许传递精确百分比(即0.003123412% )的ORACLE、SAP或MSSQL数据库相比，PRESTO不允许您使用，尽管功能非常相似，而且所有的转换都在1-100%的范围内。
有人知道解决这个问
浏览 7提问于2020-09-03得票数 2
回答已采纳
我正在使用OMNeT++ 5.1.1模拟器。我已经从bernoulli()函数中看到了一些奇怪的行为，所以我构建了一个MWE来查看发生了什么。
MWE的工作方式是创建一个具有单个节点的网络，并在t=0上设置一个定时器(self message)。当计时器关闭时，模拟以成功概率p运行一些n个Bernoulli试验，值n和p是通过我使用Register_PerRunConfigOption()宏定义的每次运行配置选项指定的。
这是我的代码：
#include <math.h>
#include <omnetpp.h>
using namespace omnetpp;
R
浏览 2提问于2017-08-25得票数 2
回答已采纳
2回答
自举二项分布
、、
我正在处理一个简单的问题，我想写一个函数，它将从一个二项分布(比如60%1和40%0)中随机抽取n个大小为q的样本，然后将每个样本的1份额保存到一个数组中(这样我以后就可以分析那个数组)。例如，如果我从一个具有60% 1和40% 0的分布中抽取了1000个样本，每个样本的大小为30，我希望输出一个长度为1000的数组，其中每一行表示该样本中1的份额。
浏览 21提问于2020-02-18得票数 1
回答已采纳
在试图通过pymc中的Bayesian方法对黑客进行攻击的过程中，我看到了以下代码：
first_coin_flips = pm.Bernoulli("first_flips", 0.5, size=N)
我尝试用以下方法将其转换为pymc3，但它只是返回一个numpy数组，而不是张量(?)：
first_coin_flips = pm.Bernoulli("first_flips", 0.5).random(size=50)
大小之所以重要，是因为后来在确定性变量中使用了它。以下是我到目前为止所掌握的全部代码：
import pymc3 as pm
impor
浏览 0提问于2018-11-07得票数 3
回答已采纳
我正在为一个分类运行gbm模型problem.Below是我的代码和输出
library(gbm)
library(caret)
set.seed(123)
train=read.csv("train.csv")
gbm_model= gbm(DV~.,
data=train,
distribution = "bernoulli",
n.trees = 9,
interaction.depth = 9,
浏览 1提问于2015-06-10得票数 0
1回答
我想实现以下函数：
def bernoulli_tensor(shape, probability) -> tf.Tensor:
...
浏览 2提问于2020-02-24得票数 0
我想知道为什么scipy的随机变量类stats.bernoulli以dtype int32为样本产生ndarray：
> stats.bernoulli.rvs(0.3, size=10)
array([0, 1, 0, 1, 0, 1, 0, 0, 0, 1]
> stats.bernoulli.rvs(0.3, size=10).dtype
dtype('int32')
对二进制结果使用32位整数值似乎效率极低。我本来希望有一个np.bool或np.int8的d型。
有人知道将int32数组作为示例生成的原因吗？
备注：我正在处理大样本(10^8)。奇怪
浏览 2提问于2017-11-17得票数 1
我试图使用查询中的计数在预置中对表进行示例。所以有效地，我试图做这样的事情。我有一个表，我想根据标准进行筛选，并使用筛选过的表的总数从另一个表中进行抽样。我不认为下面的查询是以任何形式或形式工作的，但实际上这正是我所要做的。有人能帮我吗？
WITH less_100_table AS (Select * from Table A where my_total < 100),
WITH greater_100_table AS (Select * from Table A where my_total > 100),
total_count AS (Select count(*) f
浏览 9提问于2022-10-08得票数 0
回答已采纳
当使用logistic回归时，我们使用交叉熵作为损失函数。然而，根据我的理解和https://machinelearningmastery.com/cross-entropy-for-machine-learning/，交叉熵计算如果两个或多个分布是相似的。假设分布为Bernoulli分布或Multinoulli分布。
那么，我的问题是:为什么我们总是可以在回归问题中使用交叉熵，即Bernoulli？实际值和预测值是否总是遵循这样的分布？
浏览 0提问于2021-11-08得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
