posterior0=function (n, init, prior, lh=c("gev"),...,psd,burn = 0,thin = 1)
{
ar <- list(...)
if (lh == "gev") {
nas <- !is.na(ar$data) ##判断空值
ar$data <- ar$data[nas]
if (!is.null(ar$trend))
ar$trend <- ar$trend[nas]
}
initar <- c(list(par = init, prior = prior, lh = lh), ar)
inittest <- do.call("dpost", initar)
if (is.infinite(inittest))
stop("density is zero at initial parameter values")
ar <- c(list(n = n, init = init, prior = prior, lh = lh),
ar, list(psd = psd, thin = thin, burn = burn))
do.call("gibbs", ar)
}
哪们大神给理解下,那个... ar <- list(...) (ar$trend) (ar$data) 分别是 什么么