新手问个很弱的问题：用R如何抽样啊

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我是个R新手，请各路高手指点用R 怎么对核密度函数进行MCMC抽样啊（核密度函数形式很复杂，很难积分！）急！急！急！


不甚感谢啊！！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：核密度函数 密度函数 mcmc 高手指点 函数形式 新手

x <- out$x
y <- out$y
lupost <- function(beta, x, y) {
    eta <- as.numeric(x %*% beta)
    logp <- ifelse(eta < 0, eta - log1p(exp(eta)), -log1p(exp(-eta)))
    logq <- ifelse(eta < 0, -log1p(exp(eta)), -eta - log1p(exp(-eta)))
    logl <- sum(logp[y == 1]) + sum(logq[y == 0])
    return(logl - sum(beta^2)/8)
}   
在完成上面数据以及函数的定义(它们都是mcmc算法的输入参数)之后，我们就可以利用mcmc包中的metrop()来产生随机数据，模拟模型的后验分布。也就是说，我们将要从后验分布中进行取样。
   
set.seed(42)
beta.init <- as.numeric(coefficients(out))
out <- metrop(lupost, beta.init, 1000, x = x, y = y)
names(out)