b********1
发帖数: 2861 | 1 有没有人用R 或者 Matlab 做过RJMCMC model? 我知道Winbugs有函数可以做normal
distribution的RJMCMC， 但是我现在做的Poisson 和NB model 不能用.求哪位大侠指
点在选择proposal方面的经验。多谢！ |
w*****n
发帖数: 375 | 2 还真有人用Rerversible Jumping MCMC. 能好奇问一下你是在研究什么问题需要RJMCMC
? |
g*******i
发帖数: 258 | 3 http://lemur.mcs.st-and.ac.uk/Book-website/R-codes/caprecapRJ.R
which is R code for RJMCMC for capture-recapture models
【在 b********1 的大作中提到】
: 有没有人用R 或者 Matlab 做过RJMCMC model? 我知道Winbugs有函数可以做normal
: distribution的RJMCMC， 但是我现在做的Poisson 和NB model 不能用.求哪位大侠指
: 点在选择proposal方面的经验。多谢！
|
g*******i
发帖数: 258 | 4 RJMCMC is widely used in time-varying parameter context, for example:
factor models (for the number of factors)
mixture models (for the number of mixtures)
spline regressions (for the number of node points)
However, I do have concern about its usability in practice. Some
applications I read are really trashy
RJMCMC
【在 w*****n 的大作中提到】
: 还真有人用Rerversible Jumping MCMC. 能好奇问一下你是在研究什么问题需要RJMCMC
: ?
|
b********1
发帖数: 2861 | 5 我想做mixture model+variable selection, 就是同时允许模型空间和变量空间变化。
请问你为什么如此惊讶呢？是模型太复杂还是结果不可靠？多谢！
RJMCMC
【在 w*****n 的大作中提到】
: 还真有人用Rerversible Jumping MCMC. 能好奇问一下你是在研究什么问题需要RJMCMC
: ?
|
b********1
发帖数: 2861 | 6 多谢了!我研究一下这个code.
【在 g*******i 的大作中提到】
: http://lemur.mcs.st-and.ac.uk/Book-website/R-codes/caprecapRJ.R
: which is R code for RJMCMC for capture-recapture models
|
b********1
发帖数: 2861 | 7 我觉得单纯的mixture model 或者variable selection都可以通过稍微简单的方式解决
，但是RJMCMC看起来更fancy一些，而且我想实现模型空间和变量空间都变化，暂时只
是YY啦。
我是做土木工程的，在统计方面做的fancy一点，写文章比较容易忽悠，跟专业的统计
系或者经济系的人思维模式不一样。
【在 g*******i 的大作中提到】
: RJMCMC is widely used in time-varying parameter context, for example:
: factor models (for the number of factors)
: mixture models (for the number of mixtures)
: spline regressions (for the number of node points)
: However, I do have concern about its usability in practice. Some
: applications I read are really trashy
:
: RJMCMC
|
k*****u
发帖数: 1688 | 8 支持一下 我做research的时候做过这个，用来解决mixture model，还有更fancy的，
比如birth-death process， dirichlet process用来解决mixture model
但是工作以后发现什么都用不着。只要搞各种regression就好了。
【在 b********1 的大作中提到】
: 我觉得单纯的mixture model 或者variable selection都可以通过稍微简单的方式解决
: ，但是RJMCMC看起来更fancy一些，而且我想实现模型空间和变量空间都变化，暂时只
: 是YY啦。
: 我是做土木工程的，在统计方面做的fancy一点，写文章比较容易忽悠，跟专业的统计
: 系或者经济系的人思维模式不一样。
|
k*****u
发帖数: 1688 | 9 ps 其实em算法做mixture model已经够好了
取不同的初始值然后repeat个1000次最后用平均做参数估计，好像效果已经很好了。
【在 b********1 的大作中提到】
: 我觉得单纯的mixture model 或者variable selection都可以通过稍微简单的方式解决
: ，但是RJMCMC看起来更fancy一些，而且我想实现模型空间和变量空间都变化，暂时只
: 是YY啦。
: 我是做土木工程的，在统计方面做的fancy一点，写文章比较容易忽悠，跟专业的统计
: 系或者经济系的人思维模式不一样。
|
a***g
发帖数: 2761 | 10 em可以用来做mm是因为连续变量起码可以写似然的式子
如果是树结构的mm，通常连似然函数都没有很好的形式，所以在用b-d process，这也
是为什么很多做基因进化的人用它来解决问题
实际上不是什么方法fancy用什么而是看手头的数据符合什么类型才用什么 |
h*******h
发帖数: 14 | 11 slow and bad mixing performance, esp when your model space is large...
【在 b********1 的大作中提到】
: 有没有人用R 或者 Matlab 做过RJMCMC model? 我知道Winbugs有函数可以做normal
: distribution的RJMCMC， 但是我现在做的Poisson 和NB model 不能用.求哪位大侠指
: 点在选择proposal方面的经验。多谢！
|
K*****2
发帖数: 9308 | 12 搭车求问，RJMCMC是不是连十几个变量的都很难做？
【在 h*******h 的大作中提到】
: slow and bad mixing performance, esp when your model space is large...
|
