- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在寻找一个内置函数,它可以在 MATLAB 以外的其他东西中返回样本分位数和估计的置信区间(MATLAB 的 ecdf
可以做到这一点)。
我猜 R 有这个内置功能,只是我还没有找到它。
如果您有任何独立的代码来执行此操作,您也可以在此处指向它,不过我希望找到一些包含在更大的开放代码库中的东西。
-试图摆脱 MATLAB。
最佳答案
survfit
函数可用于获取具有置信区间的生存 函数。由于它只是 1-ecdf,因此分位数之间存在直接关系。要使用它,您必须创建一个变量来表明您的每个观察结果都是完整的(未审查):
library(survival)
x <- rexp(10)
ev <- rep(1, length(x))
sf <- survfit(Surv(x,ev)~1)
输出:
>summary(sf)
Call: survfit(formula = Surv(x, ev) ~ 1)
time n.risk n.event survival std.err lower 95% CI upper 95% CI
-1.4143 10 1 0.9 0.0949 0.7320 1.000
-1.1229 9 1 0.8 0.1265 0.5868 1.000
-0.9396 8 1 0.7 0.1449 0.4665 1.000
-0.4413 7 1 0.6 0.1549 0.3617 0.995
-0.2408 6 1 0.5 0.1581 0.2690 0.929
-0.1698 5 1 0.4 0.1549 0.1872 0.855
0.0613 4 1 0.3 0.1449 0.1164 0.773
0.1983 3 1 0.2 0.1265 0.0579 0.691
0.5199 2 1 0.1 0.0949 0.0156 0.642
0.8067 1 1 0.0 NaN NA NA
事实上,survfit
会计算中位数及其置信区间,但不会计算其他分位数:
>sf
Call: survfit(formula = Surv(x, ev) ~ 1)
records n.max n.start events median 0.95LCL 0.95UCL
10.000 10.000 10.000 10.000 -0.205 -0.940 NA
计算中位数置信区间的实际工作很好地隐藏在 survival:::survmean
函数中,您可以使用它来泛化到其他分位数。
关于c# - R、python 或 Octave : empirical quantile (inverse cdf) with confidence intervals?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/3442810/
我有两个 TIMESTAMP 和一个 INTERVAL。我想知道间隔有多少次适合时间戳之间的间隙。 减去两个 TIMESTAMP 得到一个 INTERVAL。 SELECT TIMESTAMP '20
我最近在 Tardos 和 Kleinberg 的算法设计的第 4 章中阅读了有关间隔调度算法的内容。为间隔调度问题提供的解决方案是这样的: Sort the n intervals based on
这个问题在这里已经有了答案: Stop setInterval call in JavaScript (9 个回答) 关闭 4 年前。 文档指出 clearInterval() 需要在 setInt
我希望能够使用一个蓝牙适配器在两个广告数据之间切换,最佳的广告切换间隔是多少,以及要立即被 iOS 识别而没有太多延迟的广告间隔是多少? (我想在 iBeacon 和 GATT 服务广告之间切换)。现
直观上,这两个间隔代表相同的时间量。但区别在于夏令时的变化,在这种情况下,“1 天”在 Spring 可能表示“23 小时”,在秋季可能表示“25 小时”。 我用 PostgreSQL 进行了测试,这
我需要在 Oracle 数据库中生成一个值列表,其中包含以下列数据: ITEM_TYPE VARCHAR2(20) ITEM_LAST_UPDATED DATE ITE
我遇到了一种合并重叠间隔的方法。它有一个部分我不明白: Collections.sort(intervals, new Comparator() { @Override public
我正在尝试从 Analytics Vidhya 做贷款预测的机器学习练习题。当我使用随机森林分类器时,它显示: TypeError:float() argument must be a string
我想为数据帧中的每一行分配一个间隔,这样所有行就不会重叠并覆盖整个可能的范围。因此,我可以根据给定间隔内的值过滤行。 我使用过 pd.Interval,但是当我尝试“正常”过滤时它不起作用: df =
在 Joda-Time 中, Interval.contains(Interval) 的实现看起来像这样: return (thisStart <= otherStart && otherStart
我有一个使用 AngularJS 间隔功能连续运行的毫秒时间。我创建了一个暂停 和继续 按钮。当点击 Pause 按钮时,它将触发 AngularJS 的 Interval.cancel 功能,如果我
我在训练 LR 模型时使用 sklearn2pmml.preprocessing.CutTransformer 和 sklearn.preprocessing.LabelEncoder 对目标进行编码
我想用 interval = alt.selection_interval(encodings=['x']) 在只有一个图表上缩放/调整大小到我选择的间隔(从 mousedown 到 mouseup)
我正在使用kafka 0.10.1.1,并与以下3个属性混淆。 heartbeat.interval.ms session.timeout.ms max.poll.interval.ms heartb
我正在尝试保留数据库标志 30dexpf = 1 直到项目的到期日期等于或介于接下来的 7 天和 30 天之间这是我的数据库表 30dexpf | expiry --------|-------- 0
我正在 Informix 中执行以下操作以删除超过 20 秒的行。 delete from sometable where someDateColumn < (current - interval (
我在 Controller 中有一个基于 setInterval 的函数。我注意到,在第一次调用它之后,它工作得很好,但即使在路线改变后它仍然继续。之后我尝试使用 native $interval 服
我正在尝试创建一个按钮,单击该按钮将启动计时器,第二次单击将停止计时器,并将“按钮”的innerHtml更改为您停止计时器的时间。然而,它不适用于我在下面编写的代码,间隔只是继续进行,似乎每次点击都会
我是编程新手,最近一直在研究 AngularJS。 为了练习,我决定尝试制作一个简单的秒表。 从初始“时间”值 0 开始,我使用 $interval 将“时间”每 10 毫秒递增 0.01。我可以毫无
我正在尝试将间隔*[a,b]*除以保存在浮点变量数组中的*npt*点。 我需要 *a* 和 *b* 始终出现在最终数组中,并且 *npt* 可以根据我的需要进行变化。 我尝试过: delta = (b
我是一名优秀的程序员,十分优秀!