R median和ecdf()函数给出了不同的结果-为什么？
R median和ecdf()函数给出了不同的结果的原因是因为它们计算的是不同的统计量。
首先，R median函数用于计算数据集的中位数。中位数是将数据集按大小顺序排列后，处于中间位置的数值。如果数据集中的观测值有偶数个，则中位数是中间两个数的平均值。中位数是衡量数据集的集中趋势的一种方式。
而ecdf()函数用于计算数据集的经验累积分布函数（ECDF）。ECDF函数给出了对于给定的值，数据集中小于等于该值的观测值的比例。它表示了数据集中每个值的累积频率分布情况。
因此，R median函数和ecdf()函数计算的是不同的统计量，所以它们会给出不同的结果。median函数关注的是数据的集中趋势，而ecdf函数关注的是数据的分布情况。
无论是在什么场景下，不同的统计量都有其独特的用途和应用。如果您需要衡量数据的中心趋势，可以使用median函数；如果您需要了解数据的分布情况，可以使用ecdf()函数。
对于腾讯云的相关产品和产品介绍，由于要求不能提及特定的云计算品牌商，我无法提供具体的产品和链接。但是，腾讯云作为一家知名的云计算提供商，提供了各种云计算服务和解决方案，涵盖了云服务器、云数据库、人工智能、物联网等领域。您可以访问腾讯云官方网站，查看他们的产品和解决方案以获得更多相关信息。
相关·内容
我有一个有80个值的向量vec，如果我应用median(vec)函数，我会得到一个值。然而，我想做的是相反的，给出一个数字估计它所属的百分位数。我已经找到了ecdf()函数，但是我得到了不同的结果。这是一个简化的例子 > vec = c(100,150,150,150,150,150,200)
> median(vec)
# This gives the expected result
[1] 150
# However if I go the other way around, meaning I pass the value and try to return the pe
浏览 10提问于2020-04-23得票数 1
我试图在这里使用这个解决方案。
dat[, .(quant = quants, val = ecdf(dist)(quantile(dist, quants))), by = rowval]
这给出了每个行的不同结果。
然而，当我将它应用于我的数据时，总是得到相同的输出。
示例：
ecdf(gg)(quantile(gg, quants))
ecdf(ff)(quantile(ff, quants))
为什么我对ff和gg都有同样的看法？
浏览 3提问于2022-07-06得票数 -1
回答已采纳
1回答
我用ecdf()函数计算了一些随机样本的经验累积密度函数：
set.seed(0)
X = rnorm(100)
P = ecdf(X)
现在，P给出了ECDF，我们可以绘制它：
plot(P)
abline(h = 0.6, lty = 3)
我的问题是：我如何才能找到样本值x，例如P(x) = 0.6，即ECDF的0.6分位数，或者ECDF和h = 0.6的交点的x坐标。
浏览 0提问于2016-07-23得票数 6
回答已采纳
2回答
我一直在研究R中的MatchIt()包，并想知道如何计算这个包中的eCDF平均值。我已经使用了这个包中的数据，并运行了matchit包
library("MatchIt")
data("lalonde")
m.out1 <- matchit(treat ~ age + educ + race + married +
nodegree + re74 + re75, data = lalonde,
method = "nearest", distance = "g
浏览 4提问于2021-11-06得票数 1
回答已采纳
1回答
这里对问题的回答提出了一种通过应用ecdf的方法。
但是，我使用的是Hmisc包中的Ecdf()，因为它提供了一种方便的方法来绘制互补累积分布函数(Ccdf)图。(通过将what选项设置为'1-F')
默认情况下，Ecdf()执行绘图并返回包含x和y的嵌套列表。
如何提取某个x值的y值？然后把它画在原图上？
仅供参考：
> str(Ecdf(rnorm(20), lwd = 2))
List of 2
$ x: num [1:21] -1.46 -1.46 -1.18 -1.17 -1.16 ...
$ y: num [1:21] 0 0.0
浏览 6提问于2013-05-29得票数 2
回答已采纳
我正在尝试绘制数据集的这种。我尝试过使用Ecdf()
library(Hmisc)
dd<-read.table('critical.1.dat',head=F)
cdf<-Ecdf(dd$V1)
df<-hist(dd$V1)
ll<-((1-cdf$y[df$mids])/(df$density))
plot(df$mids,ll)
似乎一切都很顺利。在这个之后，我尝试了使用ecdf()
dd<-read.table('critical.1.dat',head=F)
cdf<-ecdf(dd$V1)
df<-hist
浏览 3提问于2012-06-21得票数 0
回答已采纳
我是R的新手，所以请耐心地回答我这个基本的问题。我有一个使用data.table包创建的数据集，数据。我在0到1之间创建了200个随机数，然后做了10000次，最后为每次迭代创建了一个带有描述性统计信息的数据表。我的代码是这样的：
rndm<-runif(200, min=0, max=1)
reps <- data.table(x=runif(200*10000),iter=rep(1:200,each=10000))
DATA <- reps[,list(mean=mean(rndm),median=median(rndm),sd=sd(rndm),min=min(rnd
浏览 6提问于2013-05-02得票数 4
回答已采纳
1回答
我试图从数据框架的子集中重新排序一个因素，该子集由另一个使用forcats::fct_reorder()定义的因素定义。
考虑以下数据框架df
set.seed(12)
df <- data.frame(fct1 = as.factor(rep(c("A", "B", 'C'), each = 200)),
fct2 = as.factor(rep(c("j", "k"), each = 100)),
val = c(rnorm(100, 2), # A
浏览 3提问于2020-06-10得票数 3
我用R生成了下面的图：
代码：
ecdf.bdel<-ecdf(bDE[,6]) #### cumulative distribution of a vector of numeric values
curve(1-ecdf.bdel(x),col="red",xlim=r,xlab="r2",ylab="Fraction of SVs tagged") ###plotting inverse ecdf curve
剧情和预期的一样。但是，它会在轴标记和轴标签之间以及轴标签和轴实验室之间留下巨大的空白。
有没有人能提供一些减少
浏览 10提问于2016-08-14得票数 4
回答已采纳
2回答
我想要检查一个概率是否来自经验民防框架所指定的分布。kstest给出了我认为是错误的p值，有什么问题吗？
我编写了一个测试函数来验证p值.我正在比较两个相同分布的样本数组，并检查从kstest和ks_2samp函数中得到的p值。由于零假设是正确的(分布是相同的)，那么p值必须在0,1上均匀分布，换句话说，我必须看到错误发现率等于使用的p值阈值。然而，这种情况只适用于ks_2samp函数给出的p值.
from scipy import stats
from statsmodels.distributions.empirical_distribution import ECDF
def test
浏览 6提问于2019-06-19得票数 0
回答已采纳
1回答
考虑以下序列和累积图：
x=c(0,0,0,0.5,10,1500)
qplot(x,geom='step',stat='ecdf')
这将生成一个从零的左侧开始的图形。换句话说，它看起来好像x有负值。在右侧，它在达到1500和100%后继续向右移动。所以，看起来好像有x个点大于1500。
当我手动完成整个过程时，我得到了我所期望的：
xs=sort(x)
qplot(xs,1:length(xs)/length(xs),x,geom='step')
但这似乎违背了stat='ecdf'快捷方式的全部目的。我遗漏了什么？
浏览 6提问于2013-07-17得票数 2
1回答
我正在寻找一个内置函数，返回样本分位数和估计的置信区间，而不是MATLAB (MATLAB的ecdf这样做)。
我猜R有这个内置的，只是我还没找到。
如果您有任何独立的代码可以这样做，您也可以在这里指出，虽然我希望找到一些内容，作为一个更大的开放代码库的一部分。
-Trying离开MATLAB。
浏览 1提问于2010-08-09得票数 2
回答已采纳
在这个小例子中，两个"res“变量给出了不同的结果。有人能解释一下为什么吗？我预计他们都会返回大约5。
import numpy as np
import matplotlib.pyplot as plt
dist1 = np.random.normal(100., 10., 10000)
dist2 = np.random.normal(0.05, 0.005, 10000)
res1 = dist1
res1 *= dist2
res2 = dist1 * dist2
print np.median(res1)
print np.median(res2)
# 4.9868
浏览 0提问于2019-01-31得票数 0
回答已采纳
6回答
R: ecdf优于直方图
、、、
在R中，使用ecdf我可以绘制经验累积分布函数
plot(ecdf(mydata))
使用hist，我可以绘制数据的直方图
hist(mydata)
如何在同一图中绘制直方图和ecdf？
编辑
我试着做这样的东西
浏览 1提问于2015-03-27得票数 4
1回答
我试图在数据帧中获取每一列c的ecdf()，然后将列的行r输入到该列的ecdf(c)中，以返回每个单元格对应的值ecdf(c)(r)。下面的功能起作用。
getecdf <- function(data){
# initialise result data
result <- data.frame(matrix(NA, nrow = nrow(data), ncol = ncol(data)));
rownames(result) <- rownames(data); colnames(result) <- colnames(data);
for (c
浏览 1提问于2013-11-28得票数 2
回答已采纳
1回答
下面您将找到我的python代码，用于几周前给我的类赋值，但我无法成功地调试。问题是要用FFT来找出总损失随机变量的风险值(即p%分位数)。给出了一个明确的数学过程，通过它我们可以得到一个关于总损失随机变量的离散性函数的估计。然而，我的结果是严重的，我正在犯某种错误，即使在调试我的代码几个小时之后我也找不到。
给出了总损失随机变量S，使得S=sum(X_i for i in range(N))，其中N是负二次分布的r=5, beta=.2，X_i是指数分布的theta=1。这种参数化的概率生成函数是P(z)=[1-\beta(z-1)]^{-r}。
我们被要求近似于S的分布
discretiz
浏览 5提问于2022-03-30得票数 3
回答已采纳
我想使用dplyr的case_when函数来动态地在函数之间进行选择。
例如，假设我想选择一个平均值函数，然后调用它：
average = function(.data, .avg){
f = case_when(
.avg == 'mean' ~ mean,
.avg == 'median' ~ median,
)
f(.data)
}
现在，如果我调用这个函数，就会得到以下错误：
> average(c(1, 2), 'mean')
Error in value[[1]][rep(
浏览 3提问于2021-03-25得票数 1
回答已采纳
3回答
我正在尝试计算一个数字样本("S")的平均偏差。使用"mad()“函数和一步一步地计算平均偏差时，我得到的结果是不同的。为什么？
s<- c(100,110,114,121,130,130,160)
使用"mad()“函数，我得到：
> mad(s)
[1] 13.3434
当分解公式并一步一步地执行相同的操作时，我得到：
> sum(abs(s-mean(s)))/length(s)
[1] 14.08163
为什么这些结果会有所不同？
我在输入公式时出错了吗？(这并不奇怪--我才刚刚开始学习R)。我的公式有什么问题？
或者R用于计算平均偏
浏览 2提问于2015-06-28得票数 7
我正在尝试手动计算两个随机样本的KS统计量。据我所知，KS统计量D是两个CDF之间的最大垂直偏差。但是，手动计算两个CDF之间的差异并从基数R运行ks.test会产生不同的结果。我想知道错误在哪里。 set.seed(123)
a <- rnorm(10000)
b <- rnorm(10000)
### Manual calculation
# function for calculating manually the ecdf
decdf <- function(x, baseline, treatment) ecdf(baseline)(x) - ecdf(tr
浏览 15提问于2019-04-12得票数 1
回答已采纳
我用R来计算一些数据的ecdf。我想在另一个软件中使用结果。我使用R只是为了做‘工作’，而不是为我的论文制作最终的图表。
示例代码
# Plotting the a built in sampla data
plot(cars$speed)
# Assingning the data to a new variable name
myData = cars$speed
# Calculating the edcf
myResult = ecdf(myData)
myResult
# Plotting the ecdf
plot(myResult)
输出
> # Plotting the
浏览 0提问于2014-07-27得票数 1
回答已采纳
2回答
我有一个数值变量，范围从1 (min)到5(最大值)。该值范围跨越8个不同的变量。因此，第一行如下所示：
Var1 Var2 Var3 Var4 Var5 Var6 Var7 Var8
4 4 1 4 5 4 4 1
我计算了这8个变量中每一行的中值。有时，中位数将是一个中点值，例如4.5 (因为它是偶数变量)。因此，结果行可能如下所示：
Var1 Var2 Var3 Var4 Var5 Var6 Var7 Var8 Median
1 2 3 4 5 5 5 5 4.5
当我在Medina变量中计算的每个
浏览 0提问于2018-04-05得票数 0
回答已采纳
我试着用蒙特卡洛的方法。在R中，为了找到这个问题的一个解:我有一个x~ U(0,1和Y=log(X) )。我想要的是一个估计的pdf和民防部队。问题是，我的目标是在不使用ECDF的情况下获得对CDF的估计。那么，有什么方法可以接近我的民防没有这个科曼德？理论上，我可以集成我的pdf，但我不知道它的确切形状。
为了获得这两个代码，我创建了以下R代码：
X = runif(1000) # a= 0 and b=1 default
sample = log(X)
hist(sample, xlim=c(-6,0), main="Estimated vs true
浏览 8提问于2022-12-02得票数 0
在将median()函数应用于数据帧时，我注意到一些不一致的行为。“不一致的行为”通常意味着我不理解某事，所以我希望有人愿意为我澄清这一点。
我意识到一些函数(例如，min()、max())将数据帧转换为向量，并返回整个df的相应值，而mean()和sd()则为每列返回一个值。虽然有点令人困惑，但这些行为上的差异不会造成太多问题，因为如果返回标量而不是向量，大多数代码都会崩溃。然而，median()似乎并不一致。例如：
dat <- data.frame(x=1:100, y=2:101)
median(dat)
返回一个向量：[1] 50.5 51.5
但是，有时它会崩溃：
dat2
浏览 0提问于2011-05-06得票数 11
回答已采纳
1回答
我在R中使用expss软件包进行分析，我用cro_mean_sd_n()给出了平均值和标准差。我可以调用另一个函数来查找模式吗？
和其他可能性意味着，中位数，模式，和，N在一个表中。
浏览 0提问于2018-12-06得票数 0
回答已采纳
1回答
第一次在向量上使用mad()中位数绝对偏差，大部分是4和5，还有一些NA。我是不是错了，但是为什么mad函数不会为
mad(c(4,5,5,5))
？相反，我得到了一个0。只有当向量中的所有值都相同时，我才会期望0，如果值有变化，则不会。
请给我开导一下。我使用的是R版本4.05。
浏览 10提问于2021-05-22得票数 3
我有数据和一个有效的汇总函数：
[data:] x = runif(100)
[function 1:] sapply(list("mean", "median", "sd"),
function(f) array(do.call(f, list(x)),
dimnames=list(f)))
mean median sd
0.5198142 0.5566567 0.3048372
以下代码也会运行，但不会给出dimname
浏览 5提问于2013-05-05得票数 3
回答已采纳
4回答
在Python中R的ecdf(x)(x)函数在numpy或scipy中的等效值是什么？ecdf(x)(x)基本上与以下内容相同吗？
import numpy as np
def ecdf(x):
# normalize X to sum to 1
x = x / np.sum(x)
return np.cumsum(x)
还是还需要别的东西？
编辑如何控制ecdf使用的回收箱数量
浏览 0提问于2013-04-03得票数 24
回答已采纳
1回答
我不太确定CDF (累积分布函数)和ECDF (经验累积分布函数)之间的区别，但我通常使用CDF图来观察我的数据。
我最近一直在使用R，并且拼命地试图找出如何绘制我的数据的CDF和CCDF (互补CDF)。我所能找到的就是R有ecdf，但是我不确定这是否是我要找的。绘制ECDF非常简单，如下所示：
plot.ecdf(data)
有人知道如何使用R绘制数据集的CDF和CCDF吗？
浏览 0提问于2011-06-14得票数 6
回答已采纳
我和data.table一起工作。从data.table对象中选择列很容易：
> head(data.table(mtcars)[,2])
cyl
1: 6
2: 6
3: 4
4: 6
5: 8
6: 6
但是，尝试在ecdf调用中使用此语法选择列会产生一个错误：
> ecdf(data.table(mtcars)[,2])(data.table(mtcars)[,2])
[.data.frame(x，i)中的错误:所选未定义列
有人能解释一下为什么是吗？
实际上，有一种方法是这样做：
> ecdf(data.table(mtcars)
浏览 3提问于2017-09-04得票数 2
回答已采纳
我想添加垂直段到脊线图，其直方图显示定制的分位数。
如果我用..x..映射填充颜色，我设法得到垂直段。但是我想在密度图中显示分位数。我编写了以下代码：
library(datasets)
library(ggplot2)
data("iris")
iris_lines <- data.frame(Species = c("setosa", "versicolor", "virginica"),
x0 = c(5, 5.9, 6.5))
Figure1 <- ggpl
浏览 1提问于2019-06-18得票数 1
回答已采纳
1回答
我想找到每一列的中位数，但是它不像我想要的那样工作。
1 2 3
3 2 1
2 1 5
我期待着
2 2 3
然而，对于结果，它只是给出了和错误和列的一些“和”。下面是“中间列”的代码片段。
while read -r line; do
read -a array <<< "$line"
for i in "${!array[@]}"
do
column[${i}]=${array[$i]}
((length[${i}]++))
result=${column[*]} | sort
浏览 0提问于2015-10-13得票数 1
回答已采纳
我希望替换一些使用dplyr::do的R代码，因为这个函数很快就会被废弃。我的很多工作都需要创建分层的CDF地块。在使用dply时:是否将I分层上的变量作为变量传递给结果数据框架，然后我可以轻松地使用该变量进行绘图。
我有一个替代dplyr::do使用dplyr::group_split和purrr::map_df的解决方案。但是，我在dplyr::group_split中传递的变量在结果数据帧中没有命名。这使得绘制分层数据变得困难。如何确保在dlyr::group_split中传递的变量在结果数据框架中被命名？
下面是创建使用dplyr：：do绘制数据所需的一些代码：
library(dpl
浏览 0提问于2019-07-18得票数 3
回答已采纳
2回答
这是我第一次在Stack Overflow上找到一个R问题--请原谅我没有找到任何东西的原因是因为我正在寻找的东西类型是我所不知道的(是吗？)
我想以累积频率显示数据。由于我的焦点更多地集中在分布的边缘，所以将y轴缩放到正态分布是很有帮助的。结果应该是这样的：
我读过分位数图，但老实说，如果我想保留X轴的话，我不知道如何应用它们。
我尝试了base图形和ggplot2，但都搞不清楚。因此，我目前的解决方案是，例如
plot(ecdf(trees$Volume))
或
ggplot(data=trees, aes(Volume)) + stat_ecdf()
浏览 7提问于2016-04-12得票数 1
回答已采纳
我正在尝试使用putplot来改变ggpairs图的对角线上的图。我可以一个一个地做这个，但是当我使用一个循环时，对角线元素是相同的！
library(ggplot2)
library(GGally)
p=ggpairs(iris,
columns=, c("Sepal.Length", "Sepal.Width", "Petal.Length", "Petal.Width"),
colour='Species',
lower=list(continuous='points
浏览 1提问于2014-03-07得票数 5
回答已采纳
2回答
我正在尝试使用另一个调色板，用于R.中的ggpairs库中的GGally库中的散点图(参见 )。
library(ggplot2)
library(GGally)
Works
ggplot(iris, aes(x=Sepal.Width, colour=Species)) + stat_ecdf() + scale_color_brewer(palette="Spectral")
还可以工作
ggplot <- function(...) ggplot2::ggplot(...) + scale_color_brewer(palette="Spectra
浏览 2提问于2014-03-06得票数 7
2回答
我喜欢plyr将一个数据帧分割成多个数据集，然后在每个数据集上执行相同的操作的能力。最好的部分是，它以整洁、紧凑、标签良好的表格的形式向您显示结果。我喜欢使用each()将一堆计算放到一行中。然而，我不明白为什么在ddply参数中使用摘要函数会破坏输出，并使其变得冗长且无标签。看看这里，看看我的意思。你能告诉我我哪里做错了吗？我更喜欢使用摘要。
让我们首先设置一个示例数据帧。想象一下，你有60名参与者参与了一项研究。其中20个是有趣的，20个是聪明的，20个是友好的。然后每个受试者都会得到一个分数。
type<-rep(c("funny","clever"
浏览 2提问于2011-03-10得票数 3
回答已采纳
2回答
我正在尝试实现R函数ecdf()。
我正在考虑两种情况:一种是t一维的，另一种是以t为向量的.
#First case
my.ecdf<-function(x,t) {
indicator<-ifelse(x<=t,1,0)
out<-sum(indicator)/length(x)
out
}
#Second case
my.ecdf<-function(x,t) {
out<-length(t)
for(i in 1:length(t)) {
浏览 1提问于2014-12-07得票数 0
回答已采纳
如何返回R中相同函数内的均值、中值和标准差？我所能得到的只是计算标准差的函数的最后一部分。我在想，当我返回'b‘时，通过给b分配Tail_wags，我将拥有所有的三个值。在“b”之后为函数之外的三个值添加结果，以查看应该是什么值。
Dog_biscuits <- c(0,1,2,3,4,5,6,7,8,9,10)
Tail_wags <- c(0,0,1,3,8,13,14,12,15,16,14)
dog_wags<-cbind(Dog_biscuits,Tail_wags)
dog_wags
summarystat<- function(x) {
z
浏览 4提问于2019-11-03得票数 3
回答已采纳
2回答
我想在R中绘制1-CDF。
我使用的是ggplot stat_ecdf。
g1=ggplot() +
stat_ecdf(data=data_ploting, aes(x, colour=ggg), alpha=0.8, geom= "smooth", pad = FALSE) +
theme_test ()
浏览 15提问于2022-04-12得票数 0
我想使用list在R中创建一个未计算函数的列表。我想要一个可以由以下代码生成的"myList“：
xlist = c("A", "B", "C", ..., "Z")
myList = alist(print(xlist[1]), print(xlist[2]), print(xlist[3]), ..., print(xlist[26]))
然而，只有当xlist很短时，上述方法才是可行的。如何使用一些聪明的函数生成myList？我已经试过了。
tempfun = function(x) alist(print(x))
浏览 0提问于2016-05-03得票数 1
你好，我正在构建一个R函数，并想了解为什么其中一个工作，另一个不工作。
这个不起作用
isola_por_fator_em_col <- function(data,col,fator)
{
y <- data[which(data$col==fator),]
x <- select_if(y,is.numeric)
summary(x)
}
isola_por_fator_em_col(data=desempenho_aluno_escola,col=priv,fator="privada")
Warning messa
浏览 1提问于2022-06-02得票数 0
我写信是想问你对如何解释这个案例有什么看法。我有两个向量"a“和"b”，我试图比较它们。
wilcoxon检验给出了a/b的p值为5.139217e-303，另一个选择是“更大”。现在，如果我对它们中的每一个做一个总结，我会得到以下内容
> summary(a)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.0000000 0.0001411 0.0002381 0.0002671 0.0003623 0.0012910
> summary(c)
Min. 1st Qu.
浏览 1提问于2012-09-20得票数 0
回答已采纳
1回答
我正在尝试构建一个R包，其中包含几个ecdf对象。它们是由ecdf(一些变量，如p1)创建的。我将它们放在一个列表中，并将其保存为数据文件夹中的rda文件，但是当我在包内运行函数(假设命名为b1)时，在安装包后，我得到了以下错误：
b1(zzz[1,]) (zzz is a data frame and I ran the function on one row of it)
Error in fc(p1) : could not find function ".approxfun"
fc是存储在保存列表中的ecdf函数，我在函数中调用list by data(list1)，
浏览 1提问于2013-09-13得票数 1
回答已采纳
1回答
我正在将分布拟合到给定的data.then，我已经估计了分布的参数，我还使用R中的ecdf()命令绘制了经验cdf。我该怎么做呢？ecdf()命令还能帮助我吗？
浏览 3提问于2013-10-15得票数 2
1回答
有可能区分ECDF吗？例如，以下面获得的结果为例。
set.seed(1)
a <- sort(rnorm(100))
b <- ecdf(a)
plot(b)
为了得到它的概率密度函数，我想取b的导数。
浏览 0提问于2018-07-20得票数 3
1回答
叠加箱形图和线状图出海
、、、、
我正在尝试在同一个框图上叠加一个箱形图(基于另一个变量的一系列箱形图)和该变量的中线的线形图。像下面这样的简单代码可以很好地工作。
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
dfx=pd.DataFrame({'S':np.random.randint(10,100,9)*10,'C':
['X','X','X','Y','Y','
浏览 3提问于2019-07-13得票数 1
我们可以使用以下命令创建ECDF
import numpy as np
from statsmodels.distributions.empirical_distribution import ECDF
ecdf = ECDF([3, 3, 1, 4])
，然后使用以下命令获取ECDF
ecdf(x)
但是，如果我想知道百分位数97.5%的x，该怎么办？
从http://www.statsmodels.org/stable/generated/statsmodels.distributions.empirical_distribution.ECDF.html?highlight=ecdf上看，
浏览 0提问于2017-05-23得票数 7
这个问题是来自交叉验证的migrated，因为它可以在堆栈溢出时得到回答。Migrated 15天前 我正在写一个使用加权非线性最小二乘(WNLS)方法估计GPD参数的代码。WNLS方法由两个步骤组成 步骤1：$(\hat{\xi_1}，\hat{b_1}) = arg\ \min_{(\xi，b)} \sum_{i=1}^{n} \log(1-F_n(x_i)) - log(1-G_{\xi，b}(x_i))$，这里$F_n$是ECDF，$1-G_{\xi，b}$是广义帕累托分布。 有人能告诉我如何计算R中数据"X“的EDF函数$F_n$吗？ecdf(X)(X)会计算ECDF吗？如
浏览 29提问于2021-11-15得票数 1
回答已采纳
1回答
我正在尝试使用ecdf，但我不确定我是否做得对。我的最终目的是找出什么分位数对应于一个特定的值。例如：
sample_set <- c(20, 40, 60, 80, 100)
# Now I want to get the 0.75 quantile:
quantile(x = sample_set, probs = 0.75)
#result:
75%
80
# Let's use ecdf
ecdf(x = sample_set) (80)
#result
0.8
为什么会有这种差异？我是做了一些微不足道的错误，还是取决于分位数的计算方式？
谢谢，麦克斯
浏览 2提问于2016-03-10得票数 6
回答已采纳
1回答
用条形图来绘制R中的ecdf函数，而不是用线段来表示ecdf函数是否可行？
矿石是否有另一种方法可以用y轴上的累积密度而不是频率来绘制ggplot中的累积直方图？
浏览 6提问于2017-03-29得票数 0
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
