gerry 发表于 2007-11-16 18:12

MCMC的等值线该如何画

markov chain monte carlo 得出的粒子分布如下，该如何绘制等值线呢？


是划分网格，计算各处的粒子分布密度，还是从已知的各点数值拟合曲面？

没有找到书面资料，从另一个程序中看到人家是用插值实现
但我的链并不是顺序排列， 无法直接使用interp 而且粒子数较多 （1e+5） 不知该如何处理

donkeyxu 发表于 2007-11-16 19:10

回复 #1 gerry 的帖子

有意思，模拟的是一个星系恒星的分布？
“是划分网格，计算各处的粒子分布密度”就这方法！

gerry 发表于 2007-11-16 22:34

是MCMC的链在参数空间的投影，不对应实体，
划分网格的话，其实是二维hist，但我已有图中每点对应的高度值，直接插分曲面应该是更快的方法吧

donkeyxu 发表于 2007-11-16 22:42

回复 #3 gerry 的帖子

有高度，没有密度啊，不划分网格，怎么求密度？

gerry 发表于 2007-11-16 23:18

要密度干什么，直接生成三位曲面不就能做contour了么？

donkeyxu 发表于 2007-11-17 11:43

回复 #5 gerry 的帖子

你的坐标值不是密度值，你要等密度线就要求出各点密度。
页: [1]
查看完整版本: MCMC的等值线该如何画