24小时热门版块排行榜    

查看: 4662  |  回复: 12

yoyo0201

金虫 (著名写手)

[求助] 关于概率密度函数为负数问题 已有3人参与

概率密度函数为非负,但是用ksdensity 函数得到的概率密度函数有负数,现在我想将负数去掉请问可以吗?如何操作?求指点
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nonsmooth

银虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
yoyo0201: 金币+5, ★★★★★最佳答案 2015-08-02 14:38:20
【用ksdensity 函数得到的概率密度函数有负数】说明这里的处理有问题。
学术无国界。
2楼2015-08-02 13:34:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
yoyo0201: 金币+5, ★★★很有帮助 2015-08-02 14:42:35
应该是处理有问题吧,PDF应该非负的。
知其然,知其所以然。
3楼2015-08-02 14:02:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yoyo0201

金虫 (著名写手)

引用回帖:
2楼: Originally posted by Nonsmooth at 2015-08-02 13:34:59
【用ksdensity 函数得到的概率密度函数有负数】说明这里的处理有问题。

恩,matlab用Ksdensity函数直接得到的概率密度函数就有负数,网上查了下说里面的核函数对数据样本进行了平滑处理,核函数的参数带宽会影响平滑的程度,没有查到如何处理。。求助一下
4楼2015-08-02 14:42:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yoyo0201

金虫 (著名写手)

引用回帖:
3楼: Originally posted by mygt_hit at 2015-08-02 14:02:31
应该是处理有问题吧,PDF应该非负的。

那请问下,哪里处理有问题呢?我是看网上说用这个函数直接得到概率密度函数,但是得到非负不符合我的实际,求助如何解决
5楼2015-08-02 14:43:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by yoyo0201 at 2015-08-02 14:43:30
那请问下,哪里处理有问题呢?我是看网上说用这个函数直接得到概率密度函数,但是得到非负不符合我的实际,求助如何解决...

matlab默认的核函数是高斯核,本身是非负的。通过不同高斯核的线性组合估计PDF,如果出现负值,应该是线性组合系数中出现负数了。你试试其他的带宽参数或其他核函数。PDF估计也有其他方法,也可以试试。
知其然,知其所以然。
6楼2015-08-02 14:54:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yoyo0201

金虫 (著名写手)

引用回帖:
6楼: Originally posted by mygt_hit at 2015-08-02 14:54:01
matlab默认的核函数是高斯核,本身是非负的。通过不同高斯核的线性组合估计PDF,如果出现负值,应该是线性组合系数中出现负数了。你试试其他的带宽参数或其他核函数。PDF估计也有其他方法,也可以试试。...

其他的带宽参数或者核函数不知道如何设置?试了直方图画概率密度函数图,貌似准确性没有这个好。。。请求帮助~~~急待解决的问题
7楼2015-08-02 14:58:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by yoyo0201 at 2015-08-02 14:58:05
其他的带宽参数或者核函数不知道如何设置?试了直方图画概率密度函数图,貌似准确性没有这个好。。。请求帮助~~~急待解决的问题...

你看看Ksdensity的帮助,通过属性名-属性值可以设置核函数,也可以设置其他参数。

例如:通过核函数估计CDF的。
x = [randn(30,1); 5+randn(30,1)];
xi = linspace(-10,15,201);
f = ksdensity(x,xi,'function','cdf');
plot(xi,f);

通过kernel属性设置核函数。
f = ksdensity(x,xi,'kernel','triangle');
通过width属性设置带宽。
f = ksdensity(x,xi,'width',1);
知其然,知其所以然。
8楼2015-08-02 15:12:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yoyo0201

金虫 (著名写手)

引用回帖:
8楼: Originally posted by mygt_hit at 2015-08-02 15:12:34
你看看Ksdensity的帮助,通过属性名-属性值可以设置核函数,也可以设置其他参数。

例如:通过核函数估计CDF的。
x = ;
xi = linspace(-10,15,201);
f = ksdensity(x,xi,'function','cdf');
plot(xi,f);

...

好的,我看看。非常感谢啊~~有问题再向您请教
9楼2015-08-02 15:46:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

43fd6ys

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
应该全正,计算为负应检查算法的正确性.
10楼2015-08-02 20:58:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yoyo0201 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 4/200 2026-03-22 07:20 by ilovexiaobin
[考研] 0703化学调剂 +11 妮妮ninicgb 2026-03-15 15/750 2026-03-21 19:15 by ColorlessPI
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见