已知一个密度函数，怎么用R求四分位点

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

已知一个密度函数，不属于任何特定的分布，请问怎么用R求四分位点？
像正态分布可以用qnorm函数来求，可是现在只有一个密度函数，怎么求？
请论坛里的各位大神帮忙~~

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：密度函数 正态分布 Norm ORM NOR 正态分布

有没有人啊

提供一个简单粗暴的办法，根据定义算，给个简单的例子

1. Density<-function(x){
   
2.    (1/sqrt(2*pi))*exp(-x^2/2);
   
3. }
   
4. quan<-function(y){
   
5.     t<-integrate(Density,-Inf,y)$value-1/3;
   
6. }
   
7. uniroot(quan,c(-1,0));
   
8. qnorm(1/3);

复制代码

如果lz对精度不满意的就再调一下参数好咯

用Monte Carlo取样

1. h<-function(alpha,Theta,a,b){
   
2. f<-function(x,Theta)
   
3.   cdf(x,Theta)-alpha
   
4. uniroot(f,c(a,b),n,m)->r;r$root}

复制代码

α是下分位点，θ是参数向量，cdf()是某一确定分布的分布函数，a、b是迭代区间；
从定义出发，与3L的数值积分不同。h即为所求α下分位数

想想分位点的定义，编程然后对应找到分位数

求证1加1 发表于 2014-1-31 21:01
提供一个简单粗暴的办法，根据定义算，给个简单的例子如果lz对精度不满意的就再调一下参数好咯

蟹蟹！