numEff: Compute Numerical Standard Error and Relative Numerical Efficiency
Description
numEff computes the numerical standard error for the mean of a vector of draws as well as the relative
numerical efficiency (ratio of variance of mean of this time series process relative to iid sequence).
Usage
numEff(x, m = as.integer(min(length(x), (100/sqrt(5000)) * sqrt(length(x)))))
Arguments
x
R x 1 vector of draws
m
number of lags for autocorrelations
Value
stderr
standard error of the mean of x
f
variance ratio (relative numerical efficiency)
Warning
This routine is a utility routine that does not check the
input arguments for proper dimensions and type.
Details
default for number of lags is chosen so that if R = 5000, m =100 and increases as the sqrt(R).
References
For further discussion, see Bayesian Statistics and Marketing
by Rossi, Allenby and McCulloch, Chapter 3.
http://www.perossi.org/home/bsm-1