我想知道如何从参数为p的伯努利随机变量的另一个泊松过程中模拟。
为了在区间[0,t]上模拟参数为\lambda的第一泊松过程,通常
pois = rpois(1, \lambda)
v = runif(pois, O, t)
w = sort(v)
现在,我知道我们可以将一个伯努利随机变量与第一个泊松过程的到达时间联系起来,来模拟另一个参数为p * \lambda的泊松过程,但是怎么做呢?
我正在开发一个网络模拟器,其中数据包到达和传输尝试的事件遵循泊松分布。我有一个对Knuth算法的改编:
public class Poisson {
private double λ;
private Random rand;
/** Creates a variable with a given mean. */
public Poisson(double λ) {
this.λ = λ;
rand = new Random();
}
public int next() {
double L
我正在尝试对通常符合泊松过程的事件的时间戳进行建模,其中均值=方差。我使用下面的代码实现了这一点,它从一个指数中采样(对于泊松,间隔时间计数是指数),它似乎工作得相当好。 def GenerateTimes(self):
intervals = [random.expovariate(self.CountsPerSecond) for i in range(self.nDataPoints)]
timeStamps = [0.0]
timeStamp = 0.0
for t in intervals:
问:一个有效的零充气准泊松模型能在R中拟合吗?
我知道零充气泊松和零充气负二项环都可以适用于每个pscl::zeroinfl()和glmmADMB::glmmadmb()。
然而,在pscl::zeroinfl()的帮助文件示例中,准泊松是在没有通货膨胀的情况下安装的,但在通胀演示中却被忽略了:
## data
data("bioChemists", package = "pscl")
## without inflation
## ("art ~ ." is "art ~ fem + mar + kid5 + phd + ment
在模拟一个基于泊松分布的随机点过程时,它包含1000个点,它们似乎都占据了窗口中心的一个小区域。
利用Donald逆采样算法实现了基于泊松的伪随机数发生器.
Lambda值(也称成功率)被设置为window_维度/2,并获得了这个结果(截图)
代码:
public double getPoisson(double lambda) {//250
double L = Math.exp(-lambda);
double p = 1d;
int k = 0;
do {
k++;
p *= Math.random();
我正在运行一个Fortran代码,它执行一个标记泊松集群过程的随机模拟。在实践中,事件属性(例如,发生时间)通过反转方法,即通过累积分布函数的随机采样来生成。由于泊松随机性,我希望每个生成的序列都是不同的,但事实并非如此。我猜原因是伪随机数生成器的种子在每次模拟时都是相同的。我不懂Fortran,所以我不知道如何解决这个问题。下面是伪随机数生成器的代码部分,你有什么想法吗? subroutine pseud0(r)
c generation of pseudo-random numbers
c data ir/584287/
data ir/574289/
我试图创建一个泊松模拟使用rpois()。我有两个小数位利率的分布,我想知道它们是否有泊松分布,而不是正态分布。
rpois()函数返回正整数。我希望它返回两个小数位正数。我试过以下几种方法
set.seed(123)
trialA <- rpois(1000, 13.67) # generate 1000 numbers
mean(trialA)
13.22 # Great! Close enough to 13.67
var(trialA)
13.24 # terrific! mean and variance should be the same
head(trialA, 4)
我试图在Matlab中模拟车辆到达交叉口的过程。车辆按泊松分布随机生成。
让我们说,在一个肮脏的地方,交通流量是每小时600辆车的强度。据我从理论上了解,泊松分布的λ应该是600/3600 (1小时内3600秒)。
然后我运行这个循环:
for i = 1:3600
vehicle(i) = poissrnd(600/3600);
end
有一个问题:当我计算数组vehicle中的“1”时,从来没有600个,它总是一些数字,比如567,595等等。
问题是,我是不是做错了,也就是兰巴应该有所不同吗?或者,数字永远不会相等是正常的吗?