rbprobitGibbs 帮助文档模拟输入值？

Question

当我查看上面的帮助文件时，它在 R 中给出了这个示例程序：

##
## rbprobitGibbs example
##
if(nchar(Sys.getenv("LONG_TEST")) != 0) {R=2000} else {R=10}

set.seed(66)
simbprobit=
function(X,beta) {
##  function to simulate from binary probit including x variable
y=ifelse((X%*%beta+rnorm(nrow(X)))<0,0,1)
list(X=X,y=y,beta=beta)
}

nobs=200
X=cbind(rep(1,nobs),runif(nobs),runif(nobs))
beta=c(0,1,-1)
nvar=ncol(X)
simout=simbprobit(X,beta)

Data1=list(X=simout$X,y=simout$y)
Mcmc1=list(R=R,keep=1)

out=rbprobitGibbs(Data=Data1,Mcmc=Mcmc1)

summary(out$betadraw,tvalues=beta)

if(0){
## plotting example
plot(out$betadraw,tvalues=beta)
}

那么，我的问题是，simbprobit 的代码是否只是为 rbprobitGibbs 函数准备了一些模拟值？也就是说，如果我正在计算 rbprobitGibbs 的结果，我是否可以有效地忽略 simbprobit 并简单地准备 Data、Prior、Mcmc 输入值？

Answer 1

是的，经过一番研究，我发现它是通过simbprobit模拟的随机输入，因为函数调用rnorm意味着它是一个随机变量。来自帮助文件：

"dnorm给出密度，pnorm给出分布函数，qnorm给出分位数函数，rnorm生成随机偏差。"

rbprobitGibbs 帮助文档模拟输入值？

rbprobitGibbs help documentation simulated input values?

r

bayesian