|
f***a
发帖数: 329 | 2 我在我的的电脑上死活用不了snow，后来用了下snowfall，感觉更user friendly，你
可以看看。再后来干脆不写parallel code了，反正我做的是cross-validation with
MCMC,所以直接把data切割好submit jobs到学校的cluster上算，呵呵~~~ good luck! |
|
l*********s
发帖数: 5409 | 3 netflix price 原始的贝叶斯随机矩阵分解解法：
http://citeseerx.ist.psu.edu/viewdoc/download?
doi=10.1.1.127.6198&rep=rep1&type=pdf
我的问题是
1)如果没有适当的限制，因子分析的解不是唯一。 这个对MCMC有什么影响？比如，多模态? 边际分布
弥散？
2)要进行多维的预测，怎么做呢？后验密度函数给出预测分布的均值？模？ 如果是模，怎么求经验多元
分布的模？ |
|
l*********s
发帖数: 5409 | 4 all right>
Below is the reference paper of the original Bayesian probability matrix
factorization method for the Netflix price data set(
http://citeseerx.ist.psu.edu/viewdoc/download?
doi=10.1.1.127.6198&rep=rep1&type=pdf )
Essentially, this is a Bayesian flavored factor analysis, with user
preference matrix as the factor, and latent movie feature matrix as the
factor loading.
My question are two fold:
1) Since the factorization is not unique due to rotation invariance, I
wonder what will be the... 阅读全帖 |
|
o**m
发帖数: 828 | 5 Markov chain Monte Carlo: stochastic simulation for Bayesian inference By
Dani Gamerman, Hedibert Freitas Lopes
如果有纸质书的，可以考虑二手购买
thx |
|
q**j
发帖数: 10612 | 6 为什么叫monte carlo呢？好像不是很make sense。 |
|
a****m
发帖数: 693 | 7 casino city name in France, Monte carlo, |
|
|
q**j
发帖数: 10612 | 9 why does monte carlo have any thing to do with random number generation? 叫
一个别的名字, 比如 markov chain inversion 不是更合理？ |
|
c*******o
发帖数: 8869 | 10 generate random number from a pre-defined probability distribution is monte
carlo simulation by definition, although in statistics I do not see any
simulation that is not monte carlo.... |
|
q**j
发帖数: 10612 | 11 really? traditional method is just called AR method.
monte |
|
g********r
发帖数: 8017 | 12 在美国的应该叫MCLV - markov chain las vegas
这是俺们老师说地. |
|
|
|
s******5
发帖数: 513 | 15 monte carlo markov chain, 第二个mc不是markov chain么.......... |
|
b*****n
发帖数: 685 | 16 是么？
MH-within-Gibbs是指整个MCMC是Gibbs，不过无法得到conditional posterior
distribution的显式表达所以用MH来获得conditional posterior sample。
Gibbs-within-MH是这意思吗？能否说明白点？ |
|
b*****n
发帖数: 685 | 17 呵呵，还是跟我说的一样嘛。
不过很多人做MCMC不严谨，有些人就凭想象乱作一通，根本不管detailed balance。这
个现象很不好。 |
|
d*******o
发帖数: 493 | 18 实在看不懂，你要用nested macro做mcmc吗？ |
|
p******0
发帖数: 107 | 19 两篇文章是关于MCMC的，虽然发在物理杂志上。在物理版问过却没人理。希望发在这里
有同学帮助。
我们这学校真是穷疯了，连这样的杂志都不定！
Physics Letters B
Volume 195, Issue 2, 1987, Pages 216-222
Physics Letters B
Volume 268, Issue 2, 1991, Pages 247-252
我万分感谢，有人帮我找到这两篇文章能帮我发到
a*****[email protected]
吗？ |
|
a*********r
发帖数: 139 | 20 Simulation by Ross has a very readable chapter that can help you understand
MCMC. |
|
w*********a
发帖数: 156 | 21 搜了一下现有的R packages，好像都是已经预先定义好的proposal density，而且多是
normal proposal，但是我想用uniform proposal density，因为我的随机变量的取值
范围
是在一个有限区间内的。
多谢帮忙啊~ |
|
|
|
w*********a
发帖数: 156 | 24 难倒是感觉不会太难，但是自己编的程序感觉在效率和robustness上肯定比已有的那些
package
（比如好像有个aMCMC，看介绍，好像效率高很多）会弱不少吧？？
我不知道自己编会不会导致运算时间太长，什么的。。。没啥经验
谢谢回复啊 |
|
r********0
发帖数: 65 | 25
自己写的大部分时间还是效率不是很高的。 diagnostics 都是用package做得。
但是自己的model大部分还是得自己写。 |
|
a****m
发帖数: 693 | 26 the first course in bayesian analysis，等书里面有相似的，改一改。 |
|
j********e
发帖数: 52 | 27 I recommend "Bayesian Computation With R" along with some lecture notes.
You can first read "Bayesian statistics: a concise introduction" (google it)
to get familiar with the basics. Then read the book above and do exercises
. In later chapters, you need to get additional notes on "MCMC", "
hierarchical modeling" and "regression models".
Gelmen's book is horrible and totally inappropriate for beginners. |
|
a***g
发帖数: 2761 | 28 比较虚词这是老梗了
像现在汉语在网上出现的情况是很多词或者词组都是同音的别字只要顺序正确就不影响
表达
虽然英语也有很多缩写和谐音，但是还算是成块的，比较起来真是减少了很多麻烦
前一阵和在多伦多的同学聊天还提起来那边有做mcmc的介入了一些这方面的工作，不过
就是瞎侃，也没见过文章
另外不知道这个现在汉语的机器翻译做到什么程度了 |
|
q*****q
发帖数: 158 | 29 倒，我没说你可以直接写成显性形式。你要是想自己算，就得用MCMC之类的方法去估概
率。当时一起实习的一个兄弟就是做的这个。当然可以用软件更好了。有一个专门的软
件，忘了名字。SAS里面好像也有相关的procedure。
话说我实习的时候也是干了三个月的group sequential design。我只是希望能够给你
些可行的建议。不要这么没耐心。飘过。。。 |
|
D**g
发帖数: 739 | 30 陈老师是会SAS的吧？PROC MCMC里有个例子，你何不用你的理论，给大家演示一下具体
算法，最后看看你的算法和结果和贝叶斯结果，frequentist算法有何区别。
Example 52.6 Change Point Models
Consider the data set from Bacon and Watts (1971), where is the logarithm of
the height of the stagnant surface layer and the covariate is the logarithm
of the flow rate of water. The following statements create the data set:
title 'Change Point Model';
data stagnant;
input y x @@;
ind = _n_;
datalines;
1.12 -1.39 1.12 -1.39 0.99 -1.08 1.03 -... 阅读全帖 |
|
m**s
发帖数: 346 | 31 修一门bayesian methods in computer science。
主要是编程实现一些统计上面的算法。
开始还好，最近的一个lab，让我很头大。也许给我很多时间，我能想出来，但是会太
迟了。明天就要交了，虽然是memorial day。
不是想找人帮忙做，只需要懂的人看看，点拨一下。为了表示对您的时间的尊重，我愿
意$10/hr 补偿您的时间。钱肯定不多，您如果嫌不够，可以再商量。
下面是这个lab的大概意思：
implement MCMC using Metropolis (or Metropolis-Hastings).lab给了4个bayesian
networks，要用我们写的程序针对某个query，做sampling，然后plot mixing, prior,
and posterior. 其实我自己感觉这个lab有些地方没说清楚，但是也许懂的人能够猜
出来。这五个networks是：
1）faculty evaluation. faculty evaluation scores are distributed normally
with a normal pr... 阅读全帖 |
|
g********r
发帖数: 8017 | 32 所以我劝楼主找统计大牛实验室做个薄厚。也就两年。然后师兄弟工业界学界都有，就
能拿到面试了。计算化学的编程背景，多半mcmc也熟，还是挺吸引人的。物理博士到生
统当教授也有挺成功的例子。 |
|
g********r
发帖数: 8017 | 33 所以我劝楼主找统计大牛实验室做个薄厚。也就两年。然后师兄弟工业界学界都有，就
能拿到面试了。计算化学的编程背景，多半mcmc也熟，还是挺吸引人的。物理博士到生
统当教授也有挺成功的例子。 |
|
f***a
发帖数: 329 | 34 不会必修c++吧。。。
我感觉统计系上至教授，下至学生一般没几个c++用得好的，大多根本不会用。。。
当然你搞bayesian，老是做复杂模型mcmc的话，c/c++还是有必要，不然R/matlab太慢
了。。。 |
|
q**j
发帖数: 10612 | 35 两个假定都有。但是尽管在这么多假定的情况下大部分人还是用mcmc，这个就比较奇怪
了。 |
|
F******n
发帖数: 160 | 36 第一个假定是很多人用，但只要是正太，线性模型，正常情况没人用MCMC（莫非为了教学，模拟演示？）。你在那本书或者资料，文献里看到这么用？我很好奇。
第二个假定应该很少有人用。比如最简单的，我上面说的，标准的卡尔曼滤波就没有这个假定。 |
|
F******n
发帖数: 160 | 37 我真的不知道你想要问什么答案了。。。
这样吧，你把我上面给的West和Harrison的书里第四章看看吧，尤其你感兴趣的那个“
FORECAST DISTRIBUTIONS”（第4节）。全部结果解析可解的，哪有什么MCMC?
你的那个MLE的“简单”解法也贴出来吧，否则光说，谁也没法猜出是怎么回事。 |
|
m**s
发帖数: 346 | 38 如下的bayesian network，
(1)写出node D的complete conditional equation. 不用化简，因为最终的
distribution is not recognizable.
(2) If node B is observed, and F is queried. 哪些node可以pruned.
我心里有个答案，但是不确认，如果你对你的答案很有把握，请联系我。
在线等，急！ |
|
w**********n
发帖数: 29 | 39 谢谢大家的评论。我是C++非常入门的新手，开始以为C++里面有这样的library直接用
就好了。看来没那么简单啊。。。。
因为在R里要做MCMC, 所以比较慢，所以想在C++里面写循环。但是这一块程序里面需要
用nlminb,所以才想问C++如何写。
各位若是有这方面的经验，请指教啊~~谢谢啦。 |
|
c*****a
发帖数: 808 | 40 大学修stat & econ, 美国大学太宽松，糊糊涂涂混了毕业。
未来两年想认真读ms，我不想再后悔没好好读书了。
想做银行的credit risk。一大串问题：
银行都是用什么的？generalized linear model?logistic reg?time series?survival
analysis?要特意修什么课吗？有没有推荐的书看吗？
顺便问一下这课对银行工作有帮助吗
Bayesian Statistics
Bayes Theorem, subjective probability, conjugate priors, non-informative
priors, posterior estimation, credible intervals, prediction, sensitivity
analysis, comparison to classical procedures, MCMC, Gibbs sampling,
hierarchical Bayesian analysis. Use of statistical software. |
|
|
f*******s
发帖数: 57 | 42 Sam has done some interesting work in MCMC, e.g. EE sampler. Also it's
unfair to say what he did in single molecule models are purely applied. His
work openned some new directions in this area, which has the same impact as
many of the great theoretical breakthrough. Don't be so cynical. |
|
|
w***e
发帖数: 35 | 44 不知道我的理解对不对。
过去做Bayesian Analysis的时候，可以用gibbs sampling, M/H sampling, and slice
sampling等方法针对 posterior distribution 进行采样。对于采样的结果是否符合
要求，通常会对采样的结果做一个Convergence的测试。如果证明Converge了，就认为
采样的结果可以用于进一步的分析。
的， |
|
R*******c
发帖数: 249 | 45 这个只是看是否converge吧，我想知道的是，它是否converge到所需要的分布
slice |
|
s*****9
发帖数: 108 | 46 你说的这些个常见的算法，只要converge了，得到的dist肯定是posterior dist啊，证
明都已经有了的 |
|
l*********s
发帖数: 5409 | 47 Wrong. And the real hard problem is to actually prove convergence. |
|
r********0
发帖数: 65 | 48 用 Gelman-Rubin Diagnostic看收不收敛 |
|
R*******c
发帖数: 249 | 49 可否给些文章看看呢？（关于得到的dist肯定是posterior dist) |
|
|
