问题标签 [ecdf]

我正在尝试在 matlab 中进行生存分析，并想计算几条曲线之间的对数秩测试分数。我在这里找到了一个可能的代码来进行日志排名。但根据它的描述，它只能在两组之间进行对数秩检验。如果我有 3 条或更多曲线怎么办？

我从文件中读取参数“data1”和“data2”并使用此代码绘制 cdf，但我有两个问题：

1. 使图形更平滑
2. 将 Y 轴标记为 CDF

请注意，这段代码是正确的，但我需要一些修改。

编辑：

好吧，对不起，我会尽量说得更清楚，

我有 50 个场景（我在这里随机创建），我将所有这些场景放在一个矩阵中。在我可以应用 ecdf 函数之后，这给了我 50 个 ecdf 的列表。我想从我的 50 个场景的所有这些 ecdf 中计算分位数 90 和 10 以及中位数。

这是一个基本代码：

所以我可以轻松地绘制所有这些参数，但我只想在图表上表示三个参数（Q10、Q50、Q90）。

编辑：

我找到了它的确切方法，所以我分享它，如果有时有人需要它。

您可以尝试代码，图形非常明确，并且很好地解释了我想要做什么。感谢那些试图帮助我的人！

您可以选择使用分位数和中位数绘制所有 ecdf，或者仅绘制分位数和中位数以使其更清晰。

在 R 中使用 ecdf() （或 Ecdf ）时，我得到的结数约为 500。但是，我需要一些结来反映我在原始数据集中的观察数量（在我的情况下为 300） .

如何调整 ecdf 中的结数？或者，有可能吗？

下面是R代码：

我喜欢 ggplot2 包的stat_ecdf()功能部分，我发现它对于探索数据系列非常有用。然而，这只是视觉上的，我想知道它是否可行 - 如果是的话如何 - 获取关联的表？

请查看以下可重现的示例

我有一系列随机值a <- runif(1000, 0, 10)， 和单个值：b <- 1.5。我想绘制 ecdf，并在其上加上值以查看拥有它的概率。

但是当我这样做时：

点不出现。我想我必须用曲线拟合点才能找到与之对应的 y 轴值，但我真的不知道该怎么做。如果有人可以帮助我。谢谢！

下面是我的代码。我能够在 R studio 中生成 ecdf 图，但当我将它们放在一个闪亮的应用程序中时却不行，如下所示：

服务器.R

用户界面

我究竟做错了什么？

我正在尝试在 ggplot 中绘制两条累积频率曲线，并在某个截止处遮蔽交叉。我很久没有使用 ggplot 了，所以我希望有人可以帮助我解决这个问题。

没有填充区域的情节，看起来像这样......

我使用以下代码创建的...

然后我想要的是在 y 轴上为低于 0.2 的温度遮蔽两条曲线。理想情况下，我希望看到蓝色阴影为蓝色，红色阴影为红色。然后，他们在紫色交叉的地方。

但是，我管理的最接近的是如下...... ]

我在代码中添加了以下内容。

我已经看到了一些人为正态分布密度图着色的例子，这是我改编我的代码的地方。但由于某种原因，我的盒子似乎不想与温度曲线有任何关系。

请帮忙！我敢肯定这很简单，我真的很迷茫并且尝试了一些，产生的结果不如这些令人信服。

非常感谢您看一看。

问题解决了感谢下面的帮助...

从下面运行建议的代码

给...

这几乎是我所追求的解决方案，但最后添加了一个......就像这样

我得到了我所追求的...

我再次设置数据的原因是它只填充了两个区域中最低的 20%。

十分感谢你的帮助 ：-）

我有 6 条曲线描述了以固定价值购买的门票数量的 ECDF。

现在我想插值以在它们之间制作曲线，但遵循下一个公式。例如以 10k 的价格估算 ECDF，应以 F(10k)=αF1+(1-α)F2 为指导；0<α<1 。我绘制的曲线是通过 ecdf 函数得到的。我试图做 interp1 来模拟观察，但它对我来说无法模拟新价格的观察。

所以我尝试了@rayryeng 制作的代码

max_len = max(numel(F1),numel(F2)); F1_interp = interp1(1:numel(F1), F1(:).', linspace(1,numel(F1),max_len)); F2_interp = interp1(1:numel(F2), F2(:).', linspace(1,numel(F2),max_len));

阿尔法 = (10e3 - 9.5e3) / (11e3 - 9.5e3); Fnew = alpha*F1_interp + (1-alpha)*F2_interp;

并得到下一条曲线

但是用 F1 和 F2 绘制曲线似乎并不介于两者之间或两者之间的插值。

我想用两条线创建一个 ecdf 图，我想在其中一条线上添加误差线。

我正在使用此代码

错误栏应该添加到 x 值，但它给了我这个错误：

错误：美学长度必须为 1，或与 dataProblems 长度相同：x + xError, x - xError

我不明白 - 结果是相同的长度。

编辑

我改成了问题，所以它变得更容易了——我认为真正的问题与 ECDF 图和误差线有关。以这段代码为例：

它打印错误栏，但情节完全被破坏了。