weibull

    2热度

    2回答

    我想在给定间隔内生成一些威布尔随机数。例如,来自威布尔分布的具有形状2和在区间(0,10)中的比例30的20个随机数。 rweibull R中的函数产生具有给定形状和比例值的威布尔分布的随机数。有人可以提出一种方法吗?先谢谢你。

    1热度

    1回答

    我有一些简单的事件数据时间,没有协变量。我试图将威布尔分布拟合到它。所以我有以下代码。一切都很好,直到我加载我的首字母缩写。它说“这个链包含未初始化的变量”。但我不明白。我认为Weibull dist只有2个参数,我已经指定了它们。您能否提一些建议?谢谢! model { for(i in 1 : N) { t[i] ~ dweib(r, mu)I(t.cen[

    3热度

    2回答

    我想估计3p威布尔分布的尺度,形状和阈值参数。 我到目前为止已经做的是以下几点: 指的这个帖子,Fitting a 3 parameter Weibull distribution in R 我用过的功能 EPS = sqrt(.Machine$double.eps) # "epsilon" for very small numbers llik.weibull <- function(sha

    2热度

    1回答

    使用SciPy,我试图从this question重现威布尔拟合。 import numpy as np from scipy.stats import genextreme import matplotlib.pyplot as plt data=np.array([37.50,46.79,48.30,46.04,43.40,39.25,38.49,49.51,40.38,36.98,4

    3热度

    1回答

    我有一个数据集,我想适合已知的概率分布。目的是在数据生成器中使用拟合的PDF - 这样我就可以从已知的(合适的)PDF中采样数据。数据将用于模拟目的。目前我只是从正态分布抽样,与实际数据不一致,因此模拟结果不准确。 我第一次想用下面的方法: Fitting empirical distribution to theoretical ones with Scipy (Python)? 我首先想到的是

    0热度

    1回答

    我尝试用optim()函数估计三个参数a,b0和b1。但我总是得到错误: 优化错误(par = c(1,1,1),fn = logweibull,method =“L-BFGS-B”,: L-BFGS-B需要有限值'fn “ t<-c(6,6,6,6,7,9,10,10,11,13,16,17,19,20,22,23,25,32,32,34,35,1,1,2,2,3,4,4,5,5,8,8,8,8

    2热度

    2回答

    我想创建一个QQ图来测试,如果我的数据可以由Weibull分布使用命令 建模使用数据 x =c(3.367, 0.769,0.8,1,1.2) qqplot(x,'weibull') 我不断收到带有与错误 "In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion" ,并找不出原因。这是否意味着我无法将We

    0热度

    1回答

    我想通过在R 2与一个给定的审查矢量施加到下面的数据来估计威布尔分布的极大似然参数: 数据= 9 2 11 49 7 5 3 36 30 6 62 5 3 29 29 1 13 1 24 11 9 4 7 15 11 15 1 1 1 1 1 2 6 12 12 28 14 14 57 17 4 2 3 6 21 6 16 19 28 18 19 9 59 12 3 27 8 26 19 47 6

    2热度

    1回答

    我试图找出一个随机变量事件超过特定值的概率,即pr(x> a),其中a是某个常数,通常远高于x的平均值,并且x不是任何标准的高斯分布。所以我想要拟合一些其他的概率密度函数,并把x的PDF从a到inf的积分。由于这是尖峰建模的问题,我认为这是一个极值分析问题,并且发现威布尔分布可能是合适的。 关于极值分布,威布尔分布有一个非常“不易实现”的积分,因此我想我可以从Scipy得到pdf,并做一个黎曼和。

    3热度

    1回答

    我正在研究一个涉及组合两个威布尔分布并因此创建双峰值曲线的项目。然后,我的目标是使用这个预测。我在网上搜索,我似乎无法找到任何东西,或者如果R有一个函数,允许我结合两个Weibulls。 下面显示了我用来创建两个Weibull分布的代码,我希望这两个分布组合起来构成一个概率密度函数。 curve(dweibull(x, scale=30.59898985, shape=2.27136646),fr