求各位R大神打救小弟。。。。EM算法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

MC算法

1、假定X来自Exp(1), 而Y|X来自N(0,X), 现在已知Y=4，利用MC方法求E(X|Y=4)以及Var(X|Y=4) （从后验分布里抽取10000份样本，计算样本均数和样本方差）。

贝叶斯公式：（在下方）

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：EM算法 贝叶斯公式 样本均数 EXP 贝叶斯 样本

MC(Monte Carlo)

     生成指定分布的随机数的抽样。
1. Markov Chain （马尔科夫链）
     假设 f(t) 是一个时间序列，Markov Chain是假设f(t+1)只与f(t)有关的随机过程。
     Implement in R:
03.  
04.N = 10000  
05.signal = vector(length = N)  
06.signal[1] = 0  
07.for (i in 2:N)  
08.{  
09.    # random select one offset (from [-1,1]) to signal[i-1]  
10.    signal[i] = signal[i-1] + sample(c(-1,1),1)   
11.}  
12.  
13.plot( signal,type = 'l',col = 'red')