作者AmigoSafin ()
看板Statistics
标题[问题] 已有p.d.f.求shanonn entropy
时间Wed Oct 10 04:56:15 2018
大家好～
shannon entropy的公式是 对 p.d.f.以及log(p.d.f.)积分
已知p.d.f.为
(-a/b) [y/b]^{a-1} * e^{-[y/b]^a}
并且给了a=1
於是我就将a=1带入
得到 integrate (1/b) e^{-y/b} * log (1/b*e^{-y/b}) dy
整理往下
integrate (1/b) e^{-y/b} log (e^{(-y/b)(1/b)} dy （这边很不确定1/b可以直接成为
(-y/b)的指数吗？
再往下
1/b integrate e ^{-y/b} * (-y/b)^{1/b} dy
再往下就不知道该怎麽处理
想用 uv-integrate v du
但似乎无法处理乾净
找了一下不知道R有没有可以直接给p.d.f.算H[X]的
来跟大家请益求救～～
谢谢大家！！
--
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 129.21.70.218
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1539118579.A.999.html
1F：→ Pieteacher: 知道函数公式就可以直接写R呀 令function()然後用inte 10/10 10:04
2F：→ Pieteacher: grate(fun,xmin,xmax) 10/10 10:04
3F：→ Vulpix: pdf居然是负的。 10/10 11:45
4F：推 goshfju: 你的pdf有问题吧 10/10 11:50
5F：→ goshfju: 你对下你上一篇的cdf 很显然你微分微错了 10/10 11:51
6F：推 garrrick: pdf长得像是Weibull distribution的pdf，只是多了负号 10/10 20:29
7F：推 garrrick: a=1就变成exponential distribution了模样更简洁 10/10 20:31
8F：推 garrrick: 然後看要什麽衍生的资讯就google一下 10/10 20:34
9F：→ AmigoSafin: 谢谢大家～～少了一个负号..orz 10/12 08:54