问题标签 [weibull]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
292 浏览

r - 预测 Weibull 即时危险

我正在尝试使用不同的生存模型对基线危险进行建模。例如,我可以使用 Cox 对累积基线风险进行建模,如下所示:

但是 Cox 对此并不理想,我想模拟即时危险而不是累积,所以我尝试使用 Weibull,如下所示:

但是,我无法弄清楚如何从中预测即时危险。例如,如果我尝试以下代码,我只会在每个时间点得到一个常量值 (38.18681)。

我在哪里错了?如果我不清楚,我的目标是在危险与时间的图中可视化瞬时危险。

0 投票
1 回答
66 浏览

r - 用威布尔概率填充矩阵

我在使用 R 来填充矩阵时遇到问题,该矩阵的每一行具有不同的“比例因子”。此时每列都被相同地填充。

我已经尝试了很多方法,但下面我展示了我所做的。我认为我缺少一些基本的东西!

获取 26 个比例因子的列表:

准备矩阵:

在矩阵的每一行上应用 Weibull 函数(对于每个 x 值(风速))

或者:

我无法让任何一种方法起作用。任何帮助都会非常有帮助!

0 投票
0 回答
156 浏览

r - R EHA 中的事件

我准备了一个数据集来执行事件历史分析。当我计算事件的数量时,

当我跑步时,

我明白了,

我不明白为什么事件计数与 event_info 数据集中的计数不同。

不幸的是,我无法发布数据。

有人可以阐明如何在 eha 中计算事件吗

0 投票
0 回答
54 浏览

r - 在生存分析中纳入非事件

我是生存分析的新手。我试图了解将未观察到任何数据的测量值包含到数据集中是否重要。

我有以下虚构数据:三个患者,

  1. 观察到第一位患者在 12 时患病,
  2. 第二次在时间 13 观察两次,在时间 14 观察一次以无病,并且
  3. 观察到第三位患者在时间 1 有病。

我尝试了以下两个示例

结果表明两者之间存在明显差异,我想知道是否有人可以解释为什么包括患者没有生病的观察结果很重要。谢谢!

0 投票
2 回答
2883 浏览

r - R生存分析:使用威布尔的survreg中的错误

我正在尝试使用我之前在 SAS 中制作的 Weibull 发行版来复制一些生存分析——我现在正在使用未经许可的机器工作,所以我正在使用 R(都来自 Windows)。我的(右删失)输入数据如下所示:

事件和审查员是二进制文件,分数在 150 到 450 之间变化,时间在 1 到 28 之间。输入数据集中有大约 30,000 行。

当我尝试:

我收到一条警告消息:

在 survreg.fit(X, Y, weights, offset, init = init, controlvals = control, : 用完迭代并且没有收敛,

而且没有输出。

我已经在网上(并通过这个网站)搜索了这个味精,但没有找到任何表明问题可能是什么的东西。我在 SAS 中通过(proc 逻辑和)lifereg 放置相同的数据时没有收敛问题。

0 投票
1 回答
1143 浏览

r - 如何绘制多变量survreg生成的生存曲线

这个问题是从这里的问题开始 的,因为我们想要绘制一个多变量模型。

假设我们survreg要从包发行说明中的​​示例中绘制生存函数

模型在哪里

据我了解,以下问题似乎解决了这个问题,但不提供数据集,只使用一个变量。在这个例子中,我们有两个连续变量和一个因子

0 投票
1 回答
780 浏览

r - 模拟竞争风险数据

我的目标是模拟一个可用于测试竞争风险模型的数据集。我只是尝试使用该survsim::crisk.sim函数的一个简单示例,但它不会导致我期望的结果。

我希望这些数字与beta0.ev. 任何指向我可能做错的地方或其他如何模拟竞争风险数据的建议。

为了完成:我希望模拟数据中的事件遵循每个风险不同的 Weibull 分布。我希望能够在数据中指定一个层次和集群。审查可以遵循 Weibull 或 Bernouli 分布。

0 投票
0 回答
118 浏览

r - 具有 0 个标准错误和 0 个朴素标准错误的竞争风险模型

我正在使用相关数据在专有数据集上拟合竞争风险模型。合身有效,但我得到了奇怪的标准结果。我的模型摘要中的错误和 Naive SE。我想知道是否有人遇到过类似的问题并且知道如何解决它。

0 投票
1 回答
7518 浏览

python - 根据数据确定 Weibull 参数

我想确定我的数据的Weibull 参数(即形状和比例)。

我已经尝试过这个答案提出的建议,并且我正在使用 Python 3.4。

这给了我以下输出:

我假设我的 csv 文件中的数据被读取为 x 输入值,而不是 Weibull 函数的 y 值。当我用 bin 添加第二列(或行)时,它会给出一个错误,即字符串值无法转换为浮点数。

我如何需要修改我的 csv 文件才能将其中的数据用作 Weibull 函数的 y 值?

我想我的问题可能是我不明白这一行:

这里代表什么1, 1?那么参数不应该是负数。

0 投票
1 回答
2746 浏览

r - 使用 nls 将曲线拟合到 R 中的威布尔分布

我正在尝试将此数据拟合到 weibull 分布:

我的xy变量是:

情节如下所示:

在此处输入图像描述

我正在寻找这样的东西: 拟合情节

我想给它拟合一条威布尔曲线。我在 R 中使用 nls 函数,如下所示:

这个函数总是抛出一个错误说:

所以首先我尝试了不同的起始值,但没有任何成功。我无法理解如何对起始值做出“好的”猜测。然后我使用SSweibull(x, Asym, Drop, lrc, pwr)了 selfStart 函数。现在 SSWeibull 函数需要 Asym、Drop、lrc 和 pwr 的值,我不知道这些值可能是什么。

如果有人可以帮助我弄清楚如何进行,我将不胜感激。

数据背景:我从 bugzilla 获取了一些数据,我的“y”变量是特定月份报告的错误数,“x”变量是发布后的月份数。