4.1.3版本的R语言中建立加权有向ergm模型,参照分布项设置为泊松分布,对模型分别使用MH采样函数以及Gibbs采样函数下的MCMC算法进行参数估计,迭代设置为1000次,绘制后验密度图
时间: 2024-02-06 08:11:28 浏览: 41
好的,以下是绘制后验密度图的参考代码:
```
# 载入相关包
library(ergm)
library(MCMCpack)
# 建立模型
model <- ergm(cbind(y, x) ~ edges + gwesp(0.5, fixed = TRUE),
constraints = "symmetric",
control = control.ergm(seed = 1234),
response = "Poisson")
# MH采样函数下的MCMC算法进行参数估计
set.seed(1234)
mcmc_mh <- MCMCmetrop1R(model, burnin = 100, mcmc = 1000)
# Gibbs采样函数下的MCMC算法进行参数估计
set.seed(1234)
mcmc_gibbs <- MCMCergm(model, control = control.ergm(seed = 1234),
constraints = "symmetric",
sampler = "MH",
mcmc = 1000,
burnin = 100)
# 绘制后验密度图
par(mfrow = c(1, 2))
plot(density(coef(mcmc_mh)), main = "Posterior Density Plot (MH)")
plot(density(coef(mcmc_gibbs)), main = "Posterior Density Plot (Gibbs)")
```
以上代码中,我们使用ergm包中的ergm函数建立加权有向ergm模型,参照分布项设置为泊松分布。然后,我们分别使用MH采样函数和Gibbs采样函数下的MCMC算法进行参数估计,迭代设置为1000次。最后,我们可以使用density函数绘制后验密度图。在这里,我们使用par函数将两个密度图并排放置。
相关推荐
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![exe](https://img-home.csdnimg.cn/images/20210720083343.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![docx](https://img-home.csdnimg.cn/images/20210720083331.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![pdf](https://img-home.csdnimg.cn/images/20210720083512.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)