我很努力地理解为什么这两段代码会产生截然不同的结果:
首先是一些功能:
1
现在输入代码:
n <- 5000
alpha1 <- rep(NA, n)
Y <- c(replicate(n, sample(c(-1, 1), 1, prob=c(1/2, 1/2))))
conditional_sampler <- function(x){
sample(c(-1, 1), 1, prob = c(x, 1-x))
}
t <- .4
lower <- function(t){
runif(1,(.4-t),.4)
}
upper <- function(t){
runif(1,.6,(.6+t))
}
vs。
alpha1 <- ifelse(Y==1, conditional_sampler(lower(t)), conditional_sampler(upper(t)))
任何帮助将不胜感激。