b*******g 发帖数: 513 | 1 问题是找一个有两个mode的distribution的highest density region. 有个人写了个文
章和一个R package,链接是:http://robjhyndman.com/papers/hdrgraph.htm。要用的package is "hdrcde", subroutine is hdr.den for one dimension data. 用他的example data which is faithful$eruption,有结果.用自己的data, 总有一个error:Error in bw.SJ(x) : no solution in the specified range of bandwidths。在subroutine hdr.den 设定了bandwidth, "h=0.09", 得到一个还行的结果。而且用不同的bandwidth,得到的highest density region 不同。主要是自己对于一个mode的distribution,写了一个程序,有了一个不错的结果。用自己有一个mode的数据,用hdr.den 分析有同样err | t*****9 发帖数: 19 | |
|