quantile

0热度

1回答

我在R中使用quantreg包对一组数据运行分位数回归（95％）。我想将分位数回归的斜率设置为值1.4，就像之前的分析中所做的那样，我想比较我的结果。如果在lm()这对于函数offset()是可能的，对于固定分位数（例如0.025）使用rq()，则这不起作用。该代码不会给出错误，但1.4的值对我的结果没有影响。 fit.0.025<-rq(y~offset(1.4*x),tau=0.025,

2热度

2回答

如何在R中使用texreg进行分位数回归时设置se = boot？

我正在运行分位数回归（包quantreg）并使用texreg来创建我的模型的乳胶输出。我感兴趣的是bootstrapped s.e.并在摘要的选项中设置se =“boot”，但是当我使用texreg时，我得到“n.i.d.” S.E. 如何更改该选项？下面是我在做什么： tm3 <- rq(nback ~ cara + mat + dut + e_brown + e_green + e_bl

-2热度

2回答

由两个不同的列

我有一个数据集，看起来像这样 year month age 2007 1 17 2007 1 18 2007 1 19 2007 1 30 2007 1 31 2007 2 18 2007 2 19 2007 2 30 2008 2 41 2008 2 52 2008 2 49 2008 3 23 2008 3 19 2008 3 39

4热度

1回答

R分位数函数的行为

在处理问题时，我注意到一些有趣的事情。我不知道到底发生了什么，但发生了一些我并不期望发生的事情。这可能是我犯了一个错误，但让我用一个例子开始： x <- rnorm(100) y <- x[ x > quantile(x, 0.1) ] z <- x[ x > quantile(x, c(0.1, 0.2)) ] a <- x[ x > quantile(x, c(0.1, 0.2, 0.3

4热度

1回答

未知累积函数的反函数

我正在处理数据文件，里面的观察值是随机值。在这种情况下，我不知道x的分布（我的观察）。我使用函数密度来估计密度，因为我必须应用核函数估计。 T=density(datafile[,1],bw=sj,kernel="epanechnikov") 之后，我必须整合这一点，因为我正在寻找分位数（类似于VaR，95％）。为此，我有两个选择： ecdf() quantile() 现在我有95位数

5热度

3回答

Equavelent python在MATLAB分位数命令

我想在matlab中复制我的一些代码到python。我发现matlab中的分位数函数在Python中没有“完全”对应的一个。我发现的最接近的是python的mquantiles。例如的MATLAB： quantile([ 8.60789925e-05, 1.98989354e-05 , 1.68308882e-04, 1.69379370e-04], 0.8) 给出：0.00016958

1热度

1回答

qqline没有给出我预期的线

如果我执行我的多元数据的mahalanobis距离与卡方分布相关的多元qqplot，我预计伴随的qqline是一个截距为0和斜率为1的直线。但是，如果我运行下面的代码： scores<-matrix(rnorm(100*3),nc=3) mah_scores = mahalanobis(scores,center=colMeans(scores),cov=cov(scores)) chi_sc

5热度

1回答

绘制分位数回归线

在Stata中，可以绘制分位数回归线吗？我知道一个标准的OLS回归线可以添加到散点图中，但我不清楚如何添加其他类型的回归线。如果这是不可能的，是否可以绘制我以mx + b格式指定的一条线？然后我可以从分位数回归输出中取斜率和截距项并手动绘制它们。感谢您的帮助！

2热度

2回答

R：知道表（一个或多个），如何calc下位数（S）

s是一个大阵列，只需保存桌在数据库 > table_s s 1 2 3 4 5 3000000 1 1 999999999999999999 34 如何与R中table_s calc下分位数（S）？感谢

0热度

1回答

如何在SAS中的正态概率Q-Q图上绘制残差的95％置信区间（面积）？

在'CI区'下，我的意思是平等。像在@cbeleites' answer图形上绘制或罗恩Pearson的post描述：谢谢。