Displaying 1 result from an estimated 1 matches for "outfun".
2009 Sep 06
2
question about ... passed to two different functions
...for each parameter with Monte Carlo standard errors.
apply(out$batch, 2, mean)
sqrt(apply(out$batch, 2, function(x) initseq(x)$var.con) / out$nbatch)
### Now suppose I want Monte Carlo estimates of some function of
### the parameters other than the identity function. I write a function
### outfun that does that. Also suppose I want some extra arguments
### to outfun. This example is a bit forced, but I hit on natural
### examples with a new function (not yet released) that works like
### metrop but does simulated tempering.
outfun <- function(beta, degree) {
stopifnot(is.nume...