我现在在使用这个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)
扫码加好友，拉您进群
全部版块
我的主页
收藏
