3回答
我已经用C++编写了一个模拟，它从特定的概率分布中生成(1,000,000)^2个数字，然后对它们执行一些操作。到目前为止，我已经使用了指数分布、正态分布、伽马分布、均匀分布和泊松分布。下面是其中之一的代码：
#include <boost/random.hpp>
...main...
srand(time(NULL)) ;
seed = rand();
boost::random::mt19937 igen(seed) ;
boost::random::variate_generator<boost::random::mt19937, b
浏览 7提问于2012-04-28得票数 8
回答已采纳
1回答
我对用python pymc3包进行负二项式回归的解释感到困惑。我不知道如何解释GLM中的mu和alpha。这里我有一个简单的向量，我想为自己找到NB回归模型：
# The data
y = [100,200,300,400,50,300,60,89,90,100,100]
data = {'y':y, 'x':[1]*len(y)}
basic_model = pm.Model()
with basic_model:
fml = 'y~x'
pm.glm.GLM.from_formula(formula=fml, data=da
浏览 1提问于2018-10-15得票数 0
回答已采纳
2回答
Boost inversion_chi_square_distribution的输入参数(1-1e-3，2)与Matlab的chi2inv值不同.
有人能解释我做错了什么吗？
我用了伽马反演分布。
例子：
投入：
shape param : n/2
scale : 2
p = 1 - 1e-3;
Matlab -> chi2inv(1-1e-3, 2);
Matlab -> gaminv(1-12-3, n/2, 2);
结果-> 13.8155。这两个函数的结果相同。
C++ Boost
cdf(boost::math::inverse_gamma_distribut
浏览 5提问于2021-12-19得票数 1
1回答
我有两份名单。两者都包括归一化百分比：
actual_population_distribution =0.2、0.3、0.3、0.2
sample_population_distribution =0.1、0.4、0.2、0.3
我希望将这两个列表拟合成伽马分布，然后计算返回的两个列表，以得到KL值。
我已经能得到KL了。
这是我用来计算伽马的函数：
def gamma_random_sample(data_list):
mean = np.mean(data_list)
var = np.var(data_list)
g_alpha = mean *
浏览 4提问于2019-08-05得票数 1
1回答
R中具有NA /0值的γ族glm
、、、、
我想使用一个广义线性混合效应模型。我的数据遵循伽马分布，但包含NA和0值。然而，伽马族不允许我计算这些模型，如果我有0值。有人知道解决这个问题的方法吗？
我听说glmmTMB包允许使用带有负值的伽马发行版，但我在mac上工作，而且似乎无法下载这个包。
当我尝试时，我会得到一个错误代码，上面写着"clang: error: unsupported选项'-fopenmp'“。如果你们中有人有主意的话，那就太好了。
浏览 9提问于2022-07-22得票数 0
1回答
假设我有一些遵循伽马分布的数据，我计算了伽马分布的平均值和标准差。
我还知道数据中有一些我想要排除的异常值(噪声)。
现在，
是否可以采用Mean+ 3西格玛来覆盖几乎(99.7%)所有数据(如正态分布)？我的目标是排除那些不在3西格玛之内的数据，并将其视为噪声。这有意义吗？
浏览 32提问于2019-02-01得票数 0
1回答
生成RTT值
、、
我正在编写一个Java applet，我应该能够在其中模拟两台主机之间的连接。因此，我必须随机生成数据包往返时间。
这些RTT可以从~0到无穷大，但通常在某个平均值附近振荡(即，极大或极小的值非常不可能，但也可能)。我想知道有没有人知道我是怎么做到的？
提前感谢
浏览 1提问于2011-01-12得票数 1
回答已采纳
我正在改变合成图像的信号强度。我需要信号变化在0到0.1之间，但是我需要用伽马分布来做，这样更多的信号落在.01/.02的范围内。问题是，我使用的是2010年版本的Matlab，而没有将gamrnd函数作为其库的一部分的统计工具箱。
所有的帮助都是非常感谢的。
浏览 4提问于2017-01-31得票数 2
1回答
我目前已经开始研究图形化的概率模型。我读过Avi Pfeffer写的关于Figaro概率编程语言的书“实践概率编程”。作为练习，我正在尝试从学习集学习正态分布的参数。不管怎么说，我得到的数据并不完全是我所期望的。
我定义了一个模型，其中正态分布取决于两个参数:均值是另一个正态分布，均值为50，方差为0.01。方差是k=2和theta=2的伽马分布。
我正在进行100次观察，每个观察的值都是100。我用一个重要的抽样算法来推断均值和方差。以下是代码
val mean : Element[Double] = Normal(50,0.01)
val variance: Element[Double]
浏览 8提问于2017-12-19得票数 1
2回答
Mathematica有一个四参数广义逆伽马分布：
在那页上也给出了它的PDF。是否有人在R中实现了密度、分布、分位数和抽样的功能？
我做了一个快速的开始( PDF只是那个页面上转换成R的方程式)，但是如果它已经完成了，我就不会去实现CDF和分位函数了。
是否存在计算任何分布的PDF (PDF的积分)和分位数(反演PDF)的通用函数？
注意，这不是广义逆高斯
还请注意数学页面上的“属性和关系”下拉列表，这似乎意味着它不是任何东西的特例或概括(除了逆伽玛)。
浏览 1提问于2012-06-29得票数 4
回答已采纳
如果泊松分布的均值和方差不相同，例如均值大于方差还是方差大于均值时，要选择的模型？
浏览 1提问于2017-10-07得票数 4
1回答
Python中的加权随机生成
、、、、
我想要生成一个150 k的加权随机分布，其中分布中的权重由Python中的F(v) (见Weibull分布的附图)给出。因此，根据这个方程，代码应该生成150 k值，范围从0到30 (即v)。而这150 k值可以是范围内的任何实数，而不仅仅是整数，当这个函数在这个范围内达到峰值时，位于中下区域的数的概率就会更高。注: m，c，k是常数，v是唯一需要150 k值的变量。
浏览 9提问于2022-02-22得票数 0
1回答
我正在寻找一个等价的函数，游戏在boost中。有谁知道它是否在boost库中实现了像gamfit这样的功能？
谢谢
浏览 3提问于2014-06-28得票数 0
回答已采纳
R BTYD和Forecast软件包很棒，我能够预测和预测‘客户频率模式’，参数估计等。有趣的是，这两个软件包只关心‘客户频率’，而不是‘购买量’。但我们的问题是，我们的客户来我们这里的频率是一样的，但他们从我们这里购买的数量没有以前多了。是否有需要同时考虑频率和购买金额的R包？
感谢你的忠告。
浏览 1提问于2013-12-18得票数 0
2回答
在Java中生成概率分布
、、、
我想知道如何生成概率分布。我正在研究随机最短路径问题，其中边具有关联的概率分布和每个概率的关联成本。我能够像这样生成一个(正态)分布：
0 1 2 3 4
0.15 0.2 0.18 0.22 0.25
使得通过以下中提供的答案，所有概率的和为1。现在，我需要生成其他分布，如双正态分布、对数正态分布和伽马分布。我真的很感激任何关于这些发行版和(伪)代码的澄清，最好是用Java来生成它们。
浏览 12提问于2015-11-16得票数 0
1回答
我想知道我的数据点的分布，所以我首先绘制了我的数据的直方图。我的直方图如下所示：
其次，为了使它们适合一个发行版，下面是我编写的代码：
size = 20000
x = scipy.arange(size)
# fit
param = scipy.stats.gamma.fit(y)
pdf_fitted = scipy.stats.gamma.pdf(x, *param[:-2], loc = param[-2], scale = param[-1]) * size
plt.plot(pdf_fitted, color = 'r')
# plot the histog
浏览 2提问于2015-03-23得票数 6
回答已采纳
2回答
我有一个伽马分布适合我的数据使用libary(fitdistrplus)。我需要确定一种定义x值范围的方法，它可以被“合理地”期望，类似于使用带正态分布的标准差。
例如，两个标准差与平均值之间的x值可以看作是正态分布期望值的合理范围。对于如何根据伽马分布的形状和速率参数来定义类似的期望值范围，有什么建议吗？
...maybe，类似于识别包含95%数据的x的两个值吗？
浏览 7提问于2015-10-09得票数 0
回答已采纳
我正在尝试使用以下代码从逆伽马分布中采样，修改我在网上找到的一些代码。我不是C++编码方面的专家，所以我需要你的一些解释和帮助。
#include <iostream>
#include <fstream>
#include <random>
#include "boost/random.hpp"
#include "boost/generator_iterator.hpp"
#include <boost/math/distributions.hpp>
using namespace std;
int mai
浏览 5提问于2015-12-12得票数 0
我一直在使用amt包中的现有函数random_steps。我需要从我的数据中产生一个与观测值相匹配的伽马分布，然后从这些值中进行采样，以便为我提供可能的替代方案。问题是观察值的范围是0到53，但生成的值是0到522。显然，这个函数给我的值不仅是不可信的，而且是不可能的。
我不想翻阅源代码()来寻找要修复的东西，所以我希望有人能给我一个解脱，这样我就可以把amt抛在脑后。然而，我在其他地方找不到解决方案。当然，有一种简单的方法可以使用现有的值来生成分布，然后从中采样？
这些是在数据中观察到的值(左)和由random_steps生成的值的密度图(右)。
我的数据集有200多万行，这给我带来了
浏览 0提问于2021-03-08得票数 0
我正在做乳房X光检查。我想在matlab中使用带有伽马分布的k-means聚类来分割乳房X光片。我只知道应用k均值聚类，但不知道如何使用gamma分布来实现它。
浏览 3提问于2013-03-20得票数 0
1回答
γ-γ累积分布函数
、、、
我想要计算伽玛民防的一系列数据，我有。我已经计算了α和β参数，但是我不知道如何计算R中的of (是否有类似于Matlab的gamcdf?)
我见过一些人使用fitdistr或pgamma，但我不知道如何将alpha值和beta值放在一起，或者根本不需要它们？
谢谢。
浏览 0提问于2018-05-11得票数 2
回答已采纳
有没有一个Java库可以用来将给定的数据向量拟合到Gamma分布中，以便找出它的形状和比例。就像R中的fitdistr(data,'gamma')？
我只能找到根据给定形状和比例参数计算样本分布的库。
浏览 2提问于2013-09-24得票数 3
我需要找到伽马累积分布的倒数。我知道在excel中有GAMMA.INV(probability,alpha,beta)函数可以做到这一点。我如何在R语言中实现这一点？
浏览 0提问于2012-04-14得票数 3
回答已采纳
1回答
两个随机分布的差异
、、、
日安！我有两个伽马分布，想找出它们差异的分布。使用np.random.gamma按参数生成分布，但生成的分布随时间变化很大。代码：
import numpy as np
from scipy.stats import gamma
for i in range(0, 10):
s1 = np.random.gamma(1.242619972, 0.062172619, 2000) + 0.479719122
s2 = np.random.gamma(456.1387112, 0.002811328, 2000) - 0.586076723
r_a, r_loc,
浏览 0提问于2018-09-04得票数 0
我正在尝试使用pymc来查找时间序列中的变化点。随着时间的推移，我看到的值是“转换”的概率，这个值非常小，平均为0.009，范围为0.001-0.016。
我给这两个概率一个均匀的分布，作为零和最大观测值之间的先验。
alpha = df.cnvrs.max() # Set upper uniform
center_1_c = pm.Uniform("center_1_c", 0, alpha)
center_2_c = pm.Uniform("center_2_c", 0, alpha)
day_c = pm.DiscreteUniform("da
浏览 5提问于2017-04-10得票数 0
1回答
我试图在R中创建一个GLMM，我想了解蝙蝠的出现时间如何取决于不同的因素。在这里，我把各自蝙蝠的离开和当天日落之间的时间差作为因变量(公制)。作为固定因素，我想包括不同的天气数据(公制)以及蝙蝠的繁殖状态(分类)。此外，还有应答器号码(个人识别码)作为排除蝙蝠之间个体间差异的随机因素。
我最初是用线性混合模型(package lme4)在R上工作的，但是残差的QQ图与正态分布有很大的偏差。另外，数据的直方图表示伽马分布。因此，我实现了一个具有伽玛分布的GLMM。下面是一个有一个天气参数的例子：
model <- glmer(formula = difference_in_min ~ re
浏览 4提问于2021-03-04得票数 0
回答已采纳
1回答
我想将R包中实现的发行版应用到rstan中。我知道如果直接实现的话可以使用定制的发行版： https://mc-stan.org/docs/2_18/stan-users-guide/custom-probability-functions-chapter.html 但是，我想知道是否可以使用其他包中的发行版。例如，我想使用VarianceGamma包中的方差-伽马分布。 dvg(x, vgC , sigma , theta , nu ) 如何使用dvg(x, vgC , sigma , theta , nu , log = TRUE)中已有的实现来定义对数密度函数，例如vg_lpdf()，
浏览 24提问于2020-05-03得票数 1
回答已采纳
1回答
当我加载一个外部1024x1024 RGBA32 .png (通过PaintXP或Gimp保存)时，中间有一个(64,64,64)像素的点滴，为什么会在底部返回不正确的值？--我能得到的最接近的是一个未压缩的.png (来自Gimp)，它的值为(65,66,65)，但标准图像似乎返回为(56,56,56)。
Texture2D tex = null;
byte[] fileData;
if (File.Exists(mapPath + "/map.png"))
{
fileData = File.ReadAllBytes(mapP
浏览 5提问于2017-09-02得票数 3
回答已采纳
3回答
关于颜色功能的问题
、、、
我想帮助您理解以下代码行：
return float3(pow(color.r, 2.2), pow(color.g, 2.2), pow(color.b, 2.2));
我见过一些地方称它为伽马校正，但我不知道它是否正确。
如果它真的是伽马校正，它是如何工作的？
📷
📷
经修正
我喜欢它的艺术效果，我会保留它。我问这个问题是为了更好地理解。
浏览 0提问于2021-03-21得票数 2
1回答
有人能帮我从python中的伽马分布中生成随机数吗?我试过这两种可能性，但我仍然想知道它们之间的主要区别:第一种是：
shape, scale= 0.5,1
size=(1024,10)
np.random.gamma(shape, scale, size)
第二个是：
from scipy.stats import gamma
gamma.rvs(0.5, 1, (1024,10))
我认为这两种方法都被用来产生服从伽玛分布的随机样本，那么这些语法之间的区别是什么呢？什么时候我们应该使用第一种方法和第二种方法？
浏览 2提问于2020-06-14得票数 0
回答已采纳
1回答
我希望使用键控哈希机制来创建一些令牌来表示我们的单个客户。我计划在令牌中使用客户个人信息(例如H(Firstname,Surname,DOB))的组合来防止频率分析，但是有什么方法可以量化吗？
例如，在定性上，我假设使用：
H(Firstname,Surname,DOB,FirstPurchaseDate)
将比仅使用以下内容更安全：
H(Firstname)
由于后者将提供一些通过频率分析可以识别的值，给出了足够大的样本(为此我必须开放协议，所以我假设攻击者知道这些令牌由哪些变量组成)。我想知道什么样的变量组合会使频率分析变得不切实际--对此有什么既定的技术吗？
浏览 0提问于2019-08-23得票数 0
1回答
R中Erlang分布的样本
、、、、
有两个参数:自然数k(形状参数)和实数λ(速率参数)。如何从Erlang分布中获取大小为n的随机样本，在没有外来包的情况下使用R？
浏览 2提问于2018-03-16得票数 4
回答已采纳
1回答
当您采样的数据的精度远远高于实际情况时，天真的模式实现是无用的(每个示例很可能是唯一的)。例如，在联网的机器上取样往返时间。CPU时钟的潜在精度相当高。如果你只关心精度下降到1ms左右，并且你在从Pmax到Pmin的范围内取样，那么有什么可靠的方法来测量其中“最常见的”ping呢？
浏览 1提问于2014-07-14得票数 0
回答已采纳
I具有偏正态分布随机数的均值、方差和偏度，i具有任意随机数分布的平均方差和偏斜性。
我怎么能有相同的均值和方差，但不同的偏斜度？
浏览 0提问于2011-01-10得票数 2
2回答
我得到了一个在任何R包中都找不到的概率分布，并被告知从它生成随机值。概率分布函数称为LogGamma分布，具有如下概率分布函数。
f(x) = \frac{\alpha^\beta}{\Gamma(\beta)}(\ln{x})^{\beta-1}x^{-\alpha-1}；
因为我找不到这个表达式的积分逆的显式表达式，所以我创建了一个函数，它遍历一系列值，直到到达返回期望概率的停止点为止。
LogGammaPDF <- function(x){
return(alpha^beta * log(x)^(beta - 1) * x^(-alpha-1)/gamma(beta))
浏览 4提问于2017-03-14得票数 3
1回答
Python单边KS测试
、、、
针对大量的理论概率分布，我对我的分布进行了单向的KS检验(观测到了一个从0到100之间的公共交通网格的占用情况)：
cdfs = [
"norm", #Normal (Gaussian)
"alpha", #Alpha
"anglit", #Anglit
"arcsine", #Arcsine
"beta", #Beta
"betaprime",
浏览 6提问于2019-12-20得票数 2
回答已采纳
有没有统计方法来评估2个分布的切线关系？
If 2 distribution are seperated clearly, result is bigger;
if distributed like comb, result is smaller.
import matplotlib.pyplot as plt
plt.show()
2具有不同y轴值的分布0，1
浏览 3提问于2014-01-17得票数 1
1回答
我有两张相似的图片，其中一张比另一张更亮。
我试图得到一个颜色或伽马值来处理，并发现以下Imagick方法没有结果(0)：作为最终结果，我想调整两个图像在相同的亮度。我希望是0.5或相关的值。怎么啦？
<?php
$s1=new Imagick('01.jpg');
$s1->getImageGamma (); // Output 0
?>
浏览 2提问于2015-09-23得票数 2
从双变量正态分布和学生T分布生成随机样本的最佳方法是什么？在这两种情况下，sigma都是1，均值为0-所以我真正感兴趣的唯一参数是相关性(以及学生t的自由度)。我需要在C++中获得解决方案，所以不幸的是，我不能使用MatLab或Mathematica中已经实现的函数。
浏览 0提问于2010-01-25得票数 4
回答已采纳
1回答
实际上，我正在研究GAMLSS模型(用于位置、规模和形状的通用加法模型)。我的问题是:正确地说，这些模型是GLM和线性回归的推广？如果使用带指数族分布的GAMLSS和σ上的同一性链接函数，则得到GLM？作为特例，如果我使用正态分布和两个Idenity链接函数，我得到一个线性回归？
感谢每个人。
浏览 9提问于2022-02-05得票数 0
我正在寻找一个简单的伽马校正公式，灰度图像值在0到255之间。
假设我的屏幕的伽马值是2.2 (这是一个LCD屏幕，所以我可能需要用一个更复杂的过程来估计它，但假设我的屏幕表现得很好)。
下列哪一个公式是正确的？
Corrected = 255 * (Image/255).^2.2
或
Corrected = 255 * (Image/255).^(1/2.2)
(这些注定是MATLAB的代码，但我希望即使不是MATLAB的人也能理解)
我在网上四处看看，但发现这两种公式都在流传。我怀疑(2)是正确的，我的困惑是因为我倾向于把“伽马值”称为实际伽马值的反比，但我真的很感
浏览 5提问于2013-05-13得票数 48
回答已采纳
我需要估计一个截断的伽马分布参数(形状，比例)。
但是，我只知道数据的均值和标准差。dev.我不知道数据集。
给定均值和标准差。dev.对于来自截断伽马分布的数据集，如何找到分布参数的形状和比例？
我知道MLE可能对解决这个问题很有用。但是，它们依赖于对整个数据集的了解。
任何帮助都将不胜感激。
浏览 1提问于2014-03-18得票数 0
2回答
直方图上的泊松分布
、、、、
这个程序应该是把一个泊松分布放在一个绘图图上。但当我运行它时，什么也没有出现。我试着重新排列代码，但它告诉我x和y不在同一维度...
有没有人能弄明白？
from numpy import *
from scipy.special import gamma
import pylab
import matplotlib.pyplot
mu = 70
counts = range(45,101,1)
for counts in range(46,96,1):
p = exp(-mu) * mu**counts / gamma(counts +1)
pylab.plot(counts,p)
p
浏览 0提问于2012-10-29得票数 0
2回答
我一直在研究gamlss软件包来拟合半参数模型，并在ga()函数中发现了一些奇怪的东西。即使模型被指定为具有伽马分布，并使用REML拟合，模型的输出也是高斯的，用GCV拟合。
例子：
library(mgcv)
library(gamlss)
library(gamlss.add)
data(rent)
ga3 <- gam(R~s(Fl)+s(A), method="REML", data=rent, family=Gamma(log))
gn3 <- gamlss(R~ga(~s(Fl)+s(A), method="REML"), data=re
浏览 1提问于2022-04-13得票数 1
1回答
伽马分布
、
让X~gamma(2,1)。
查找：
P(X>=2)
(X<=10)
我在努力想办法解决这个问题。我知道你必须集成一些函数，从2到无穷，从0到10，但是不知道要集成什么。
浏览 1提问于2017-03-08得票数 0
回答已采纳
5回答
有人能帮我在python中拟合伽马分布吗？嗯，我有一些数据:X和Y坐标，我想找出符合这个分布的gamma参数……在中，事实证明fit方法确实存在，但我不知道如何使用它:s.首先，参数"data“必须采用哪种格式，以及如何提供第二个参数(参数)，因为这就是我要查找的内容。
浏览 3提问于2010-05-24得票数 32
2回答
我有一份数字清单。以下是一些基本的统计数据：
N > 1000
Max: 9.24
Min: 0.00955
Mean: 1.84932
Median: 0.97696
似乎数据是正确的，即许多小数字和一些非常大的数字。
我想找一个分布来概括这些数字。我认为正态分布，伽玛分布和拉普拉斯分布看起来都是可能的。我如何确定哪个发行版是最好的？
浏览 0提问于2014-08-10得票数 0
这里使用ggplot2和函数stat_function()编写了这段代码。我用它来创建一个样本正态分布，平均值为100，标准差为15 (IQ)。
ggplot(data.frame(x = c(0, 200)), aes(x)) +
stat_function(fun = dnorm,
geom = "area", fill = '#1e466e', alpha = 0.5,
args = list(
mean = 10
浏览 2提问于2022-01-11得票数 0
1回答
PyMC中的复合模型
、、、、
我试图使用PyMC 2.3来获得参数的估计值。所谓“复合”，我指的是一个随机变量，它跟随一个分布，其参数是另一个随机变量。(“嵌套”或“层次结构”有时用于指这种情况，但我认为它们不那么具体，在这种情况下会产生更多的混乱)。
让我们做一个具体的例子。“实”数据是由一个复合分布生成的，它是一个具有指数分布的参数的Poisson分布。在“普通”中，数据生成如下：
import numpy as np
from scipy.stats import distributions
np.random.seed(3) # for repeatability
nsamples = 1000
tau_tru
浏览 0提问于2014-10-13得票数 0
回答已采纳
1回答
关于声调和伽马校正，我没有什么疑问。
声调映射操作符是一个函数，其输出是0,1之间的值吗？如果不是这种情况，并且由于伽马编码需要一个值在范围0,1中，那么我是否需要在声调映射之后对值进行规范化？
想象一个漫射水平面，反射率为(0.5，0.5)，在无遮挡的均匀天空下(1,1,1) cd/m2。从这个表面来的图像像素上的入射亮度值为(0.5,0.5,0.5)。现在，我看到的大多数代码在保存到文件之前的这个阶段应用了伽马编码(将pxl值提高到0.45)。但是在观看图像时，颜色是扭曲的，并且比(186,186,186)颜色值要亮得多。如果我不应用伽马编码，那么图像文件将显示预期的颜色，即(127,12
浏览 0提问于2020-04-30得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
