请问大家有用ARS（adaptive rejection sampling）这个package的同学吗？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我现在在使用这个package的时候遇到了一个问题，就是在设定lower bound之后R就crash了，不知道是这个package的问题还是我的方程书写的有问题。在Mac和PC上都测试过未果，非常着急。麻烦大家了。下面是我写的两个公式：

func_alpha <- function(alpha){

                prod1 <- 0

                for (i in 1:101){

                                prod1 <- prod1 + log((1 + exp(alpha * theta_spl - beta_spl[m])))

                }

                alpha * sum(theta_spl * sample[, m]) - prod1

        }

        func_alpha_prima <- function(alpha){

                prod2 <- 0

                for (i in 1:101){

                                prod2 <- prod2 + (theta_spl * exp(alpha * theta_spl - beta_spl[m])) / (1 + exp(alpha * theta_spl - beta_spl[m]))

                        }

                sum(theta_spl * sample[, m]) - prod2

        }

        ars(1, func_alpha, func_alpha_prima, ub = T, xub = 10) * rnorm(1,mu2,sigma2)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：rejection Sampling Adaptive Adaptiv package package

呃，你的lb在哪里...