15 条回复
请问楼主解决问题了没?我也碰到了同样的问题 |
这个问题谁能解释下 不胜感激:handshake :handshake |
同求呀,大师帮忙解决一下吧 |
这个应该有点帮助 |
2013-12-26 14:22 上传
点击文件名下载附件
88.84 KB, 下载次数: 19940
数理统计常用命令
同求指导……楼主后来是怎么解决这个问题的呢? |
本帖最后由 maple1314168 于 2014-5-6 19:53 编辑 可以使用插值spline等命令获取。不保证最短的区间长度。例如:x = [randn(30,1); 5+randn(30,1)];[f,xi] = ksdensity(x,'function','cdf');F=[spline(f,xi,0.05) spline(f,xi,0.995)] 可以改变两端取最短长度。 或者在命令中加入你想要的点,不插值。 最简单可用prctile命令,直接用数据估计。 |
maple1314168 发表于 2014-5-6 19:12 前辈您好,我想问问如果想要找最短区间长度的置信区间,应该怎样编程实现?先谢谢您了! |
小呆a1 发表于 2016-3-29 11:02 用fitdist结合fmincon可以顺利解决问题。 如: x = [randn(30,1); 5+randn(30,1)]; pd = fitdist(x,'Kernel'); t=fmincon(@(x)icdf(pd,x)-icdf(pd,x-0.95),0.95,[],[],[],[],[0.95],[1]) %0.95是你设定的置信概率 t= 0.9744 [icdf(pd,t-0.95), icdf(pd,t)] %最短置信区间 ans =[-3.5774 8.2325] |
maple1314168 发表于 2016-3-30 17:06 谢谢您的耐心解答!真的很感激!!我想问如果我已知一个分布的概率密度函数和累积分布函数,可以怎么做?也可以这样做吗? |
maple1314168 发表于 2016-3-30 17:06 请问 这里t=0.9744是什么值呢 如果我换了x数据 这个t要怎么修改 |
guaishou970210 发表于 2022-6-29 16:24 t是计算出来的。t=fmincon(@(x)icdf(pd,x)... 计算最短区间长度。 |
Powered by Discuz! X3.4
© 2001-2022